# An Analytical Method for Determining the Load Distribution of Single-Column Multibolt Connection

Authors: Nirut Konkong; Kitjapat Phuvoravan
Journal: Advances in Civil Engineering
 (2017)

---

## Abstract

The purpose of this research was to investigate the effect of geometric variables on the bolt load distributions of a cold-formed steel bolt connection. The study was conducted using an experimental test, finite element analysis, and an analytical method. The experimental study was performed using single-lap shear testing of a concentrically loaded bolt connection fabricated from G550 cold-formed steel. Finite element analysis with shell elements was used to model the cold-formed steel plate while solid elements were used to model the bolt fastener for the purpose of studying the structural behavior of the bolt connections. Material nonlinearities, contact problems, and a geometric nonlinearity procedure were used to predict the failure behavior of the bolt connections. The analytical method was generated using the spring model. The bolt-plate interaction stiffness was newly proposed which was verified by the experiment and finite element model. It was applied to examine the effect of geometric variables on the single-column multibolt connection. The effects were studied of varying bolt diameter, plate thickness, and the plate thickness ratio (t2/t1) on the bolt load distribution. The results of the parametric study showed that the t2/t1 ratio controlled the efficiency of the bolt load distribution more than the other parameters studied.


---

## Body

## 1. Introduction

A bolt connection is generally used in cold-formed steel structures. It can result in a high-stress concentration at contact regions, especially with thin cold-formed steel members. The shearing of cold-formed steel, bearing on material in front of the bolt, tearing of the cold-formed steel in the net section, and shearing of the bolt are the main failure modes of a single bolt connection. In multibolt connections, the failure modes are identified by numerical and analytical methods, with the results showing combined bearing failure and net-section failure modes [1, 2]. The failure behavior of a multibolt connection is considered in two steps as (a) the load distribution between bolts, which is the process for considering the stress bypass in the multibolt connection, and (b) local stresses analysis as the bearing stress and net-section stress around a critical hole, which is used to investigate ensuring a safety margin of connection [3–5]. Practically, the lap shear connection is used for tension in a structural member such as the diagonal and tension chord truss element. It is designed using a multiple bolt connection with even bolt load distribution. The bolts in a row in the loading direction are assumed to have an equal distribution of load among the bolts. In an actual lap connection, the behavior of each bolt in a row is different due to an unequal load distribution that also results from the stiffness of the single bolt deformation and plate deformation between bolts as shown in Figure 1, where the evolution of the internal forces versus the applied force (P) is demonstrated with a five-bolt connection in which δ and δb-i are the connection elongation and bolt deformation.Figure 1
Stages of force distribution in bolt connections [6].In the first step, the connection evolves elastic behavior with the outer bolts being more loaded than the inner ones. In the second step, the outer bolt rows reach their elastic limit strength but the inner bolt rows are still in the elastic stage. In the third step, the plastic behavior is distributed to the inner bolt rows. In the fourth step, the plastic behavior is developed in all the bolt rows. However, the limitation of the cold-formed ductility between the bolt-plate in bearing may be reached before the connection develops full plastic strength. These failures result in a lower average bearing stress value per bolt hole of the multibolt connection. In this case, the connection would not develop the ultimate bearing capacity as recommended by the design method [6, 7]. Bolt load prediction is a complex procedure involving the establishment of the bolt load distribution. The stiffness components (the axial stiffness of the plate and the interaction of the bolt-plate) must be included in the investigation of bolt load distribution. However, the interaction stiffness of the bolt-plate is not easily defined in a closed-form solution because the analysis requires knowledge of the influenced parameters such as the shear deformation of the bolt, bending deformation of the bolt, bearing deformation of the plate, bearing deformation of the bolt, and deformation of the plate between bolts. Thus, the equation of the bolt-plate interaction stiffness in industry today is derived by using a semiempirical approach based on the testing data [8–10]. The weakness of these empirical equations is that they ignore the shear deformation assumption in the bolt component which is dominated by the short bolt length.This study describes the method of analysis for determining the load distribution among bolts in a single-column multibolt connection. A bolt-plate interaction stiffness equation is proposed which is verified by the experimental test and finite element analysis (FEA). The experimental investigation was performed using lap shear bolt connection testing and the connections were fabricated from cold-formed steel.


## 2. Material Properties and Model Description

