A General Solution for Eccentric Loads on

Weld Groups
G. DONALD BRANDT

Eccentric loads on weld groups traditionally have been 5. The weld element which will reach its ultimate defoi
analyzed by elastic methods in which the individual effects mation first is the one for which the ratio of A m a x di
of an axial load (applied through the centroid of the weld vided by the radius to the instantaneous center is th
group) and a pure moment (in the plane of the weld group) smallest. It is assumed that deformations vary linearl
were combined.1 In the 7th Edition Manual of Steel Con- with distance from the instantaneous center.
struction2 of the American Institute of Steel Construction,
6. Consistent deformations (A) at all other weld elemen
as in earlier editions, such a procedure was used to create
are then found from
tables which could be used to determine allowable eccentric
loads on selected weld groups.
These tables were criticized as having non-uniform
factors of safety when compared to the actual ultimate loads 7. The following parameters are then calculated for eac
which the welds could support. An "ultimate strength weld element:
method'' was proposed3 in which the resulting force per 10 + B;
unit of length of each weld element is calculated from R\i\i,i -
0.92 + 0.06030,-
R = Ruh(\ - e-^Y M z = 75^0.01 m

and the total resisting forces and moments are predicted by X,. = o.4e°- 014 <^
summing the elemental forces. In the 8th Edition Manual, 4
such a procedure was used to create tables similar to those Ri = *uit,,-(l - e-^)^
in earlier editions, but with different numerical coeffi- 8. By statics, calculate the corresponding applied force ai
cients. moment which hold the weld forces in equilibrium.
Butler, Pal, and Kulak 3 describe the ultimate strength
method in some detail. For the sake of completeness in this
report, this procedure is summarized as follows:
1. Choose an instantaneous center of rotation. (Fig. 1)
2. Assume that the resisting force on any weld element acts
perpendicularly to a radius connecting that element to
the instantaneous center.
3. Calculate the angle, 0, between the elemental force and
the axis of the weld element. (Angle 6 is expressed in I.C.
degrees.)
4. Determine the ultimate deformation which can occur
on each weld element from
A m a x = 0.225(0 + 5)-°-<7

G. Donald Brandt is Professor and Chairman, Dept. of Civil

Engineering, The City College, City University of New
York.
Figure 1

150

ENGINEERING JOURNAL / AMERICAN INSTITUTE OF STEEL CONSTRUCTION

 It should be noted that the formulas given above corre- Engineers familiar with FORTRAN should have no
spond to a V4-in.weld made using E60 electrodes. The R{ particular difficulty in understanding the program after
values are pounds per linear inch of weld. The final forces reading the following discussion.
and moments are ultimate values.
Tide 5 described the way in which these values were DISCUSSION OF THE COMPUTER PROGRAM
converted to allowable load coefficients for the AISC tables.
In summary, this requires: Computations for the Elastic Solution—The length of
each weld element, W, is calculated, and using these dis-
1. Converting from V4-in. weld to Vi6-in. weld by dividing crete elements, the centroid is located, the polar moment
by 4.
of inertia, / , is calculated, and the moment, M 0 , of the
2. Converting from E60 electrode to E70 electrode by applied unit force about the centroid is calculated.
multiplying by 70/60. A mapping function called FACTOR is used to trans-
form forces to distances and it is equal to / divided by the
3. Introducing a safety factor conforming to AISC Speci- product of the total length of weld and moment M0.
fication Sect. 1.5.3 by multiplying by 0.3. The instantaneous center is located by adding to the X
4. Checking the shear stress on the most highly stressed and Y coordinates of the centroid the quantities — Py X
weld element and, if it exceeds 21.0 ksi, reducing the FACTOR and +PX X FACTOR, respectively.
load by the ratio of 21.0 divided by that shear stress. The radius vector, D to the center of each weld element
is calculated, and the largest one noted. The allowable
The papers by Tide 5 and by Butler et al 3 indicate that moment about the instantaneous center of all elemental
solutions have been obtained using computer programs. forces is XWD2/Dm3LX times the allowable force per inch
Butler states that the method is general, but shows details of weld. The allowable (elastic) load is that moment divided
for only a C-shaped weld subjected to loading parallel to by the moment of the applied unit force about the instan-
a principal axis. Tide shows results for a pair of parallel taneous center. This load and the identification number of
line welds and a C-shaped weld with the loads again par- the critical weld element are printed.
allel to principal axes. The AISC tables include a pair of
parallel line welds, rectangular box welds, C-shaped welds, Computations for the First Approximate Ultimate
and L-shaped welds; loads are parallel to a principal axis Load—The instantaneous center located above is used as
in the first three types and parallel to a leg of the L in the the first trial center. Calculations follow the procedure
last. attributed to Butler3 described earlier in this paper. (Angle
In a previous paper, 6 the author showed how rapid so- 0 is determined as follows: Form the dot product of the
lutions could be obtained for any eccentrically loaded bolt radius vector and the weld element vector. Divide this by
groups. The same method can be extended to weld groups. the product of the magnitude of the two vectors; this quo-
In essence, the method involves: tient is the cosine of the angle between the vectors. Its
complement is angle 6.)
1. Directly finding the instantaneous center corresponding The (first approximate) ultimate load is the moment
to elastic behavior of any weld arrangement for any about the instantaneous center of all the elemental forces
eccentric load. (Hi times Wt times D{) divided by the moment about the
instantaneous center of the applied unit load. This ultimate
2. Directly determining the elastic solution for the maxi-
load, and the identification number of the critical weld
mum permissible load.
(having the largest force per inch) are printed, along with
3. Directly determining an approximate value for the the coordinates of the current instantaneous center.
ultimate load.
Iterating to Improve the Approximate Ultimate
4. Iterating to improve the approximate value. Load—When the instantaneous center is correctly located,
not only will the sum of the moments about that center be
5. Using the same procedure described by Tide to convert
zero, but the vector sum of the forces will be zero as
the ultimate load to an allowable load consistent with
well.
the AISC tables.
For the ultimate strength case, the X component of the
Inasmuch as welds are continuous, it is necessary to force on each weld element due to the unit load can be found
discretize into a finite number of weld elements. A mod- from
erately large number of discrete elements is required if
reasonable accuracy is to be achieved, so the rather long
procedure for each weld element described above probably A Puit
makes the procedure too laborious for manual calculations. and similarly for the Y component.
Accordingly, the procedure is given here as a FORTRAN The sum of all the Rx values plus Px will not be zero
computer program. (See Appendix.) unless the instantaneous center is correct. The unbalanced

151
THIRD QUARTER / 1982
force is Fxx. Similarly for the Y direction. The vector sum
of Fxx and Fyy is F> the unbalanced force. Fxx , Fyy and F
are printed when the approximate ultimate load is
printed.
It has been found that components of a desirable shift 30"

from the previous trial center can usually be predicted from

the same formulas used to locate the elastic instantaneous
center. Thus, if x\ and y\ represent the coordinates of the
instantaneous center for which the unbalanced force
components are Fx and Fyy , the coordinates of the next
center should be:
+
x2 = xx - (Fyy X FACTOR)
y2 = yx + (FxxX FACTOR)
New values for the radius vectors to all weld elements
i"
are now calculated and the ultimate strength solution is
repeated. In most cases, the unbalanced force decreases
rapidly, and the ultimate load stabilizes after a few itera- Figure 2
tions. In the program, when the unbalanced force is less
than 1% of the applied unit force, the solution is declared signal that all data has been given. (Format for all three
found. types of data cards is A2, 8X, 5F10.3. A fourth form of data
The ultimate load is converted to permissible load on a is discussed later.)
Vi6-in. weld made using E70 electrodes in the manner at- The program converts weld descriptions of either type
tributed to Tide 5 described earlier. This value is into weld elements. These are numbered consecutively and
printed. printed, giving their identification numbers, center coor-
It should be observed in using the program that, as in the dinates, and X and Y projected lengths. The LOAD data
case of AISC tables, the permissible load on the weld group is also printed.
is obtained by multiplying the permissible load given by the
program by the number of sixteenths in the weld size. Using the Program—In all the examples, the maximum
Furthermore, if the electrodes used are other than E70 permitted elastic load, Pe, the ultimate load, Pu, and the
(having Fv = 21 ksi), the permissible load should be mul- maximum permitted load, Pp, derived from the ultimate
tiplied by the allowable shear divided by 21 ksi. Unlike the load, are listed.
AISC tables, no multiplication by any length should Example 1 (see Fig. 2):
made.
Data cards were prepared using the LINE and LOAD
Input—In creating the program, provision was initially descriptions. Actual computer output is shown in Fig. 3 for
made for describing the weld configuration in two different a case where each line is subdivided into 10 equal ele-
ways. (Both of these have been retained.) ments.
1. If the first columns on a data card contain the word A reasonable question arises in every problem regarding
LINE, the description of a weld consists of the coordi- the number of subdivisions required to get good accuracy.
nates of its starting point and ending point, and the For the weld geometry and load of Example 1, computer
number of equal segments into which the weld is to be runs were made dividing each line into 2, 5,10, 20, and 40
subdivided. segments. The results are shown in Table 1.
2. If the first columns on the data card contain the word
Table 1
ELEMENT, the description of a weld element consists
of the coordinates of its center and the length of its X and Subdivisions 2 5 10 20 40
Y projections on the coordinate axes. 1.539 1.253 1.151 1.100 1.075
Pe
If the first columns on the data card contain the word Pu 26.304 24.783 24.741 24.593 24.499
LOAD, the description of a normalized load (magnitude
Pp 1.565 1.472 1.469 1.460 1.455
= 1) consists of its X and Y projections on the coordinate
axes, and the coordinates of a point on its line of action.
Observe that there may be as many data cards of LINE Values obtained from the 7th Edition AISC Manual (for
and/or ELEMENT type as required to describe the entire the elastic solution) and from the 8th Edition AISC Manual
weld. Such cards may be in any order or any mixture. The (for the ultimate strength solution) are, respectively, Pe =
last data card must be a LOAD card, since this is also a 1.05 and P p = 1.45.

152
ENGINEERING JOURNAL / AMERICAN INSTITUTE OF STEEL CONSTRUCTION
DATA FOR DESCRIPTION OF WELDS AND LOAD

LINE C,.COOO C..COOO 0 . 0000 10.0000 10.0

LINE 1. . 0 0 0 0 C,. 0 0 0 0 1 . COOO 10,0000 10.0
LOAD 0,. 0 0 0 0 - 1 .. 0 0 0 0 30. 5000 0.0000

NUMBER OF WELD ELEMENTSi = 20

COORDINATES AND PROJECTICNS
1 0.0000 0.5000 0.0000 1.0000
2 0.0000 1.5000 0.0000 l.COOO
3 0.0000 2.5000 0.0000 1.0000
4 0.0000 3.5CC0 0.0000 1.0000
5 0,0000 4.5000 0.0000 1.0000
6 oJoooo 5.5000 0.0000 1.0000
7 0.0000 6.5000 0.0000 1.0000
8 0. cooo 7.5000 0.0000 1.0000
9 o.coco 8.5000 0.0000 1.0000
10 0.0000 9. 5000 0.0000 I.0000
11 l.COOO 0.5000 0.0000 1.0000
12 1.COOO 1.5000 0.0000 1.0000
13 1.0000 2.5000 0.0000 I.0000
14 1 . cooo 3.5000 0.0000 1.0000
15 1 . cooo 4.5000 0.0000 1.0000
16 1.000 0 5 . 5C00 0.0000 1.0000
17 1.0000 6.5000 0.0000 1.0000
18 1.0000 7.50G0 0.0000 1.0000
19 1 . cooo 8.5000 0.0000 1.0000
20 1 . cooo 9*5000 0.0000 l.CCOO

PX= 0.CGOO,PY= -1.0000, POLX 30.5000,POLY = 0.0000

ELASTIC VALUE FOR MAXIMUM LOAD IS 1.151

MULTIPLY T H I S BY NUMBER OF SIXTEENTHS
AND BY ALLOWABLE KSI / 2 l . , 0
C R I T I C A L ELEMENT I S NUMBERI 11

AT T R I A L CENTER NO. 1
X= 0.2139 Y= 5.0000
FX = -0.0000, FY = 0.3313
X
F = 0.3313
APPROXIMATE ULTIMATE LOAD = 24.748
CRITICAL ELEMENT IS NUMBER I
ON WHICH ULTIMATE FORCE PER INCH IS 15.6361
Figure 3. Computer printout for Example 1

153
THIRD QUARTER / 1982
 AT TRIAL CENTER NO. 2
X= 0.30 8 7 Y* 5.0000
FX = -0.0000 FY « -0.1052
F = 0.1052
APPROXIMATE ULTIMATE LOAO * .24.743
C R I T I C A L ELEMENT I S NUMBER I
ON WHICH ULTIMATE FORCE PER INCH I S 15.6238

AT TRIAL CENTER NO. 3

X= 0.2786 Y= 5.0000
FX = -0.0000 FY = 0.0344
F = 0.0344
APPROXIMATE ULTIMATE LOAD * 24.740
C R I T I C A L ELEMENT I S NUMBER 1
ON WHICH ULTIMATE FORCE PER INCH I S 15.6277

AT TRIAL CENTER NO. 4

X= 0.2884 Y= 5.0000
• FX - -0.0000 FY = -0.0112
F = 0.0112
APPROXIMATE ULTIMATE LOAO = 24.741
C R I T I C A L ELEMENT I S NUMBER 10
ON WHICH ULTIMATE FORCE PER INCH I S 15.6264

AT T R I A L CENTER NO. 5
X= 0 . 2 8 52 Y= 5 . 0 000
FX * -0.0000 FY * 0.0036
F = 0.003 6
APPROXIMATE ULTIMATE LCAO = 24.741
C R I T I C A L ELEMENT I S NUMBER 1
ON WHICH ULTIMATE FORCE PER INCH IS 15.6269

INSTANTANEOUS CENTER HAS BEEN LOCATED

MAXIMUM PERMISSIBLE LOAD I S 1.469
MULTIPLY T H I S BY NUMBER OF SIXTEENTHS
AND BY ALLOWABLE K S I / 2 1 . 0
Fig. 3 (cont'd)

154

ENGINEERING JOURNAL / AMERICAN INSTITUTE OF STEEL CONSTRUCTION

Example 2:
Use the same weld configuration and load as in Example 1
1. Intuition indicates that more subdivisions should give L 3" 4" 10"
better results, but there are two effects. One is simply the
closer approximation to continuity achieved by using I

smaller elements. The other is the fact that the maximum CO 4

stress occurs on the outside edge of a weld element rather
than at its center. Obviously, more subdivisions of uniform
size means that the .edge of an element and its center will
be closer together. It is, however, possible to divorce the two vb X 4
effects by using non-uniform weld elements. For example, 3 X
the top and bottom 0.001 in. of each of the welds in Ex-
ample 2 are entered as individual elements, and the re- <N
mainder subdivided as in Example 1. The new results are —• x
shown in Table 2.
3" ] 4"
Table 2
Subdivisions 2 5 10 20 40
Figure 4
Pe 0.797 1.009 1.039 1.047 1.049

23.867 23.977 24.307 24.374 24.388

AUTOMATIC in place of LINE. When this was done for
Pu
both line welds, the results were as shown in Table 3.
Pp 1.417 1.423 1.443 1.447 1.448
Table 3
Subdivisions 2 5 10 20 40
It can be observed that the solutions stabilize with fewer
subdivisions, and that the approach to the true solution is Pe 0.800 1.010 1.040 1.047 1.049
from the safe side rather than the unconservative side. The Pu 23.919 23.997 24.317 24.378 24.391
same behavior has been observed in a dozen other test cases
Pp 1.420 1.424 1.444 1.447 1.448
for welds of different configurations.

Example 5 {see Fig. 4):

Example 3:
Welds need not be vertical or horizontal, or even symmet-
Use the same weld configuration and load as in Example
rical. The applied load can be applied in any direction.
1. It should obviously make very little difference if the small
Consider the bracket shown.
element welds are introduced in addition to (or even over-
Making use of the AUTOMATIC feature, and selecting
lapping) the original welds. Therefore, if the critical weld
the number of subdivisions for each line so that "real" el-
element was unknown before the first analysis, only two
ements were about 1-in. long, the data looked like this:
additional cards need be prepared and added to the previous
data for Example 1 to achieve the precision of Example 2.
AUTOMATIC 7.0 11.0 3.0 11.0 4.0
In Example 3, the data is:
AUTOMATIC 3.0 11.0 0.0 8.0 4.0
AUTOMATIC 0.0 8.0 0.0 2.0 6.0
ELEMENT 0. 10. 0. 0.001 AUTOMATIC 0.0 2.0 3.0 0.0 4.0
ELEMENT 1. 10. 0. 0.001 AUTOMATIC 3.0 0.0 7.0 0.0 4.0
LINE 0. 0. 0. 10. 10 LOAD 0.6 -0.8 17.0 9.0
LINE 1. 0. 1. 10. 10
LOAD 0. -1. 30.5 0. From the program, the value of Pe is 3.869; since the
weld size is 5/i6, the allowable elastic load is 5 X 3.869 =
The results are Pe = 1.040, Pu = 24.312, Pp = 1.444. 19.345 kips. The value of Pp is 5.452; the permissible load
is 5 X 5.452 = 27.260 kips.
Example 4: REFERENCES
Use the same weld configuration and load as in Example 1. Salmon, C. G. and J. E. Johnson Steel Structures—Design
1. In order to avoid the extra effort involved in specifying and Behavior Second Edition, Harper and Row, 1980, Chap.
data for additional elements, an automatic insertion of an 5.
extra "fake" element 1/1000th as long as the subdivision 2. Manual of Steel Construction Seventh Edition, American
size can be made at each end of any line by using the word Institute of Steel Construction, 1970.

155
THIRD QUARTER / 1982
3. Butler, L. J., S. Pal, and G L. Kulak. Eccentrically Loaded Design Tables Engineering Journal, American Institute of
Welded Connections Journal of the Structural Division, Steel Construction, Vol. 17, No. 4, 4th Quarter 1980.
ASCE, Vol. 98, ST5, May 1972, pp. 989-1005. 6. Brandt, G. Donald Rapid Determination of Ultimate
4. Manual of Steel Construction Eighth Edition, American Strength of Eccentrically Loaded Bolt Groups Engineering
Institute of Steel Construction, 1980. Journal, American Institute of Steel Construction, Vol. 19, No.
5. Tide, R. H. R. Eccentrically Loaded Weld Groups—AISC 2, 2nd Quarter 1982.

APPENDIX
$ JOB L I S T , K P = 2 9 , T IME=6C0,PAGES=90
$IBFT
C PROGRAM TO LOCATE I N SiTANTANEQUS ' CENTERS OF WELD GROUPS
C AND OBTAIN ELASTIC Af <HQ ULTIMATE RESISTANCES
1 REAL J , M, M P , MPCG, I X , I Y , MU, LAMBDA
2 INTEGER EL EM, TYPE
3 INTEGER AUTO
4 DIMENSION TEMP{5)
5 DIMENSION X C 1 0 0 ) , Y{ 1 0 0 ) , D X ( I O O ) , D Y < 1 0 0 ) r D ( 1 0 0 > ,
1 T H E T A ( I O O ) , R ( 1 0 0>1) , W ( 1 0 0 ) , W X U O O ) * WY( 1 0 0 )
6 FALLOT = 0 * 9 2 8 1
7 200 CONTINUE
GENERALIZED INPUT
8 WRITE i 6 , 6 0 0 0 )
9 6 0 0 0 FORMAT i 1 H 1 , f DATA FOR DESCRIPTION OF WELDS AND LOAD » , / I
10 I = 0
li 300 CONTINUE
12 READ ( 5 , 5CC ) •Et
TYPI ( TEMPCJJ),JJ= 1,5 )
13 500 FORMAT { A 2 , 8X, 5F 1 0 . 3 )
1* DATA LINE / • L I 1 / ELEM / «EL« / , LOAD / f L C /
15 CATA AUTO / ' A U ' /
16 I F ( TYPE . E C . LINE GO TO 3 0 1
17 I F C TYPE . E Q . AUTO GO TO 3 0 1
18 I F ( TYPE - E Q . ELEM GO TO 302
19 I F ( TYPE • EQ. LOAD GO TO 303
20 WRITE ( 6 , 6 004 ) YPE, ( T E M P { J J > , J J - 1 , 5 )
21 6004 FORMAT M H O , 'UNABLE 0 DECODE THIS CARD ', /,
10X, A2 8 X , 5 F 1 0 . 3 , / , • TYPE NOT STANDARD ' , / , / , / )
22 GO TO 300
23 301 CCNTINUE
24 WRITEI6. 6001 ) i TEMPIJJ),JJ =U 5 )
25 6001 FORMAT { 1H , f LINE«, 6 X , 4 F 1 0 . 4 , F 8 . 1 I.
26 IF ( TYPE .EG. LINE ) GO TO 306
27 WRITE ( 6, 6005 )
28 6005 FORMAT i 1H , f SHORT WE LDS ADDED AUTOMATICALLY TO LINE ABOVE1,/,/)
29 306 CONTINUE
CCNVERT TEMP INTO L I N E DATA
30 SX = T E M P ( l )
31 SY * TEMPC2)
32 EX = T E M P I 3 )
33 EY « T E M P I 4 )
34 NSEG = TEMPJ5)
XP1, YP1 = PROJECTED LENGTH OF 1 ELEMENT
35 XP1 = ( EX - SX ) / N SEG
36 YP1 * i EY - SY ) / N SEG
XBACK, YBACK -= COORD I NATE OF P O I N T H A L F S T E P BEYOND STAfU
37 XBACK = SX - X P 1 / 2 .
38 YBACK = SY - Y P 1 / 2 .
39 IF (TYPE .EQ. LINEI GO T O 3 0 4
40 FAKEX XP1 / 1000.
41 FAKEY YP1 / 1000.
42 I * I 1
43 WX( I ) FAKEX
44 WYU) FAKEY
45 X ( I ) = SX
46 Y U i = SY

