Khan-Yan2020 Article FiniteElementAnalysisOnSeismic
Khan-Yan2020 Article FiniteElementAnalysisOnSeismic
Khan-Yan2020 Article FiniteElementAnalysisOnSeismic
Received: 1 September 2019 / Accepted: 29 January 2020 / Published online: 5 February 2020
© Korean Society of Steel Construction 2020
Abstract
In prefabricated modular steel building, joint plays an important role in the structural behaviour and transfer of load among
connected modules. In this paper, a new type of bolted joint with tenon-gusset plate as horizontal connection and long beam
bolts as vertical connection has been proposed. The nonlinear structural behaviour of the novel joint was analysed against
static and quasi-static cyclic loadings by using finite element (FE) software ABAQUS. Numerical study was carried out to
predict the lateral load carrying capacity, bending moment capacity, rotational stiffness, failure modes and seismic behaviour
of the joint. Connection was simplified with nonlinear spring connector and 3-D beam elements. Parametric study was carried
out on the length of column tenon and gap between adjacent modular units. The results revealed that against lateral loads,
the gap was generated between the upper and lower parts of joint which resulted in the bearing failure of floor beam (FB).
The length of column tenon has obvious while gap between modules showed marginal effect on load carrying capacity and
structural behaviour of joint. The simplified joint accurately mimicked load–displacement curve and structural behaviour of
detailed joint. Then, FE analysis was justified by analysing test results of innovative connection in the reference. With these
evidences, the veracity of the FE analysis for studying nonlinear behaviour of novel joint was confirmed.
Keywords Prefabricated modular steel construction · Finite element analysis · Novel joint · Non-linear analysis · Load–
displacement curve · Parametric study
List of Symbols MPa Mega-Pascal = N∕mm2
P Pretension of bolt Pa Pascal = N∕m2
Ae Effective area of bolt 𝜈 Poisson’s ratio
ftv Tensile strength of bolt fy Yield strength
Δyield Lateral yield displacement fu Ultimate strength
Δ Lateral displacement 𝜀 Strain
Lx Distance of neutral axis from top of column kN 103 N
∅ Rotation P Lateral load
Es Modulus of elasticity m meter
𝜇 Friction coefficient mm millimetre
PTest Ultimate lateral load capacity of test specimen
Abbreviations
PFE Ultimate lateral load capacity of FE model
Cov Coefficient of variation
PFMS Prefabricated modular steel
MSB Modular steel building
* Jia‑Bao Yan
[email protected] HSS Hollow structural section
FE Finite element
Kashan Khan
[email protected] CB Ceiling beam
FB Floor beam
1
School of Civil Engineering, Tianjin University, FS Floor stringer
Tianjin 300350, People’s Republic of China CS Ceiling stringer
2
School of Civil Engineering/Key Laboratory of Coast Civil CP Cover plate
Structure Safety of Ministry of Education, Tianjin University,
Tianjin 300350, People’s Republic of China
13
Vol:.(1234567890)
International Journal of Steel Structures (2020) 20(3):752–765 753
13
754 International Journal of Steel Structures (2020) 20(3):752–765
HSS FS
HSS column
HSS CB B
HSS CS
HSS beam
13
International Journal of Steel Structures (2020) 20(3):752–765 755
connection is inserted in column tenon of upper part and is lower parts of the connection. As the gap between FB and
bolted with column. Each part of the connection supports GP keeps increasing, the FB keeps bending and faces bear-
ceiling beams (CBs) and FBs with beam tenons and bolts. ing stresses against beam bolts. With the gradual increase of
Except top story, the lower part of joint mainly consists of uplift force, the movement of column against tenon increases
two solid column tenons; upper solid tenon which is inserted which exerts shearing force on the column bolts.
in the column tenon of upper part and lower solid tenon is From literature review, it has found that connection in
inserted in lower column. The GP in exterior corner con- PFMS building is the most critical part in load distribution.
nection only provides sufficient gap between the upper and Therefore, structural behaviour of the modular units depends
lower modules. But in the exterior and interior middle con- on the connection system and overall behaviour of PFMS
nection as shown in the Fig. 1b, the GP has two functions; building.
one is to provide vertical gap between the upper and lower
modules while second is to horizontally connect adjacent
modular units. 3 Experimental Studies on the MSB
As shown in Fig. 2a, when bending moment as uplift Connection
force is applied on the top of the column, the disturbance is
produced and gap starts generating between the upper and Chen et al. (2017b) reported experimental works to study
working mechanism, load carrying capacity and seismic
behaviour of the innovative joint in MSB. Monotonic and
Bending force
cyclic loading analyses were performed on six specimens
Bolts in tension in Structural engineering laboratory of Tianjin University.
13
756 International Journal of Steel Structures (2020) 20(3):752–765
reinforced specimens. By column end loading method, upper behaviour, material’s yield strength “ fy ”, ultimate strength
column was supposed free end constrained condition and “ fu ” and strain values “𝜀” for corresponding strengths were
100t lateral displacement was applied on the upper end of defined according to the values used in test as enlisted in
column, while the lower column and beams ends were pin- Table 1 (Chen et al. 2017b).
constrained. Quasi-static loading pattern adopted in the test,
was taken from Chinese standard, “Regulation of seismic
4.3 Formation of FE Model of Novel Joint
test method (JGJ101-96)”. For comparative study, two axial
force ratios (AFRs) “0.2” and “0.1” were used in the test.
As shown in Fig. 1b, simplifications in the FE modelling
were assumed for accurate representation of connection with
3.2 Material Properties
FE techniques. Accordingly, the bolt nut and head was mod-
elled circular and together with bolt shank while the threads
The HSS columns, FBs, CBs, cover plates (CPs) on beams,
on the bolt shank and nuts were neglected. The gap between
GP, long beam bolts and stiffeners according to Chinese
the bolt shank and bolt hole was neglected to avoid mod-
standards were made of Q345B steel whereas, the plug-in
elling complications. The flanges of beams and CPs were
device was made from ZG35 cast steel. The columns and
modelled together to reduce contact surfaces. The chamfered
beams were welded by groove welding whereas stiffeners
edges and cylindrical solid tenon were modelled as regular
by fillet welding with wire ER50-6. The material properties
geometries for the ease of mesh and to make the contact
used in the tests are enlisted in the Table 1.
surfaces coincident.
13
International Journal of Steel Structures (2020) 20(3):752–765 757
4.3.1.2 Simplified FE Model In simplified FE models of both other. To find accurate value, detailed study was performed on
exterior and interior joints, columns and beams were replaced three different values of friction coefficient i.e., 0.2, 0.3 and
by 3-D beam elements that were fined meshed (similar with 0.4. The FE analysis with friction coefficient of 0.2 resulted
detailed model) with 2-node linear beam in space element in lowering the load carrying capacity and 0.4 increasing than
type (B31) while connection were replaced by spring connec- experimental findings. Therefore, friction coefficient of 0.3
tor elements. was chosen for contact surfaces except bolts that were mod-
elled frictionless in FE analysis. In the exterior and interior
4.3.2 Interactions middle connections, the adjacent columns and joints of the
adjacent modular units were simulated with hard contact for
In detailed FE models, the contact between columns and con- accurate sharing of pressure.
nection, beams and connection, columns and bolts and beams
and bolts were simulated with surface to surface (standard) 4.3.3 Boundary and Loading Conditions
with finite sliding and “hard contact” in the normal direction,
while in the tangential direction as “penalty friction formula- In the experiment, movements were restrained in all direc-
tion”. The hard contact formulation interaction in ABAQUS tions at bottom end of the lower column, X and Z direction at
states that two interacting surfaces will share pressure if they top end of the upper column, while at CB and FB movement
contact each other, whilst if they are separated no pressure in Y direction and rotation in Z direction was restrained. In
will be transferred between the surfaces. In ABAQUS the ABAQUS/CAE, horizontal displacement controlled loading
penalty friction formulation allows slip between two or more of 100 mm and axial compression force calculated from AFR
contacting surfaces (bolts with beams, plates) relative to each was applied on the top end of the column (Chen et al. 2017b).
other and the frictional coefficient determines the friction Using Eq. 1, pretension force was calculated for bolt load.
force among interacting surfaces. This type of “hard contact
0.9 × 0.9 × 0.9
and penalty friction formulation” allows the relative motion P= Ae ftv (1)
1.2
of the components and avoids them from penetrating in each
150 50
120
0
P(kN)
P(kN)
P(kN)
100
80 -150 -100 -50 0 50 100 150
-50
50
40 S2 Test Results
S1 Test Results -100
S2 FE Results QS1 Test Results
S1 FE Results 0 QS1 FE Results
0
0 50 100 150 -150
0 30 60 90 120
∆(mm) ∆(mm) ∆(mm)
(a) Specimen S1 (b) Specimen S2 (c) Specimen QS1
150 200
220
100 150
160
100
50 100
50
P(kN)
0 40
P(kN)
P(kN)
13
758 International Journal of Steel Structures (2020) 20(3):752–765
where, P = Pretension force; Ae = Effective area of bolt; 4.4.1 Validations of FE Model Against Uniaxial
ftv = Tensile strength of bolt = 180 MPa. and Quasi‑Static Loading Test
13
International Journal of Steel Structures (2020) 20(3):752–765 759
Table 3 Comparison of ultimate load of test specimens and FE 5.1.2 Exterior Corner Connection
Sp. no Loading method Ultimate load Ultimate load PTest/PFE
(Test)/PTest (FE)/PFE(kN) Figure 6a shows the lateral load carrying capacity of the
(kN) exterior corner connection. As load started increasing, gap
was generated between FB and GP, due to which FB experi-
S1 Static 112 118 0.95
enced bearing stresses against beam bolts and the yielding of
S2 Static 183 162 1.13
FB started at lateral displacement of 49.5 mm when lateral
QS1 QS +ve 81 79 1.02
load was 88.9 kN. It was found that FB faced bearing failure
QS −ve − 102 − 118 0.86
at ultimate lateral load and displacement of 119.773 kN and
QS1 QS +ve 81 79 1.02
90.3206 mm.
QS −ve − 102 − 118 0.86
QS2 QS +ve 118 121 0.97
5.1.2.1 Simplification of Exterior Corner Connection By
QS −ve − 137 − 127 1.07
using rotational stiffness obtained from FE analysis of
QS3 QS +ve 161 122 1.32
detailed model of joint, the novel connection was simplified
QS −ve − 183 − 162 1.13
with spring connector and 3-D beam elements. Figure 6c
QS4 QS +ve 141 129 1.09
and d compare the von Mises stress and load–displacement
QS −ve − 168 − 169 0.99
curves of detailed and simplified FE models which justify
Mean 1.05
that simplified model accurately predicted the load carrying
Cov 0.12
capacity and structural behaviour of detailed model of joint.
fracture in beams, local buckling and inward deformation 5.1.2.2 Classification of Novel Joint Based on the strength
in columns, tear in columns and gap formations showed by and stiffness calculation techniques provided by Choi and
the FE models are in well accordance with the test findings. Kim (2015) in the classification table enlisted in Euroc-
Table 3 compares the FE predicted ultimate load capaci- ode (EN 1993-1-1 2005), the novel joint was classified as
ties with the test specimens. It can be found that average test semi-rigid/partial strength joint.
to FE ultimate load capacity prediction is 1.05 and coeffi-
cient of variation (Cov) is 0.12 that are in well accordance
with the accuracy of FE analysis. The ratios of ultimate 5.1.3 Exterior Middle Connection
load capacities higher than 1.0 indicates that FE averagely
slightly underestimated whilst lower than 1.0 indicates that As shown in Fig. 7a, at application of horizontal mono-
FE averagely slight overestimated the ultimate load carry- tonic displacement and axial compression loading on top
ing capacities. The above validations assure the accuracy of mid of two columns, the gap was generated between FB
developed FE models to simulate basic structural behaviour and GP at lateral displacement of 53.4614 mm, and lateral
of novel joint. load of 250.925 kN, which resulted in generation of bear-
ing stresses and bearing failure of FB.
5 Numerical Studies on the Novel Joint 5.1.3.1 Simplification of Exterior Middle Connec‑
tion Detailed exterior middle connection was simplified
5.1 Geometry of the Joint by using same technique of spring connector and 3-D
beam elements. And it can be seen from Fig. 7b and c that
Figure 1b shows main components of novel joint that mainly simplified model has the tendency to accurately predict
consists of three components i.e., upper component, lower the load carrying capacity and structural behaviour of the
component and GP. detailed model with capacity ratio of 1.05.
Figure 5a shows the dimensions of the components of con- 5.2.1 Length of Column Tenon
nection. Three long beam bolts are inserted vertically and
each bolt is inserted in upper and lower column. The small As tenon is main component of connection so studying
scale specimens of approximately 2/3rd of the full scale effect of length of column tenon on the overall load car-
dimensions of the structural components (according to test rying capacity and structural behaviour of connection is
specimens) are modelled. The vertical gap of 50 mm is con- important. For the reason, FE analysis was carried out on
sidered between upper and lower modules for GP. four models with different lengths of column tenon such
13
760 International Journal of Steel Structures (2020) 20(3):752–765
as 100, 200, 300 and 400 mm denoted as L1 and 420, 520, for rigid connection. It has been found that model 1 shows
620 and 720 mm as L2, respectively as shown in Fig. 5b. column failure whereas, models 2, 3, and 4 show bearing
L1 is the length of the column tenon that supports column failure of FB as shown in Fig. 8. The increase in length of
on both upper and lower part, whereas, L2 is the length of column tenon resulted in obvious local web buckling of both
intermediate column tenon that rests in upper part of joint FB and CB but decrease in column damage because of least
13
International Journal of Steel Structures (2020) 20(3):752–765 761
P(kN)
100 88.901kN
FE P-Delta
50 Ultimate
Yield
0
0 50 100 150 200
Δ(mm)
Von Mises stress Load-displacement curve
(a) Behaviour against monotonic loading
100
P(kN)
50
Detailed FE Model
Simplified FE Model
0
0 50 100 150
∆ (mm)
(d) Comparison of load-displacement curves of detailed and simplified FE models
13
762 International Journal of Steel Structures (2020) 20(3):752–765
P(kN)
200
FE P-Delta
100
Yield
Ultimate
0
0 60 120 180
∆ (mm)
400
300
P(kN)
200
0
0 50 100 150
∆ (mm)
(c) Comparison of load-displacement curves of detailed and simplified FE models
column rotation. It can be seen that with increase in length 5.2.2 Gap Between Adjacent Modular Units
of column tenon, the gap between the column and tenon is
reduced which minimizes the risk of column failure at bolted In MSB, joint connects eight modular units together. So,
area while FB and CB fails prior to the failure of column as in order to study the effect of modular units on structural
seen in models 2, 3 and 4. In the model 1, the beams were behaviour of adjacent modular units, 31 models with vary-
not failed but the column failed as a result of which the ing gap from no space to 30 mm were studied under lateral
structure failed prematurely. monotonic loading. Although with increase in gap between
Moreover, from Fig. 10, it can be seen that as the length modular units, a slight rise in stresses was observed but the
of column tenon increases, initial stiffness and load carry- failure mode was still remained same due to bearing failure
ing capacity increase obviously but increment in capacity of FB as shown in Fig. 9. Figure 10 shows load–displace-
followed decreasing pattern from 50% (between model 1 ment curves of only four models (with variation of gap of
having length of 100 mm and model 2 having length of 200 10 mm) which showed marginal increase in capacity up to
mm) to 14% (between model 3 having length of 300 mm and 3.4% between model “having no gap” and the model “with
4 having length of 400 mm). gap of 30 mm”.
13
International Journal of Steel Structures (2020) 20(3):752–765 763
(a) (b)
(c) (d)
(a) (b)
(c) (d)
5.3 Quasi‑Static Analysis connection under cyclic loading loses its strength gradu-
ally and showed good ductile behaviour with proper energy
The quasi-static analysis using ATC-24 loading protocol dissipation (Fig. 6b).
(Krawinkler 1992), was carried out with limited number
of cycles to observe the seismic behaviour of the exterior
corner connection. Figure 6b shows the lateral load–dis-
placement capacity that is in good accordance with the
monotonic loading results. It can be clearly observed that
13
764 International Journal of Steel Structures (2020) 20(3):752–765
13
International Journal of Steel Structures (2020) 20(3):752–765 765
Chen, Z., Liu, J., & Yu, Y. (2017a). Experimental study on interior Lawson, M., Ogden, R., & Goodier, C. (2014). Design in modular con-
connections in modular steel buildings. Engineering Structures, struction. Boca Raton: Taylor & Francis Group, ed., CRC Press.
147, 625–638. Lawson, R. M., Ogden, R. G., & Bergin, R. (2012). Application of
Chen, Z., Liu, J., Yu, Y., Zhou, C., & Yan, R. (2017b). Experimental modular construction in high-rise buildings. Journal of Architec-
study of an innovative modular steel building connection. Journal tural Engineering, 18(2), 148–154.
of Constructional Steel Research, 139, 69–82. Lawson, R. M., & Richards, J. (2010). Modular design for high-rise
Choi, K., & Kim, H. (2015). An analytical study on rotational capac- buildings. Proceedings of the Institution of Civil Engineers—
ity of beam-column joints in unit modular frames. International Structures and Buildings, 163(3), 151–164.
Journal of Civil, Structural, Construction and Architectural Engi- Liu, X. C., Pu, S. H., Zhang, A. L., & Zhan, X. X. (2017). Performance
neering, 9(2), 100–103. analysis and design of bolted connections in modularized prefab-
Ding, Y., Deng, E. F., Zong, L., Dai, X. M., Lou, N., & Chen, Y. ricated steel structures. Journal of Constructional Steel Research,
(2017). Cyclic tests on corrugated steel plate shear walls with 133, 360–373.
openings in modularized-constructions. Journal of Construc- Liu, X. C., Xu, A. X., Zhang, A. L., Ni, Z., Wang, H. X., & Wu, L.
tional Steel Research, 138, 675–691. (2015). Static and seismic experiment for welded joints in modu-
EN 1993-1-1. (2005). Eurocode 3: Design of steel structures - Part larized prefabricated steel structure. Journal of Constructional
1-1: General rules and rules for buildings. Eurocode 3. Steel Research, 112, 183–195.
Fathieh, A., & Mercan, O. (2016). Seismic evaluation of modular Liu, X. C., Zhan, X. X., Pu, S. H., Zhang, A. L., & Xu, L. (2018).
steel buildings. Engineering Structures, 122, 83–92. Seismic performance study on slipping bolted truss-to-column
Hwan Doh, J., Ho, N. M., Miller, D., Peters, T., Carlson, D., & Lai, P. connections in modularized prefabricated steel structures. Engi-
(2017). Steel bracket connection on modular buildings. Journal neering Structures, 163, 241–254.
of Steel Structures & Construction, 02(02), 0437–2472. Yan, J.-B. (2015). Finite element analysis on steel–concrete–steel sand-
Kamali, M., & Hewage, K. (2016). Life cycle performance of mod- wich beams. Materials and Structures, 48(6), 1645–1667.
ular buildings: A critical review. Renewable and Sustainable
Energy Reviews, 62, 1171–1183. Publisher’s Note Springer Nature remains neutral with regard to
Krawinkler, H. (1992). ATC-24 Guidelines for cyclic seismic testing jurisdictional claims in published maps and institutional affiliations.
of components of steel structures. California: Redwood City.
13