The laboratory test was performed using G550 cold-formed steel with a plate thickness of 1.0 mm for the bolt connection test. G550 is a hot-dipped zinc-coated structural steel and has a minimum yield strength of 550 MPa. The material properties for G550 steels were specified by following ASTM Standards A370-07 [11] as shown in Table 1. Following the research objective and scope, the possibility of bolting shear failure was eliminated by using M5 Grade 8.8 steel bolts conforming to ASTM A325M [12] as shown in Table 2.Table 1
Material properties of cold-formed steel, G550-1T.
Test number
SD
COV
CT-1
CT-2
CT-3
CT-4
CT-5
Mean
T (mm)
1.045
1.045
1.044
1.045
1.045
1.04
0.0004
0.04
W (mm)
12.64
12.65
12.59
12.58
12.52
12.59
0.0522
0.41
F
y (MPa)
608.47
605.77
605.90
617.50
606.46
608.82
4.9715
0.81
F
u (MPa)
627.00
622.50
625.50
630.00
620.00
625.00
3.8890
0.62
E (GPa)
213.51
213.49
213.49
213.58
213.58
213.53
44.1757
0.02
F
u
/
F
y
1.03
1.02
1.03
1.02
1.02
1.02
0.0054
0.53
Remark. Fy is the yield stress, Fu is the ultimate stress, E is the modulus of elasticity, T is the plate thickness, W is the width of coupon test specimen, COV is the coefficient of variation, and SD is the standard deviation.Table 2
Material properties of bolts.
Test number
SD
COV
BT-1
BT-2
BT-3
BT-4
BT-5
Mean
d (m)
4.93
4.97
4.97
4.97
4.94
4.95
0.019
0.393
F
y (MPa)
930.00
949.00
927.00
976.00
872.00
930.80
38.232
4.107
F
u (MPa)
963
930
963
997
909
952.40
33.908
3.560
E (GPa)
204.08
204.08
204.09
204.10
204.09
204.08
7.385
0.003
Remark. d is the bolt diameter.The bolt connection specimen and the hole size are shown in Figure2. The bolts and nuts were assembled with integral washers that were tightened by hand, which allowed for slip in the connection after the initial loading.Figure 2
Bolt connection specimens.
(a)
Bolt connection geometries
(b)
Hole sizeThe bolt connection specimens were tested using a universal testing machine (UTM) in the structural laboratory of Kasetsart University, Bangkok, Thailand, as shown in Figure3(a). A typical bolted shear connection test setup is shown in Figure 3(b). The load and the axial deformation were measured to describe the connection failure behavior. The load was applied to the specimen using displacement control with a speed of 1 mm/minute. A linear variable differential transformer (LVDT) was installed to measure the elongation of a single bolt connection over a distance of 200 mm from the center line. The load was gradually applied until either the specimen failed or the specimen entered into plastic deformation. Once the applied load began to drop, even though the specimen continued to elongate, the test was stopped.Figure 3
Bolted connection test setup.
(a)
Test apparatus
(b)
Test specimen


## 3. Finite Element Model and Validation

In the present study, the ANSYS finite element package [13] was used to predict the failure behavior of the bolt connection as shown in Figure 4. Nonlinear finite element analysis was used to compare and validate the results of the single bolt connection test. The nonlinear geometry, material nonlinearity, and contact problem were modeled in the finite element model. The cold-formed steel was modeled using eight nodes with six degrees of freedom at each node shell elements (SHELL281) as shown in Figure 4(b).Figure 4
Geometry, element type, and material model of FEA.
(a)
Geometry model
(b)
SHELL281
(c)
SOLID186
(d)
Contact element
(e)
Pure penalty algorithm
(f)
Mesh and boundary condition
(g)
Actual stress and plastic strainThree-dimensional, 20-node solid elements (SOLID186) with three degrees of freedom at each node were used to model the steel bolts and the washers, as shown in Figure4(c). Contact between the steel plates, steel bolts, washers, and wall of the holes was defined as the contact elements (CONTA173 and TARGE170). The pure penalty contact algorithm was used to solve the constrained optimization contact problem as shown in Figures 4(d)-4(e). The mesh size was controlled by the aspect ratios (length-to-width ratio) which were selected to be close to 1.0 for all components. Around the bolt holes, fine mesh was utilized in order to transfer the stress from the bolt to the cold-formed steel plate. The nodes along the end plates were restrained with geometric boundary conditions as pinned and roller support. The uniform load was applied along the end plate line B with the displacement control as shown in Figure 4(f). The static friction between the steel plates was neglected due to the smooth coating surfaces. The actual stress-strain values from the cold-formed steel tensile testing were converted to the true stress (σtrue) and strain (σtrue) which imposed a multilinear isotropic hardening model as shown in Figure 4(g). It was converted using(1)σtrue=σtest1+εtest,εtrue=ln⁡1+εtest.The large deflection option was activated to achieve nonlinear geometry behavior. By turning on this option, iterative solving and updating of the stiffness matrix based on the incremental nodal displacements at each equilibrium iteration were performed into the nonlinear geometry analysis. Finally, the Newton-Raphson algorithm was activated to obtain the solution after each displacement increment.