156
ENGINEERING JOURNAL / AMERICAN INSTITUTE OF STEEL CONSTRUCTION
47 304 CCNTINUE
48 00 1001 . K * I t NSEG
49 1 = 1 + 1
50 wx( n = X P I
51 W Y d ) •= YP1
52 X ( I ) = XBACK • K * X P I
53 Y d ) « YBACK • K * YP1
54 1001 CCNTiNUE
55 I F ( TYPE . E Q . LINE I GO TU 305
56 1 * 1 * 1
57 W X ( I ) = FAKEX
58 W Y I U = FAKEY
59 x m * EX
60 Y d ) = EY
61 305 CONTINUE
62 GO TO 3 0 0
63 302 CCNTINUE
64 WRITE ( 6 t, 6002
6002 )) ( TTEEM
MPP((JJJJ)) ,, JJJ
J = 1I , 4 )
65 6 0 0 2 FORMAT ( iH i •ELEMENT* t 3 X , 4 F 1 0 „ 4 )
C CONVERT TEMP TO ELEMENT DATA
66 1 * 1 * 1
67 X U ) = TEMPdi
68 Y d ) = TEMPI2)
69 WX( I ) = T E M P O )
70 W Y ( I ) = TEMP(4)
71 GO TO 3 0 0
72 303 CONTINUE
73 WRITE ( 6 , 6003 ) ( TEMP(JJ) , J J = U
74 6 0 0 3 FCRMAT ( 1H , * L0A0 f , 6X , 4F10.4 )
I CCNVERT TEMP TO LOAD DATA
75 PX * TEMP(1)
76 PY « TEMP(2)
77 PGLX = TEMP(3)
78 POLY = T E M P U )
79 N = I
80 WRITE ( 6 , 6 0 2 ) N , d , X { I ) , YC I ) , W X d ) , WY I ! ) , I=IiN )
81 602 FCRMAT ( 1 H 1 , • NUMBER OF WELO ELEMENTS = • , 1 6 , / ,
f
1 • COORDINATES AND PROJECTIONS , /,
2 (IX, 14, 4F10*4) )
82 WRITE ( 6 , 6 0 3 ) PX, P Y , PGLX, POLY
83 603 FORMAT I 1 H 0 , # P X = f , F l 0 . 4 , • , PY= • , F 1 0 . 4 ,
f
1 , POLX = ' , F 1 0 , 4 , V,PGLY = f , F 1 0 . 4 )
C CALCULATE LENGTH OF EACH WELD ELEMENT
84 DO 101 I * It N
85 W(U = SORT < VtXi I) * • 2 + WYU) ** 2 )
86 101 CCNTINUE
C LOCATE CENTER OF GRAVITY OF WELD GROUP
87 SUMX » , 0 .
88 SUMY = 0 .
89 SUMW = 0 .
90 DC 1 0 2 1 = 1,N
91 SUMX = SUMX + X ( I ) * W d I
92 SUMY * SUMY «- Y ( I ) * W ( I )
93 SUMW - SUMW> W( I )
94 102 CONTINUE
95 XCG = SUMX / SUMW
96 YCG = SUMY / SUMW
C CALCULATE MOMENT OF P ABOUT CG
97 MPCG = PY * (POLX - XCG) - PX * (POLY - YCG)
C CALCULATE* J
98 J = 0.
99 DO 1 0 3 I = 1,N
100 DX(I) * X( I) - XCG
101 DY( I ) * Yd) - YCG