## 4. Single Bolt Connection Test Results

The load capacity and failure mode of single bolt connection tests are shown in Figure5(a). The deformation shape from the experimental tests and the FEA showed that the specimens were subjected to bearing failure as shown in Figure 5(b). The ratio of the ultimate loads between the experiment and the FEA (Pn,test/Pn,FEA) is shown in Table 3 and this indicates that the ultimate load capacity of the experimental results was in good correlation with the FEA.Table 3
Summary of the lap shear connection results.
Test
Loading capacity,Pn,test (N)
Failure mode
P
n
,
t
e
s
t
/
P
n
,
F
E
A
S-1
10,812
Bearing
0.97
S-2
10,660
Bearing
0.95
S-3
10,707
Bearing
0.96
S-4
10,800
Bearing
0.97
S-5
10,758
Bearing
0.96
Mean
10,747
0.96
FEA
11,167
BearingFigure 5
Load-deformation curves and failure mode of single bolted connection.
(a)
Load-deformation curves
(b)
Failure modeIn order to identify the connection failure modes, the Von-Mises stress contour was normalized by the ultimate stress of the cold-formed steel material as shown in Figure6.Figure 6
Stress normalization and strain distribution of bolt connection.Hart-Smith [14] and He and Wang [2] have suggested a method for bearing failure prediction of the shear connection. The monitoring strain point at C, D, and E is zero when the connection suffers from bearing failure, and the monitoring stress normalization at point C is equal to 1 for net-section failure. The stress-strain analysis of the FEA results showed that the lap shear bolt connection suffered from bearing failure.


## 5. Bolt-Plate Interaction Stiffness Analysis