157
THIRD QUARTER / 1982
 102 IX = W U ) / 12. * W Y H ) ** 2 «• W U J * 9Y < I ) ** 2
103 IY = W( i) / 12* * W X U ) ** 2 * W(I) * 0X(I) ** 2
104 J = J + IX + IY
105 103 CONTINUE
106 FACTOR = J / C SUMW * MPCG )
107 XIC = - P Y * FACTOR + XCG
108 Y I C = *PX * FACTOR «• YCG
109 ITER = 0
110 FPREV * i.E+12
111 201 CONTINUE
RECALCULATE DX, DY, D ANO DELTA FROM TRIAL CENTER
112 DO 104 1=1,N
113 DX<I) = X(I) - XIC
114 DY<Ii « YCI) - YIC
C CALCULATE LENGTH OF RADIUS VECTOR
115 D(I) = SQRT { DX(I)**2 * D Y U J * * 2 I
116 104 CCNTINUE
C CALCULATE MCMENT OF P
117 MP = PY * ( P G L X - X I C ) - P X * ( P O L Y - Y I C )
118 I F ( ITER . G E . 1 I GO TO 2 0 4
C GET THE ELASTIC SCLJTION
119 BIGD = 0 *
120 SUMWD2 = 0 .
121 DC 107 I * 1 . N
122 lfiQ(l) .LE. BIGD ) GO TO 2 0 3
123 BIGD = Oil)
124 IBIGD » I
125 203 CONTINUE
126 SUMWD2 * SUMWD2 f M ( I ) * O U ) * * 2
127 107 CCNTINUE
128 ELPMAX^ABS i SUMWD2 * FALLOW / ( MP * BIGD ) I
129 WRITE ( 6 , 6 0 4 1 ELPMAX, I B I G D
130 604 FORMAT i 1 H 0 , »ELASTIC VALUE FOR MAXIMUM LOAD I S % F 1 0 . 3 t
1 10X» • MULTIPLY T H I S BY NUMBER OF SIXTEENTHS • , / ,
2 1 0 X , f AND BY ALLOWABLE KSI / 2 1 . 0 ' , / »
3 1H , ' C R I T I C A L ELEMENT I S NUMBER f , 14, / )
131 204 CONTINUE
132 DLDMIN ^ l . E + 1 2
133 DO 1 0 5 I = 1,N
C CALCULATE ANGLE BETWEEN RADIUS VECTOR AND WELD AXIS
134 PHl = i A R C G S U B S U W X U ) * D X C I ) • W Y U ) * D Y U ) ) /
1 4W(I) * D M ) ) ) ) ) * 180. / 3.1415927
C CALCULATE ANGLE BETWEEN FORCE AND WELC A X I S
135 T H E T A ( I ) = 9 0 . - PHI
£ CALCULATE PRELIMINARY DELTA VALUES
136 DELTA = 0 . 2 2 5 * ( THETACI) • 5 . 0 ) * * ( - 0 . 4 7 )
137 DELD « DELTA / D i l i
C FIND SMALLEST DELTA/D
138 I F { DELD .GE. DLDMIN ) GO TO 105
139 DLDMIN = OELD
140 DMAX « D ( I )
141 BIGDEL = DELTA
142 105 CONTINUE
143 BIGR - 0.
144 SUMM = 0 .
C CALCULATE REVISED D E L T A , MU, LAMBDA, RULT, R, M
145 DO 1 0 6 I * 1, N
146 DELTA = BIGDEL * D ( I ) / CMAX
147 MU = 7 5 . * EXP ( 0 . 0 1 1 4 * T H E T A ( I ) )
148 LAMBDA » 0 . 4 * EXP C 0 . 0 1 4 6 * T H E T A U ) )