In multibolt connections, the bolts in different rows carry different amounts of the load which depend on the bolt diameter, plate stiffness, and other geometry configurations [15, 16]. The connection stiffness is fundamental to understanding the load transfer and the load distribution of the bolt connection. In accordance with the typical bolt, the deformation of the steel bolt, illustrated in Figure 7, was simplified as a short beam. The concept of bolt flexibility (Cb) prediction can be defined by(2)Cb=1kb=δbF,where kb is the stiffness of the steel bolt,F is the applied force, and δb is the bolt deformation.Figure 7
Deformation of single bolt connection.An analytical beam theory with the shear deformation and rotational bending effects included is the theory known as Timoshenko beam theory [17]. For Timoshenko beams, plane cross-sections will rotate due to shear forces. The rotation (θ) and deflection of a beam (u) are written by(3a)θx=∫0xMsEIds+θ0,(3b)ux=∫0xθsds-∫0xVsκAGds+u0,where V0, M0, θ0, and u0 are constants of integration which can be determined using the boundary conditions for the particular problem, A is area cross-section of the beam, G is a shear modulus of the beam, and κ is the shear coefficient due to the nonuniform distribution of shear stress along the cross-section. Figure 8 shows the assumption of the load components of a steel bolt which was used to analyze the bolt-plate interaction stiffness. Point A, point B, and point C are the reference points on the bolt that refer to the middle of the plate thickness and the bolt length.Figure 8
Forces in bolt.The bolt tolerated the bearing force, shearing force, and bending moment initiated by an eccentric load from the steel plates. The statically indeterminate system was divided in Figure9.Figure 9
Statically indeterminate system.A cantilever beam with partially uniform loads and end moments idealizes the load components of a steel bolt. The superposition analysis method was used to analyze the bending deformation and shear deformation for a cantilever beam carrying a uniform load of intensity over part of the span as shown in Figures10 and 11.Figure 10
Bending deformation of cantilever beam.Figure 11
Shear deformation of cantilever beam.The bending moment was introduced by eccentric force from the plate. It was derived by sum of the rotation of the neutral axis at the end of the bolt as shown in Figure12. This rotation can be rewritten by(4)ϕ=ϕ1+ϕ2+ϕ3=θ1+γ1+θ2+γ2+θ3+γ3=0,where θ1, θ2, and θ3 are the bending rotation components due to the loads ω1t1, ω2t2, and M, respectively. γ1, γ2, and γ3 are shear rotation due to the loads ω1t1, ω2t2, and M, respectively.Figure 12
Rotation of the bolt.The bending rotation and shear rotation components can be written as (5a)-(5b) which refer to Figures 10-11.(5a)θ1=-ω1t136EbIb,θ2=ω2L3-t136EbIb,θ3=MBLEbIb,(5b)γ1=0,γ2=0,γ3=0,where Eb is the modulus of elasticity of the bolt, Ib is the moment of inertia of the bolt, t1 and t2 are the thickness of the plate, ω1 and ω2 are the loads on the bolt, and L is the thickness of the connection. The force in the bolt (Fb) can be presented by(6)Fb=ω1t1=ω2t2.The relationship of the uniformed load in the connection can be written by(7)ω1=ω2t2t1.Thus, (4) can be rewritten as (8)-ω1t136EbIb+ω2L3-t136EbIb+MBLEbIb=0.The bending moment in the bolt can be rewritten as(9)MB=-ω2t1+t23-t13-t12t26t1+t2.The connection stiffness calculation was divided into two main stiffness components being the bolt stiffness (kb) and the contact stiffness (kc) as shown in Figure 13. The bolt stiffness was analyzed by considering the flexural stiffness (kflex,b) and the shear stiffness (kshe,b) of the bolt. The contact stiffness was the interaction stiffness between the plate and the bolt which is referred to as the plate bearing stiffness (kbea,p) and the bolt bearing stiffness (kbea,b).Figure 13
Stiffness components of single bolt connection.The bolt-plate interaction stiffness (kbp) is associated with the bending, shearing, and bearing stiffness of the assembly parts. It can be written as the spring components using (10a)-(10b): (10a)1kbp=1kb+1kcwith(10b)1kb=1kflex,b+1kshe,b,1kc=1kbea,p+1kbea,b.In the bending stiffness component, the bolt deformation was analyzed using (3b) which ignored the shear deformation term as written by (11). θ0 and u0 were set equal to zero for the cantilever beam boundary condition.(11)ux=∫0xθsds+u0=∫0x∫0xMsEIds+θ0ds+u0.The flexural deformation of the bolt was calculated by sum of deformation at pointA and point B as shown in Figure 14.Figure 14
Flexural deformation of the bolt.The flexural deformations of the bolt at pointA (δflex,A) can be written as(12)δflex,A=δflex,A,ω1t1+δflex,A,ω2t2+δflex,A,M,where δflex,A,ω1t1, δflex,A,ω2t2, and δflex,A,M are the flexural deformation at point A due to ω1t1, ω2t2, and M, respectively, which are defined by(13a)δflex,A,ω1t1=-ω1t1324EbIb6a2-4ax+x2,(13b)δflex,A,ω2t2=ω2bx212EbIb3L+3a-2x,(13c)δflex,A,M=Mx22EbIb,(13d)a=t1,b=t2,L=t1+t2,x=t12.The deformations of the bolt at point B (δflex,B) can be written by (14)δflex,B=δflex,B,ω1t1+δflex,B,ω2t2+δflex,B,M,where δflex,B,ω1t1, δflex,B,ω2t2, and δflex,B,M are the flexural deformations at point B due to ω1t1, ω2t2, and M, respectively, which are defined as(15a)δflex,B,ω1t1=-ω1a324EbIb4x-a,(15b)δflex,B,ω2t2=ω224EbIbx4-4Lx3+6L2x2-4a3x+a4,(15c)δflex,B,M=Mx22EbIb,(15d)x=t1+t22.The total bolt flexural deformation due to ω1t1, ω2t2, and M was defined using (16)δflex,b=δflex,B-δflex,A=1128t1+t23t12+13ttt2+3t22EbIbω2t2,where Fb=ω2t2. Thus, the total bolt flexural deformation can be rewritten as (17)δflex,b=3t12+16t22t2+16t1t22+3t23128EbIbFb.The bolt bending flexibility (Cflex,b) is defined as the ratio of the bending deformation due to the force of the bolt which can be calculated using (18)Cflex,b=1kflex,b=δflex,bFb=3t12+16t22t2+16t1t22+3t23128EbIb.In the shear stiffness component, the bolt shear deformation is analyzed using (3b) in which the flexural deformation term is ignored as written by (19). u0 is set equal to zero for the cantilever beam boundary condition.(19)ux=-∫0xVsκAGds+u0.The Timoshenko beam theory [17] proposed that the shear stress is uniform over the cross-section. Actually, it is not uniform over the cross-section. Thus, the shear coefficient (κ) was introduced as a correction factor to allow the nonuniform shear strain to be expressed as a constant. Hutchinson [18] proposed the shear coefficient equation of a circular cross-section using (20)κ=61+ν27+12ν+4ν2,where ν is the Poisson ratio. For ν=0.3 for a steel bolt, the shear coefficient is 0.93. The shear deformation of the bolt was calculated by sum of deformation at point A and point B as shown in Figure 15.Figure 15
Shear deformation of the bolt.The deformations of the bolt at pointA (δs,A) can be written using (21)δs,A=δs,A,ω1t1+δs,A,ω2t2+δs,A,M,where ⁡δs,A,ω1t1, δs,A,ω2t2, and δs,A,M are the shear deformations at point A due to ω1t1, ω2t2, and M, respectively, which are defined using (22a)δs,A,ω1t1=∫0xγ1xdx=-3ω1t128κGbAb,(22b)δs,A,ω2t2=∫0xγ3xdx=ω2t1t22κGbAb,(22c)δs,A,M=0,(22d)x=t12,where Ab is the cross-sectional area of the bolt and Gb is the shear modulus of the bolt. The deformations of the bolt at point B (δB) can be written using (23)δs,B=δs,B,ω1t1+δs,B,ω2t2+δs,B,M,where δs,B,ω1t1, δs,B,ω2t2, and δs,B,M are the shear deformation at point B due to ω1t1, ω2t2, and M, respectively, which are defined by (24a)δs,B,ω1t1=∫0aγ1xdx+∫0xγ2xdx=-ω1t122κGbAb,(24b)δs,B,ω2t2=∫0aγ3xdx+∫0xγ4xdx=ω2t28t1+3t28GbAb,(24c)δs,B,M=0,(24d)x=t22.The total bolt shear deformation due to ω1t1, ω2t2, and M was defined using (25)δs=δs,B-δs,A=38t1+t2κAbGbω2t2,where ω1=ω2t2/t1. Thus, the total bolt shear deformation can be rewritten by (26)δs=38t1+t2κAbGbω2t2.The bolt shear flexibility (Csd,b) is defined as the ratio of shear deformation due to the force of the bolt which is calculated using (27)Csd,b=1ksd,b=δsFb=38t1+t2κAbGb.To calculate the bearing stiffness of the bolt and plate, the bearing stiffness of connection was calculated by using the schematic diagram shown in Figure 16. The bearing stress (σbr) was defined by the applied load (P) divided by the hole diameter and the plate thickness which is written as (28a). The bearing strain (εbr) was defined by the plate bearing displacement divided by the hole diameter which is written using (28b):(28a)σbea=Pdt,(28b)εbea=δbea,pd.The bearing stiffness of bolt hole was calculated by using linear stress-strain relationship as presented using (29a)σbea,p=Epεbea,p,(29b)kbea,p=Pδbea,p=tdEpd=tEp.Similarly, the bearing stiffness of the steel bolt was calculated in the same way using the same manner which was written as (30a)σbea,b=Ebεbea,b,(30b)kbea,b=Pδbea,b=tdEbd=tEb.Thus, the bolt bearing flexibility (cbea,b) and the plate bearing flexibility (cbea,p) can be calculated using (31)cbea,b=1kbea,b=1t1Eb+1t2Eb,(32)cbea,p=1kbea,p=1t1Ep,1+1t2Ep,2,where Ep1 and Ep2 are the modulus of elasticity of plate 1 and plate 2, respectively. Finally, the bolt-plate interaction stiffness was newly proposed by the equilibrium of interaction force between the bolt and the plate which can be written using (33)1kbp=3t12+16t22t2+16t1t22+3t23128EbIb+38t1+t2κAbGb+1t1Eb+1t2Eb+1t1Ep1+1t2Ep2.Figure 16
Schematic diagram of bearing deformation in the plate and bolt.
(a)
Bearing deformation in the plate
(b)
Bearing deformation in the boltA common method for modeling a connection stiffness association of the single bolt connection is represented by the spring model as shown in Figure17.Figure 17
The spring model of single bolt connection.
(a)
Connection model
(b)
Spring model
(c)
Plate stiffnessThe matrix equation represented a single bolt connection which can be written as(34)00F3=kbp+kp,2-kbp0-kbpkbp+kp,1-kp,10-kp,1kp,1u1u2u3.The plate stiffness was calculated in two zones as shown in Figure17(c). The plate zone represented the plate stiffness in front of the bolt. The bolt zone represented the plate stiffness alongside the bolt. The stiffness of the plate was found considering a plate subjected to uniform tensile load which is presented in Table 4.Table 4
Plate stiffness.
Stiffness
Stiffness equation
Value (N/mm)
1
k
p
,
1
1
k
1
-
1
+
1
k
1
-
2
=
1
ω
1
t
1
E
p
,
1
/
l
1
-
0.5
d
h
+
1
ω
1
-
d
t
1
E
p
,
1
/
0.5
d
h
1/63692
1
k
p
,
2
1
k
2
-
1
+
1
k
2
-
2
=
1
ω
2
-
d
t
2
E
p
,
2
/
0.5
d
h
+
1
ω
2
t
2
E
p
,
2
/
l
2
-
0.5
d
h
1/63692
1
k
b
p
Equation (33)
1/50660
Remark. kp,1 and kp,2 are the plate stiffness and Ep,1 and Ep,2 are the modulus of elasticity of plate 1 and plate 2.The stiffness, load vector, and the node displacement matrix have been implemented in the Maplesoft program [19]. The load-deformation curve of node 3 was compared with a single-shear bolt connection test result which is shown in Figure 18(a). The results show that the analytical method was in good correlation with the FEA and the experimental results. Figure 18(b) shows the fastener stiffness prediction by several known semiempirical methods comparing with the proposed equation. The Grumman equation [8] is an empirically derived fastener stiffness formula that was presented by the Grumman Aerospace Corporation as given in (35)1kbp=t2+t12Ebd+3.721t1Ep1+1t2Ep2.The fitting equation to the load-displacement curves of the single steel bolt connection type was proposed by Huth [9] which is presented in (36)1kbp=t2+t12d2/331t1Ep1+1t2Ep2+12t1Eb+12t2Eb.Figure 18
Comparison of the load-deformation curve of the proposed spring model.
(a)
Comparison with tests and FEA
(b)
Comparison with other equationsThe Swift equation [10] is a flexibility empirical equation that was established based on testing data of a single-shear metallic connection as given by (37)1kbp=1.6660.5dEL+ETplate-1EL+ETplate-2+0.861EL+ETt1plate-1+1EL+ETt2plate-2.The load comparison results showed that the stiffness values from the Grumman equation [8] and the Swift equation [10] were conservative, which does not accurately reflect the actual connection strength, especially for a single bolt connection. The Huth equation [9] matched the testing data but the shear deformation of the bolt and bending deformation of the bolt were ignored.


## 6. Bolt Load Distribution Analysis

An application of bolt-plate interaction stiffness analysis for the multibolt connection is presented by the triple bolt lap connection as in Figure19. The structural model was assembled using the 1 mm thickness of cold-formed steel and a 5 mm diameter steel bolt. The analytical application investigated the load applied distribution among the bolt rows due to relative deformation between the top plate (plate 1) and bottom plate (plate 2).Figure 19
Multibolt connection analysis.The matrix equation representing the multibolt spring model in Figure20 can be written as (38), and the stiffness components are presented in Table 5. (38)000000F=k1+k2+k3-k2-k30000-k2k2+k40-k4000-k30k3+k5+k6-k5-k6000-k4-k5k4+k5+k70-k7000-k60k6+k8-k80000-k7-k8k7+k8+k9-k900000-k9k9u1u2u3u4u5u6u7.The nodal deformation was calculated by inverting the stiffness matrix and multiplying by the load vector. The carried bolt load (P) was calculated using (39)P1=K8u6-u5,P2=K5u4-u3,P3=K2u2-u1.To validate the analytical model, the finite element model was generated for comparison with the analytical result as shown in Figure 21(a). An analysis was performed in which material failure was not involved in the study. Thus, the small applied load was modeled using a finite element model to investigate the load distribution behavior. The validation of the analytical model results is presented in Figure 21(b). The connection load and load distribution indicated by the analytical model were in good agreement with the FEA. The load distribution results in Figure 21(c) have been expressed as a percentage of the total connection load (F). It can be seen that the load was shared unevenly between each bolt. The uneven load distribution in the triple bolt connection was analyzed using the parametric study.Table 5
Calculation of stiffness.
Stiffness equation
Stiffness component
1
k
1
=
1
k
1,1
+
1
k
1,2
k
1,1
=
ω
2
t
2
E
p
2
l
d
,
k
1,2
=
ω
2
-
d
h
,
3
t
2
E
p
2
0.5
d
h
,
3
1
k
3
=
1
k
3,1
+
1
k
3,2
+
1
k
3,3
k
3,1
=
ω
2
t
2
E
p
2
l
c
,
k
3,2
=
ω
2
-
d
h
,
3
t
2
E
p
2
0.5
d
h
,
3
,
k
3,3
=
ω
2
-
d
h
,
2
t
2
E
p
2
0.5
d
h
,
2
1
k
4
=
1
k
4,1
+
1
k
4,2
+
1
k
4,3
k
4,1
=
ω
1
t
1
E
p
1
l
c
,
k
4,2
=
ω
1
-
d
h
,
3
t
1
E
p
1
0.5
d
h
,
3
,
k
4,3
=
ω
1
-
d
h
,
2
t
1
E
p
1
0.5
d
h
,
2
1
k
6
=
1
k
6,1
+
1
k
6,2
+
1
k
6,3
k
6,1
=
ω
2
t
2
E
p
2
l
b
,
k
6,2
=
ω
2
-
d
h
,
2
t
2
E
p
2
0.5
d
h
,
2
,
k
6,3
=
ω
2
-
d
h
,
1
t
2
E
p
2
0.5
d
h
,
1
1
k
7
=
1
k
7,1
+
1
k
7,2
+
1
k
7,3
k
7,1
=
ω
1
t
1
E
p
1
l
b
,
k
7,2
=
ω
1
-
d
h
,
2
t
1
E
p
1
0.5
d
h
,
2
,
k
7,3
=
ω
1
-
d
h
,
1
t
1
E
p
1
0.5
d
h
,
1
1
k
9
=
1
k
9,1
+
1
k
9,2
k
9,1
=
ω
1
t
1
E
p
1
l
a
,
k
9,2
=
ω
1
-
d
h
,
1
t
1
E
p
1
0.5
d
h
,
1
1
k
2
,
1
k
5
,
1
k
8
Equation (33)
Remark. k1 and k9 are the stiffness of the plate from the bolt to the ends of plates, k3, k4, k6, and k7 are the stiffness of the plate between bolts, and k2, k5, and k8 are the bolt-plate interaction stiffness.Figure 20
The multibolt connection stiffness.
(a)
Multibolt connection spring model
(b)
Plate stiffnessFigure 21
Connection load and load distribution of analytical model and FEA.
(a)
Finite element model of triple bolt connection
(b)
Connection load
(c)
Load distribution