158

ENGINEERING JOURNAL / AMERICAN INSTITUTE OF STEEL CONSTRUCTION

149 RILT * I 1 0 . + THETAC I I ) / I 0 . 9 2 * 0 , 0 6 0 3 * T H E T A ( I ) )
150 Ril) * RULT * I 1 . - EXP i -MU * DELTA ) ) * * LAMBDA
C CALCULATE LARGEST R
151 IF ( R i l l . L E . BIGR J GO TO 2 0 2
i52 BIGR = Rill
153 JBIG = I
154 2 02 CONTINUE
155 M * R ( I ) * D ( I ) * W<It
156 SUMM=SUMM + M
157 106 CCNTINUE
C CALCULATE RUNIT DUE TC U N I T LOADt CALCULATE PULT
158 RUNIT * - M P / SUMM
159 PULT = ABS ( SUMM / MP )
Z CALCULATE RX/P AND R Y / P
160 SUMRX * 0 ,
161 SUMRY « C.
162 DO 1 0 8 I*lfN
163 RX = - D Y I I J / D I D * R ( I ) * RUNIT * M( I I
164 SUMRX « SUMRX + RX
165 RY = 0X< I-> / 0 ( 1 ) * R U I * PUNIT * U ( I )
166 SUMRY = SUMRY • RY
167 108 CCNTINUE
168 FX = PX + SUMRX
169 FY = PY + SUPRY
170 F = SQRTi ?X**2 • FY**2 )
171 ITER = ITER * 1
172 WRITE ( 6 , 6 0 1 1 I T E R , X I C , Y I C , F X , F Y , F , PULT, J B I G , BIGR
173 601 FORMAT I 1 H 0 , • AT TRIAL CENTER NO. «, 14, / ,
1 i O X t «X* f , F 1 0 . 4 , • Y= • , F10.4, /,
2 2 0 X , *FX * f , F 1 0 . 4 , f FY - . * , F 1 3 . 4 , / ,
3 3 OX, «F * f , F10.4, /
4 2 U X , 'APPROXIMATE ULTIMATE LOAD = • , F 1 0 . 3 , /,
5 1 0 X , ' C R I T I C A L ELEMENT I S NUMBER • , 14, / ,
6 1 0 X , «CN WHICH ULTIMATE FORCE PER INCH IS • » F 1 0 . 4 , / !
174 I F C ABSiFI *LE. 0 . 0 1 } GO TO 2 06
175 I F { ITER . G E . 30 ) GO TC 2C8
C TRY NEW CENTER
176 I F i F . L T . FPREV ) GO TO 2 0 5
177 FACTOR = FACTOR / 2 .
178 205 CONTINUE
179 FPREV = F
180 XIC « X I C - FY * FACTOR
181 Y I C « Y I C * FX * FACTOR
182 GC TO 2 0 *
183 206 CONTINUE
C CCNVERT PULT TO 1 / 1 6 WELD OF E70 ELECTRODE
184 PMAX = PULT / 4 . * 7 0 . / 6 0 .
c INTRODUCE SAFETY FACTOR
185 PMAX a o . 3 * PMAX
186 R L I M = FALLOW / 0 . 3 * 4 . 0 * 6 0 . / 70
187 I F i BIGR . L E . RLIM ) GO TO 207
188 PMAX = PMAX * RLIM / BIGR
189 207 CONTINUE
190 WRITE ( 6 , 6 0 5 ) PMAX
191 60 5 FCRMAT ( 1 H 0 , •INSTANTANEOUS CENTER HAS BEEN LOCATED • , /,
1 1 0 X , MAXIMUM PERMISSIBLE LOAD* IS • , F 1 0 . 3 , / ,
2 1 0 X , •MULTIPLY T H I S BY NUMBER OF SIXTEENTHS f , / ,
3 1 0 X , f AND BY ALLOWABLE KSI / 2 1 . 0 • , / )
192 GO TO 2 0 0
193 208 WRITE (6„606> ITER
194 606 FCRMAT { 1 H 0 , •ITERATIONS = •# 15,f...TRIALS TERMINATED•• I
195 GO T3 2 0 0
196 END

$DATA

159

THIRD QUARTER / 1982