## 7. The Parametric Study of Multibolt Connection

The analytical spring model as in the previous section was used to examine the effects of various connection parameters. The parametric study calculated the load distribution of the bolt connection with variations in the bolt diameters, plate thickness (t1=t2), and the plate thickness ratio (t2/t1). The geometry model and the load condition as in the previous section were adapted in the parametric study. The efficiency (ψ) of the load distribution of multibolt connection was proposed shown in (40)ψ=nPmin∑i=1nPi,fori=1,2,3,…,n,where Pi is the bolt load in the i-th row, Pmin is the smallest load in the bolts, and n is the number of rows in the connection. Sensitivity analysis investigated these potential changes and their impacts on inferences to be obtained from the analytical model. The conception of sensitivity analysis is presented in Figure 22. The slope of the normalized efficiency graph with respect to changes in the geometric parameter was used to measure the sensitivity level. A steeper graph indicates a higher sensitivity level while a gentler slope indicates a lower sensitivity level. In each case study, the efficiency graph was normalized by the reference value. The reference cases involved bolt diameter of 8 mm, plate thickness of 1.0 mm, and t2/t1 of 3.0 which are used to analyze the sensitivity of variation in the bolt diameter, plate thickness, and t2/t1 sensitivity, respectively.Figure 22
Concept of the sensitivity analysis.Variation in the bolt diameter introduced efficiency of load distribution shown in Table6 and Figure 23(a). Increasing the bolt diameter had a significant effect on the load distributions in a multibolt connection. The change in the bolt diameter affected the proportion of the bolt-plate interaction stiffness between the bolt and plates which resulted in an even load distribution.Table 6
Bolt load distribution for variation in bolt diameter.
d (mm)
Bolt load distribution (%)
Efficiency,ψ (%)
Bolt 1
Bolt 2
Bolt 3
4
34.99
30.03
34.99
90.09
6
34.53
30.95
34.53
92.84
8
34.27
31.46
34.27
94.37
10
34.11
31.78
34.11
95.33
12
34.00
32.00
34.00
95.99
Remark. t1 and t2=1 mm and l1 and l4=100 mm.Figure 23
Bolt load distribution and sensitivity analysis results of varying bolt diameters.
(a)
Bolt load distribution
(b)
Parameter sensitivityThe efficiency of bolt load distribution had small changes with increased plate thickness as shown in Table7 and Figure 24(a). The bolt-plate interaction stiffness and axial stiffness of the plates were related to the plated thickness, as a change in the plate thickness produced an equal change in the bolt-plate interaction stiffness and the axial stiffness which slightly changed the efficiency of load distribution.Table 7
Bolt load distribution for variation in plate thickness.
t
1 and t2 (mm)
Bolt load distribution (%)
Efficiency,ψ (%)
Bolt 1
Bolt 2
Bolt 3
0.40
34.75
30.51
34.75
91.52
0.60
34.72
30.56
34.72
91.69
1.00
34.67
30.67
34.67
92.01
1.40
34.59
30.83
34.44
92.62
1.60
34.48
31.05
34.48
93.14
Remark. l1 and l4=100 mm, l2 and l3=15 mm, ω=30 mm, and d=5 mm.Figure 24
Bolt load distribution and sensitivity analysis results of varying plate thickness.
(a)
Bolt load distribution
(b)
Parameter sensitivityThe efficiency of load distributions was changed suddenly with changes to thet2/t1 ratio as shown in Table 8 and Figure 25(a).Table 8
Bolt load distribution for variation in plate thickness ratio (t2/t1).
t
2
/
t
1
Bolt load distribution (%)
Efficiency,ψ (%)
Bolt 1
Bolt 2
Bolt 3
1
34.72
30.56
34.72
91.69
2
38.14
30.65
31.21
91.95
3
39.67
30.75
29.59
88.76
4
40.39
30.86
28.75
86.25
5
40.69
30.98
28.33
84.98
Remark. t1=1 mm, l1 and l4=100 mm, l2 and l3=15 mm, and d=5 mm.Figure 25
Bolt load distribution and sensitivity analysis results of the plate thickness ratio (t2/t1).
(a)
Bolt load distribution
(b)
Parameter sensitivityThe efficiency of load distribution decreased when thet2/t1 ratio was greater than 1, but the carried load for the first bolt row increasingly changed. Once the t2/t1 ratio was greater than 1 (t2/t1≫1), the behavior of plate 2 was established as a rigid bar as shown in Figure 26. The steel bolt was constrained by the axial stiffness of plate 2 which was simulated as a semi-fixed support condition. Then, bolt 1 could not transfer the load before bearing failure or shearing of the bolt failure occurred.Figure 26
Effect of thet2/t1 ratio more than 1 (≫1) to the steel bolts.From the parametric study results, it was clear that varying the geometries influenced the connection stiffness, which significantly affected the load distribution. The sensitivity analysis results of the multibolt connection from variation in the bolt diameter, plate thickness, andt2/t1 are shown in Figures 23(b), 24(b), and 25(b). The sensitivity value was 0.0605, 0.0136, and −0.0646 for the bolt diameter, plate thickness, and t2/t1 sensitivity, respectively. The variation in t2/t1 showed high sensitivity with the triple bolt connection, with sudden decay when t2/t1>1. Thus, differences in the plate thickness should be emphasized in the connection design.


## 8. Conclusion

The bolt load distribution efficiency of a cold-formed steel, multibolt connection was presented. The analytical precedence with a spring model was used to simplify the stiffness of the bolt connection. A new bolt-plate interaction stiffness equation was proposed which was verified by the FEA and experimental tests. The spring model with bolt-plate interaction stiffness was used to estimate the bolt load distribution of the multibolt connection. The results showed that the analytical model was accurate and had good correlation with the FEA. Furthermore, the analytical procedure was used to examine the efficiency of load distribution for the multibolt connection and these results showed an uneven sharing of the load among bolts with one of them carrying the major component of the load. Thus, the stress distribution of the major load carrying hole resulted in high-stress concentration which caused weakness in the mechanically fastened connection. Moreover, increasingt2/t1 had a greater influence on the load distribution and efficiency of a multibolt connection than changing the bolt diameter and plate thickness. The results of the study indicate that the practical design procedures with assuming equal load distribution among the bolts may not be conservative because some components may be failed before the connection develops the full plastic strength. Finally, the proposed procedure can be used together with the conventional design of the bolt connection to get insight into the more realistic load distribution.



---

*Source: 1912724-2017-07-06.xml*