UILU-ENG-74-2029

I ENGINEERING STUDIES
RAt RESEARCH SERIES N.O. 413

E
A C PUlER PR GRA
F RINELASlIC RESP NSE
/ C FRA ES T EARTH UAKES

By
SHUNSUKE OTANI

A Report on a Research
Project Sponsored by
THE NATIONAL SCIENCE FOUNDATION
Research Applied to National Needs Program,
Research Grant GI 29934

UNIVERSITY OF ILLINOIS
at URBANA-CHAMPAIGN
URBANA, ILLINOIS
November, 1974
 SAKE

A COMPUTER PROGRAM
FOR INELASTIC ANALYSIS
OF Ric FRAMES TO EARTHQUAKE

by

Shunsuke Otani

A Report on a Research Project Sponsored by

THE NATIONAL SCIENCE FOUNDATION,

Research Applied to National Needs Program,
Research Grant GI-29934

University of Illinois
at Urbana-Champaign
Urbana, Illinois
November, 1974
 ii

ABSTRACT

A special purpose computer program was developed to analyze an

inelastic reinforced concrete plane frame structure subjected to an

intense earthquake motion in one horizontal direction. A reinforced

concrete structure may consist of more than one regular rectangular

unbraced frame of arbitrary numbers of stories and bays.

The function, usage and background of the program are described

in this report with an example solution.

 iii

ACKNOWLEDGMENT

This study was part of a research project, entitled "The Evaluation

of Structural Damage Caused by the San Fernando Earthquake (1971) toward

the Development of Improved Methods for Earthquake-Resistant Design

(principal investigator Professor Mete A. Sozen)." The project was

sponsored by the Research Applied to National Needs program of the

National Science Foundation under Grant GI-29934. The study was

carried out at the Department of Civil Engineering, University of

Illinois at Urbana-Champaign.

The writer wishes to thank Dr. Mete A. Sozen, Professor of Civil

Engineering at the University of Illinois at Urbana-Champaign, for his

constant encouragement, generous support and critical comments. With-

out his encouragement, this report could not have been made in this

form.

The IBM 360/75 computer system of the Office of Computing Service

at the University of Illinois at Urbana-Champaign was used during

programming. Computer subroutine programs (DSINV and DMFSD) in the

IBM Scientific Subroutine Package were used in the program for

inverting a symmetric matrix. A computer subroutine program (SYBAN),

developed by Dr. Robert E. Miller, Professor of Theoretical and Applied

Mechanics, was used to solve a system of linear equations.

 iv

TABLE OF CONTENTS

Page
Chapter

1. INTRODUCTION. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 1
1.1 Obj ect . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 1
1.2 Review of Previous Work ......... 0 2

2. ASSUMPTIONS OF ANALYSIS. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2.1 Assumptions about Structures ...................... 6
2.2 Assumptions about Member and Joint .................. 7
2.3 Assumptions in Numerical Integration ........... 8

3. STRUCTURAL CONFIGURATION .................... 0 10

3.1 Regular Frame. 0 0 0 10
3.2 Irregular Frame ............................. 11
3.3 Treatment of Hinged and Free Joints ........ 12

4. MEMBER CHARACTERISTICS .............................. 13

4.1 Elastic Line Element.................................. 13
4.2 Inelastic Flexural Rotational Springs ................ 14
4.3 Rigid Zones ................................ "........ 16
4.4 Inelastic Joint Rotational Spring ............. 16
5. ANALYSIS METHOD...... . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18

6. DATA PREPARATION.. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 21
6.1 Structural Cards...................................... 22
6.2 Member Property Cards. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 27
6.3 Component Frame Cards ............................. 32
6.4 Earthquake Cards...................................... 37

7. DESCRIPTION OF COMPUTER OUTPUT ......................... 44

7.1 Print Output.......................................... 44
7.2 Output on Magnetic Tape .............................. 46

LIST OF REFERENCES..... . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 53

FIGURES. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 58

APPENDIX

A. EVALUATION OF MEMBER PROPERTIES ........................ 70

A.l Flexural Properties ................................... 70
A.2 Properties of Flexural Elements ....................... 71
A.3 Properties of Joint Springs ................ 73
 v

APPENDIX

B. HYSTERESIS RULES FOR FLEXURAL AND

JOINT SPRINGS............................................. 76
B.l Takeda's Hysteresis System............................. 78
B.2 Simplified Takeda Hysteresis System ................... 83

C }.1El.1BER STIFFNESS MATRIX.... . 90

D. STRUCTURAL STIFFNESS MATRIX . 97

E. MASS AND DAMPING MATRICES. 104

F GRAVITY EFFECT .. v. . . 106

G. NUMERICAL INTEGRATION OF EQUATION

OF MOTION . ., .. ., ....... 0 ..... 0 " " a II 0 iii e 110

H. MEMBER RESPONSE AND UNBALANCED FORCES 112

I. LISTING OF FORTRAN IV PROGRAM . 120

J EXAMPLE PROBLEM.............................................. 137

 1. INTRODUCTION

1.1 Object

A special purpose computer program was developed in order to under-

stand the realistic, inelastic behavior of a multistory reinforced

concrete frame structure to an intense earthquake motion in one horizontal

direction.

A structure on a rigid foundation may consist of more than one

regular rectangular unbraced plane frame of arbitrary numbers of bays

and stories.

An equivalent spring model was used in the program to simulate

inelastic flexural deformation of a member. The analytical model

recognizes stiffness changes caused by cracking of the concrete, yielding

of the reinforcement, and stress reversals. Hysteresis rules for the

flexural behavior of a reinforced concrete member under load reversals

were adopted after Takeda et al (1970).

In addition to inelastic flexural deformation, a deformation caused

by bond slip of the longitudinal reinforcement within a joint core was

considered.

The behavior of a structure depends on the behavior of individual

structural elements under a loading condition. In order for the

structure to behave sound and ductile, each member must be designed

so that premature failure or excessive stiffness deterioration,

caused by shear and axial loads in the member or by bond between the

concrete and reinforcement, must be properly prevented up to a

1
 2

probable maximum deformation level during an intense earthquake motion.

I

This program provides the maximum deformation that might occur in

a frame member during a given earthquake motion, within limitations

defined by the assumptions of the program.

This report was prepared as a user's guide to the program (code

name: SAKE). A program user is assumed to be liable for the applica-

tion and results of the program.

1.2 Review of Previous Work

A considerable amount of work has been carried out on the development

of analytical models and numerical methods for the studies ,of the behavior

of buildings during strong earthquakes, especially after the development

of high-speed electronic digital computers.

A review in this section is limited to those works which were

related to the development of nonlinear dynamic analysis methods with a

special emphasis on inelastic analysis of plane reinforced concrete frame

structures to earthquake motions.

When a high-speed digital computer was not readily available to

those who worked in the fields of structural dynamics and earthquake

engineering, a graphical solution was used in computing nonlinear dynamic

response of structures (Tanabashi, 1956). Or as an extension of a

classical normal mode approach for elastic structures, inelastic

response of a simple structure was computed (Bleich and Salvadori, 1955,

and DiMaggio, 1958) on the basis of normal modes at different stages

of plasticity in the structure. Application of these methods was

naturally limited to a simple structure.

 3

With advent of a digital computer, numerical methods (Newmark, 1959)

were developed to solve a set of second-order ordinary differential

equations: the equation of motion. Multistory structures were

idealized as a shear-beam system with bilinear stiffness characteristics.

The behavior of analytical models was studied under earthquake excitations

(Sheth, 1959, and Penzien, 1960).

A numerical procedure was developed by Berg and Daddepo (1960)

to analyze a rectangular plane frame with elasto-plastic members to

earthquake motion. Heidebrecht et ale (1963, 1964) developed a

numerical method to analyze plane frame structures with elasto-plastic

members under dynamic excitations, considering the effect of gravity

loads applied in girder members.

During the third World Conference on Earthquake Engineering in

New Zealand, two mathematical models were presented for the inelastic

analysis of frame structures. In the first model by Wen and Janssen

(1965), a frame member was divided into short segments along its length.
A bilinear moment-curvature relation was assigned to each short

segment. The application of this model was limited to a relatively

simple structure because a relatively large size computer program was

required for the analysis.

In the second model by Clough, Benuska and Wilson (1965), a frame

member was divided into two imaginary parallel elements: an elasto-

plastic element to represent a yielding characteristic, and a fully

elastic element to represent strain hardening behavior. It was found

difficult for this model to simulate continuously varying stiffness

(Giberson, 1967).
 4

An equivalent spring model was proposed by Giberson (1967); the

model consisted of a linearly elastic member with two equivalent nonlinear

springs at the member ends. This model may be generalization of a simple

elasto-plastic member model, which had two rigid-plastic springs at

the ends of an elastic .member. The deformation of an equivalent spring

must reflect the distribution of inelastic deformation along a member,

which, in fact, depends on the distribution of loads. Suko and Adams

(1971) determined the spring stiffness for a force distribution defined

by the location of the point of contraflexure at the initial elastic

stage.

Takizawa (1973) assumed a distribution of flexural rigidity, along

a frame member, the amplitude and distribution shape which varied as a

function of member end moments.

Otani and Sozen (1972) proposed to determine the stiffness of a

frame member as, a function of stress-history of the member and location

of inflection point.

The major difficulty in analyzing a reinforced concrete frame

dynamically and statically is that: (1) The inelastic stiffness of

a member section changes continuously with its stress history; and (2)

the inelastic stiffness distributes along a member length in accordance

with a member stress distribution.

The computer code developed in this report employed an equivalent

spring model for a frame member.

A brief review of force-deflection rules applied to an analytical

model may be useful. An elasto-plastic model was dominantly used in

inelastic analysis (Penzien, 1960, and Berg, 1961) because of its

 5

mathematical simplicity and an observed elasto-plastic stress-strain

relationship of the steel under constantly increasing loads. Then

a bilinear model became popular to simulate strain hardening character-

istics and the Bauschinger effect in steel member. Reinforced concrete

members was also idealized by the same type of hysteresis models as

steel members (Walpole and Shepherd, 1969 a and b).

Jennings (1963) adopted the "Ramberg-Osgood function" to simulate

a hysteretic behavior of a steel structure. This model was supported

by experimental observations.

Clough and Johnston (1966) proposed a "degrading stiffness"

model to simulate the behavior of reinforced concrete members under load

reversals. This model may be the first that differentiated the behavior

of reinforced concrete structures from that of steel structures in a

dynamic analysis of a structure.

Takeda, Sozen and Nielsen (1970) proposed a more complicated hysteresis

model on the basis of experimental observation of a number of reinforced

concrete members under flexural load reversals with or without axial

loads. The dominant failure mode of those test specimens was in

flexure. This set of hysteresis rules was referred to as "Takeda

hysteresis rules" in this report, and has been successful in reproducing

the behavior of reinforced concrete frame structures observed in dynamic

tests to failure (Gulkan and Sozen, 1971, and Otani, 1974).

Takeda hysteresis rules were adopted in this report to simulate

moment-rotation relations of equivalent springs.

There were more works, which stimulated the progress in the

inelastic analysis of reinforced concrete frame structures to earthquakes.

A limited number of them are shown in the list of references.

 6

2. ASSUMPTIONS OF ANALYSIS

An engineering analysis is generally carried out in the following

three steps; (a) formulation of a mathematical problem through a set of

engineering assumptions, (b) solution of the mathematical problem by an

exact or approximate method, and (c) evaluation of the solution through

the same set of assumptions.

A set of assumptions were required to simplify a complicated

engineering problem into a simple, solvable mathematical problem. At

the same time, the set of assumptions will clarify the limitations and

applicabilities of the solution.

An inelastic analysis program of a reinforced concrete frame

structure to dynamic base disturbance was developed. The solution of

the method described in this report is limited by the following

assumptions.

2.1 Assumptions about Structures

(1) A reinforced concrete structure consisted of a set of rectangular

regular plane frames. The height of floor levels was the same in any

constituent frames. The analysis was limited in the direction parallel

to the planes of frames (Fig. 2.1).

(2) Every member in the structure was treated as a massless line

member represented by its centroidal axis. The centroidal axis of

columns on a column line must be straight and continuous from the first

story to the top. The centroidal axis of beams at a floor level must

be straight and continuous from the first bay to the last in a frame.
 7
(3) The analysis was limited to the small deformation which

retained the initial configuration. However, the effect of gravity

loads (p - ~ effect) on story shear and overturning moment through the

lateral story displacement was considered.

(4) The idealized frames were assumed to be fixed at their

base on an infinitely rigid foundation.

(5) All mass at a floor level was assumed to be concentrated

outside the structure, and linked to the floor levels by rigid truss

elements. No rotational inertia was considered at the structural

joints.

(6) Axial deformation of a member (beam or column) was assumed to

be zero. Or to be more specific, all joints of any frames at a floor

level were assumed to displace the same amplitude in the horizontal

direction, and vertical displacements of any joints in the structure were

assumed to be zero.

(7) Damping forces were assumed to act at floor levels in the horizontal

direction. The amplitude of a damping force was proportional to the hori-

zontal velocity of the floor level relative to the base of the structure.

The damping matrix was assumed to be the sum of a part proportional to

the constant structural mass matrix, and a part proportional to the

variable structural stiffness matrix.

(8) Ground motion was considered only in one horizontal direction,

parallel to the planes of frames.

(9) Compatibility of deformations and equilibrium of forces at a

structural joint were assumed to be satisfied

2.2 Assumptions about Member and Joint

(1) Axial and shear deformations of a member were ignored. A

 8

member was assumed to have infinite strength and stiffness to shear and

axial loads.

(2) Inelastic flexural deformation of a member was assumed to be

concentrated and approximated by the deformation of two nonlinear

rotational springs placed at the ends of the flexible part of the

member (Fig. 2.2).

(3) The moment-rotation relationship of a "flexural rotational

spring" was assumed, under any loading conditions, to be defined by

Takeda hysteresis rules (Takeda et ale 1970) on the basis of a force-

deflection curve (a backbone curve) under monotonically increasing

loading. The backbone curve of a flexural spring was considered to be

symmetric with respect to the origin of the relationship.

(4) The force-deflection relationship was assumed to be constant

irrespective of dynamic axial loads, which varied with time in a column

member due to overturning moments.

(5) A beam-column joint core, defined as the part common to both

the beam and the column, was assumed to be infinitely rigid.

(6) Rotation at a member end due to bond slip of longitudinal

reinforcement within a joint core was simulated by the deformation of

an inelastic rotational spring placed outside of the joint core. The

moment-rotation relationship of a "joint (rotational) spring" was

assumed to follow simplified Takeda hysteresis rules on the basis of

symmetric bilinear backbone curve. The hysteresis rules did not assume

bond deterioration along longitudinal reinforcement within a joint core.

2.3 Assumptions in Numerical Integration

(1) Linear variation of acceleration was assumed in the numerical

integration of the equation of motion during a short time interval.

 9

(2) structural characteristics (mass, stiffness and damping)

were assumed to be constant during the short time interval.

 10

3. STRUCTURAL CONFIGURATION

A structure may consist of more than one regular rectangular

unbraced frame. Each component frame can have different numbers of

stories and bays. However, the floor levels, at which floor weights

were assigned, must be the same for all component frames.

Figure 3.1 shows an example of a structure, which consisted of

four frames. Frames 1 and 4 had two stories and one bay, while Frames

2 and 3 had three stories and three bays. This structure was analyzed

by the program as four frames which were connected in series at the

floor levels with rigid truss elements (Fig. 3.2). The order of four

frames was not relevant in the program. If Frames 1 and 4, Frames 2

and 3 were identical, it is sufficient to consider a structure with

Frames 1 and 2, and with one half of floor weights of the original

structure.

A total floor weight was considered at each floor level. The

distribution of a floor weight among component frames or columns, or

along girders was not considered. A floor weight was used in this

program to define the horizontal mass inertia at a floor level and

the P - 6 effect.

A program user must provide story heights and floor weights in

order to define a structure.

3.1 Regular Frame

Each substructure of a structure must b~ a regular rectangular

empty frame. In other words, the centroidal axis of each column on a

column line must be straight and continuous from the first story to
 11

the top in a frame. The centroidal axis of each girder at a floor level

must be straight and continuous from the first bay to the last.

A program user has only to provide bay widths in order to define

the configuration (dimensions) of a frame, because story heights were

assumed to be common among component frames of a structure.

3.2 Irregular Frame

With a special treatment, the program is capable of handling a case

in which either a column or beam line is not continuous in a frame.

A fictitious truss member with two mechanical hinges at both ends

can be inserted to the location to make the column or beam line continuous.

The fictitious member must be considered as a regular member in preparing

data to the program.

A program user must be reminded of the two assumptions employed in

the analysis (a) vertical displacements at any joints were considered

to be zero by ignoring axial deformation of columns, and (b) horizontal

displacements at any joints at a floor level were assumed to be identical

by ignoring in-plane deformation of girders.

Consequently, a fictitious member did not contribute the stiffness

of a structure. The fictitious members were required in the program in

order to systemize the bookkeeping of member and joint numbers.

A program user should exercise cautions in applying the program to

an irregular frame structure. The program should not be used for a

irregular frame unless the effect of the irregularity is small enough

so that the assumptions stated above would be able to approximate the

behavior of the structure.

 12

Two extreme examples are given in Fig. 3.3 to illustrate misusages

of the program. The behavior is obviously different among the two

structures with and without fictitious members in these two examples.

3.3 Treatment of Hinged and Free Joints

If all member ends of a joint have mechanical hinges, a special

treatment is required in this analysis program due to a numerical

difficulty (singular stiffness matrix). The program does not have the

capability to eliminate this condition. The same difficulty occurs in

the analysis of a cantilever column if a mechanical hinge is specified

at the free end.

In order to avoid this problem, at least one member end must be

free of a mechanical hinge at any joint.

If three member ends around a joint have mechanical hinges, then

the joint behaves as a hinge no matter if the last member end has a

mechanical hinge. Therefore, the mechanical hinge is not necessary at

the fourth member end, and should not be specified in order to avoid

the numerical difficulty of the program (Fig. 3.4).

As for the top of a cantilever column, the program automatically

assumes zero moment at the free end, and the special treatment to assign

a mechanical hinge is not necessary.

 13

4. MEMBER CHARACTERISTICS

A member was defined in the program as a structural element that

connected two adjacent structural joints of a frame (Fig. 2.2). Axial

and shear deformations of a member were ignored, and only flexural

(elastic and inelastic) deformation was considered. A part at a joint

common to both the beam and column was assumed to be infinitely rigid.

However, bond slip of longitudinal reinforcement within a joint core

was taken into consideration in terms of a member end rotation.

A member (horizontal beam or vertical column) was assumed to consist

of (Fig. 2.2).

(1) a flexible elastic prismatic line element,

(2) a nonlinear flexural rotational spring at each end of the

elastic element,

(3) a nonlinear joint rotational spring at each end of the elastic

element, and

(4) two rigid line segments outside the elastic line element.

Individual members can have different stiffness properties in a

frame. The functions of individual elements are described in the following

sections.

4.1 Elastic Line Element

A flexural deformation of a member, that was calculated on the basis

of initial elastic stiffness, was called "elastic" flexural deformation.

The rest of flexural deformation was called "inelastic" flexural deforma-

tion in this report.

 14
An elastic prismatic line element was used to simulate the elastic

flexural deformation of a member. The properties of an elastic line

element can be defined by a flexural rigidity (EI) and its length.

4.2 Inelastic Flexural Rotational Springs

Inelastic flexural deformation of a member was approximated by

the deformations of two inelastic flexural rotational springs, placed

at the ends of the flexible line element. The spring characteristics

were defined by an inelastic moment-rotation relationship of the member,

as a simply supported member, under a symmetric double curvature moment

distribution (Fig. 4.1). The moment-rotation relations of a spring

under a monotonously incleasing load is calculated in the program

(Appendix A) based on a moment-curvature relation provided by a program

user.

The approximation, to use two flexural springs, is not the best way

to simulate the inelastic flexural behavior of a member because the

rotation at an end of a member becomes independent of the existing moment

at the other end. Apparently, this is not true in an elastic beam.

However, this approximation could be reasonable, if the location of

an inflection point of a member does not shift extensively and if the

spring characteristics were to be determined on the basis of an average

location of the inflection point. The average location of an inflection

point is generally not known prior to an analysis. The initial location

of inflection point was suggested (Suko et aI, 1971) to replace the

average location. However, the initial location of an inflection point

is generally quite senstive to a loading condition, and may not be

representative throughout the analysis.

 15

Therefore, the program assumed the average location of the inflection

point of a member to be at the mid-point of the flexible line element

in evaluating a moment-rotation relationship of a inelastic flexural

spring.

This assumption about the location of inflection point is obviously

not realistic for shear wall members of a shear wall-frame structure

(Fig. 4.2), in which the inflection point may reside outside a shear
wall member. Therefore, the use of the program to a shear wall-frame

structure should be discouraged.

The inflection point does not locate close to the midspan in a

column near the top or bottom story of a high-rise frame structure

(Fig. 4.3), or in a beam close to the exterior columns. In either cases,

if member stresses are small, the behavior of the member is elastic. If

the member stresses become larger, the inelastic flexural springs start

to function. As the member is stressed further, then a member end moment

reaches to a limiting (yielding) value. Additional shear must be

carried by the member by increasing mo~ent at the other member end.

Consequently, the inflection point tends to shift inward closer to the

midspan of the member. Therefore, the assumption about the location of

inflection point may be reasonable for a frame structure.

Some effort was made to include a shifting location of inflection

point in the analysis. However, a numerical solution yielded unstable

when the inflection point moved large amount for a small change in member

end moments (Otani, 1972). Some experience was required to control the

unstable numerical conditions. Therefore, the previous program was not

suited for a general user. The program has been revised on the basis
 16
of the above mentioned simple assumption to avoid the numerical unstable

conditions.

The inelastic moment-rotation'relationship of a flexural spring

under load reversals was assumed to follow Takeda hysteresis rules

(Appendix B.l, and Takeda et aI, 1970). The hysteresis rules were

developed on the basis of observed behavior of many reinforced concrete

members under static lateral load reversals with light to medium

amount of axial loads. Severe stiffness reduction due to shear failure

nor bond deterioration, neither of which is obviously not desirable in

a reinforced concrete structure, was not considered in the hysteresis

rules.

The major characteristics of the hysteresis rules are (1) the stiff-

ness changes with a stress (strain) history, and (2) some energy can be

dissipated even from a small amplitude stress reversal.

Finally, it should be noted that the actual location of inflection

point of a member varies in the analysis by ,this program in accordance

with the distribution of forces and stiffness in a structure.

4.3 Rigid Zones

A beam-column joint core, defined as the part common to both the

beam and the column, was assumed to be infinitely rigid in flexure.

Therefore, a part between the center of a joint to the column face of a

beam member, or to the beam face of a column member, was considered to

be rigid in the program.

The length of a rigid zone can be varied for individual members.

4.4 Inelastic Joint Rotational Spring

It has been pointed out by many researchers (for example, Bertero

and Bresler 1969, Ismail and Jirsa 1972) that the stiffness of a
 17
reinforced concrete member was influenced by the deterioration of bond

between concrete and steel within a joint core due to large load

reversals. Therefore, the rotation caused by bond slip of tensile reinforc-

ing steel along its embedded length in a joint core was simulated by the

deformation of an inelastic rotational spring, although a joint core was

assumed to be rigid in flexure.

Excessive bond deterioration may prevent the longitudinal reinforce-

ment from developing its strength until after a large deformation is

taken place by the bond slip, which is definitely not desirable in a

structure. Therefore, the program did not consider such a case.

Morita et al (1973) proposed hysteresis rules for a steel stress-

bond slip relation under load reversals. The proposed rules were

characterized by a very steep unloading slope after a peak stress was

reached, and by a wide flat slip range after a stress was reversed.

However, Morita's hysteresis rules were not available when the

program was originally developed, nor were available any other explicit

hysteresis rules about the bond slip. Therefore, the program adopted

simplified Takeda hysteresis rules (Appendix B.2) with a bilinear

backbone curve to simulate a moment-rotation relation of an inelastic

joint rotational spring.

The joint rotational spring can also be used to simulate mechanical

hinge at an end of the flexible portion of a member.

 18

5. ANALYSIS METHOD

A nonlinear dynamic analysis of an idealized structure is generally

carried out on the basis of the following four relations:

(1) nonlinear relations between forces and resultant displacements;

(2) compatibility relation of displacements at a structural joint;

(3) equilibrium relation of forces (inertia, damping, resistance,

and external forces) at a structural joint; and

(4) differential relation among time response functions (accelera-

tions, velocities and displacements).

The solution is not available in a closed form because of the nonlinear

nature of stiffness characteristics and differential relation among time

response functions.

Nonlinear stiffness relation was approximated by a piecewise linear

relations by considering a small increment of force and displacsnent.

The numerical method in this report employed two assumptions to

linearize the relations among time response functionSe The first

assumption was that the change in structural properties during a short

time interval of response calculation was small enough to be neglected.

In other words, incremental forces and displacements were assumed to

be linearly related by instantaneous stiffness during a short period,

and incremental damping forces was proportional to incremental velocites.

The second assumption was that acceleration changed linearly with

time during a short period. This assumption provided the linear relations
 19

among displacements, velocities and accelerations, and has been known

to be reasonable in solving numerically the equation of motion of a

linearly elastic system (Newmark, 1959), if a sufficiently short time

interval is chosen in a step-by-step numerical integration procedure.

The procedure of analyzing a nonlinear structure can be summarized

as follows (Fig. 5.1):

(1) Linearization of a nonlinear member end force-displacement

relation by an instantaneous stiffness matrix for a small increment of

force and displacement (Appendix C).

(2) Formation of an instantaneous structural stiffness matrix on

the basis of the compatibility relation of displacements at a joint

and the equilibrium of external forces at a joint with resisting

structural forces (Appendix D). The structural stiffness matrix

related incremental joint displacements to incremental joint forces.

(3) Formulation of the equation of motion on the basis of the

equilibrium condition of all forces at each joint. Inertia forces

(proportional to absolute acceleration, Appendix E), damping forces

(proportional to velocity, Appendix E), resisting forces (proportional

to displacement), gravity effect forces (p - 'D. effect terms, Appendix F)

and external forces were set to be in equilibrium. The equation of

motion is an ordinary differential equation of the second order, and

is not solvable in a closed form.

(4) Numerical integration of the ordinary differential equation

(Appendix G). Various numerical methods are available and have been

used to solve the quation of motion of structures (Newmark, 1959, Goel,

1967, Clough et aI, 1965, Jennings, 1963, Wilson, 1973, and others).
 20

The solution of the equation of motion provided incremental displacements

at each structural joint.

(5) Transformation of the incremental displacement at structural

joints to incremental displacements at member ends through the compati-

bility relations.

(6) Computation of incremental member forces as a product of an

instantaneous member stiffness matrix and the incremental member dis-

placements. The resultant member forces were calculated as the sum of

the previous and incremental member forces.

If calculated member forces and displacements satisfy the force-

deflection relationship of the member, the solution of a mathematical

problem is completed within an error range defined by the accuracy of

numerical integration assumptions and of numerical calculations involved.

Calculated member forces and displacements may not satisfy the force-

deflection relationship because of a nonlinear nature of the relationship

(Appendix H.). The calculated deflection and force must be modified to

fit a specified relationship. A modification of either a member force

or displacement will violate either the equilibrium or the compatibility

condition at a structural joint. This program corrects member force to

fit a given relation at the calculated deflection. The mbdified force

was considered at the joint as an external force in the equation of motion

in the next step of the step-by-step response computation.

 21

6. DATA PREPARATION

The data cards for the analysis of a structure can be grouped into

the following four sets:

(1) structural cards,

(2) typical member cards,

(3) component frame cards, and

(4) earthquake cards.

Data cards must be prepared and arranged in a fixed order, using a

fixed format. The basic rules are as follows:

(1) All input values must be defined in a consistent unit system.

The unit system can be arbitrarily chosen by a program user.

(2) An integer number must be punched right-justified in a

specified area (usually 5 columns) of a card. If an integer is to be

punched in columns 1 through 5 of a card, number 3 in the fourth column

indicates number 30.

(3) A floating (decimal) number must be punched within a specified

area (usually 10 columns). The decimal point can be placed at any

location within the area.

The content and format of individual data cards are described card

by card in the following sections. "Column" is referred to as column

number of a data card. Each data card has 80 columns. "Notation" is

a symbol used in the program to represent the value in a specified

column area. "Entry" is the description of the value.

 22

6.1 Structural Cards

A set of data within this group is needed to define the overall

structure, or to provide common constants among members and frames.

Five cards must be provided in this group for a structure less than

or equal to eight stories high. Two additional cards must be provided

for each eight stories in excess of the first eight stories.

Figure 6.1 is provided for reference.

6.1.1 Structure Name Card (20A4)*

The structure name card (one card) can be used to identify a

structure, problem, solution, or anything. Only one card is allowed.

The information in this card will. be printed in the first page of output.

COLUMN NOTATION ENTRY

1-80 STNAME Anything to identify a structure. Any

alphameric (alphabet or numeric) symbols

are allowed in this card.

6.1.2 Unit Card (3A4)

The unit card (one card) is to identify units of length, force and

time employed in the preparation of data. All data must be provided

in this unit system. The information is not used in the program, which

does not have capability to change the units. The calculated results

will be given in the same unit system.

*FORTRAN IV input format used in the program

 23

COLUMN NOTATION ENTRY

1-4 UNIT(l) Unit of length (alphameric)

5-8 UNIT(2) Unit of force (alphameric)

9-12 UNIT(3) Unit of time (alphameric)

13-80 (blank)

6.1.3 Control Information Card (515, 5X, 4FIO.0)

The control information card (one card) defines constants to be used

in a structure. A structure may consist of more than one frame.

Members in a structure can be grouped by common dimensions (lengths

of rigid and flexible zones), or by common flexural characteristics

(a moment of inertia and moment-curvature relationship)) or by common

joint rotational characteristics (moment-rotation relationship). Each

set of common properties is referred to as a typical member. dimension,

a typical flexural element, or a typical joint spring.

The damping matrix of a structure was assumed to be made of a part

proportional to a structural mass matrix and a part proportional to a

structural stiffness matrix (Appendix E).

With this background information, the contents of the control

information card can be described belowg

COLUMN NOTATION ENTRY

1-5 NSTORY Maximum number of stories in any of the constitu-

ent frames of an entire structure. Positive

integer. Limitation in the program is 10.

6-10 NROWS Number of constituent frames in a structure.

 24

COLUMN NOTATION ENTRY

Positive integer. Limitation in the program

is 5.

11-15 NSTIFO Number of different typical member dimensions

in a structure. The actual values of NSTIFO

typical member dimensions are to be given in

a set of cards described in 6.2.1. Positive

integer. Limitation in the program is 20.

16-20 NSTIFl Number of different typical flexural elements

in a structure. The actual values of NSTIFl

typical flexural element properties are to be

given in a set of cards described in 6.2.2.

Positive integer. Limitation in the program

is 20.

21-25 NSTIF2 Number of different typical joint rotational

springs in a structure. The actual values of

NSTIF2 typical joint spring properties are

to be given in a set of cards described in

6.2.3. Positive integer. Limitation in the

program is 20.

26-30 (blank)

31-40 AE Modulus of elasticity of concrete. Positive

decimal.

41-50 G Value of gravity acceleration. Positive decimal.

Computed acceleration response will be normalized

by this number in output.

 25

COLUMN NOTATION ENTRY

51-60 ALPH A proportionality constant of a part of a

damping matrix, which is proportional to a

mass matrix. For a linear structure, the

value can be determined by

in which hI and h2 = damping factors for the

first two modes, and WI and w = circular

2
frequencies of the first two modes. Zero.

or positive decimal.

61-70 BETA A proportionality constant of a part of a

damping matrix, which is proportional to a

stiffness matrix. For a linear structure,

the value can be determined by

Zero or positive decimal.

71-80 (blank)

6.1.4 Weight Card (8FIO.0)

Weights at each floor level of an entire structure must be given

sequentially from the first to the top floor using 10 columns of a data

card. If one card is filled (the number of stories greater than 8) ,

the following card must be used. Positive decimal. The number of weight

values must be equal to NSTORY.

 26

For a three story structure, the card looks as follows:

COLUMN NOTATION ENTRY

1-10 WT(l) Total weight at the first floor level.

11-20 WT( 2) Total weight at the second floor level.

21-30 WT(3) Total weight at the third floor level.

31-80 (blank)

6.1.5 Inter-story Height Card (8FIO.0)

Inter-story height is a dimension between two adjacent floor levels.

Inter-story height of each story must be given sequentially from the first

to the top story, using successive 10 columns of a data card for each

story. If one card is filled (the number of stories is greater than 8),

the following card should be used. The number of values provided in the

inte~story card must be equal to the number of stories of the structure

(NSTORY). Positive decimal. The inter-story heights must be the same

in any constituent frames of a structure.

For a three story structure, the card should look as follows:

COLUMN NOTATION ENTRY

1-10 HT(l) Inter-story height of the first story.

11-20 HT(2) Inter-story height of the second story.

21-30 HT(3) Inter-story height of the third story.

31-80 (blank)
 27

6.2 Member Property Cards

Many members in a structure often share common dimensions, and

stiffness properties. A set of data within this group define all

different member properties in an entire structure.

Member properties were divided into three groups: (1) member

dimensions (lengths of flexible and rigid zones), (2) flexural properties

of the flexible part of a member, and (3) stiffness properties of a

joint rotational spring.

6.2.1 Member Dimer:sion Cards (15, 5X, 3FIO.0)

Each member can have a flexible portion and rigid portions at its

ends. The total length (sum of lengths of flexible and rigid portions)

of a member must be equal to one of inter-story heights or bay widths

of a frame.

The total number of member dimension cards must be equal to NSTIFO

of the control information card in 6.1.3. The information needed here

can be used for both vertical (column) and horizontal (beam) members.

Figure 6.2 is provided for reference.

COLUMN NOTATION ENTRY

1-5 I Typical member dimension identification number.

Positive integer from 1 to NSTIFO. Dimensions

of any member in a structure must be specified

by this number in 6.3.3.

6-10 (blank)
 28

COLUMN NOTATION ENTRY

11-20 BA(I) Length of a rigid portion of left end of a beam,

or top end of a column. Zero or positive

decimal.

21-30 AL(I) Length of a flexible portion of a member.

Positive decimal.

31-40 BB(I) Length of a rigid portion of right end of a

beam, or top end of a column. Zero or positive

decimal.

41-80 (blank)

6.2.2 Flexural Property Cards (15, 5X, 6FIO.0)

Flexural properties of a typical member section should be provided

in a flexural property card. The total number of cards provided here

must be equal to the number (NSTIF1) of typical flexural elements defined

in the control information card in 6.1.3. stiffness properties of a

flexural rotational spring are computed on the basis of this information.

The moment-curvature relationship of a member section under monoton-

ously increasing load was assumed to be symmetric with respect to the

origin of the relationship.

The flexural properties (Fig. 6.3) referred here consist of

(1) moment of inertia of uncracked transformed section about the

centroid of the section;

(2) moment CM about the centroid of the section at the tensile

cracking of the concrete under the existing constant axial load;

 29

(3) moment YM about the centroid of the section and curvature

YP at the yielding of tensile reinforcement under the existing constant

axial load; and

(4) moment UM about the centroid of the section and curvature UP

at the crushing of compressive concrete under the existing constant

axial load.

Curvature CP at the tensile cracking of concrete will be found from

cracking moment, modulus of elasticity and moment of inertia. The

moment-curvature relationship of a section was assumed to be piece-wise

linear (Fig. 6.4). If the behavior of a member is to be limited to

elastic one, a large 'cracking moment should be specifIed.

COLUMN NOTATION ENTRY

1-5 J Typical flexural element identification number.

Positive integer from 1 to NSTIF1. Flexural

properties of any member will be specified

by this number.

6-10 (blank)

11-20 AI(J) Moment of inertia of uncracked transformed

section about the centroid of the section.

Positive Decimal.

21-30 CM(J) Moment about the centroid of a section at the

tensile cracking of concrete under an existing

constant axial load. Zero or positive decimal.

 30

COLUMN NOTATION ENTRY

31-40 YM(J) Moment about the centroid of a section at the

yielding of tensile reinforcement under an

existing axial load. YM(J) > CM(J). Positive

decimal.

41-50 UM(J) Moment about the centroid of a section at the

crushing of compressive concrete under an

existing axial load. UM(J) > YM(J). Positive

decimal.

51-60 YP(J) Curvature of a section at bending moment

YM(J). Positive decimal.

61-70 UP(J) Curvature of a section at bending moment

UM(J). UP(J) > yp(J). Positive decimal.

71-80 (blank)

The program did not consider moment UM to be the failure moment of

a section, but used point (UM, Up) to be a point on the last slope of

piece-wise linear moment-curvature curve after yielding.

6.2.3 Joint Spring Cards (15, 5X, 4FlO.0)

A joint spring simulated the bond slip of tensile reinforcement in

a joint core. The stiffness characteristics of a joint spring were

defined by yield moment VY and rotation WY, ultimate moment VU and

rotation WU. If loads were to be reversed in a joint spring, stiffness

 31

characteristics were assumed to follow simplified Takeda hysteresis

rules (Appendix B.2).

The joint spring can be used to simulate a mechanical hinge at the

end of a flexible portion of a member.

The properties of each typical joint spring must be provided in

a card. The total number of cards provided here must be equal to the

number (NSTIF2) of typical joint springs defined in the control information

card in 6.1.3.

If tensile reinforcement was not provided with sufficient development

length in a joint core, or if an extensive bond deterioration is suceptible,

the hysteretic behavior of bond slip of a tensile reinforcement will be

quite different from the way it was assumed in this program. The program

is not recommended in such a case. The behavior under an extensive bond

deterioration has not been fully understood to formulate general realistic

hysteresis rules.

COLUMN NOTATION ENTRY

1-5 K Typical joint spring identification number.

Positive interger from 1 to NSTIF2. Properties

of a joint spring will be referred to by this

number.

6-10 (blank)

11-20 VY(K) Moment at the break-point of an idealized

bilinear back-bone curve (Fig. 6.5), if positive.

If zero, the program assumes no joint

spring at the location.

 32

COLUMN NOTATION ENTRY

If negative, a mechanical spring is

assumed at the location instead of a spring.

21-30 VU(K) Moment at a point on the second slope of an

idealized bilinear curve. This value is

ignored if a mechanical hinge or no spring

is specified (VY(K) ~ 0.0).

31-40 WY(K) Rotation at the break-point of an idealized

bilinear back-bone curve (Fig. 6.5). This value

is ignored if a mechanical hinge or no spring

is specified (VY(K) ~ 0.0).

41-50 WU(K) Rotation at a point (m~ment = VU(K)) on the

second slope of an idealized bilinear back-bone

curve. This value is ignored if a mechanical

hinge or no spring is specified (VY(K) ~ 0.0).

51-80 (blank)

6.3 Component Frame Cards

A structure may consist of more than one frame. The properties of

an individual component frame are defined in the component frame cards.

A set of cards described in this section must be grouped for each

frame, and placed one group after another. The number of groups must be

equal to the number (NROWS) of constituent frames in the control infor-

mation card (6.1.3).

 33

6.3.1 Frame Control Card (215)

A frame control card defines the numbers of stories and bays of

a frame.

COLUMN NOTATION ENTRY

1-5 MSTORY Number of stories of IROW-th frame. Positive

(IROW)
integer, but less than or equal to the number

(NSTORY) of stories of a structure in the

control information card (6.1.3) .

6-10 MBAYS Number of bays of IROW-th frame. Zero or

( IROW),
positive integer. Limitation in the program

is 10.

11-80 (blank)

Limitation in the numbers of stories and bays of a frame is as

follows due to dimensioning of variables in the program.

Number of Stories Number of Bays

10 4
9 5
8 5
7 6
6 7
5 8
4 9
3 10
2 10
1 10
 34

6.3.2 Bay Width Card (8FIO.0)

If the number of bays of a frame in question is zero (for example,

an independent wall, or a cantilever column), this card should be omitted.

Bay width is a dimension between two adjacent column lines. Bay

widths must be given sequentially from left to right, using 10 columns

of a data card for each bay. If one card is filled (more than 8 bays in

a frame), the following card should be used.

For a three-bay frame, the card looks as follows:

COLUMN NOTATION ENTRY

1-10 WIDTH(l) Width of the first bay from left in a frame.

Positive decimal.

11-20 WIDTH(2) Width of the second bay from left in a frame.

Positive decimal.

21-30 WIDTH(3) Width of the third bay from left in a frame.

Positive decimal.

31-80 (blank)

6.3.3 stiffness Pointer Cards (615)

A stiffness pointer card related a member in a frame with typical

properties (typical member dimensions, typical flexural elements, and

typical joint springs).

The members in a frame must be numbered (Fig. 6.6) in a fixed

manner in the program in order to take advantage of regular nature of

a rectangular frame.
 35
(1) Beam members must be numbered first from left to right at a

given floor level, and from the top down to the first floor level,

(2) Column members must .be then numbered from left to right at

a given story, and from the top to the first story, and

(3) Member numbers must start from 1 in each frame, being sequential

and continuous for beams and columns.

The number (NE) of members in a frame is given as

NE = MSTORY(IROW) * (MBAYS(IROW) * 2 + 1)

and the number of ~tiffness pointer cards of a frame must be equal to

the number of members for each frame.

A member can have one set of dimensions (lengths of flexible and

rigid parts), two sets of flexural properties for two member ends, and

two sets of joint spring properties for two member end springs.

COLUMN NOTATION ENTRY

1-5 L Member identification number in a frame.

Positive integer, but less than or equal to

the number (NE) or members in a frame.

6-10 ISTFO(L) Typical dimension identification number. The

dimensions specified under this number in

6.2.1 are to be used for this member. The

number must be positive, and less than or

equal to the number NSTIFO of typical member

dimensions of the control information card

(6.1.3).
 36

COLUMN NOTATION ENTRY

11-15 ISTFlA(L) Typical flexural element identification number.

The flexural properties specified under this

number in 6.2.2 are to be used in the left

half part of a beam member, or the top half

part of a column member. Positive integer,

but less than or equal to the number (NSTIFl)

of typical flexural elements of the control

information card (6.1.3).

16-20 ISTFIB(L) Typical flexural element identification number.

The flexural properties specified under this

number in 6.2.2 are to be used in the right

half part of a beam member, or the 'bottom

half part of a column member. Positive

integer, but less than or equal to the number

(NSTIFl) of typical flexural elements of the

control information card (6.1.3).

21-25 ISTF2A(L) Typical joint spring identification number.

The spring properties specified under this

number in 6.2.3 are to be used in the left

end joint spring of a beam member, or the top

end joint spring of a 'column member. Positive

integer, but less than or equal to the number

(NSTIF2) of typical joint springs of the control

information card (6.1.3).

 37

COLUMN NOTATION ENTRY

26-30 ISTF2B(L) Typical joint spring identification number.

The spring properties specified under this

number in 6.2.3 are to be used in the right

end joint spring of a beam member, or the

bottom end joint spring of a column member.

Positive integer, but less than or equal to

the number (NSTF2) of typical joint springs

of the control information card (6.1.3).

31-80 (blank)

6.4 Earthquake Cards

An earthquake acceleration history must be given from data cards.

The record can be either accelerations at a constant interval, or pairs

of time and acceleration.

6.4.1 Earthquake Name Card (20A4)

The earthquake name card can be used to identify an earthquake

record, its magnitude, length and scale factors. Only one card is

allowed. The information in this card will be printed to identify

the earthquake.

COLUMN NOTATION ENTRY

1-80 EQN~ME Anything to identify the earthquake record

to be analyzed. Any alphameric symbols are

allowed in this card.

 38

6.4.2 Response Control Card (415, 5FIO.0)

The response control card contains the information about the

earthquake record data cards, and response calculation.

The response of a structure is calculated at a small time increment,

the value of which must be much smaller than a value required by the

convergence criteria of numerical integration in order to reduce the

effect of unbalanced forces (Appendix H). It is recommended to study

the effect of the size of time increment on the maximum response of a

structure. If the maximum response values are greatly affected by the

choice of a time interval size, the size may be too large. A value of

approximately 1/50 of the shortest elastic period of a structure may

be used as a conservative estimate.

An earthquake acceleration record at a uniform time interval or at

a random time interval with time coordinates can be accepted in this

program. The data format of earthquake record cards can be of any form

as long as the format does not change from a card to card.

COLUMN NOTATION ENTRY

1-5 IESIZE Number of acceleration data in a single earthquake

record card. Positive integer, but less than

or equal to 10.

6-10 NSTART The first (NSTA~T-l) points of acceleration

record data will be ignored in the response

calculation. If the number is less than or

equal to 1, the response calculation will be

made from the first acceleration history point.

 39

COLUMN NOTATION ENTRY

11-15 NSKIP Number of calculated response points that will

not be used to check the maximum response

values and stored on the output tape. As

the time interval of response calculation

is very small, it is not recommended to store

every calculated response value.

If the number is less than 1, the program

will assume number 1 for NSKIP, and store

every point.

16-20 KFLAG Indicator to describe if the earthquake

acceleration history is given by acceleration

coordinates at a constant interval (KFLAG=l),

or by pairs of time and acceleration coordin-

ates (KFLAG = 0). Number must be either 0 or 1.

21-30 DT Time increment for rrwmerical integration of

the equation of motion (response calculation).

A value of approximately 1/50 of the shortest

elastic period of a structure may be used as

a conservative value.

31-40 ASCACE Scale factor to be multiplied to the accelera-

tion data in order to modify the intensity of

an earthquake motion, or to change the unit

of accelerations to the unit system used in

an analysis.
 40

COLU:MIif NOTATION ENTRY

41-50 TSCALE Scale factor to be multiplied to the time

coordinates given in earthquake record cards.

Generally this value is 1.0.

51-60 ZERO Value to be subtracted from the original

acceleration coordinates given in earthquake

record card before it can be multiplied by

ASCALE. This is to shift the base line of

the input acceleration waveform. Generally,

this value is 0.0.

61-70 DELT Time interval of earthquake acceleration record,

if KFLAG=l. In other words, if earthquake

acceleration coordinates are to be given at a

uniform time interval, this interval must be

given here. If earthquake time and acceleration

coordinates are given (KFLAG=O), this value

should be ignored.

71-80 (blank)

6.4.3 Earthquake Format Card (20A4)

This card describes the format of earthquake record cards, and

should have ~ither of the following forms:

(1) KFLAG=O: variable time interval

(n (Fn . n , Fn4 . n ) ,I)

l 2 3 5
 41
(2) KFLAG=l: constant time interval

in which

n = number of acceleration data in a card, positive integer.

l
This number should be equal to IBSIZE given in the response

control card (6.4.2).

n = positive integer constant specifying the number of columns
2
to be used for a time coordinate in an earthquake record

card.

n3 = positive integer constant specifying the number of decimal

places of fractional portion within column n for a time co-

2
ordinate.

n4 = positive integer constant specifying the number of columns

to be used for an acceleration coordinate in an earthquake

record card.

n5 = positive integer constant specifying the number of decimal

places of fractional portion within column n4 for an

acceleration coordinate.

COLUMN NOTATION ENTRY

1-80 FMT Format of earthquake record cards, which is

of the form

or

in which a blank is not allowed in the parentheses.

 42
6.4.4 Earthquake Cards (FMT)

Earthquake acceleration history must be given in this set of cards,

defined either by acceleration coordinates at a uniform interval or by

time and acceleration coordinates at an arbitrary interval. In the

latter case, a time coordinate must be given prior to an acceleration

coordinate.

The format of earthquake card data must be consistent with the one

given in earthquake format card (6.4.3). The unit of acceleration must be

consistent with the unit system used throughout the analysis and given

in the unit card (6.1.2). If the unit of acceleration in earthquake

cards is different from the one used in the analysis, ASCALE of the

response control card (6.4.2) must be used to change the unit.

If an acceleration value is not given at the right time of response

calculation, the acceleration value is determined by the linear inter-

polation of two adjacent points.

Suppose the earthquake format card (6.4.3) is of the following

form (KB'LAG = 0),

(3 (F6.3, F12.4),/)
then the earthquake card must be as follows:

COLUNrN NOTATION ENTRY

1-6 TC(l) Time at the first acceleration point given

in this card. Positive decimal. Three columns

to the right of the decimal point.

 43

COLUMN NOTATION ENTRY

7-18 AC(l) Acceleration at the first point given in this

card. Decimal. Four columns to the right

of the decimal point.

19-24 TC(2) Time at the second point given in this card.

TC(2) > TC(l).

25-36 AC( 2) Acceleration at the second point given in

this card.

37-42 TC(3) Time at the third point given in this card.

TC ( 3 ) > TC ( 2) .

43-54 AC(3) Acceleration at the third point given in

this card.

55-80 (blank)

A blank card must be placed at the end of the whole earthquake

record cards to terminate the program.

 44

7. DESCRIPTION OF COMPUTER OUTPUT

The computer program prints input and output information about a

structure and an earthquake, and has an option to write calculated

results on a magnetic tape.

7.1 Print Output

All the input information is identified and interpreted in the

computer program, and then printed on the output sheet.

A user of the program is urged to examine the printed information

carefully so that a structure defined by the user and that interpreted

by the program may be the same.

Common and frequent errors have been caused by a wrong data format

~mployed in data cards. A user may often be able to find an error by

checking the printed information.

7.1.1 Numbering System in a Frame

It should be necessary to describe the joint numbering system of a

frame in the program. Joints of a frame were numbered as follows (Fig.

7.1) :
(1) starting with 1 for each frame,
(2) sequentially from left to right at a given floor level,

(3) continuously from the top to the first floor, and

(4) without numbering joints at the base.

The last joint number of an n-story m-bay frame must be equal to

n(m + 1).
 45

Members of a frame were also numbered in the same manner:

(1) Beams were numbered first, starting with 1 for each frame,

sequentially from left to right at a given floor level, and from the

top to the first floor, and

(2) columns were numbered, continuous with beam numbers, sequentially

from left to right at a given story, and from the top to the first story.

The last beam number of an n-story m-bay frame must be equal to n x m,

and the last column number must be equal to n(m + m + 1).

For each member (beam or column), member ends were called A and B

(Fig. 7.2). The A-end of a member means the left-end of a beam member,

or the top end of a ~olumn. The B-end of a member in~icates the right

end of a beam member, or the bottom end of a column member.

7.1.2 Maximum structural Response

The calculated maximum response is printed at the end of the

response computation. The maximum structural response values printed

include

(i) maximum and minimum (largest value in the negative direction)

displacements at floor levels relative to the base,

(ii) maximum and minimum absolution accelerations at floor levels,

and

(iii) maximum and minimum base shears and overturning moments.

The values of acceleration response were normalized by the gravity

acceleration, and given in g.

7.1.3 Maximum Member Response

The maximum response values of each frame member were printed in

terms of moments and rotations in flexural and joint springs placed

at the ends of a flexible portion of a member.

 46
A printed rotation of a flexural rotational spring consists of an

elastic and an inelastic rotation at an end of a unit-length simply

supported member due to a moment at the same end (Fig. 7.3). The value

may be interpreted as a rotation at an end of a simply supported member

by multiplying the member length (Appendix A.2).

Rotations of a unit-length simply supported member at tensile cracking

of the concrete, tensile yielding of the longitudinal reinforcement, and

crushing of the concrete at the extreme compressive fiber are printed

as rotation characteristics of flexural springs (elastic deformation

included) for typical flexural springs. Therefore, a "ductility factor"

of a member end can be calculated as a ratio of maximum response rotation

to the rotation at tensile yielding of the reinforcement.

Maximum and minimum rotations of each joint spring are also printed

with the corresponding moments.

Clockwise rotations and moments (Fig. 7.4) were considered to be

positive throughout the program.

7.2 Output on Magnetic Tape

The program is capable of writing calculated response values on a

magnetic tape if alphabetical symbols "c" are removed from the first

column of the following cards.

Those cards are:

(1) Twelve consecutive cards in MAIN program.

C WRITE (21) (STNAME(I),I=1,20),(EQNAME(I) ,I~1,20),

(UNIT(I) ,1=1,3),

C 1 NSTORY,NROWS,NMMBR,NNODE,G,DT,NSKIP,
C 2 NSTIFO,NSTIF1,NSTIF2,
C 3 (SM(I),I=l,NSTORY),(HT(I),I=l,NSTORY),
 47

C 4 (MBAYS(I),I=l,NROWS),(MSTORY(I), I=l,NROWS),
C 5 (AL(I),I=l,NSTIFO),
C 6 (BA(I),I=l,NSTIFO),(BB(I),I=l,NSTIFO),
C 7 (CM(I),I=1,NSTIF1) ,(YM(I),I=1,NSTIF1),
C 8 (DC(r) ,I=1,NSTIF1) ,(DY(I) ,I=1,NSTIF1),
C 9 (VY(I) ,I=1,NSTIF2),(WY(I) ,I=1,NSTIF2) ,
C 1 (EI(I) ,I=l,NMMBR),(ISTFO(I) ,I=l,NMMBR),
C 2 (ISTF1(I),I=1,NMMB2) ,(ISTF2(I),I=1,NMMB2)

(2) Four consecutive cards in STORE4 program.

C WRITE (21) GR2,(ACC(I),I=1,NSTORY),(VEL(I) ,I=l,NSTORY) ,

C 1 (DSP(I), 1=1, NSTORY)
C 2 (F2( I) ,I=l ,NMMB2) ,( D2( I) ,1=1 ,NMMB2) ,
C 3 (W2(I),I=1,NMMB2)

(3) Fifteen consecutive cards in MXRsp4 program.

C GR2=0.0
C DO 40 I=l,NSTORY
C ACC(I)=O.ODO
C VEL(I)=O.ODO
C 40 DSP(I)=O.ODO
C DO 50 I=1,NMMB2
C F2(r)=O.o
C D2(I)=0.0
C 50 W2(I)=0.0'
C DO 60 J=l,lO
C 60 WRITE (21) GR2, (ACC (I) ,1=1 ,NSTORY) , (VEL( I) ,1=1 ,NSTORY) ,
C 1 (DSP(I) ,I=l,NSTORY) ,
C 2 (F2(I) ,I=1,NMMB2) ,(D2(I) ,I=1,NMMB2),
C 3 (W2(I),I=1,NMMB2)
c REWIND 21

In order to read the content of the tape, a user has to have a

knowledte about prograw~ing. The content of the tape will be described

in the following section to aid a user make a short program for printing

or plotting stress or strain history of some specific members.

7.2.1. Basic Information

The program writes on a magnetic tape the basic characteristics of

the structure in the main program. The information that is written is

as follows:
 48

STNAME: a string of alphameric symbols to identify a structure

(6.1.1) .
EQNA.ME: a string of alphameric symbols to identify an earthquake

record (6.4.1).
UNIT a string of alphameric symbols to identify the units of

length, force and time (6.1.2).

NSTORY: integer to designate the number of stories of a structure

(6.1.3).
NROWS integer to designate the number of constituent frames in

a structure (6.1.3)

NMMBR integer to show the total number of members in a structure

(sum of numbers of members of each constituent frame).

NNODE integer to show the total number of joints in a structure

(sum of numbers of joints of each constituent frame without

counting joints along the base line).

G decimal number to give a value of gravity acceleration in

the unit system used (6.1.3).

DT decimal number to give a time interval of numerical integration

(6.42).
NSKIP integer number to indicate how many points were skipped along

time axis in recording calculated response values. The time

interval of the response points stored on the tape will be

~T*NSKIP (6.4.2).
NSTIFO: integer to show the number of different typical member

dimensions (6.1.3).
NSTIF1: integer to show the number of different typical flexural

elements (6.1.3).
NSTIF2: integer to show the number of different typical joint springs

(6.1.3).
SM array of decimal numbers that describe the values of mass

at each floor level from the first to the top floor level

(6.1.4).
HT array of decimal numbers that describe the interstory height

from the' first to the top story (6.1.5).

MBAYS array of integers to denote the number of bays of each

constituent frame (6.3.1).

MSTORY: array of integers to denote the number of stories of each

constituent frame (6.3.1).

AL array of decimal numbers that represent the length of the

flexible portion of different typical member dimensions (6.2.1).

BA array of decimal numbers that represent ratios of a rigid

zone lengths at the left end of beams or at the top of columns

to the flexible length (AL) of typical member dimensions.

BB array of decimal numbers that represent ratios of a rigid zone

length at the right end of beams or at the bottom of columns

to the flexible length (AL) of typical member dimensions.

eM array of decimal numbers to represent moments at tensile

cracking of the concrete of typical flexural elements.

YM array of decimal numbers to represent moments at yielding of

the tensile reinforcement of typical flexural elements.

 50

DC array of decimal numbers that represent rotations at tensile

cracking at an end of a unit length simple beam without moment

CM applied at the same end. The section consists of a typical

flexural element. In order to find a rotation of a simply

supported beam, value DC must be mutliplied by the beam

length.

DY array of decimal numbers to represent rotations at yielding

of tensile reinforcement of a typical flexural element at an

end of a unit length simple beam with moment YM applied at

the same end. Value DY must be multiplied by the beam length

in order to find a rotation of an arbitrary length beam.

VY array of decimal numbers to represent moments at a break point

(Fig. 6.5) of typical joint springs.

WY array of decimal numbers to represent rotations at a break

point (Fig. 6.5) of typical joint springs.

EI array of decimal numbers to represent an elastic flexibility

constant (L/6EI) of each member.

ISTFO array of integers that point an identification number of

typical member dimensions to be used for a member.

ISTFI array of integers that point a pair of identification numbers

of typical flexural elements to be used for a member.

ISTF2 array of integers that point a pair of identification numbers

of typical joint springs to be used for a member.

7.2.2 ResponseHistory

Response values at a time interval of DT*NSKIP are written in

subroutine STORE4. These values that are written on a magnetic tape

 51

are as follows:

GR2 earthquake acceleration at the base of the structure in the

units used in the analysis.

ACe array of decimal numbers which represent accelerations at floor

levels relative to the base. The values are given from the

first to the top floor level. The units of acceleration is

consistent with the units used in the analysis.

VEL array of decimal numbers which represent velocities at floor

levels relative to the base. The values are given from the

first to the top floor level in the units consistent with those

used in the analysis.

DSP array of decimal numbers which represent displacements of floor

levels relative to the base. The values are given from the

first to the top floor level.

F2 array of decimal numbers that represent a pair of bending

moments at the ends of flexible part of a member. The flexural

spring and joint spring share the same moment F2 because the

two springs locate in the same place. The values are given for

all members in the structure continuously.

D2 array of decimal numbers that represent total flexural rotations

at the ends of the flexible portion of a member. The rotations

given here are not direct rotations at each member end, but

are those flexural rotations developed at the end of a unit-

length simply supported member due to bending moment F2 applied

at the same end, including elastic deformation.

W2 array of decimal numbers that represent a pair of rotations

in the joint springs at the ends of the flexible part of a

member due to bending moment F2.

 52

The total number of response history points will be printed in the

last page. After the response calculation, ten record blocks of response

history numbers of null values will be written on the magnetic tape in

subroutine MXRsp4.
 53

LIST OF REFERENCES

1. Anderson, J. C., and V. V. Bertero, "Seismic Behavior of Multistory

Frames Designed by Different Philosophies," Earthquake Engineering
Research Center, University of California at Berkeley, Report
No. EERC-69-11, October 1969.

2. Anderson, J. C. and R. P. Gupta, "Earthquake Resistant Design of

Unbraced Frames," Journal of the Structural Division, ASCE, Vol.
98, No. STll, November, 1972, pp. 2523-2539.

3. Aoyama, H. and T. Sugano, "A Generalized Inelastic Analysis of

Reinforced ,Concrete Structures Based on the Tests of Members,"
Recent Researches of Structural Mechanics, Uno Shoten, Tokyo,
Japan, 1968, pp. 15-30.

4. Bathe, K. J., and E. L. Wilson," Linear and Nonlinear Earthquake

Analysis of Complex Structures," Proceedings, Fifth World Conference
on Earthquake Engineering, Rome, Italy, June, 1973, Session 5B,
Paper No. 224.

5. Berg, G. V. and D. A. DaDeppo, "Dynamic Analysis of Elasto-Plastic

Structures," Journal of the Engineering Mechanics Division, ASCE,
Vol. 86, No. EM2, April 1960, pp.35-58.

6. Berg, G. V., "Response of Multi-Story Structures to Earthquake,"

Journal of the Engineering Mechanics Division, ASCE, Vol. 87,
No. EM2, April 1961, pp. 1-16.

7. Bertero, V. and B. Bresler, "Seismic Behavior of Reinforced Concrete

Frame Structures," Proceedings, Fourth World Conference on Earthquake
Engineering, Santiago, Chile, January, 1969, Vol. I, Session B2,
pp. 109-124.

8. Bleich, H. H., and M. G. Salvadori, "Impulsive Motion of Elasto-

Plastic Beams," Transactions, ASCE, Vol. 120, 1955, p. 499.

9. Clough, R. W., K. L. Benuska, and E. L. Wilson, "Inelastic Earthquake

Response of Tall Buildings," Proceedings, Third World Conference
on Earthquake Engineering, New Zealand, 1965, Vol. II, Session
II, ppm 68-89.

10. Clough, R. W., K. L. Benuska, and T. Y. Lin & Associates, FHA

Study of "Seismic Design Criteria for High-Rise Building~(HUD
TS-3), Federal Housing Administration, Washington, D.C., August
1966.
 54

11. Clough, R. W. and S. B. Johnston, "Effect of Stiffness Degradation

on Earthquake Ductility Requirements, II Proceedings, Japan
Earthquake Engineering Symposium, Tokyo, October 1966, pp. 227-32.

12. Clough, R., and K. L. Benuska, "Nonlinear Earthquake Behavior of

Tall Buildings," Journal of the Engineering Mechanics Division,
ASCE, Vol. 93, No. EM3, June, 1967, pp. 129-146.

13. DiMaggio, F. L., "Dynamic Elasto-Plastic Response of Rigid Frames,"

Journal of the Engineering Mechanics, ASCE, Vol. 84, No. EM3,
Proc. Paper 1693, July, 1958.

14. Giberson, M. F., "The Response of Nonlinear Multi-Story Structures

Subjected to Earthquake Excitation," thesis submitted to the
California Institute of Technology at Pasadena in partial fulfill-
ment of the requirements of the degree of Doctor of Philosophy,
1967.

15. Gael, S. C., "Inelastic Behavior of Multistory Building Frames

Subjected to Earthquake Motions," thesis submitted to the Univer-
sity of Michigan at Ann Arbor in partial fulfillment of the require-
ments of the degree of Doctor of Philosophy, 1967.

16. Gael, S. C., "P-6 and Axial Column Deformation in Aseismic Frames,"
Journal of the Structural Division, ASCE, Vol. 95, No. ST8,
August, 1969, pp. 1693-1712.

17. Gulkan, P. and M. A. Sozen, "Response and Energy-Dissipation of

Reinforced Concrete Frames Subjected to Strong Base Motions,"
Civil Engineering Studies, Structural Research Series No. 377,
University of Illinois, Urbana, May 1971.

18. Gurn, B. P., and A. C. Heidebrecht, "Factors Influencing the

Inelastic Response of Multi-story Frames Subjected to Strong
Motion Earthquakes," Proceedings, Fourth World Conference on
Earthquake Engineering, Santiago, Chile, January 1969, Session
A4, pp. 1-14.

19. Hays, C. 0., and H. Matlock, "Nonlinear Discrete Element Analysis

of Frames," Journal of the Structural Division, ASCE, Vol. 99,
No. STIO, October, 1973, pp. 2011-2030.

20. Heidebrecht, A. C., J. F. Fleming, and S. L. Lee, "Dynamic Analysis

of Inelastic Multi-Degree Systems," Journal of the Engineering
Mechanics Division, ASCE, Vol. 89, No. EM6, December, 1963,
pp. 193-215.

21. Heidebrecht, A. C., S. L. Lee, and J. F. Fleming, "Dynamic

Analysis of Elasto-Plastic Frames," Journal of the Structural
Division, ASCE, Vol. 90, No. ST2, April, 1964, pp. 315-343.
 55

22. Imbeault, F. A. and N. N. Nielsen, "Effect of Degrading Stiffness

on the Response of Multistory Frames Subjected to Earthquakes,"
Proceedings, Fifth World Conference on Earthquake Engineering,
Rome, Italy, June, 1973, Session 4B, Paper No. 165.

23. Ismail, M. A. F. and J. o. Jirsa, "Bond Deterioration in Reinforced

Concrete Subject to Low Cycle Loads," Journal, ACI, Vol. 69,
No.6, June, 1972, pp. 334-343.

24. Ismail, M. A. F. and J. o. Jirsa, "Behavior of Anchored Bars under

Low Cycle Overloads Producing Inelastic Strains," Journal, ACI,
Vol. 69, No.7, July, 1972, pp. 433-438.

25. Jennings, P. C., "Response of Simple Yielding Structures to Earthquake

Excitations," thesis submitted to the California Institute of
Technology at Pasadena in partial fulfillment of the requirements
of the degree of Doctor of Philosophy, 1963.

26. Jennings, P. C. and R. Husid, "Collapse of Yielding Structures during

Earthquakes," Journal of the Engineering Mechanics Division, ASCE,
Vol. 94, No. EM5, October 1968, pp.

27. Kroenke, W. C., M. J. Gutzwiller, and R. H. Lee, "Finite Element

for Reinforced Concrete Frame Study," Journal of the Structural
Division, ASCE, Vol. 99, No. ST7, July, 1973, pp. 1371-1390.

28. Lionberger, S. R. and W. Wiever, Jr., "Dynamic Response of Frames

with Nonrigid Connections," Journal of the Engineering Mechanics
Division, ASCE, Vol. 95, No. EMl, February 1969, pp. 95-114.

29. Meyer, I. L., J. B. Hoerner, and M. Zayed, "Inelastic Dynamic Analysis

of a 60 Story Building," Proceedings, Fifth World Conference on
Earthquake Engineering, Rome, Italy, June, 1973, Session 4B,
Paper No. 165.

30. Morita, S., and T. Kaku, "Local Bond Stress-Slip Relationship under
Repeated Loading," A paper presented in RILEM, Theme IV, 1973.

31. Newmark, N. M., "A Method of Computation for Structural Dynamics,"

Journal of the Engineering Mechanics Division, ASCE, Vol. 85,
No. EM3, July 1959, pp. 67-94.

32. Otani, S., and M. A. Sozen, "Behavior of Multistory Reinforced

Concrete Frames during Earthquakes," Civil Engineering Studies,
Structural Research Sereis No. 392, University of Illinois at
Urbana-Champaign, November, 1972.

33. Otani, S., "Inelastic Analysis of Ric Frame Structures," Journal

of the Structural Division, ASCE, Vol. 100, No. ST7, July,
1974, pp. 1433-1449.
34. Padilla-Mora, R., W. C. Schnobrich, "Non-linear Response of Framed
Structures to Two-Dimensional Earthquake Motion," Civil
Engineering Studies, Structural Research Series No. 408,
University of Illinois at Urbana-Champaign, Urbana, Illinois,
July, 1974.

35. Pekau, o. A., and R. Green, "Analysis and Collapse Behavior of

Inelastic Frame Structures," Computer & Structures, Vol. 4,
pp. 337-362, Pergamon Press 1974.

36. Penzien, J., "Dynamic Response of Elasto-Plastic Frames," Journal

of the Structural Division, ASCE, Vol. 86, No. ST7, July 1960,
pp. 81-94.

37. Penzien, J., "Elasto-Plastic Response of Idealized Multistory Struc-

tures Subjected to a Strong Motion Earthquake," Proceedings,
Second World Conference on Earthquake Engineering, Tokyo and
Kyoto, Japan, 1960, Vol. II, Session II, pp. 739-60.

38. Powell, G. H., and A. Kanaan, "General Purpose Computer Program

for Inelastic Dynamic Response of Plane Structure," Proceedings,
Fifth World Conference on Earthquake Engineering, Rome, Italy,
June, 1973, Session 6B, Paper No. 272.

39. Sheth, R. M., "Effect of Inelastic Action on the Response of Simple

Structures to Earthquake Motions," M. S. Thesis, University of
Illinois, Urbana, 1959.

40. Spencer, R. A., "The Nonlinear Response of a Multistory Prestressed

Concrete Structure to Earthquake Excitation," Proceedings, Fourth
World Conference on Earthquake Engineering, Santiago, Chile,
January 1969, Session A4, pp. 139-154.

41. Suko, M. and P. F. Adams, "Dynamic Analysis of Multibay Multistory

Frames," Journal of the Structural Division, ASCE, Vol. 97,
No. STIO, October 1971, pp. 2519-33.

42. Sun, C-K., G. V. Berg, R. D. Hanson, "Gravity Effect on Single-

Degree Inelastic System," Journal of the Engineering Mechanics
Division, ASCE, Vol. 99, No. EMl, February, 1973, pp. 183-200.

43. Takeda, T., "Study of the Load-Deflection Characteristics of

Reinforced Concrete Beams Subjected to Alternating Loads," Trans-
actions, Architectural Institute of Japan, Vol. 76, September
1962, p. 98.

44. Takeda, T., M. A. Sozen, and N. N. Nielsen, "Reinforced Concrete

Response to Simulated Earthquakes," Journal of the Structural
Division, ASCE, Vol. 96, No. ST12, December 1970, pp. 2557-73.

45. Tanabashi, R., "Studies on the Nonlinear Vibrations of Structures

Subjected to Destructive Earthquakes," Proceedings, World Con-
ference on Earthquake Engineering, Berkeley, California, June
1956, Paper No.6.
 57

46. Umemura, H., H. Aoyama, and H. Takizawa, "Analysis of the Behavior

of Reinforced Concrete Structure during Strong Earthquakes Based
on Empirical Estimation of Inelastic Restoring Force Character-
istic of Members," Proceedings, Fifth World Conference on
Earthquake Engineering" Rome, Italy, June, 1973, Session 6B,
Paper No. 275.

47. Walpole, W. R., and R. Shepherd, "Elasto-Plastic Seismic Response

of Reinforced Concrete Frame," Journal of the Structural Division,
ASCE, Vol. 95, No. STIO, October, 1969, pp. 2031-2055.

48. Walpole, W. R., and R. Shepherd, "The Inelastic Response of a

Steel Frame," Proceedings, Fourth World Conference on Earthquake
Engineering, Santiago, Chile, January 1969, Session A4, pp. 195-
204.

49. Wen, R. K., and J. G. Janssen, "Dynamic Analysis of Elasto-Ineastic

Frames," P::t;'oceedings, Third World Conference on Earthquake
Engineering, New Zealand, January, 1965, Vol. II,pp. 713.;...729.

50. Yoshioka, K., T,. Takeda, and K. Nakagawa, "Inelastic Earthquake

Response of Reinforced Concrete Buildings," Proceedings,
Fifth World Conference on Earthquake Engineering, Rome, Italy,
June, 1973, Session 6B, Paper No. 271.
 58

W6

W5

W
4

W3

tJ
"2

WI

Fig. 2.1 Idealized Structure

/ Joint Rotational Spring

~ ~ Inelastic Flexural Spring

Elastic Prismatic Element

Fig. 2.2 Idealized Member

 59

Frame 1

Frame 2

Frame 3

Frame 4

(a) Plan View

~)

Frames 2 and 3 Frames 1 and 4

(b) Elevation

Fig. 3.1 Example of Structure

 M3

M2

Ml

Frame 1 Frame 4 Frame 2 Frame 3

(a) Structural Model

0\
o
1/2 M3

1/2 M2

1/2 Ml

Frame 1 Frame 2

(b) Equivalent Model

Fig. 3.2 Analytical Model of Structure

 61

(a) Irregular Column

(b) Irregular Beams

Fig. 3.3 Irregular Frames

Fig. 3.4 Hinged Joint

 (a) Original Flexible Member

(b) One-half Flexible Member

(c) Moment Distribution

Inelastic Rotation

(d) Curvature Distribution

Fig. 4.1 Moment-Rotation Relation of Flexural Spring

 63

...

Fig. 4.2 Moment Distribution in Shear Wall Member

(a) Exterior Beam (b) First Story Column

Fig. 4.3 Moment Distribution in Frame Member

 64

Start

Nonlinear Member Stiffness

(Linearization by Small Time Increment)

Instantaneous Linear
Member Stiffness

(Compatibility)

Instantaneous Linear
Structural Stiffness Mass Inertia
Force
Damping Force
(Equilibrium) External Force

Equation of Motion

(Nonlinear)

(Numerical Integration Assumption)

Linear Equation about Incre-

mental Joint Displacement

(Solution of Linear Equation)

Incremental Member End

Displacement

(Instantaneous Member Stiffness)

~ ____~:T_ncremental Member End

Force

Fig. 5.1 Flow Chart of Nonlinear Structural Analysis

Fifth Story Height
L
1
Fourth Story Height

Third Story Height

Second Story Height

First Story Height

F.irst Bay Second Bay Third Bay'

Fig. 6.1 Structural Definition

&
IV
I'

,- r-- - 1""-

II I
g~
Joint Core
I I
n
J I ~~ I I
i
I

I '-1 I j
,

I I I
I I I
I
,- b- I I ,- ..... _-

BA AL BB
p

Fig. 6.2 Member Definition

 66

Fig. 6.3 Flexural Properties

U
M
u

M
y

+.J
~
(lJ
g
~
M
c

Curvature

Fig. 6.4 Idealized Moment-Curvature Diagram

 VU

Idealized
+> Curve
~
Q)
S
0
~

Vy

WY wu
Rotation

Fig. 6.5 Idealized Backbone Curve for Joint Spring

1 2 3

10 11 12 13

4 5 6

14 15 16 17

7 8 9

18 19 20 21

Fig. 6.6 Member Numbering

 68

1 2 3 4

6 1 8
5

10 11 12
9

Fig. 1.1 Joint Numbering

A-end

A-end B-end

(b) Beam

B-end

(a) Column

Fig. 1.2 Member End Notation

 69

Unit-Length

Rotation

Fig. 7.3 Rotation of Flexural Spring

(b) Beam

(a) Column

Gif. 7.4 Sign Convention of Moment

 70

APPENDIX A. EVALUATION OF MEMBER PROPERTIES

This appendix describes methods of determining the stiffness constants

of flexural elements (A.I) and joint rotational springs (A.3). A program

user is not restricted by the methods described here in determining

input constants. However, the moment-rotation relationship of a flexural

element (A.2) will be calculated in the program based on the moment-

curvature relationship provided by the user.

A.I Flexural Properties

The primary moment-curvature curve, which is defined as a moment-

curvature curve for a constantly increasing load at a given axial load,

provides a good index to understand the behavior of a section. The

flexural behavior of a member can be determined by the moment-curvature

relationship.

The numerical values of bending moments and curvatures can be computed

from the geometry of the section, the amount of existing axial load, the

properties of concrete and reinforcing steel, on the basis of Bernoulli's

Hypothesis, which assumes a linear strain distribution across the depth

of the section. The primary moment-curvature curve may be idealized by

a piecewise linear relationship (Fig. 6.4) with break-points represented

by moments and curvatures at tensile cracking of the concrete, yielding

of the tensile reinforcement, and crushing of the concrete at the extreme

compressive fiber.

Numerical values of moments and curvatures at cracking, yielding

and ultimate conditions can be computed as described in various text

 71

books on reinforced concrete design (for example, Winter and Nilson,

Ferguson, others), and are used as input information to the program.

A.2 Properties of Flexural Elements

The moment-rotation relation of a flexural spring element was cal-

culated in the program based on the idealized moment-curvature curve.

A flexible part of a member was assumed to have uniform section

properties from an end to the center.

The rotation at an end of a flexural element depends on a moment

distribution. A linear moment distribution was assumed with the point

of contraflexure to be at the mid-point. In this cas~, the moment-

rotation relationship at an end of a flexible element becomes equivalent

to that of a one-half length simply supported member (Fig. 4.1). The

numerical values of moments and rotations at cracking, yielding and

ultimate conditions were expressed as follows:

(1) Cracking of tensile concrete:

8 =.Q,.M 13EI (A.l)

c c

in which

EI = initial elastic flexural rigidity,

1 = length of a simply supported member,

M = moment at cracking of tensile concrete.

c

( 2) Yielding of tensile reinforc ement :

8 = ,Q,[(1_,,3) + ,,2 ]/3 (A.2)

y y c
in which

" = Mc 1My
 72

c = curvature at at cracking of tensile concrete,

y = curvature at yielding of tensile reinforcement, and

M = moment at yielding of tensile reinforcement.

y
(3) Crushing of compressive concrete:

in which

Al = Mc 1Mu

A2 = My 1Mu

M - M
a =( u ~) (~)
u y M
y

M
u
= moment at crushing of compressive concrete.

The moment-rotation relationship was idealized by a piece-wise

linear curve with break points (8 , M ), (8 , M ) and (8 , M ).

c c y y u u
It should be noted that the rotations above are proportional to

the length of a simply supported member. Therefore, it is sufficient

to prepare the moment-rotation relationships only for unit-length elements

with a typical section properties (flexural elements).

Takeda hysteresis rules (Appendix B.l) were used in the program to

define the moment-rotation relationship of unit-length flexural elements

(including elastic deformation) under load reversals. The deformation

of a unit-length flexural element must be multiplied by the length of a

simply supported member in order to get a realistic value. By the same

token, the stiffness of a unit-length flexural element must be divided

by the length of a simply supported member in order to formulate the

stiffness of a frame member.

 73

The stiffness of an inelastic flexural rotational spring was defined

in the program by subtracting the contribution of elastic stiffness from

the stiffness of a half-length simply supported flexural element in a

following manner:

I I I
K =K Ke
(A.4)
sp t

in which

K
sp = stiffness of an inelastic flexural spring element,

K
t = stiffness of a half-length simply supported flexural

element including elastic and inelastic d'eformation,

and

K
e
= stiffness of an elastic half-length simply supported

member (3EI/Q,).

A.3 Properties of Joint Springs

In addition to the flexural deformation of members in a structure,

rotation due to the slip of the tensile reinforcement along its embedded

length was considered at the ends of a member, and was simulated by a

rotation of a joint (rotational) spring.

The moment-rotation relationship of a joint spring depends on various

parameters such as anchorage length of tensile and compressive reinforce-

ment, bond stress distribution and bond strength between the steel and

concrete, amount of reinforcement in the joint core, and others.

The following simple-minded method may be used in evaluating the

stiffness of a joint rotational spring. However, a program user should

not be restricted to the following method. More elaborate methods may

be necessary to evaluate a realistic behavior of bond slip in a joint core.

 74

Let us consider the moment-rotation relationship of a joint spring

under constantly increasing bending moment, assuming the failure will

not occur in the anchorage zone.

If the tensile reinforcement has sufficient development length in

a joint core, and if the bond stress is assumed to be constant along the

development length, the elongation of the tensile reinforcement can be

expressed as

e = ~ . f~ / 8 Es u

in which

e = elongation of the tensile reinforcement

d .= diameter of the tensile reinforcement

b
E = modulus of elasticity of steel
s
u = average bond stress

f
s
= stress in the tensile reinforcement at the joint spring.

If the compressive reinforcement does not slip, the rotation due to the

slip can be evaluated as

(A.6)

in which

d = depth of the tensile reinforcement

d~ = depth of the compressive reinforcement

If a stress in the tensile reinforcement is assumed to be proportional

to the moment, in other words,

(A.7)
 75

then

8 (A.8)

This parabotic moment-rotation relationship may be idealized by

a bilinear relationship (Fig. 6.5) to determine a break-point (VY and

WY), and an ultimate point (VU and WU).

 APPENDIX B. HYSTERESIS RULES FOR FLEXURAL
AND JOINT SPRINGS

Assumed stiffness characteristics have an important role in

determining an inelastic dynamic response of an analytical model.

Experimental observation of reinforced concrete members under load re-

versals has revealed that a bilinear hysteresis expression is not

representative nor realistic. Reinforced concrete members tend to

show stiffness degradation as a function of stress history.

Two hysteretic systems were adopted from Takeda et al (1970) in

this program: (1) the original Takeda hysteresis system (Fig. B.l)

with a tri-linear backbone curve to simulate inelastic flexural behavior,

and (2) a simplified version of the Takeda hysteresis systmfr (Fig. B.2)

with a bilinear primary curve to simulate inelastic behavior of a joint

spring.

The following definitions and notations were used to simplify the

d~scription of the hysteretic rules.

loading = an increase, without change in sign, of the

absolute value of force

unloading = a decrease, without change in sign, of the

absolute value of force

load reversal= change in sign of force with respect to the

one in the last load step

K = stiffness to be used in the next load increment

o = origin of the primary curve

C = cracking point on the primary curve

Y = yielding point on the primary curve

 77

U = ultimate point on the primary curve

P = point at which current calculation begins

U.
l
= point at which unloading begins

X.
l
= load reversal point or intersection point of

the force-deflection curve with the displacement

axis.

S.
l
= slope
F(A) = force at point A on the force-deflection curve

n(A) = deflection at point A on the force-deflection

curve

SCAB) = slope of line segment AB

At each load step, the force active is assumed to be positive in describ":"'

ing the hysteresis rules. An apostrophe is used to indicate the

corresponding point in negative region of the force~deflection diagram.

In the following section, the rules of the hysteresis system are

described in terms of the definition given above. An example is

provided below to help in the interpretation of the rules.

Example

Rule 3: loading on the primary curve after yielding

3.1 loading: K = S(YU), go to rule 3.

3.2 unloading: unloading point = Um ,

K = Sl' go to rule 4.
 78

Interpretation

Rule 3 governs if loading starts on the primary curve after yield-

ing has occurred.

If the absolute value of load continues to increase in the same

direction (3.1 loading), the path of the force-deflection curve follows

the primary curve for that increment with the stiffness (K) defined by
the slope of segment YU (slope defined by yield point Y and ultimate

point U). For the next increment, rule 3 continues to govern.

If the absolute value of load is decreased (3.2 unlaoding) from the

point U , the unloading slope is defined by the product of the slope

m

(C'Y) and the ratio of the yield deflection to the maximum deflection

reached in either direction. Calculation for that increment uses the

slope Slo The next increment is governed by rule 4.

B.l Takeda's Hysteresis System

Moment-rotation relationship of a unit-length simply supported

flexural element (Appendix A.2) was defined by a set of Takeda's

hysteresis rules. A piece-wise linear primary (backbone) curve was

defined by two break points at cracking of tensile concrete and yielding

of longitudinal reinforcement. The primary curve was assumed to be

symmetric about its origin.

The hysteresis rules were developed on the basis of experimental

observation of reinforced concrete members (beams and columns) under

load reversals. The hysteresis rules may simulate the flexural behavior

of reinforced concrete members (beams and columns) with light to medium

amount of axial loads. The failure mode of a member must be dominantly

flexure. Failure or extensive damage caused by shear or bond deterioration

was not considered in the hysteresis rules.

 79

The following set of rules were used to determine the stiffness

of a frame member at each load level (Fig. B.l).

Rule 1: elastic stage

1.1 loading

1.1.1 F(P) ~ F(C) K:::: S(OC), go to rule 1.

1.1.2 F(P) > F(C) K :::: S(CY), go to rule 2.

1.2 unloading and load reversal K:::: S(OC), go to rule 1.

Rule 2: loading on the primary curve up to yielding

2.1 loading

2.1.1 F(P) .$ F(Y) K :::: S(CY), go to rule 2.

2.1.2 F(P) > F(Y) K :::: S(YU), go to rule 3.

2.2 unloading : unloading point :::: U , S. :::: S(PC'), K ::::

Sl' go
m 1
to rule 5.

Rule 3: loading on the primary curve after yielding

3.1 loading: K:::: sCYU), go to rule 3.

3.2 unloading: unloading point:::: U ,

m

S :::: S(C'Y) *1.D(Y)/max{D(U ) ,D(U' )}

1 mm

K :::: Sl' go to rule 4.

Rule 4: unloading from point U on the primary curve after yielding

m
4.1 loading

4.1.1 F(P) ~ F(U ) K :::: Sl' go to rule 4.

m

4.1a2 F(P) > FeU ) K :::: S (YU) .~. go to rule 3,

m

4.2 unloading: K :::: Sl' go to rule 4.

 80

4.3 load reversal

4.3.1 uncracked in the negative range

K = Sl' go to rule 15.

4.3.2 otherwise: load reversal point =

= S(X u ),
o m

K = S2' go to rule 6.

Exception in 4.3.2

If F(Y) > F(U ) and S(X Y) > S(X U )

moo m
then S2 = S(Xo Y), Um ,= Y, K = 8 , go to rule 6.
2

Rule 5: ~nloading from point U on the primary curve before yielding

m

5.1 loading

5.1.1 F(P) ~ F(U ) K = Sl' go to rule 5.

m
5.1.2 F(P) > F(U ) K = S(CY), go to rule 2.
m
5.2 unloading: K = Sl' go to rule 5.

5.3 load reversal

5.3.1 uncracked in the negative range

go to rule 14.

5.3.2 otherwise the same as 4.3.2.

Rule 6: loading toward point U on the primary curve

m
6.1 loading

6.1.1 F(P) ~ F(U )

m
K = 8(XoUm), go to rule 6.

6.1.2 F(P) > F(U ) : the same as 2.1.

m

6.2 unloading: unloading point = Uo' K = Sl' go to rule 7.

 81

Rule 7: unloading from point U after rule 6.

o
7.1 loading

7.1.1 F(P) ~ F(U )

o K = Sl' go to rule 7.

7.1.2 F(P) > F(U )

o
K= s(x oUm ), go to rule 6.

7.2 unloading: K = Sl' go to rule 7.

7.3 load reversal load reversal point = Xl~ ~ = 8(XI Um),

go to rule 8.

Rule 8: loading toward point U on the primary curve

m
8.1 loading

8.1.1 F(P) ~ F(U )

m

8.1.2 F(P) > F(U ) the same as 2.1

m

8.2 unloading: unloading point = Ul' K = 81 , go to rule 9.

Rule 9: unloading from point U after rule 8

l
9.1 loading

9.1.1 F(P) ~ F(U

l
) K = 81 , go to rule 9.

9.1.2 F(P) > F(U ) the same as 8.1

l
9.2 unloading: K = 81 , go to rule 9.

9.3 load reversal load reversal point = X2 ' K = S(X 2Uo ),

go to rule 10.

Rule lO:loading toward point U

o
10.1 loading

10.1.1 F(P) ~ F(U )

o
K = 8(X 2Uo )' go to rule 10.

10.1.2 F(P) > F(U ) the same as 6.1

o
10.2 unloading unloading point = U2 ' K = 81 , go to rule 11.

Rule 11: unloading from point U after rule 10

2
11.1 loading
 82

11.1.1 F(P) ~ F(U )

2
K = Sl' go to rule 11.

11.1.2 F(P) > F(U ) the same as 10.1

2
11.2 unloading : K = Sl' 'go to rule 11.

11.3 load reversal load reversal point = X3 ' K = S(X3Ul ),

go to rule 12.

Rule 12: loading toward point U

l

12.1 loading

12.1.1 F(P) ~ F(U )

l
K = S(X3 Ul ), go to rule 12.

12.1.2 F(P) > F(U ) the same as 8.1 ..

l
12.2 unloading: unloading point = U3 , K = Sl' go to rule 13.

Rule 13: .unloading from point U after rule 12

3
13.1 loading

13.1.1 F(P) ~ F(U )

3
K = Sl' go to rule 13.

13.1.2 F(P) > F(U ) the same as 12.1.

3
13.2 unloading: K = Sl' go to rule 13.

13.3 load reversal : the same as 9.3 ..

Rule 14: loading in the uncracked direction after cracking in the other

direction

14.1 loading

14.1.1 F(P) ~ F(C) K = Sl' go to rule 14.

14.1.2 F(P) > F( C) K = S(CY), go to rule 2.

14.2 unloading : K = Sl' go to rule 14.

14.3 load reversal K = Sl' go to rule 5

 83

Rule 15: loading in the uncracked direction after yielding in the

other direction

15.1 loading

15.1.1 F(P) ~ F(C) K = Sl' go to rule 15.

15.1.2 F(P) > F(C) let Q be point on the curve for

rule 15 and F(Q) = F(C), then

K = S(QY), go to rule 16.

15.2 unloading: K = Sl' go to rule 15.

15.3 load reversal

15;3.1 F(P) ~ F(U )

m
K = Sl' go to rule 4.
15.3.2 F(P) > F(U )
m
K = S(YU), go to rule 3.

Rule 16: loading toward the yield point after rule 15

16.1 loading

16.1.1F(P) ~ F(Y) K = S(QY), go to rule 16.

16.1.2 F(P) > F(Y) K = S(YU), go to rule 3.

16.2 unloading: U
m
= Y, Uo = P, S2 = S(QY), K = Sl'
intersection of YQ and deflection

axis = Xo , go to rule 7.
B.2 Simplified Takeda Hysteresis System

Behavior of the longitudinal reinforcement in a joint core

becomes a popular target in reinforced concrete studies. However, our

knowledge has not advanced far enough to formulate the hysteresis behavior

of bond slip between concrete and rein~orcement.

Based on an assumption that a sound structure must be free of bond

deterioration within a joint core, Takeda's hysteresis system was simpli-

fied and used to simulate a member end rotation due to bond slip of the

reinforcement within a joint core.

 84

The simplified Takeda's hysteresis rules have a pilinear primary

curve, and are similar to the one proposed by Clough and Johnston

(1966). The simplified model includes more rules for small-amplitude

load reversals.

The following set of rules were used to determine the stiffness of

a rotational spring at each load step (Fig. B.2).

Rule 1: elastic stage

1.1 loading

1.1.1 F(P) ~ F(Y) K = 8(OY), go to rule 1.

1.1.2 F(P) > F(Y) K = 8(YU), go to rule 2.

1.2 unloading and load reversal: K = 8(OY), go to rule 1.

Rule 2: loading on the primary curve after yielding

2.1 loading : K = 8(YU), go to rule 2.

2.2 unloading: unloading point = Um , K = 8(OY), go to rule 3.

Rule 3: unloading from point U on the primary curve

m
3.1 loading

3.1.1 F(P) ~ F(Y) K = 8(OY), go to rule 3.

3.1.2 F(P) > F(Y) K = 8(YU), go to rule 2.

3.2 unloading: K = 8(OY), go to rule 3.

3.3 load reversal load reversal point = Xo , K = 8(X0 U*),

m
go to rule 4.
Rule 4: loading toward point U on the primary curve
m
4.1 loading

4.1.1 F(P) ~ F(U )

m
K = 8(XoUm), go to rule 4.

4.1.2 F(P) > F(U )

m
K = 8(YU) , go to rulE; 2.

4.2 unloading: unloading point = U0' K = S( OY), go to rule 5.

*If U in the new direction is not defined yet, U is taken as yielding

m m
point Y.
Rule 5: unloading from point U after rule
o
4
5.1 loading
5.1.1 F(P) ~ F(U')
o
K = 8(OY), go to rule 5.

5.1.2 F(P) > F(U o ) the same as 4.1.

5.2 unloading: K = 8(OY), go to rule 5

5.3 load reversal load reversal point = Xl' K = 8(Xl Um),

go to rule 6.

Rule 6: loading toward point U* on the primary curve

m

6.1 loading
6.1.1 F(P) ~ F(Um ) K = 8 (XlU ), go to rule 6.
m
6.1.2 F(P) > F(Um) K = 8(YU), go to rule 2.

6.2 unloading: unloading point = Ul ' K = 8(OY), go to rule 7.

Rule 7~ unloading from point U after rule 6

l
7.1 loading
7.1.1 F(P) ~ F(UI ) K = 8(OY) , go to rule 7.

7.1.2 F(P) > F(Ul ) the same as 6.1.

7.2 unloading: K = 8(OY), go to rule 7.

7.3 load reversal load reversal point = X2 , K = 8(X2Uo )'

go to rule 8.

Rule 8: loading toward point U

o
8.1 loading
8.1.1 F(P) ~ F(U o ) K = 8(X2Uo ), go to rule 8.

8.1.2 F(P) > F(U o ) the same as 4.1.

8.2 unloading: unloading point = U2 ' K = 8(OY), go to rule 9.

*U in load level 6 is of the sign opposite to the sign of U in load

m m
level 4.
 86

Rule 9: unloading from point U after rule 8

2
9.1 loading

91.1 F(P) ~ F(U )

2
K = S(OY), go to rule 9.

91.2 F(P) > F(U ) the same as 8.1.

2
9.2 unloading: K = S(OY), go to rule 9.

9.3 load reversal load reversal point = X3 ' K = S(X3Ul)'

go to rule 10.

Rule 10: loading toward point U

l

10.1 loading

10.1.1 F(P) ~ F(U )

l
K = S(X3Ul ), go to rule 10.

10.1.2 F(P) > F(U ) the same as 6.1

l
10.2 unloading: unloading point = U3 , K = S(OY), go to rule 11.

Rule 11: unloading from point U after rule 10

3
11.1 loading

11.1.1 F(P) ~ F(U )

3
K = S(OY), go to rule 11.

11.1.2 F(P) > F(U ) the same as 10.1

3
11.2 unloading: K = S(OY), go to rule 11

11.3 load reversal load reversal point = X , K

2
= S(X2Uo )'
go to rule 8.
 IS,(

CD Rule
Rule 2 Q)
Y
(,J

@ Rule 5
....
0
LL
@ Rule l.i.;

('~ Rule '1

L-

'G':
\t", Rule 3
0 Hule Lt,
Rule 6 Deflection

(D Rule
(%~i Hule 2
Urn
~) Rule 3 Q)
Q
Rule 4 ....
0
u...
rJ;) Rule is
() Rule 16
(f) Rule J
~

@) Rul~ 4
Rule 6 Oaf laction

(b)

.-----+"
u m
(!)

Takeda's Hysteresis Rules

 88

CD: 4
@: 6
@: 7
@): 8
@~ 9
@: 10
(!): 11
12
13
New

0: Um

Fig. B.1 (Contd) Takeda!s Hysteresis Rules

 89

CD: Rule
(g): Rule 2 C
(.)
@: Rule 3 r..
&
@: Rule 4
@: Rule 2
: New Rule 3
(i): New Rule 4
@: Iiule 5
@: Rule 6

CD: 3
@: Rule 4
(.)
-
@; Rule 5 0
LA-
@): Rule 6
@: Rule 7
@: Rule a
(1): Rule 9
@: Rule 10
@: 11
: New 8

(b)

Fig. B.2 Simpl ified Takeda1s Hysteresis Rules

 90

APPENDIX C. MEl\1BER STIFFNESS MATRIX

A member (horizontal beam or vertical column) consisted of a

flexible elastic line element, two nonlinear rotational springs each

at an end of the flexible elastic line element, and two rigid zones

outside of the rotational springs as shown in Fig. 2.2.

If the portion of a member between two rigid zones is simply sup-

ported, and external bending moments are applied at the supports

(Fig. C.l), then the rotation at support A~ is the sum of the rotations

of two rotational springs (instantaneous flexibility constants of fA

and jA ) and a rotation at the end of the flexible elastic line element

(flexural rigidity of EI and length ~)o Therefore, in an incremental

form,

(C.l.a)

By the same token

(C.l.b)

If the relation is written in a matrix form,

~eAl _ ~'/3EI + fA + jA -i'/6EI

(C.2)
~eBf - -~'/6EI ~'/3EI + fB + jB ~~
 91

This relation can be easily rewritten by inverting the flexibility

matrix into a form:

= (C.3)

With this information, let us consider a simply supported member

with two external bending moment applied at the ends (Fig. C.2).

Incremental moments ~A and ~MB at the ends of a member are related to

incremental moments ~MA and ~MB at the rotational springs by the

following matrix relation because of a linear moment distribution along

the member:

(c.4)

in which

BA : length ratio of a rigid zone at A-end to the flexible elastic

portion,

BE length ratio of a rigid zone at B-end to the flexible elastic

portion.

Incremental rotations ~8A and ~8B at the ends of a member and

, ,
those 68 and ~8B inside rigid zones are related in the following form:
A

I \ / \

68"
A

=
1 + BA BB ~8A
l (C.5)
6 e" BA 1 + BB ~eB
B J
 92

By combining Eqs. C.3, c.4 and C.5, a stiffness relation of a

member can be formulated for member end moments and rotations in a

following form:

= (c.6)

For a beam member, in which a member axis is always horizontal

even after deformation and coincides with the horizontal global coordinate

axis, no further transformation is necessary to express local member

end rotations and moments in the global coordinates. However, a member

axis of a column is rotated by lateral column end displacements (Fig. C.3).

Some transformation is required to express global member end (joint)

displacements in the local member coordinate. Transformation matrix [C]

can be expressed as

[c] = l: o

1
-l/~

-l/~

in which ~ = total length of a column.

The member end rotations as a simply supported member can be

expressed by lateral displacements U and UB' and rotations G and

A A
G at structural joints A and B in the global coordinates, using
B
transformation matrix [C],

~
A
~e
B
= [C] (c.B)
U
~eB A
U
B
 93

Magnitudes of bending moment at the ends of a column in the local

coordinates are the same in the global coordinates. Lateral forces

P and P at the ends of a colUmn can be calculated from bending moments

A B
MA and MB" In an incremental form, the transformation of member end

forces in the local member coordinates to those in the global structural

coordinates can be written in the following form:

~1
~B
- [C]T (C.9)
~PA ~~

~PB

in which [C]T is transform of matrix [C].

If Eqs. c.6, c.8 and C.9 are combined, a 4 x 4 member stiffness

matrix of a column in the global coordinates can be formulated, and

the force-displacement relation will be in an incremental form

~ kll kl2 k
l3 k14 ~8
A
~B k2l k22 k23 k24 ~8
B

~PA
= -----T----- ~UA
(C.IO)
k3l k32 k33 k34
I
~PB k4l k42 I k43 k44 ~UB

A two-by-two submatrix in the upper left corner of a column stiffness

matrix is of the same form as a beam member stiffness matrix in Eq. c.6.
The program considers the cases when one or both joint springs act

as mechanical hinges. The elements of the stiffness matrix in Eq. C.3

can be expressed as follows.

 94

(i) Fixed-hinged member (M - 0)

B

(C.ll)

(ii) Hinged-fixed member (M - 0)

A

k ...
11 = k12 = k21 = 0
(C.12)
k22 = 1/ ('Q,'" /3EI + fB + jB)

(iii) Hinged-Hinged member (M

A = MB 0)

kll = k12 = k21 = k22 = 0 (C.13)

The treatment will be the same in Eq. C.3 through Eq. C.9.
 95

Fig. C.l Moment-Rotation Relation of

Flexible Portion of Member

s -"
B

Fig. C.2 Boment-Rotation Relation

of Member
 96

A B
8
B =
eB

VA

I
I (~)

I ~A

I eA
I
I
I
I
I
I
I eB
II 8
B
I
I U
B

Fig. C.3 Displacements at Joints and Member Ends

 97

APPENDIX D. STRUCTURAL STIFFNESS MATRIX

From the assumptions stated in Chapter 2, bending moments and

horizontal forces were considered as generalized forces of a structure

at each joint, and rotations and horizontal displacements as generalized

displacements. The global coordinate system consists of a horizontal

axis parallel to the base line of a frame, and a clockwise rotational

axis (Fig. D.l). A structural stiffness matrix was formulated for each

constituent frame separately.

Force-deflection relations of members in the global coordinates were

formulated in Appendix C. In order to construct force-deflection relations

at all structural joints in a frame, two assumptions were used:

(l) resultant forces acting in a joint must be equal to the sum of

all forces at member ends around the joint (equilibrium condition), and

(2) member end displacements aro~nd a joint must be equal to the

joint displacements (compatibility condition).

If all member end forces at a joint are added and expressed in terms

of member end displacements which can be replaced by structural joint

displacements on the basis of the compatibility condition, force-deflection

relations at structural joints are formulated for an entire frame (Fig.

D-2)

The joint force-deflection relations can be further reduced to a more

compact form by an assumption that all joints in a floor level displace

the same amount in the horizontal direction. A lateral force was computed
 98
as the sum of horizontal forces at all joints at a floor level, and

horizontal displacements of the joints were replaced by a lateral

displacement of the floor level.

Therefore, force-deflection relations of a frame can be expressed

in a matrix form,

(_~!L
t ~8
_j
in which

{~p} = vector which contains incremental horizontal forces

at each floor level,

{~} = vector which contains incremental bending moments at

each joint,

{~U} = vector which contains incremental horizontal

displacements at each floor level,

{~8} = vector which contains incremental rotation at each

joint,

[K ] = symmetric submatrix of size, number of stories by

ll
number of stories,

[K
21
] = submatrix of size, number of joints by number of

stories,

[K ]
12
= submatrix of size, number of stories by number of

joints, and transpose of matrix [K ],

21
[K
22
] = symmetric banded submatrix of size, number of joints

by number of joints.

As a part of an effort to reduce the computer program size, symmetric

matrix [K ] was stored only the upper triangle portion of the original
ll
matrix (Fig. D.3). Symmetric banded matrix [K ] was stored only the
22
 99

upper banded portion of the original matrix. Banded matrix [K ] was

21
stored only the banded part of the original matrix.

Equation D.I can be expanded using the submatrices,

(D.2)

As rotational inertia was not considered in the program, bending

moment at any joint must be zero. However, unbalanced forces {6M}

and {6P} were created at a' structural joint by the correction of

member end forces to fit a specified force-deflection ,relation at the

end of the previous time interval of response calculation.. If unbalanced

forces {6P} and moments {6M} were included, Eq. D.2 was expressed as

(D. 3)

Solving the second equation of Eq. D. 3 for 68 ,

(D. 4)

Substituting Eq. D.4 into the first equation of Eq. D.3,

in which

(D.6)

Matrix [K] is sometimes called a "reduced" structural stiffness matrix,

or simply a structural matrix.

 100

If a structure consists of more than one frame, the structural stiff-

ness of the entire structure can be expressed as the sum of reduced

structural stiffness matrices of component frames.

Incremental displacements {~U} were computed as a numerical solution

of the equation of motion (Appendix G). Incremental joint rotations

{~8} were calculated by Eq. D.4 from known {~U} and {iM}.
 101

P,U

Fig. D.l Global Coordinates

M1 + M2 + M3 + M4 =0

Fig. D.2 Equilibrium and Compatibility at a Joint

 102

Story Number a.
In

a
(Symmetric) nn

n = number of stories

Story Number
~

level i
r
,

~ I
/
t
joints at
,I level (i+
(m+l)(n-l-i) + 1

I
_.J_ joints at
I
I I level (i)
I I
__ L_'
I joints at
I level (i-l) (m+l)(n+2-i)
I

m = number of bays

Fig. D.3 Storage of Stiffness Matrix

 103

Main
Joint Number Diagonal

(Symmetric)

Band Width =m + 2

1 n 2 (m+l)

m+ 2 2(m+l)

[.,. ., . (m+l)(n+l-i)
I~

-=<

tl

.1

1 m

(d)

Fig. D.3 (Cont'd) Storage of Stiffness Matrix

 104

APPENDIX E. MASS AND DAMPING MATRICES

Weight of a floor level of an entire structure was assumed to be

lumped and attached to the floor level by a truss element as shown in

Fig. 2.1. Therefore, the weight does not create any stress in frame

members under a static condition. The floor weight was used to define

a horizontal mass inertia and the P-6 effect.

Displacements in the vertical direction was not considered in this

program. Hence, mass was not considered in the vertical direction.

As a member in this analysis was assumed to be a massless line element,

a rotational inertia was not assigned to a structural joint.

A mass matrix of an entire structure was reduced to a diagonal

matrix of size, number of stories by number of stories. Only diagonal

elements of the structural mass matrix were stored in the program.

The word "damping" has been used in structural dynamics with damping

recognized as a resisting force as well as a source of energy dissipation.

The mechanism of damping has not been clearly understood quantitatively.

A viscous type damping has been adopted in structural dynamics probably

because of its mathematical simplicity in the expression in addition to

its resisting and energy dissipating capability.

If the mechanism and source of damping are clearly known and quantita-

tively definable, damping should be considered at the member (or material)

level in a structural analysis. However, a precise quantitative value

 105

of damping is seldom available. Therefore, damping was considered at the

structural level in this program.

The expression should be simple and easy to solve because the

complicated mathematical model is not justifiable due to the unknown

mechanism and magnitude of damping. Consequently, viscous type damping

was adopted, with a damping matrix partially proportional to a constant

structural mass matrix, and partially proportional to a variable

instantaneous structural stiffness matrix.

[C] ::::: a[M] + (3 [K]

[C] ::::: instantaneous damping matrix,

[M] ::::: mass matrix,

[K] ::::: instantaneous stiffness matrix, and

a,S ::::: constants.

The amplitude of a damping matrix must be determined to have relatively

small modal damping factors (say 0.01 to 0.03) at the elastic stage of

the structure, because a large amount of energy will be dissipated through

inelastic hysteresis of structural members.

 106

APPENDIX F. GRAVITY EFFECT

The effect of gravity loads (commonly called P- 6 effect) on

story shear and overturning moment through the lateral story displace-

ment was considered in the program.

The gravity effect in unbraced multistory frames was studied by

S. C. Gael (1969) on both elastic and inelastic response to earthquake

motions. The effect was found to be significant in the elastic response,

and insignificant in the inelastic response for his example structures.

However, the inelastic response studies of single degree of freedom

systems by P. C. Jennings and R. Husid (1968), C. K. Sun, G. V. Berg, and

R. D. Hanson (1973) indicated the significant influence of the gravity
effect on the collapse of the systems.

When the gravity effect is included in the biaxial inelastic

response study by R. Padilla-Mora and W. C. Schnobrich (1974), the dis-

placements of systems increased considerably in a biaxial response

compared with a. uniaxial response.

The gravity effect is generally replaced by equivalent lateral

loads which will yield the same magnitude of overturning moment. Let

us consider a free body diagram (Fig. F.l) of a story of a structure.

Forces acting at the top and the bottom of the story are shear force V.,
l

overturning moments M. and M. l' and gravity load P . Relative dis-

l l- l

placement 6. of the two adjacent floor levels is

l

6.l = X.l - X.l - 1 (F.l)

 107

Shear force V. at story i can be calculated from the equilibrium

l

of moment about floor level (i-I),

V.l = (M.l + M. l)/h. + P.~./h.

l- l l l l
(F.2)

in which h.
l
= inter story height. Gravity load P. at floor level i
l

is the sum of all weights at and above floor level i. The first term of

the right hand side of Eq. F.2 represents shear force normally used in

the conventional analysis, and is equal to the sum of total horizontal

forces action at and above floor level i. The second term is an additional

shear due to the gravity effect. The total shear force V. must be resisted
l

by the structure.

Let us replace the gravity effect term Vi of story shear by

fictitious horizontal force F. acting at floor level i. From the equili-

l

brium of horizontal force at level i, (Fig. F.2),

= (X.l x.l - l)P./h.

l l
X. )P. I/h.
l l+ l (F. 3)

- - ___
Pl+ l X P'+l + ...l:.
+ (_l_ P.
h l'+l h h
i+l i+l i

Therefore, a vector .{F} of fictitious horizontal force can be expressed

as follows:

{F} = [K ] {x}
p
(F.4)

in which [K ]
p
= tridiagonal constant matrix with main diagonal element of
 108

i-th row to be (P. l/h. 1 + P./h.), the left

l+ l+ l l

off-diagonal element (-P./h.), and the right

l l

off-diagonal element (-Pi+l/h +l ).

i

{x} = displacement vector

Matrix [K ] can be considered as a sort of stiffness matrix, and some-

p
times is called as a "geometric" stiffness matrix, which may be added

to an instantaneous stiffness matrix in an analysis.

 109

.--v.
1

P.
1

Fig. F.l The P~6 Effect

F.
1

Fig. F.2 Equivalent Lateral Force for the

P-6 Effect
 110

APPENDIX G. NUMERICAL INTEGRATION OF EQUATION OF MOTION

The equation of motion can be solved numerically assuming the

properties of the structure does not change during a short time

interval. Based on the instantaneous stiffness and damping matrices,

the equation can be expressed in an incremental form,

[M] { ~i} + [C] {~x } + [ K] { x } = - [M] {y} - {~F} (G.l)

in which

[M] = diagonal mass matrix;

[C] = instantaneous damping matrix;

[K] = instantaneous stiffness matrix including the P - ~

effect terms;

{ ~x} = displacement vector at floor levels relative to

the base;

{~y} = acceleration vector of the base;

{~F} = unbalanced force at the beginning of the present

time interval.

Linear variation of acceleration over a short time interval ~t

was assumed in solving the equation of motion by a numerical method.

With this assumption, the incremental acceleration {~x} and velocity {~x}

were expressed in terms of incremental displacement {~x} and the previous

step acceleration {x}, velocity {x} and displacement {x}.

(G. 2)
2
{~x} = 6{~x}/~t - 6{x}/~t - 3{x}
 111

Equation G.l combined with Eq. G.2 can be solved for the incremental

displacement {~x} in the form

{~x} = [A]-l {B} (G. 3)

in which

[A] = (6/~t2) [M] + (3/~t) [c] + [K]

..
{B} = [M] [(6/~t) {:~:} + 3 {x} {~y} ] - {~F} (G.4)

..
+[c] [3{:~:} + (~t/2) {x} ]

The incremental velocity and acceleration were calculated from Eq. G.2.
 112

APPENDIX H. MEMBER RESPONSE AND UNBALANCED FORCES

Response of a structure could be computed by numerical integration

of the equation of motion in an incremental form as described in Appendix

G in terms of horizontal displacements at floor levels. These horizontal

displacements were assumed to be the same in any constituent frames.

Rotations at the joints of a frame in the global coordinates could

be calculated from Eq. D.4 in Appendix D, based on instantaneous stiffness,

horizontal floor displacements, and known unbalanced moments at the begin-

ning of current time interval.

Member end rotations of a beam in its local coordinates could be

taken as rotations at the corresponding structural joints without trans-

formation. Rotations of a column (Fig. C.3) must be computed by Eq. c.8

in Appendix C from floor displacements and structural joint rotations.

Incremental rotations inside the two rigid zones as a simply supported

member (Fig. C.l) could be determined by Eq. C.5 in Appendix C. Moments

at the ends of the flexible portion of a member were calculated from

Eq. C.3. Rotations of each spring element or flexible line element could

be calculated from the moment and instantaneous linear flexibility relation

of each element. The new values of element force and displacement were

determined by adding the current increments to the element response

values at the beginning of current time interval.

The calculated element force and deflection may not satisfy a specified

force-deflection relationship because a nonlinear force-deflection

relationship of a member element was approximated by an instantaneous

 113

linear stiffness relation of incremental force and deflection (Fig. H.l)

on the assumption that the characteristics of member elements might not

change drastically during a short time period.

There are three ways to cope with the discrepancy between a calcu-

lated and a specified force-deflection relationship.

The first measure is a rigorous method, which involves a "lengthy"

iteration process of recalculating element response values, within the

same time interval of numerical integration, until calculated force and

deflection converge to a point on the specified force-deflection curve

of an element (Fig. H.2). Each iteration step requires the formulation

of a structural stiffness matrix for modified element'stiffnesses, numerical

integration of the equation of motion within the same interval, and the

calculation of element forces and defections. The amount of computation

depends on the speed of convergence.

The second measure is an approximate method, which ignores the

discrepancy between the calculate and specified force-deflection relations.

In other words, the force-deflection curve of an element is used only to

define the stiffness of a member based on either calculated force or

displacement (Fig. H.3). This method has an obvious weak point; calculated

force or displacement may not be developed in an element as shown in

Fig. H.3.

Third measure, which is an alternative of the second measure, is to

modify either calculated force or displacement or both to fit a specified

force-deflection relationship (Fig. H.4), although the resultant force

or deflection may violate the equilibrium or compatibility condition in

a member or at a structural joint. In order to restore the equilibrium

condition at the point, an unbalanced amount of force must be applied.

 114

Or in order to make the deflection continuous, the modified amount of

deflection must be given to the point. The restoration of either equilibrium

or compatibility condition at a structural joint will change the dis-

placement configuration of an entire structure. As the displacements of a

structure is related to velocities and accelerations through the numerical

integration assumptions, a change in displacement will induce a change in

velocities and accelerations, and complicates the correcting procedure.

Therefore, the restoration of equilibrium or compatibility condition will

be carried over to the next time interval.

The discrepancy between calculated and specified force-deflection

relationships may be reduced in the second and third measures by using a

very short time increment in the numerical integration of the equation

of motion. However, it was found (Goel, 1967) that the gain in closeness
of calculated force and deflection to the specified curve by reducing the

time step interval was discouragingly slow as compared with the increase

in the computation time involved.

This program employs a. principle described as the third measure to

correct the discrepancy. The continuity of displacements were maintained

at a structural joint. A member was subjected to member end displacements,

which were compatible with the displacements at corresponding structural

joints. For a given pair of member end displacements, the distribution

of bending moment was determined such that bending moment distributed

linearly within a member, and that the calculated forces and displacements

of the constituent elements of a member fitted the specified relations.

The procedure of finding bending moment at the member ends may deserve

more detailed description. The bending moment at the ends of the flexible
 115

portion of a member could be computed from the instantaneous elastic

stiffness of a member, a pair of member end displacements, and geometry

of a member (Eqs. C.3 and C.5}. Then, the computed moment on the two

springs was modified (Fig. H.5) as follows. Computed moment MI was

first modified to value M2 so as to fit a specified moment-rotation

curve based on computed rotation eIj of the joint spring. Modified moment

M2 was used in a flexural rotational spring to determine rotation e~

on the basis of the instantaneous stiffness. Moment M2 in the flexural

rotational spring was, if necessary, modified to value M3 in order to

fit the specified moment-rotation curve of the flexural spring based

on rotation e~.
The equilibrium to be satisfied in the two springs, previous moment

M2 in the joint spring must be modified to M3" If a hysteresis rule

was to be applied to determine rotation e~ at moment M , some cautions

3
must be exercised because of the fact that the hysteresis rules depended

on a stress history. If the moment in the joint spring took a stress

path from MO to M3 by way of M , then the stage of stress M3 was unloading

2
j
from peak stress M0 and rotation e~ must be determined on an unloading
c J

curve. On the other hand, if the moment increased from MO to M3 directly,

the stage of stress M3 was loading, and rotation e~ should be determined

on a loading curve. Obviously, the latter stress path must be considered

in this case. Therefore, the check points created by a path from MO to

M2 must be completely revised and replaced by a new set of check points

in the hysteresis rules.

The sum of modified rotations from an elastic flexible element, a

joint spring and a flexural spring would be smaller than a given rotation.
 116

Therefore, the difference between the given rotation and the existing

rotation was considered as a new incremental rotation, and the pro-

cedure was repeated until the difference became negligibly small.

However, as the moment-rotation relations of the two springs were

piece-wise linear, not more than two cycles of iteration were necessary

within a member.

Member end moments thus modified do not satisfy the equilibrium

condition at a structural joint. The unbalanced moments thus created

at structural joints were treated as external forces in the response

computation at the beginning of the next time interval.

If the intensity of unbalanced forces is equivalent or larger than

that of incremental member forces, a solution will be disturbed by the

unbalanced forces. On the other hand, if the step time increment size is

taken to be sufficiently small (much less than a value required by the

convergence criteria of numerical integration), the effect of unbalanced

force can be reduced.

 119

Flexural Spring

Joint Spring

. I
, I
I
1J I
--f--

I
'I
I i ........... _____ ..........
I 71
11 2 J
I
l . . . . . . . . . . . . . . . . . . . . . . . . . . . _. . . . . s-

t=t
o

eJ ej ej ef ef ef
0 3 1 0 2 1
Joint Spring Flexural Spring

Fig. H.5 Force-Deflection Correction of Joint and

Flexural Springs
 120

APPENDIX I. LISTING OF FORTRAN IV PROGRAM

The FORTRAN IV program code is listed in the following in this

appendix. The total length of the program is approximately 132 k bytes.

The program does not require a tape, disk or drum during execution.

However, computed response values may be written on a magnetic tape with a

minor modification as described in Section 7.2.

In order to analyze a larger structure, the program can be expanded

easily by changing the size of arrays in the program in accordance with

the comment in the main program.

LISTING OF THE PROGRM~

5AKE0005 C ~T3(.' UNLDAOIN~ CCEFFICIENT SAKE0310

( SAi<f:: SMEOOIO C (Q~MCN Ir~NSN/AL(NSTICOI,BA(NSTlcOI.AB(NSTIFOI SAI<E0315
c 5AKEJ015 C AU*) FLEXIFILF. LE"'GTH r:!F ME"IBER 1.' SAKE0320
c 5AKE0020 RA,., L~NGTH PATIO OF RIGID ZONE TO FLEXIBLE lDNE AT A-END SAKE03Z5
c INELAST!C RE~PONSE ANALYSIS PROGPAM $AKE0025 c PBI*I LENGTH RATIO DF PIGID ZONE T~ FLEXIBLf ZlNF AT R-ENC SAKE0330
c FOP REINF~RCED C0NCPETE FPAMF STRUCTURES. SAKEO::J30 c CQ~MON IELEMNT/EIINMM8R),ESTI~~MBR*3).FO(NMMB~.21,FR(NMMBR$21. SAKEC335
c 5AKE0035 c F 2 NMMBI< II< 2) .021 IIMMRR *ZI , W2 (NMM8R*Z) I STFO (NM"IBR I, SAl( (0340
c BY CT ANI. S. <;AKE0040 c ISTFIINMMBR*2I,ISTF2(NMMBR*2I,lVL(NMMBR*ZI. SAKE03~5
C 'IN ,..ARCH 26, 1913 SIlKE0045 C KVL(~~MBR*2) SAKE0350
C MODIFIED eN NO VEMBEP. 4. 1 q14 SAI(E0050 c EII*) LIbEl OF MEMBER ,.) SAI<F0355
f. SAKE0055 C EST'''', STIFFNESS MATRIX ISYMMETRIC MODEl OF FLEXIBLE PART SAI<E0360
C DATA FRO,.. C,IIRDS SAI<EOC60 c OF ME"IBER SAI<[0365
C 20A4 STNAMEI*),*=1,20 SAI<E0065 C FC'.I INELASTIC FLEXIBILITY OF FLEXURAL ELEMENT 1*' SAKE0310
C 311 1, lINIT(>lrI,$=1,3 SAI<F0010 c ~R(*' ROTATION FLEXIBILITY OF ROTATIQNAL SPRING '$) SAI<E0315
C 515.5)t,4FIO.(\ NSTQRV,~ROWS9NSTIFO.NSTIFl,NSTIF2.AE.G. SAI<E~"15 c F2(*) CURRENT MOMENT AT THE FLEXIBLE ENO OF ~=~BERS SAKE0380
C ALPH.~ETA SAI<E0080 C r2(*1 CURRENT DISPLACE~ENT SAKE0385
C AFllO.O WT($).$=l,NSTORV SAI<EOOA5 C W2(*1 ROTATICN OF THE SPRING SAKEf.) 390
C 8 F 11).11 HT(~I.$=l,NSTORV SAKEOOqO C ISTFOI*1 STIFFNESS POINTER DIMENSIONS' OF MEMBER *) SAKE03q5
r I5,5X.3F10.0 I,8AII hAL I hElBU) .I=l.NSTIFO) SAKE0095 C ISTFl(~1 STIFFNESS POINTER (FLEXURAL ELEMENTI OF MEMBER (*) SAKE0400
c I5~5X.bFIO.O I, AI I ,CM( I I. Y"II I I .UHf II. VP I I.UP I. SAI<EOI.OO C ISTFZ(*1 STIFFNESS POINTER (ROTATIONAL SPRING) OF MEMBER ,.) SAKE0405
c I=l,~STIFll SAI<E0105 ( LVL(*) STEP POINTER IN HYSTERESIS RULE SAKE0410
C I5.5X.4FIO.0 (Y.VVlrl.VUIII.WY(II.WUIII.I=1.NSTIF2) SAI<EOHO C KVLC) LEVEL POINTER SAKE0415
c 2IS (MSTORV(I~OW).MAAVS(IROW. 5AKEn115 C CQ~MQN IHYSTR/FROINMMBP*ZJ.FRIINMMBR*2',FR2lNMMBR*21.FR3(NMMBR*21,SAI<E0420
c RFlO.C (Wln .... 1 I1.I=l."IBAYS(IR'JWII. SAKE0120 C DROINMMBR*21,OR1INMMBR*21.DR2(NMMBR*21.0R3INMMBR*2),SAI<E0425
c 16I5 II.ISTFOII),ISTFlAIII,ISTFIB(II. SAI<E0125 C XDOINMMBR*ZI.XD1(NMMBR*21,XD2(NMMBR*2I,XD3INMMBR*21,SAKE0430
c: ISTF2AIII,ISTFZeIII,I=I.NFI,IROW=1.NRQWSI SAKE0130 C VRO(NM~BR*21.VP1(NMMBR*2).VR2(NMMBR*2),VP3INMMAR*2'.SAI<E0435
C ZOll4 ECNA~E(~I.*=1.20 SA~E0135 C WRO(NMMBR*21.WRl(N"IMBR*2I,WR2(NM~BR*2).WR3(NMMBR*2).SAKE0440
C 4Y5,bF10.0 IBSIZE.NSTART,NSKIP,KFLAG.DT.ASCALE.TSCALE,SAKE01~0 C XWO(~~~BP*21,XWl(NMMBR*2I,XW2(NMM6R*2),XW3(NMMBR*2.,SAKE0445 ~
C ZERQ,DELT SAI<E0145 C FDUNfi4MBR*2D.FD2(NMMBR*I I.FWUNMMBP.*21. SAKE0450 I\)
C 20jl4 F~T(~I =1.20 5AI<E0150 C FRMINMMBR*41.CRMINMMBR*4I,VRM(NI4MBR*41,WR"IINMMBR*41 SAKE0455 ~
C F"'T TC,-),AC*'.$=l.IBSIIE (I<FLAG=O) SAI<E0155 C FRO(.) UNLOADING FORCE LEVEL FPOM STEP 6 SAKE0460
C OR AC(*I.=l.IBSIIE (I<FLAG=lJ S~KE0160 C FRll*1 UNLOADING FORCE LEVEL FROM STEP 8 SAKE0465
C SAKE0165 C FR2(.' UNLOADING FORCE LEVEL FROM STEP 10 SAKE0410
C SAKEO 110 C FP3'.' UNLOADING FORCE LEVEL FROM STEP 12 SAKE0415
c CQMMON ICONST/IFLAG,NMMBR,NROWS,NNODE.NSTORV. SAKE0115 C nRO{"'l UNLOADING DISPLACEMEfiT FROM STEP 6 SAKE0480
r: MBAYSINROWSI,MSTORYINROWSI SAKE0180 C DRU"') UNLOADING DISPLACEMENT F~Of'l STEP {I SAKE0485
c IFLAG CONDITION FLAG SAKE0185 C DR2(*1 UNLOADING DISPLACEMENT FROM STEP 10 SAKE0490
c ~MMBR TOTAL NUMBER OF MEMBERS SAI<E0190 C I)R31 *1 UNU1ADING DISPLACEMENT FR.O~ STEP 1.2 SAKEO~95
r:: NROWS NUMBER OF PARALLEL STRUCTURES SAKF 0195 C XQO(*I INTERCEPT OF DISPLACEMENT AXIS WITH STEP 6 SAKE0500
c NN~DE T1TAL NUMBER OF ~ODES SAI<E0200 C Xfl(*l iNTERCEPT OF DISPLACEMENT AXIS WITH STEP 8 SA~f0505
c NSTOPY NUMBER OF STORIES SAKE0205 C XD2(*1 INTERC~PT OF DISPLACEMENT AXIS WITH ST~P 10 SAKE0510
C M~AYS(*I NUM~EP ~F BAYS CF *-TH STRUCTURE SAKEOZIO C XD3(*' rNT~RCEPT OF DISPLACE~ENT AXIS WITH STEP 1Z SAI<E0515
C MSTORVC*j NUMBER OF STORIES CF *-TH STRUCTURE SAI<[0115 C VRO(*I UNLOADING FORCE LEVEL FRO~ STEP 4 SAKE0520
C COMMCN ICURVE/DCINSTIFll.DVINSTIFll.CtJiINSTIFl,.YM(NSTIF1,. SAI<E0220 C VRl(*) UNLOADING FORCE LEVEL FROM STEP 6 SAKE0525
C FDCINSTIF11,FDY(NSTIF11.FDU{NSTIFll,~TO(NSTIF11, SAKEC225 VP2(*1 UNLOAOING FOPCE LEVEL FRC~ ST=P 8 SAKE0530
BT1(NSTIFll.~TZ(NSTIFll.VV[NSTIF2I,WV(NSTIF21. SAKE0230 c VP3C*j UNLOADING FORCE lEVEL FROM STEP 10 ~AKE0535
C ~SPINSTIF2I,PSP(NSTIF2),8T3(NSTIF21 SAI(E0235 c WROI*1 UNLOADING ROTATION lEVEL FROM STEP 4 SAKE0540
C DC($I (RACKING DISPLACE"'ENT SAKE024r) WRl(~1 UNLOADING ROTATION LEVEL FROM STFP 6 5AI<E0545
r: !:IV(") VIELDING DISPLACEMENT SAKE0245 c WR21*1 UNLOADING RCTAT!ON LFV=L FRO~ STEP B SAI<E0550
c 0'("'1 CRACKING MOM~NT SAKEOZ50 c WP3l) UNLOADING ROTATION LEVEL FRCM STEP 10 SA~E0555
c V""., VIELDING MCMENT SAI(EC255 C XWO(*' INTERCEPT OF POTATION AXIS WITH STEP 3 SAI<E05bO
c FrCI*1 DISPLACEMENT FLEXIBILITV BEFORE CRACKING SAKE0260 'C XWl(*1 !NTERCEPT OF ROTATION AXIS WITH STEP 5 SAKEC565
c FDV~ I DISPLACEMENT FLEXIBILITY BEFORE YIELDING SAKF0265 C XW2(*) INTERCEPT QF FOTATION AXIS WITH STEP 1 SAKF0510
c FrUI"" nISPLACE"IENT FLEXIBILITV AFTER VIELDING SAKE0210 C XW3C*) I~TERCEPT OF POTATION AXIS WITH STFP q SAKE0515
C PTO(*I UNLOADING COEFFICIENT SAKE0215 C FR~(-I UNLOADING FQRCE LEVEL ON PRIMARY CURVE SAI<E058C
c FlT1('" I UNLOADING COEFFICIE~T SAI<EOZ80 C OPMC.' UNLOADING DISPLACEM~NT eN PRIMARY CURVE SAKF0585
c pT?(>t1 UNLOADING COEFFICIE~T SAKE0285 C VP~(*1 UNLOADING FORCE LFVFl ON PRIMARV CURVE FOR SPRING SAKE05QO
c VY("'I YIELDING F~RCE ~At<'=0290 C WP~($) UNLOADING ROTATION ON PRI"IARY CURVE FOR SPRING SAKEC595
c WVI*' ~~TATYJN AT YIELDING SIlK[02q5 C FDl('" I IJ!SPl/ICEMENT FlEl<IBILITY AT UNLOADING AFHR YIELDING SAt<E0600
C ~SPlllrl ELASTIC FLEXI8ILITY CONSTaNT SAI<F0300 C FOZI.' DISPLACEMENT LOACING FLEXIBILITY AT STEP b SAKE0605
c PSP(*I PLASTIC FLEXIBILITV C~NSTANT SAKE0305 FWll~1 UNLOADING STIFFNESS AFTER VIEL~ING OF JQINT SPPING SAI<E0610
 (~~w:~ IRCNST/Alo~.R~TA.crl,C(2.D~T.DElA.G,GF2 SAKE0615 rJME~5IrN ~(1250I,R(6Z51.SI551.FMT(2Q).TCllll.A(1221 SAKE!'lnC
alPH Cl=FF!CIENT TO O[TER~!~~ DAMPING
~ATRIX IA(=alPF*AMI SAKE0621" ~1UAl~ PP~(ISIQN B.S SAKE0925
~cTA C~~FFI(IENT TO ~fTERMIN: rAMPING MATRIX (A(=BETA*A~1 SAKFOb25 "OUIVAlENr~ (TC(1).~C(1211 SAKF.0930
rrl (6.0/CT+3.0*AlP~)/DT SAK~0630 C SAKF;0935
C ((2 1.0+3.0$BE1A/DT SAt<E0635 PEAr (5.50?) STNA~E SAI<E0940
C rOT TIME I~TERVAL OF ~UMERICAl INTEGRATION SAt<E0640 REAQ (5.50ZI UNIT SAI(E0945
nFlA !NCREMENTAl GRCUNC ACCElEqATIQN S IIKF. 0645 R"AD (5,500) NSTORy,NROWS,NSTIFj,NSTIFl,NSTIF2,AE,G.ALPH,BETA SAKEC9S0
c G EARTH GRAVITY SAKE0650 WRITE 16,6001 SAKE0955
c (P2 CURRENT GROUND ACCELERATION SAKr=0655 WRITE (6,6021 STNAME.NSTORY,NRCWS.NSTIFO,NSTIFl,NSTIF2,AE,G.ALPH, SAKf0960
C CO~~r:N ISTOT/DSP(NSTOPY,VELINSTORYI,AC(INSTORYi,RSFINSTOPY" SAK!=0660 1 BETA, UNIT SAKE0965
C DDe(~STQRVi.S~(NSTORV~(NSTQRY+1'/2).R8TINNODEI. SAKE0665 NP~YS=lO SAI<E0970
( BK(NNODE*(NSTOR'+lll.HT(~STORYI,SH(NSTORYI, SAI(E0670 N1=1+NBA~S SAKE0975
r SMNSTORYI .~~(NSTapVI.WT(NST~RYI SAKE0675 NZ=Nl+NSTIFl SAKE0980
C DSP(qj DISPLACEMENT RESPONSE AT I~)-TH L~VEL SAI<E0680 ,N3=N2+NSTIFl SAKE0985
C VEL'-) VELOCITY RESPONSE AT (~)-TH lEVEL SAKE0685 N4=N3+NSTIFl SAKE0990
f ACC*) R~lATIVE ACCLERATION RESPONSE AT (~)-TH lEVEL SAI<E0690 N5=N4+~STIe1 SAKE0995
C RSFI*I P~SIDUAl FORCE AT '*)-TH lEV~L SA~E0695 N6=N5+NSTIFl 5AKE1000
( DPF(*, ~AMPING F~RCE AT (*)-T~ LEVEL SAt<f0100 f\l7=N6+NSTJl'" 1 SA"El~05
C ROTIAI ROTATICN OF JOI~T 1$1 SAKE0705 NB=N7+NSTIF2 5AI<EI010
r S~(.I REDUCFD STRUCTURAL STIFFNESS MATRIX SAKE01l0 CALL STOA T 4 (A(1"A(N1I,AIN21,AIN3D.AIN4I,A(N5B.AIN6I,AIN7t,A(N8I,SAKE1015
C ~KC*I TRANSFrRMATION ~ATRIX OF STQRV DISPLACEMENT TO JOINT SAK~01l5 1 NSTIF~.NSTIF1.NSTIF2.AE.GI SAKE1020
( HT(*' RELATIVE STORY HEIGHT OF STORY '*1 SIII<"0170 C SAKFl025
(. S~I*) MASS TIMES ~cIGHT OF LEVEL 1*1 FROM THE SASE SAKf0125 NCHANl=NSTORY+NSTORY+3 SAKE! 030
C SM,) MASS OF LEVEL ~) SAK~()73:> NMMB2=~~MBP+NM~8R SAKEI035
( W~(' TryTAL WEIG~T AenVE LEVEL (*) SJII(E0735 NOS=(NSTORY+ll*NSTORY/Z SAI<EI040
r WT,-j P-DELT~ EeFECT TERM (=WH(*I/HT(*11 SAKE 0740 Ml=l+NSTORY SAKEll)45
r cOMMeN ITITlE/STNAME201.EC~A~~(201.UNIT(31 SAI(EC745 MZ=Ml+NST~P~ SAKEI050
(
STNAME(*J NAM~ OF THe STRUCTUPE TQ ~E ANALVZeD
ECNAMEC$, NA~E ~F T~E EART~QUAKE REC~RO
SAKE0150 c INFOR~ATI~N ABOUT EARTHQUAKE RECORD SAKEl055
SAKE0155 READ (5,5021 EQNAME SAKfl060 f-!
C UNIT (*l UNITS Of lENGTH.FCRCE AND TIME SAI(E0160 RfA~ (5.S04) IBSIZE,NSTART.NSKIP,KFLAG.DT,ASCALE,TSCALE,ZERO,DELT SAKEl065 I\)
r SIIKE0765 REAn (5.502) F~T SAKElOl'O I\)
C DIMENSla~ ~(lZ*~MMBR+NSTORY*4.61. SAKf0770 IF (NSTART.LE.01 NSTART=1 SAKEl075
c 8 (NSTORY>I<Z+NNJOE) SAKE077S IF (NSKIP.LE.O) NSKID=l SAKEl080
C OR (MS*(MS+l'/2+MS*3*'MB+1D+MS*(MB+lDt'M8+2', SAKE0780 IF (ASCILE.EC.O.O) ASCALE=(.O SAKElOa5
C SIISTORV+l)$NSTOP.Y/ZB SA~E0185 IF (TSCALE.EQ.O.O TSCALE=i.C SAKEl090
r. MS NUMBER OF STORY CF A FRAME SAKF.0790 IF (KFlAG.l~.O KFlAG=O SAKEI095
C MA NUMBER OF BAYS OF A FRAME SAKE0795 IF (KFlAG.GF.11 KFLAG=l SAKEHQO
r SAKE0800 I=KFL~G 5AKEl105
CC !'liMe N ICONSTfIFLAG.NMM8R,NRQWS.NNODE.NSTORY. SAKE0805 KFltiG=KFLAG+l SAI(EH10
1 MBAVSi5I,MSTORY51 SAKf0810 DDT=IJBLE (Or) SAKH1l5
cr""~CN ICURVE/DCZO).DV(20'.C~(201.YMI20I, SAKEOR15 WPITF (6,6041 EQNA"'E.F~T SAKE1l20
1 FOC(20I,FOVIZO)tFDU(20).~TOI20).BTl(ZOI.BT2(20). SAKE0820 WPITf {6.6061 Y8SIZE,NSTART.NSKIP,I.DT,ASCALE.TSCAL~,ZERO.DELT SAKE1125
2 VY(201.WY(ZOJ.ESP(201.PSP(20D.BT3(20J SAKE0825 N1=1+NCH~Nl SAKFl130
CO""M[ll\! IDMNSN/AlZO BA20j.P~(20' SAI(E0830 1IJ2=Nl+IIJCHANL SAI(H135
CQ"IMCN IELC.MNT/EIIIOOI.ESTI300i.FDIZOo,FRI2001 .F2(2001.DZ(2001. SAKE0835 fl!3=NZ+~~"'82 SAKF.1l40
1 WZI200I,ISTFOlIOOI,ISTFIIZOOI.ISTF212001. SAKE0840 N4="'3+fV"''''~2 SAKEl145
Z LVLI200J.KVlI2CC, SAKE0845 N5=N4+f\M~1'l2 SAI(~1150
COP>lIdCN I~YSTP/e~0(2001.FR1(200).FR2(200i.FR3(2001. SAKEOB50 \\I6=1II5+I'\M"'A2 SAKE1l55
I. DRO(ZrOI.OR1(zrCI.OR2(200I,DR31200), SAKECB55 N7=N6+NP-'Mfl2 SAI<E1l60
2 X~0(200'.KDI(2CO).XDZI200).XD3200). SAI<E0860 CAll INITl4 fA ( 1 I , A (I'll 11 A HI 2 I AI N3 I A IN 4 I , A ! I'll 5) A I'll 6 A I I'll 71 , SAKH165
3 VRO(ZOOI.VRIlZOOI.VR2(ZOO).VR3IZ00I, SAKF.:0965 1 NCI-lANL .~MMB21 SAKE1l70
4 WROlzonl.WRl(ZQ~J.WR2(2~QI.WR3(200). SAKEOR7r. C WRITE BASIC II'PJR~ATIGN Cf\ Ie DEI/ICE 12U SAKEll75
5 X~0(200).KW1(2CO'.XWZ(200).XW3(2001. SAI(E0615 C WP n e (21) (STNAME( I I.! =1,201, (r:QNAMEI I 1.1= 1,20). (UNIT( I ,1=1.31, SAKEHAO
6 FD1(2001.FD2(200D.FWl(200). SAI<E 08RO ( I NS:ORY.NROWS.f\"'M~P.I'\N~DE,G.nT.NSKIP, SAKEH85
1 FR~140C).ORM(4~~.VRMI4001.WRM(4001 SAKE0885 C 2 NS T IFC.NSTl e l,NSTIF2, SAKHl<110
cc~~rN IRCNST/ALPH.BETA,CCI.CCZ.DOT.~ELA.G.GR2 SAt<E0890 3 ISM(YI.I=l,NSTIJRYI.IHTIII,I=l,NSTORYI, SAKEH95
C
OOU~lE PRECISION ALPH.BETA.CC1.CC2.00T,DElA SAI(EOfl'~5 4 (MPAYS(I).I=l.~~WSI,("'STCRV(II,I=l.NROWSI, SAKElZOO
C
CCMMrN ISTOTfDSPIIOI.VELIIOt,ACCIIOI,RscIIOI.DPF(10),SK(551. SAK e 0900 C 5 (Al(II,I=l,NSTIFC). SAKE1Z05
1 RQT(601.R~(6601.~T(10"SH(101.SM(10I,WH(101.WT(lOI SAKE0905 C 6 (AA(lI.I=l.NSTtFOI'(BBH 1,I=I,NSTIFOI, SAKEl210
D~UPlE PRECISI0N DSP.VEL.ACC.RSF.DPF.SK.R1T.6K SAKfC910 1 (CM(!).I=l,NSTIFI'.iYM(II.I=l.NSTIFll. SAKElZ15
C
CO~~0~ ITITLE/STNAME(201.EONA~E(20I,UNIT{31
SAKE0915 C 8 CfC(YI.I=l,NSTIF1),(DY(II.!=I.I\!STIFU, SAKEIZ20
 9 (VVITI,I=l,NSTIF21.IWV(I I,I=1.~STlc21, SAKE 1225 WRITf {6,61CI JST=P SAKE1530
1 (~III).l=l.~~~RR).(ISTFOII).I=l.NM~ARI, S AI< "tHO ST~p SAKE1535
2 IISTFl( I I.I=t,f\;MME\ZI,( rSTFZ(! I. !=l,NW~B21 St-KE1235 500 F!.'P"'AT {5I5,5X,4 F IO.01 SAt<:E1540
WI' IT E (6,6';,,) SAKElZ4n 502 F0R~AT (20A41 SAKE1545
ISKIP=IIJWtP SAKE1245 504 FOR~j\T (4I5,6F10.01 SAKEl550
JSHP=O SAKEIZ50 600 F'lPMAT (lHl,/lI,lHO,lOX, SAKE1555
NSTFP=O SAKEIZ55 1 'N1NLINFAR PIC FRA~E RESPONSE PROGRAM (SAKE",III, SAKEl560
KBSIlf=IBSIZE+l SAKElZ60 2 IHO,TZ5,'DEVFLOPED AT UNIVERSITY OF ILLINOIS AT URBANAo,I,SAKEl565
ACIKPSIlEI=ZFR:J SAKE1265 ::l IH ,T35,'RV 'JTANI, So',,, SAKEI57(1
Tf(KPSTIEI=O.O SAKEl270 4 IH .135,'1N MARCH Z6. 1973',/. SAKEl515
TT=O.O SAKF.IZ75 5 IHO,TZ5,'SPCNSCRED BY NATIONAL SCIENCE FOUNDATION',I, SAKE15BO
r.~2=O.0 SAKElZ80 6 14 .T35,'GRANT NUMBER GI-Z9934',1111 SAKE1585
G1 Te(Z4,16I,KFLAG SAI<EIZ85 602 FQRMAT (lHO,lOX, 'PROBLEM',TZO,ZOA4,111. SAKEl590
16 TCI11=0.0 SAKflZ90 1 IHO.lOX,NUf'lBER OF STORIES' ,T55,IS'/, SAKEl595
D0 18 I=Z.KBSIIE SAKI"1Z95 2 IH .10X.'NUMBER OF CCNSTITUENT FRA~ES'.T55.15,1, SAKE1600
18 TC(I'=TC(Y-ID+DELT SAI(E1300 3 1M .10X,'NUf'lBER OF TYPICAL MEMBER DIMENSI1NS',T55,I5,I, SAKFl605
ZO AC(ll=AC(KRSIlEI SAKf1305 4 IH ,lOX,'NUMBER OF TYPICAL FLEXURAL ELEMENTS',T55.15,1. SAKEl610
P!ClID 15,FMTI IACI),I=Z,KBSIlEI SAKE1'310 5 IH ,lnX,'NUMBER OF TYPICAL R~TATIONAL SPRINGS',T55.15,I, SAKE161S
WRTTE (n,FMT) (AC(II.I=Z,KBSIZEI SAKE1315 6 IHO,lOX,YOUNG'S MOQULUS',T40,lPEZO.3,I, SAKE16Z0
DfLT=O.O SAKE 13Z0 7 IHO,lOX,'EARTH GRAVITY',T40,E20.3,1. SAK~16Z5
na ZZ i=Z.KBSIIE SAKE1325 8 IHO.IOX,'COFFFICIENT ALPH fOR DAMPING MATRIX',T40,020.3.I.SAI(E1630
Z2 DELT=DcLT+APS(ACIIII SAt(El330 9 14 ,lOX,'CQEFFICIENi BETA FO~ DAMPING ~ATRIX.T40,DZO.3,I.SAKEI635
IF OELTI 210,Z10.Z8 SAK!:1335 1 1 4 0,lOX,'UNIT OF LEI\GTH' ,T56.A4'/, SAKE1640
24 TCIl'=TC(KBSIlfl SAI<E1340 2 IH .10X,'UNIT ~F FORCE',T56,A4,I, SAKE1645
A(lll=AC(KRSIlF) SA K[1345 3 IH ,lOX,'UNIT OF TI~E.T56,A4.111 SAKE1650
RFAD (5.F~TI (TCIII.ACIII,I=Z,KBSIlEI SAKE1350 604 fORMAT (lHl.l/lII.lHO.lOX,'PASE "10TIOIII',TZ5,ZOA4,I, SAKE1655
WRITE (6.FMTI nCIII,AC(II.I=Z.I(RSIIt::1 SAKEl355 Z IHO,10X,8FORMAT'.TZ~,ZCA4) SAKE1660
Z8 IF (INSTEP+IRSIlEI.GE.~STARTI GQ TQ 30 SAKE1360 606 FORMAT (lHO.IOX,'INPtT CA~A'.I. SAKE1665
NSTEP=NS1~P+IBSIIE SAK=1365 1 IHO,15X,'NUMBER OF OATA POINTS PER CARO',T55,I5,1. SAKE1670 f-'
GO TO (24.201, KFLAG SAKE1370 Z IH ,15X,"START POINT OF THE CALCULATIQN',T55.15.I, SAKE1675 I\)
30 CO 100 J=Z,KBSIlE SAKF 1375 3 IH .15X,'NSKIP FOR RECCROINGG,TS5,I5,I, SAKE16BO LA)
NSTEP=NSTEP+l SAK!:1380 4 1H ,15X. 'FLAG FOR OUAKE DATA CARDS' ,T55.I5", SAKEl685
IF (NSTEP.LE.NSTARTI G1 TO lCO SAI<El31l5 5 IH .15X,'TIME INTERVAL OF COMPUTITION',T45.1PElS.3,I, SAKE1690
DELTAT=(TCIJI-TCIJ-lll*TSCAlE SAI(El39C 6 IH ,15X,"SCAlE FACTeR FOR DATI",T45,EI5.3.I, SIKEl6QS
IF (OELTAT! Z10.100,40 SAKE139S 1 IH ,15X,'SCALE FACTOR FOR TIME',T45,E15.3,I, SAKEl100
40 GP=(~C(J-II-IEROI$~SCAlE SAKEl400 R 1H ,15X,"BASE LINE S~IFT OF EO. ACC.',T45,E15.3,I, SAKE1705
SLOPE=(~C(J)-ACCJ-l)*ASCALE/CELTAT SAKE140S 9 IH .15X,'T!ME INTER\I~L OF QUAKE RECORO,T45,El5,;3.1ln SAKEl1l0
TI=-TT SAKEl410 60S FOR~AT (lHl,I,lHO.lOX,"EATHQUAKE RECORD',I!II SAKf1715
iT= TT+DELT AT SAKF 1415 610 F':RMAT (lHO.lII,lHO,lOX, SAKEl120
50 TF (DT.GT.TT) GO TO 100 SAI(E14ZC 1 'T8TAl NUMBER OF PESFC~SE STEPS ',14,1111, SIKE17Z5
TT=TT-DT SAKE1425 2 IHO,40X,'ALL DATA PReCESSFD') SAKE1730
T I=TJ+OT SAKE'l430 END SAKE 1735
(:Rl=GR+SlOPE$TI SAI(E1435
nIOUI=OBL~ (GRI-GR Z t SAKH440
GRZ=f,Rl SAKF1445
C RESPONSF CALCULATION SAKEl4S0
IF IIFLAC.EQ.OI GO T0 60 SAKEl455
C~Ll STRST4 IS,B.NOSI SAKE 1460
6') CALL RESPN4 IS.fHlI,811"U.8(plZI,NDSI SAKEl465
ISKIP=ISKIP+l SAKE1470 SUBRCUTINE STCAT4 (WIDTH,AI,CP,VP,UP,DU,UM,VU,WU,NSTIFO,NSTIF1, SAKE1775
IF (ISKIP.LT.NSKIP) GQ TO 50 SAK::1475 1 NSTIFZ,AE,GI SAKEl7BO
ISKIP=O SAKE1480 C
JSTEP=JSTEP+l SAKFl785
SAKE1485 C PFAD STRlCTURAL INFCRMATICN SAKE1790
CALL STORI04 (A(11.AINll.A(NZI.AIN3I,AIN4I,ACN5"AIN61.AIN71. SAI(El49Q C SAKEl 795
1 NCHANL,NMMB21 SAKFl495 C CP"~CN ICONST/IFLAG.NMMBR,NRCWS,NNOOE.NSTORY, SAKFl800
GO TCo 5C SAI(E1500 1 MAAVSI5I,MST~RV(51 SAKEl805
1 '')() CCNTINUE SAKE1505 COMplCf\J ICUR\lF/f1CI201 .DV(ZO) ,CV"ZOI,YMIZOI,
r,n Tr (24,ZO),KFLAG
SAKE1810
SAKEl510 1 FDC(ZOI,FDV(ZOI,FDUIZOI,BTOIZOI.BT1(ZOI,BTZIZOI, SAKElB15
C ENO 0F THE EARTHOUAKE RECCRO SAKE1515 Z VVI20I,WVIZOI,ESPIZOI,PSP(ZClI.BT3(ZI)I
Z1,) Cill L MXRS P4 (A ( 11 , II ( NU ,~I N2 II I N31 , A( N4 Ii. Ai N5 I AI N6 AI N7 ) , SAKE18Z0
SAKFl520 CO'-'MrN I(lMNSN/AL(201."IA( ZO) ,RAIZOI SAKEl8Z5
1 ~CHANL,~~plPZ) SAKE1525 Cr'lPo1Mr.N 1Ft Eplf\T / '= I ( 100 I, EST( 300 "Fe I 200 I FP ( ZOO I F Z( 200 I, DZ( ZOO I SAKFlI130
 1 W2 P () (; & 1ST F C1 ( 1 00 , I S TF 1 200 , ! STF2 ( 2 00 I SAKElfl35 YI31=C~IIl/u~(II SAKE2140
2 LVl(200),KVlI2001 SIIKflA40 YI41=Y"" IlIUM! I) SAKF2145
COMMON ISTDT/DsollOI.VEl(lO ACC(10i,RSF(lOI,DPFllOI,SKI551, SAKf1845 Y(51=YI31"YI3I SAKE2150
1 RIJT(60I,I3tt:1660),HHIOJ,SHIIO),SM(l0I,WHIIOI,WTIIOI SAI(El850 Y(61=VI31*V(51 SAKE2155
DOUBLE PRECISION CSP,VfL,ACC.RSF.DPF,SK,POT,BK SAI(HB55 Yn I = III 0- Y(4) I II< YI 4 HI 1 .0- Y4 I 1/ Y 2 I SAKE2160
CQ~~CN ITITlE/STNAME2~1.~QNA~E(ZOI.UNIT(31 SAKEl860 YI81=YPIII/V(41 SAI<EZ165
COMMON IwI)R~/T.J.K.l.NB.NE.KA.KB,IROW.INODE,JNJDE,ISTORY.JSTORY, SAI(El 865 DUII1=Q.16667*1!2eO+Y(4ll*V(7)+V!4)*(leO+YI411-Ylb)-Y(6I*Y(B~ SAI<E2110
1 V101.Xl,YL,X.Z.LLI51 SAKFl810 1 +0.3333*Y(51*CPIII SAKE2115
DlAAENSIO~ WIOTHII),AIINSTIFlJ.CP(NSTIFll.YPINSTIFll. SAKEIB15 FOCI I 1=IDU( I I-~VI I II/fUM( I l-Y~ H)) 511.KE2180
I UPINSTIFII.DUINSTIFll,UMINSTIF1,.VU(NSTIF2I,WU(NSTJF2J SAKH880 36 WilITF (6.6261 T.r.CII),DYIII.NI(I1 SAKEZ185
C SAI(E1885 C 5AKE2190
PEdn (5.500 (WT(I.I:l.~STCR~1 SAKEl890 WI< IT E 6. 6 3 0 SAKE2195
PFAC (5,500) (HTIII,I=I.NST1PYI SAI<EI895 ~~ 40 J=1.NSTI~2 5AI<EZ200
W~YTE 16,6001 UNIT(ZI,UNIT(ll SAKE1900 READ (5.5021 I.VVIIJ.VUIII.I>JV(II,WU(ll SAKEZ20'5
WPIT~ (6.602) II.WTIII,HTII.I-I,NSTCRYI SAKEl905 40 WP IT r: I 6 6'3 2 I I. VY( 1) , VU II I .;, V u wU( H 5AI<EZZI0
WH(NSTaRYI=~T(NST0PYI SAKE1910 ~0 60 I=1.NSTI F 2 5AKE2Z15
IF (NSTCRY.EQ.1J GO TO 12 SAI<E1915 8T3(1'=0.0 SI11(E2220
DO 10 I=2,NST~RV SAKf1920 IF (VV(III 50,52,56 S;l.I<E2225
J=~STORV+1-1 SAI<E1925 C ~EC~A~ICAl HINGE SAKE2230
10 WH(J'=WH(J+1)+WT(JI SAK!;:l930 50 ESPII)=-i.O SAKEZ235
12 Xl=OeO SAKEl935 PSP( 1.=-1.0 SAKE2Z40
~~ 14 I=l.~STaRV SAI<E1940 GO TO 54 S.II.I<E2245
SM(I)=WTII/G SAI<E1945 C N~ SPRING SAKE2250
WT(I'=WH(II/HT(II SAKEl950 52 ESP(Il=O.O 51\I(E2255
Xl=Xl+HT(II SAKEl955 PSP( Y1=0.0 SAKE2260
14 SHIi=S~iJ~Xl SAI<H960 54 vvn 1=1.0::' 24 511.1<1:2265
WRITf (6,610) SAKEl965 vuo 1=2.0E '24 5At<EZ210
DO 20 X=l.NSTIFO SAKEl910 GO TO 60 5'\I<EZ215
Rf4D 45,5021 J.BA(J8.Al(JJ.BB(JD f--'
SAKEi915 56 fSP(li"'Il/~(II/VV(n SAI<E2280 I\.)
WRITE (6,6121 J.8A(J).AliJI.RBIJD SAKEl'9I'1Q PS P n,. WU I )-WY I II VU I I-VV Ii 5AKE2285 +:-
IF (RA(J'elT.O.OI GO TO 16 SAKE1985 60 CONTINUE 5AKE2290
IF (~B(J'GlTeOeOJ GO TO 16 SAI<El990 C SAI<EZ295
IF (Al(JII 16.16,18 SAKH995 NMNlBP:::O 5AI<E2300
16 WRITE 6.6141 SAKE200('! NNODE=D SAKE2305
STOP SAK!:2005 n~ 100 IROW=l.NROWS 5AKE2310
18 A~Ji=BA(JD/AliJI SIII<1:2010 READ 15.5041 MST~RYIIROWI.MBAY5(IROWI 5A1(E2315
2~ RPJ,=BBJD/AlJI SAKE2015 IF (MBAVS(IROWI.lE.OI ~BAYSIIROWI=O 5AKEZ320
C SAKE2020 WRITE 6,640) IROW. MSTORY IPOI.I.MBAVS( IPO"l1 5AKE2325
WRITE (6.6201 SAKf2025 IF IMST'JRY(IPOWI.GT.OI CO Tr 74 SAKf2330
VO 30 J=l.NSTIFl SAKI:2030 WPTTF (6.66<:) SAKE2335
~FAC 5.50Z hAHU.C'HU II.WHlhVPU I.UPUI SAI<E2035 STOP SAKEZ340
30 WPITE (6,6221 I.IIID.CM( 1"~II.U~lll.YP(I'.UPIII SAKEZ040 74 IF (MSTCRY(IRQW).lE.NSTORVI GC TO 16 SAKE2345
WQITf 16.62'4) SAI<E2045 WRITf (6,662> SI11<E2350
~n 36 !=I.NSTIFl <;AI<E2050 STOP SAKEZ355
BTOi II=Co5 SAKE2!)55 16 N8=MRAYS(IR[~I~MSTORVIIROWI SAI<Z360
~n! 1-100 SAI(E2060 NF=NB+NS+MsynRVIIROW. SAKE2365
BT 2CTI=lloO SAKI:Z065 NNODE=NNGn~.NA+MSTOPY(IR~WI SAI<EZ310
U(II=AIHD*AE SAI<EZ'J7n. IF (~BAYSIRnWI.lE.OI GO T( 17 SAI<E2315
IF (AY(I'eGTeOeOD GQ TO 32 SAI<E2015 J="'IBAVS! IRrW) 5AI<E2380
CMIII=O.O 5AI<I:Z080 PEJlD 15.500' (WIDTHlIhY=l.JI SAI(E2385
32 If ICIVl(I'.GT .. OeO. GO TO 34 SAI<EZnS WRITE 16.642) UNITel} SAKE239C
A H I 1= YM (l DIYP II 5.11.1(1:2090 WRITE (6.6441 (I.WIDT~(II.I=l.J) SI\I(E2395
34 CP(TI=C~II'/AI{II SAI<E2095 17 l~~~~BR+~MMBR-l SAI<E2400
rc (! I=CP( I 13.0 SMEZlOO WRITE 16,6461 5AI(E2405
FnCII&=0.333333/AYIY! SAKE2105 D(1 110 J=l.NE SAI(E2410
V(31=CM(II/YM(II SAKE 2' 11 0 READ (5,5041 I,llCKI.K=l,51 SAKE2415
V!41=Y!3''''Y(31 SAKE2115 KA=l+r+I SAKE24Z0
Y 51 =VB ''''YI4o) SAI<E2120 KB=I<A+l 5AI<EZ4Z5
DYII)=0.3333*1.O-Y(5*VP(II+YI4J*CPllli SAKEZ1Z5 ISTFO(N~MBR+I'=lllll SAflE2430
FDV(tl=(DY(I)-DC(III/(YM(I)-C~IIII SAKE2130 ISTFl(KAI=lU21 SAKE2435
V(21=!1.0-U~(II/YM(I)/(1.o-UP(II/YP(III SAKEZ135 ISTFUK~)=lU31 5AKI:'2440
 ISTI=2(I<A)=LLU.' 511KE2445 IH ,TZ, 'ELEMENT' .T1Z, llIiERTIA' .TZZ, 'CRACKING', SAKE2750
tSTF2( KA l=lU51 SIIKE2450 4 T33,'YIELDING',T44,'ULTIMATE',T55,'YIELDING', SAKE2155
IF II.GT.N~I GO TO 78 SAKE245,) ') T66. 8 ULTIMATE',/1 SAKE2160
C BFAMf.'EMBE\R SAKE:?460 62Z FORMAT (IH ,I4.3X.lPQEll.3./1 SAKEZ765
INO':lE=!+I-ll/"'l8AYS( IROWI 5AKE2465 6Z4 FORMAT (IH ,11.l~O.lOX. SAKEZ170
JNrJrF=!NCDE+l SI\KE2470 1 'FLEXURAL ROTATION AS A UNIT-LENGTH SIMPLE MEMBER',I. SAKE2175
IST~RY=MSTORVIROWJ-(I-11/MPAYSCIRoW) 5AKf2475 Z 1M ,lOX.'(ELASTIC DEFORMATICN INCLUDED)',II, SAKE2180
JSTCR Y= I SFIRY SAKE2480 3 IHO.T2.EL~MENT,TI6,CRACKING.T31.YIELDING. SAKE2165
GO TO 80 SAKE2485 4 T46,'lJL TIMATE' ,n SAKEZ790
C C~LUMN MEMBER 5AKE24QO bZ6 FORMAT (1M ,I5.3X,IP5EI5.3./J SAKEZ195
78 TNODE=I-NB SAKE2495 630 FrRMAT IIHl,/I.IHO,10X.'CHARACTERISTICS OF TYPICAL JOINT SPRINGS',SAKEZ800
JNODF=INCDE+~BAYS(IROW+1 SAKE2500 I II.l~O.T17.SPRING'.T24,'~CMENT'.T49,ROTATION.I. SAKE2805
ISTCQY=MSTQRYIR~W)-'INODE-11/(MBAYS(IROW+I) 5111<E2505 Z IH ,T16,'YIELDING',TZ9,'ULTIMATE', SAKEZ810
J9>TCRY=! ST'JIlY-I SAKEZ510 3 T42,YIELDING',T55,'ULTIMATE',/1 SAI<E2815
80 WRITF 46.6481 I,ll(K).K=I.5'.INODE,JNODE,ISToRY,JSTiJRY SAI{E2515 632 FORMAT (1M .IIO.IP8EI3.3,/1 SAKE2820
L=NMMeR SAKEZ5Z0 640 FORMAT IIHI.II,l~ ,10X.'FRAME',IZ.I, SAKE2825
WR IT!: 6,650 J SAKEZ525 I IHO.15X,'NUM~ER OF STORIES',T35.IZ.I, SAKE2830
DO 90 I=I.NIE SAKE2530 2 IH .15X,'NU~PEP OF 8AY5,T35.IZ,/' SAKEZ835
L=L+l SAKE2535 642 FOR~AT (IHO.III.IHO.IOX.'STRUCTlJRAL DATA ALONG BAY',I, SIIKEZ840
K=ISTFO(U SAKE2540 I IHO. T9.BAY.TZO,~IDTH'.I, SAKEZ845
K/l=L+L-I SAKE2545 2 IH .TZO,"',A4.,,../' SAKE2850
1C8=KA+1 SAKF2550 644 FORMAT (IH TQ.I3,TI2.1PIE15.3,/1 SAKEZ855
KA=ISTFHKA' SAKEZ555 646 FQP~AT (IH~.III.IHO,IOX.'Mf~8ER IDENTIFICATION',II. SAKEZ860
ICfI=IST"'HKIU SAKE2560 1 IHO,T17.'FLEXURAl ELEMENT',T34,'JOINT SPRING', SAKE2865
IF (I.GTeNS GO TO 86 SAKE2565 2 T49,'JOINT NUM~ER.T64.'STORY NUMBR",I, SAI(E2810
C BEAM \ll!EIl4~E~ SA'<F.2510 3 IHO. T2.'MEMBER'.TIO,'OIMENSN',TIB,'(A-ENDI(B-ENDa'. SAKEZ875
J=I- 1-1& IMeus IP.mn ,lIrM8AYS( !ROW) SAI<EZ515 4 133,' A-END H 8-EM)", T48, , A-END) B-END', SAI(EZ880 f--J
Xl=WIOTHlJD SIIKE2580 5 T63,'(A-ENO)(A-ENO)'./D SAKE2885 l\)
GO TO 88 $AKEZ585 64B FOR~AT (IH .I4.I7,2X,4118,I7' 5AKE2890 \Jl
C CCLUMN I'4IEMBER SAKEZ590 650 FORMAT (IH .III.IHO.IOX.ME~BER PROPERTIES',/I. SAKEZ895
86 J=MSTQRY(IPO~-I-NB'-ll/MBAYS(IROWI+l' SAKEZ595 1 IH .TIO.'STIFFNESS,T23,DFlEXIBLE'.T39,'RIGID ZONES', SAKE2900
XL=HT( J SAKE2600 2 T63,"TOTAL LENGT~'.I. 5AI(E2'905
88 YL=~L(K).(1.O.8A(KD.BB(K~1 SAKEZ605 3 IH ,T2.'MEMBEP',TI2,GCEI",TZ3.olENGTH'.T35,'IA-END", SAKlE2CJlO
X=AUIO*SA(K) SAKEZ610 4 T41,'(B-ENDI'.T5B.'STR. DMN.'.T10,'MEM. DMN',/J SAKE2CJ15
Z=AUI(.*BIHIO SAKE2615 652 FORMAT (lH ,Y4.1)(,lP8EIl.3./' SAKE2CJ20
f:J(L)=O.5.(AIKA)~AI(KB). SAKEl620 654 FORMAT (IHO.III,lH ,IOX.'''$$$ ERROR ",1, SAKlEl9Z5
WRITE (6,6521 1.EICLI,AL(KI.X.l.XL.YL SAKE26l5 1 IH .20X,'DIMENSrON SPECIFICATION IS WRONG') SAI(E2CJ30
XL=ABS(XL-VLI/(XL+Yll SAKEZ630 660 FORMAT (IHO.III.IH .10X,~*$** ERROR **.*.',1. SAKE2935
IF XL.LT.O.OOOID GO TO ~o SAKEZ635 1 IH.20X,'MSTORV SPECIFIED IS NOT POSITIVE'.I' SAKE2'940
WRITE (6,6541 SAKE2640 662 FORMAT (IHO./II.IH ,IOX,'****~ ERRO~ ***.*',1, SAKlE2945
STrtP SAKE26'1>5 1 IH ,ZOX,'MSTOPY SPECIFIED IS GREATER THAN NSTORY',/) SAKEZ950
C fIe. )=1I6EI SAKE2650 END SAKEZ955
90 EICl'=Al(Kl/lEI(li$6.0) SAKEl655
100 NMtlI~R=NII'''BP.NIE SAKE2660
RFTURI\I SAKE2665
500 FORMAT 8FIO.IIB SAKE2610
502 FOIH4.H I5,5)(,1HOsO) SAKEZ615
504 FOP~AT H6I5) SAKE2680
60~ FORMAT HHO.II,lHO.IOX,[lATA ALONG STORY HEIGHP.I. SAKE2685
L IMO,TIO,STORy,TZO,"WEIGHT",T35,"HEIGHT',I, SAKE26CJO SURRCUTI~E INITL4 (C~X,DMN,FX.FN,DX.DN.WX.WN.NCHANL,NMMBZ) SAKH995
Z II-1.TZO,'(I,.V Q',T35."C'"V.,P,/1 SAI(E2695 C SAKE300C
602 FORfAAT (IH 5X.15,lPEI5.3, 1"15.3) SAKEl100 r INITIAL CONOITIONS SAKE3005
610 Ff1IHJlAT (lHI./,/,lHO.I0)(.DI~EN5ION OF TYPICAL ELEMENT",II. SAKE2705 C SAKE3010
1 IHO. T4. NUMBER '11, e-END'. T30.' ("II DOLE', SAKE2HO CO"'MON IC~NST/IFLAG,~M~BR.NRC~S.NNODE.NSTGRY, SAKE3015
2 T'I>3. (S-END''./J SAKE2715 1 MBAVS(51.MSTORY(SI SAKE3020
612 FORMaT IH , 3X.IZ,5X.IP8E13.J./) SIIKF27Z0 COMMCN I(URVE/DC(20J.DY(ZOI.~!Y(ZOI.YMIZOI. SAKE3025
614 I:ORMAT (lHO./II,lH .10X,.*~*. ERROR **.*.',1. SAKE2125 1 FDCIZOJ,FDY(ZO),FDU(20I,BTOIZOI.8TIIZOl,BTZ(201. SAKE3030
1 IH .ZOX,"DIMEIIISION GIVEN IS NOT POSITIVE'I SAKE2730 Z VY(ZOI,WY(20I,ESP(lOI.PSP(201,BT3(ZOI SAKE3035
620 FCR~AT IH1.11.1HO,1I.0X, SAKE2735 (OMMeN IELEMNT/~I(1001.~STI300).FD(ZOOI,FR(ZOOi.F2(ZOO.D2IZ001, SAKE3040
1 'MOMENT CURVATURE PELATION CF FLEXURAL ELEMENTS.,!II. SA!(E2740 I W2(2001,YSTFOIICOI,!STFl(200I,YSTI:Z(:?00i. SAKE301t5
Z U' TIZ. "MOMENT' "T34, '!YOMENT", T61,' ClJRVATUR!:' ,/" SAKE2745 Z LVL(ZOOI.KVL(ZOCI SAKE3050
 (OMMeN IRCNST/ALPH.RETA.CCl.CC2.rOT.O~lA.G.GR2 SAKF30'55 cn~MCN IELE~NT/=I(lOO},ESTI3001.FDI20~I.FR{2001.FZIZOO).DZ(ZO01, SAKE3360
O'JUP.lc PPFCISYiJ"I AlPH.BETA.CCl.CC2.0'n.OEU\ SAKE3060 I WZ!ZOOI.ISTI=C(1001,ISTF1(ZOOI,ISTFZ(ZOOI. SAKE3365
cr~MrN I~T9T/OSP(10J.Vfl(10).ACC(lOI.PSF(iOI.OPFIIOI.SKI55). SAKf3065 Z lVl(ZOOI.KVL(200) SAKE3310
1 R'lT(60D .~K (660)> .... T1 OJ. SHllO). SM(lOI,WH HO I.WT (10' SAKD070 Or~E~SI~N V(41 SAKE3375
DOUBLE PPECISION DSP8VEl.ACC.RSF.DPF,SK.~OT.BK SAKE 3015 C SAKE3380
r.OMMC~ IWORK/X!41,I.K.IE.~F.KA.KB SAKE3t)8(} LE=YSTFOIIEI SAKF3385
DI"'ENSION DMX(NCHANl) .DMN1NCHANU .FX(NMMAZ) .FN(NMMB2. SAKE30A5 IF (I=RiKAI.GE.O.OI GO TO 10 SAKE3390
1 DXNMMBZ.CN(NMMBZI.~X(NMMB2).WN{NMMBZI SAKE3090 e STIFFNESS MATRIX QF HINGED-FIXED ~EMBER SAKF3395
SAKE3095 EST(f4E 1=0.0 SAKE3400
YFlAG=3 SAKBIOO F.ST(MF+IJ=O.O SAKE3405
CCi=(6.0DO/DDT+3.0DO.ALPHI/OOT 5AKF3105 EST(MI=+Z)=1.0/EI(IEI+EI(IE!+FDK81*Al(lEI+FR(~B&1 SAKE3410
CC2=1.ODC+3.0DO$BETA/DOT SAKE311C C STIFFNESS MATRIX OF HINGED-HINGED ~E~BER SAI<IE3415
01] 10 K=1.~MMB2 SA1<E3115 IF (F~KR).lT.O.OI EST!f4[+2J=0.0 SAKE3420
lVUKI=l SAKE3110 GQ TO 30 SAKE3425
I=DIK.=O.O SAKE3125 10 IF (FRIKBI.fE.O.OI GO Te ZO SAKE3430
1'2 lit' 1=000 SAKIC3130 C STIFFNESS MATRIX O~ FIXED-HINGED ~E~BER SAKB435
n2(K)=0.0 SAKE3135 fSTC~EI=l.O/(EI(IEJ+EIIE)+FC(KA)*AL(lEI+FR(KA SAK'=3440
KVUK '=1 C;AI(E3140 EST( i"'E+11=O.O SAKlE3445
1= IS1F100 SAKE3145 ~ST(ME+2)=0.O SAKE3450
FRCK'=ESPCI) SAKE3150 rc 30
Ge SAKE3455
11) 14200=0.0 SAKE3155 C FlFXIBILITY MATRIX OF FIXEr-FIXED MEMBER SAKE3460
KA=-1 SAKEH60 20 V(21=-EHI<:' SAKB465
ME=-2 SAKE3165 Y(I'=-YI21-VC21+FDIKAi*Al(lEI*O.5+FR(KAI SAKE3470
Dn 12 IE=I,NMMBR SAKf3170 V(31=-Y!7'-V(2!+fD!KBJ*AlllEI*C.5+FR(KB) SAI<IE3475
I(A=KA+2 SAKE3175 c INVERSION OF FLEXiBILITy ~ATRIX SAKE34f110
Kf'\=KA+l SAKE3180 Y(41=VI11*V(31-Y(2'*VI11 SAKE3485
Mf=I\IIE+3 SAKF3U5 ESTI~EI= Y(31/Y(41 SAKE34QO
12 CALL STIfF4 X.IE.KA.KB.MED SAI<E3190 EST(~E+11=-V(2IfV(41 SAKE3495 /---'
DO 20 Y=l,NNOOE SAKE3l<~5 EST(~E+2)= V(11/Y(41 SAKE3500 I\)
?O ROTC 1)=0.000 SAKE3200 30 RETURN SAKE3505 0\
DO 30 I=1.NS10RY SAKE3205 EIIID SAKE3510
DPFI I D=C .. noo SAI<E3210
DS!)( I )=OeOOO SAKE3115
vnn D=OcODO SAKE3210
30 ~CC U=O"OOO SA1<F.32Z5
c !NITIALIl~lION Of MAXI~UM RESPONSE ON 10 DEVISE (llD SAKE3230
DO 40 I=l,NCHA~L SAKE3l35
I)MXU ':000 SAKE3240
40 01\111'1(1)=0.0 51\1<1"3145 SUPPQUTINE STRST4 (R,S9~DSD SAKE3550
DO 50 I=1.~MM82 SAKE3150 c SAI(n555
FXI I )=0 cO 5AKI=3255 r. FCRMATION OF STORY STIFFNESS MATRIX SI*) SAKE3560
FN(II=OcO SAKE3260 c I=CRMATION OF ~ACK SUBSTITUTICN MATRIX 8~$1 5111<:3565
OX(IiI=lJ"t.:! SAI<F3265 C SAKf3570
(INn 1=000 SAI(E3270 r.OM~ON ICCNSTfIFLAG.N~MBR.NROWS.NNODE.NSTORYt SAKE3575
"'X(1)=0.0 SAKE3275 1 MRAYS(51.MST1RYI51 SAKE3580
50 WN(YB=JoC 5AKE32BO IDMNSN/Al(20).BA(lOI,RB(201
C'M~CN SAKE3585
R!:TURfI! SAKE3285 COMMON IElcMNT/EI(100I.EST(300I,FD42001.FR(2001 .Fl(200),02(2001. SAI(E3590
Fill 0 5A1<E329f) 1 Wl(1001.ISTFO(1001.IST F 1(2001.ISTFZ(200), SAKE3595
2 LVLlOO.KVlI100. SAI(E3600
C[MMCN IRCN S1 I AU'H. BETA oCC 1.CC 1.ODT. DELA. G, GR2 SAKE3605
9QUBL~ PRECISION AlPH.AETA.CCl.CC2.DOT.OElA SAI<E3610
COM~CN ISTOT/OSPlOi,Vflll01.ACCI10i.RSF(lOI,DPFllOI,SK(55', SAI<E3615
1 ROH601.9KI660l.HT(101.SHllO),SMllOI.IiHI10i.WHlOI SAKE3610
OOUBlE PRJ;:C I SI"N OSP .VEl ,Ace.R SF .DPF. SK. ROT. BK SAI<E3625
(OMMeN 1~~RK/FA.fB.X.I.J.I<,l,II.JJ.KK.Ll.MM.NN.IE.IN.IS.JS.KA.KB, SAKE3630
SUBROUTIN~ STIFF4 V.IE,~A.KB.~EI SAI<E3330 1 MB.MS.NB.~M.NE.~f4~,NNN.KKK.NS1.NS2.NS3.IROW SAKE3635
c SAI<f3335 OIM~NSICN Slll.R(NDS) SAI<E364(l
c FnR~ATI[N STIFFNESS MATRIX OF A MEMBER WITH EQUIVALENT
~F ~ SAKE3340 rnURLE PPECISION R.S,X.FA.FB SAKE3645
C INFlASTI( SPRING AND FIXED POINT OF INFLECTION .. SAKE3345 C SAKlE3650
r SAKE3351) 00 10 1=1. r-WS SAI<E3655
e~M~rN fOMNSN/AL(20).RA(ZOD.RR(20D SAKF3355 10 SK (l1=O.OCO SAKF3660
c P-~ELTAEFFECT TEP~S SAt<E3665 KA=KA+2 SAKF3910
IF IN5TQPY.EQ.l1 G~ TG 18 SAKE3670 R!r,ID f.~C M1MENTS SAKE3975
S K ( 1 1= 51( 11) -CBL E (WT ( 1) HIT( 2 I I SAI([3675 F 11= C FlL F. ( ( 1. a +0,11 (L c: I I "'F 2 (KA) "PA (LE I'" F 2 ( KA+ 11 ) SAK"'39BO
p: INSTORY .. EQ.2) GC TC 16 SAKE368(1 F B =0 P LE (AB ( L F I *F 2 (K IIH ( 1. O+P 8 ( LE I I F 2 (I< A+ 11 I SAKE39~5
J=l SAKE3685 IF (IF.GT.~RI GO TO 48 SAKE39QO
I<=NST"RV-l SAKE3690 If\J=IE+(I~-lI/Mo,+II SAKE3995
DO 14 1=2,1( SAKE3695 AK(I~I=RK(!N)+FA SAKE4000
J=J+I SAKE3100 AW(IN+l'~AK{IN+ll+FB SAKf4005
SKIJ-l'=SK(J-ll+CBLE(WT(III SAKE3705 GO TC 50 SAKE4010
14 SKIJ1=SK(JI-DBLE(WTIID+WTII+11 SAKE3710 48 IN=IE-N8+II SAI<E40 15
16 SK (NDS-lI=SIC (NDS-1)+CBLEOIT( NSTOR VI I SAKE3715 IS=MS-IIE-NR-ll/(MB+ll SAKE4020
18 SK(NDSI=SK(~DSI-CBLECWT(NSTORYII SAKE3120 AKIIN)=~K(I~I+FA SAKE4025
P-DELTA EFFECT FnRCE SAKE3125 PSFIISl=RSF(ISI-(FA+FBI/DBLE(HT(ISII SAKE4030
CALL SPRDG (SK,DSP,RSF,NSTORV,ll SAKf3130 IF IIS.LF.ll GO TO 50 SAKE4035
c DAMPING FORCE AND I~ERTIA FCRCE AT EACH ST~RV LEVEL SAKf3735 !N=II'\+MB+l SAI<E4040
00 20 I=l,NSTQPV SAKE3140 AK(IN)=RKIINI+FA SAKE4045
20 RSF(I1=RSF(I1+DPF(II SAKE3145 RSF ( IS-l'=R SF ( IS-II +( FA +FB ID~ LEi HTI I SI I 5AKE4050
1 +DBlECSM(II)$(ACC(II+DBLEIGRZI-DELA+ALPH$VEL(Iil SAKE3150 50 CCNT INUE SAI<E4055
KKK=O SAKE3155 C PK(NII.,~S+1)=S(NN,NN).$1-11*BK(NN.MS+11 SAKE4060
MII'M=O SAKE3760 I=KI(K+l SAKElt065
"INN=/) SAKE3165 1T="'5+1 SAI<E4070
DO 100 IROW=l,NROWS SAKE3110 .IJ=F"P+2 SAKE4075
MS=MSTORV (I POlO SAKE3715 I'" (MS.LE.il JJ=2 SAI<E4080
MB=IIo'BAYS(IROkl SAKE3180 "'~=1+II.Sl"NS2 SAKE4085
Ll=3*INlB+l1 SAKE3785 CAll SVBAN (SCI'!MI.BK( 1I.I\N,JJrIII SAKE4090
NII.=PlS*'MB+lI SAKE3790 COP~ECTICN OF ROTATI1N DUE TO UNBALANCED MOMENTS SAK;E4095
KK=NN*MS SAKE3195 00 60 I=l.~N 5AKE4100
NSl="'S* iMS+ll 112 SAKE3800 60 R~TINNN+I)=ROTNNN+lt-BK(KKK+KK+I) SAKE4105 f--J
NSZ=PlS'4tLL SAKE3805 C PEDUCED STIFFNESS MATRIX SAKE4UO f\)
NS3=NN*(MB+ZI SAKE3IHO L=O SAKlE4U5 --.:]
IF (PlS.LE.ll NS3=NN*Z SAKE381'5 DO 10 JS"l.MS SAKE'H20
r Fr.RMATYON OF STIFFNESS MATRIX SAKE3A20 JJ=KKK+""*JS-MB-l SAKElt125
IY=l+NS1 SAKE3825 I I=NSI-LL SAKIE4130
JJ=II+NS2 SAKE3830 00 68 I S= 1, J S SAKE4135
CALL SKFPM4 (S11,SIII.S(JJI.HT,PlS,~B.PlMM.NSl.NS2.NS31 SAKE3835 II:II+LL SAKE4140
c RECONSTRUCTION OF MATRIX BK(NN,MSI FROM COMPACT FORM USING VACASAKE3840 JJ=JJ-MB-l SAKE41lt5
c AR'=A IN AK(*' SAKE3845 L=L+l SAKE4150
II=KKK+l SAK=3850 X=O.ODO SAKE4155
JJ=I<KK+ICK+f\J1\ SAKE3855 00 66 I=I.LL SAKE4160
DC 30 1=1 I,JJ 5AKF386(\ MF"=JJ+I-KKK SAKE4165
30 ~K(II=O.ODO SAI<E3865 I~ (MM.LE.OI GO TO 66 SAKE4170
L=II.S1 SAKE3810 IF "'M.GT.KK) GO TO 66 SAI<E4115
II=KKK-MB-l SAKE3875 X=X+S(II+II*~KIJJ+I) SAKE4180
('11'1 40 IS=l,MS SAKE3880 66 CONTINUE SAKE41B5
II=II+N~-(MB+ll SAKE3885 68 SK(LD=SKCL}+S(LI-X SAKE4190
00 4(1 l=l,LL SAKE3890 70 CrNTTNUE SAKE4195
l =L+ 1 SAKE3895 C CORREC~ICN F~RCc: DUE TO U~BAlANCED MCMENTS SAKE4200
/oIM=II+I-KKK SAKF3900 II=KKK+KK+NN-IIo'~-l SAKE420r;
IF (MM.LE.C) GO TO 40 SAKE3905 L=NSl SAKf4Z10
IF (MM.GT.KKI GO TO 40 SAKF3910 00 90 IS= 1 ,MS SAKE4215
RKI II+Y I=BK II+! )+S III SAKE3915 X=O.OCO SAKE4220
4 n CCNTINUE SAKE3920 I I =1 1-MB-l SAI<E4225
CALCULATICN OF UNBALANCED FORCES AT STORY LEVELS SAKE3925 on ~8 I-=l.LL SAKE4230
!I=I<KK+I<f( SAK!:3930 l=L+l SAKE4235
f-!A="'S*IIfI'J SAKE3935 tJPI=II+I-I<KK-KK SAKE4240
N"1= I "'P+MB+ 11"'jIIS SAKE3940 IF (MM.LE.el GO TO BB SME4245
N!:,="''''''' SAK"'3945 IF (MM.GT.NNI GO TO 88 SAKE4250
KA="'I"M+,..MM-l SAKF3950 X=X+S(LI*BK(IJ+II SAI(E4255
on ')0 IC:=l.NM <;AK0955 BB CrNTYNUE SAKE4260
NE=NE+l SAKE3960 90 PSF(ISI=~SF(ISI-X SAKE4265
lE=ISTFOINEl SAKE3965 Kf(K=KKK+KK SAKE4270
 ~~M=~MM+~S~(MB+M~+lD SAKE4Z15 32 If\J=If-N8 SAKF45B('
100 "~'''=NN''l+~N SAI<E42BO IN= IN+NE /lYS+l SAKE4585
IFtA!;=/') SAI<.f4285 IS=~ST~RV-(I~-ll/(NPAVS+ll SAKE45'W
T=O SAKF42CJO Xll)=!Yll'+(ZJlfHTIISI SAKE45Q5
r~ 200 KA=l.NSTORV SlhKE4Z9S X(ZI=(V(ZI+Y(3)/HTIIS) SAI<E4600
on 20') Kf\=l.K/l. SAK(4300 X(3)=(Xll'-+X(ZII!HT(ISI SAK'=4605
1=1+1 SAKE430'5 I I = IS * ( ! S+ 1 1/2 SAKE4610
II=CCZ*SK (I I
1<<< SAKE4310 AI( ! I I =AI< (I ! I +C Bl E ( X D I ) SAKE4615
ZOO IF (KAofC.K~1 R(II=R(II+CCl~CBlE(SM(KA1 SAKE43IS IF {IS.LE.ll GO T~ 34 SAKE4620
CAll DSYNV (R,NSTnRVI SAKE4320 II=I I-I SAI<E4625
Pt:;TURN SAKElt325 AKdItl=AK(II)-CelEIX(3)1 SAKE4630
EN!' SAKf4330 H=II-IS+l SAKE4635
AI< 0 I 1= AK (I n +0 Bl E( X 3 II SAI<E46ItO
34 II=(NBAYS+l'*14*IS-NSTORY-21 SAI<E4645
Bl<iII+I~I=BI(III+IN)-CBlE(Xll SAI(E4650
IF IIS.lE.l) GO TO 36 SAKE4655
SI(III+JNI=RKlII+JNI-CBlEIXIZII S'&'KE4660
II=II-4*INBAVS+l) SAKE4665
BK I I + I N1= BI( II I. IN I .DBl E Xi 111 SAKE4670
SU~RCUTINE SKFRM4 (AK.BK,DK.HT.NSTORV.NRAVS,NNN,Nl,N2,N31 SAKE4370 BK(II+J~I=RKIIr+JNI+DRlFX(ZII SAKE4675
C SAKE4315 36 OKIYNI=QK(INI+DBlEIYI1 SAKE4680
C FnRMATIC~ nf STIFFNESS MATRIX CF THE STRUCTURE SAKE4380 IF (IN.GT.NNOI')El GO TO 40 SAKE4685
r. SAKf4385 II=IJN-INI*N~ODf+IN SAKE4bq()
C)~~ON IOMNSN/Al(ZOI,BA(ZOl,B8(201 SAKE4?1'lO OK I I 1= DK I I I )+I')SlE I V( 211 SAKE46CJ5
cnMMnN fElEMNTfEIIIOOI,ESTI300I,FDIZOOI.FR(ZOO.FZ(ZOO.,DZI2001. SAKE43CJ5 DK(JNI=DK(J~l+OfllE(V(311 SIII<F.4100
I WZI2001.ISTFOCIOO),ISTFl(ZOO),ISTF2IZ001. SAKE4400 40 CONTINUE SAKElt105
2 lVl(ZOO).KVl(ZOOI SAKE4405 RETURN SAKE'4HO
DIME"SIC~ AKIN1).BK(NZI.DK(N31.HTINSTORVl,X(41.V(41 SAKE4410 END SAKHHS
DOUBLE PRECISION AK.BK.DK SAKE4415 I-'
C AK(*1 STIFFNESS MATRIX (NSTORV$NSTORVI IN SYMMETRIC MODE SAKE4420 f\)
C BKI.) STIFFNfSS MATRIX (NNrCE*NSTORVI IN A COMPACT FORM SAI<"44Z5 (
C DKC.) STIFFNESS MATRIX INNCDE$NNOOEI IN SYMMETRIC BANDED FCRM SAKE4430
Dr) 10 IE = 1 N1 SAKE4lt35
10IHqIE)=0.ODO SAKE4440
Dr) ZO IE=l.NZ SAKf4445
ZO BIIF)=O.ODO SAKE4450 SU~RCUTINE RESPN4 IS,DD,DV.DR.NDSI SAKE415'5
DO 30 IE= 1.N3 SAKE4455 C SAKElt160
30 l'KIIEI=O.OOO SAI(E4460 C SOLUTION OF THE SECCND ORDER DIFFERENTIAL EQUATION OF MOTION SA91:E4165
N~=IIlSTORV"'Ne IlYS SAKE4465 C USING lI~~~P ACCELERATION MET~cn. SAKE'4110
NNOOE=III!'I-+NSTORV SAKE4410 C SAI<E4115
Nf=lIIfI+NWIDE SAK{=4415 CO~MON !CONST/IFl~G,N~MBP.NPCWS,NNODE.NSTORV. SAI(E4180
f'~=/IlNN 5AKE4480 1 MqAVS(51.MST~RY(51 SAI<E4785
KA=NNN+NNN-l SAK!:4485 CC~Mr~ !CURVE/DC!201.DV(ZO),C~(ZOI.YM(znl. SAKE41CJO
!4E=KA+""'N-l SAKE4490 1 FOC(201.FDY(ZOI.FDU(ZO),BTO(20),BTI(201.BTZ(ZOI. SAKE47Q5
(Ie 40 IF=l.NE 5AKE4495 Z VY(201.WY(ZOI.ESP(ZOI.PSP(201,BT31201 51OKE4800
NN=NN-+l SAKIC4500 CaMMGN IDMNSN/Al(20I,BA(201.PBIZOI SAKE4805
!( /I=K /I+Z
SAKE4'50'5 CCM~CN !~l~MNT/EI(lOOI.ESTI300I,FD(2001.FP(2001 ,F2(2001.02l200I, SAKE4810
1<8=1<1\+1 SAK[4510 1 W2(ZOOI,ISTFO(100I,ISTF1(ZOOI.ISTF212001. SAKE4815
MF=~E+3 SAK!:4'H5 2 lVlIZOO).KVl(ZOCD SAKE4820
I.F= ISTFO(NN I SAKE4520 COMMCN IPCNST/AlPH.8ETA.CCl.CC2.DDT.DElA,G.GR2 SAKE4825
Ylli=FST(MEI~(1.O-+BA(lE)!+EST~E+11$eA(lEI SAI<E4525 D1UBlE PPEClSIJN AlPH.BETA,CCI.CC2.DDT.DElA SAKE4830
Y(ZI=FST(ME+11*1.O+BAllE)I+EST~F+Z-8AqlEI SAKE4530 COMMON !~TOlfDSP(101.Vfl(101.ACC(10,RSF(lOI.DPFI101.SK(551. SAKElt835
VI3D=FST(~f)$PB(lEI+EST(ME+ll$(1.O+8R(LEII SAKE4535 1 R~T(601.B~lb60I,~T(lOI.SHIIO).SM!lOI.WH{101.WTlOI SAKI:4840
Yi4i=~STI~E+l)*P~(lFI+EST(ME+Z)$(1.0+RB{lEI I 5AKE'454C DOURlE PRECISION OSP.VEl.ACC.PSF.OPF,SK,POT,BK SAKE4845
V{I'=Vll*ll.O+AA(lEJ'+Y(ZI$fA(lEI SAKE4545 COMMON IW0RK/OA.OR,FA.FR.RA.RP.Rl.I.IE,IN,IS,KA,KB.KK.MB.MS,ME. SAKE4850
Y(21=Y(31*(loO+BA(lEll+V(41.eAllEI SAKE4550 1 NE.NB.NN.~~.~~~.N~N,IST.IPOW.JFlAG.KFlAG.ll.lE.X(41 SAKE4855
Y{31=Y(3)*8B(lEI+Y(41*ll.C+PB(lEII SME45">., DIMENSIO~ SINOSB,OO(NSTORVl,rVINSTQRYI.OR(NNQOEI SIOKE4B60
IF lIE.GT.NP) GO TO 32 SAKF4560 rOUBLE PQECtSlr)N S,DD,DR,DV
TN=IE.IIE-ll!NBAYS SAI<E4865
SAKE45b5 C SAKE487(,
IN=IN+l SAKE4570 00 10 I=l.NSTORY
r;o TC 36 SA'<E4815
SAKF4575 10 OD (It=RETI\"'O.ODO*VElIIl+0.5D(}lltACClIl*DDT) SAI(E4880
 117

Calculated by linear
instantaneous stiffness

True Value

Rotation

Fig. H.l Force-Deflection Relation

of Member Element

Rotation

Fig. H.2 Iteration Method to Satisfy

Nonlinear Force-Deflection
Relation
 118

Calculated
Response

Specified Relation

Rotation

Fig. H.3 _Stiffness Based on Force-Deflection

Rotation

Unbalanced Force

Calculated
Response

Specified Relation

Rotation

Fig. H.4 Unbalanced Force to Meet Specified

Force-Deflection Relation
 CALL SPRDG (SK,DD.DV,NSTQRY,ll SAKE4B85 MFMPEP FNO RoTATIC~ SAKE519C
("10 20 I==l,NSTCRV SAI<E48'W 70 flL=AL(Lfl.O.5 SAI<F5195
20 OV (I )=CAL!: (SM( I 11*16.0r')0*VEU I 1/00T+3. OOO"'ACC( I) SAKE4S95 IF (Fql~AI.LT.O.OI BL=ALILEI SAKE5200
1 +ALPH*(3.0DO$VELIII+C.500$ACC(II*90TI-9ELAI+DVtI)-RSFIII SAKf:'4900 IF (FR(KB).LT.C.OI RL=ALIlEI SAKE5205
C CAlCULATY~N CF INCREMENTAL STORY DISPLACEMENTS 00($1 SAKE4905 KFLAC:=O SAKF5210
CALL SPRDG (s.OV,On,NSTORY,l) SAKE4910 AO KFLAG=KFLAG+1 SAKE5215
C CALCULATE ST"RV RESPQNSES SAI(E4915 IF IKFLAG.GE.41 GO TO 90 SAKE5220
DO 30 I=l.~STORV SAI(E4920 !FLAG=O SAKF.5225
("1VIII=3.0DO*(DDIII/DDT-VEllI',-0.5DO*ACCIIl*ODT SAK1'4925 IST=ISTFUKA' SAKE5230
ACC(I'=ACCII1+6.0DO*(OOIY)/rOT-VELIIII/DOT-3.0DO$ACCIII SAI(E4930 ~A=(D2(KA'-FCC(IST.F2IKAIJ.Bl+W2IKA'+(2.0.F2(KAI-F2(KB)$EIINEI SAKE5235
VELCII=VFL(I)+DV(II SAKE4935 IF IA~S(CA-RAI.LT.0.01*CC(ISTI*Bl) RA=DA SAI<E5240
30 [\ S P I I I == 0 S P I )i-D 0 I I SAKE4940 IST=I5TFl(KAI SAKE5245
DAMPING FORCE THAT IS PROPORTIONAL TO STIFFNESS SAKF4945 RB=ID2IKB1-FDC(IST)*F2(KBll*Bl+W2IKBI+(f2IKBI*2.0-F2(KAII*EIINEI SAKE5250
CALL SPRDG ISK,DV,RSf.NSTORV.ll SAI(E4950 IF CABSICA-RAI.lT.0.01.CC(IST)*ALI RB=DB SAKE5255
DO 40 I=l.NSTORV SAKE4955 C INC~E~E~TAL PIlO~ENT IFA,FBD AT FLEXleLE ENDS SAKE5260
DPF I I=CP"'( I I+SETA*RSF( I) SAKE4960 FA=EST(M~I$4DA-PAI+EST(ME+l)*([B-RBI SAKE5265
40 RSF(I1=O.OOO SAK~4965 fA=EST(~f~11*(OA-RAI+ESTIME+21.(DB-RBI SAKE5210
CALCULATION Of JOINT RESPCNSE SAKF4970 ME~BER STRESS-STP.AIN RELAiICNSHIP SAKE5215
KA=-1 SAKE4915 CALL HVSTRS (FA,KA.IFLAGI SAKE5280
PIlF=-2 SAKE4980 CALL HVSTRS IFA,KB.IFLAGI SAKF5285
NE=O SAl< E4985 JFLAG=JFLAG+IFLAG SAI<E5290
~1I4"'=0 SAKE4990 IF (IFLAGI 90,90,88 SAKE5295
fI'''IN=O SAKE4995 ~EMB~R STIFFNESS MATRIX SAI<{:5300
JFLAG=O SAKE5000 88 CALL STIFF4 (X,NE,KA.KR.MEI SAKE5305
C FOR AN INDIVIDUAL FRAME STRUCTURE SAKE5005 GO Te 80 SAKE53l0
DO 100 IROW=l.NROWS SAKE5010 90 CONTINUE SAKES3l5
I14S=I"ST'JRV(IIlOW) SAKI:5015 fdMM==MWoII+KK SAKE5320
I14B="'BAVSIIIH'WI SAKE5020 100 I\NI\=fI,NN+NN SAKE5325 I--'
f\)
NB=MB*MS SAKE5025 TFLAG=JFlAG SAKE5330
NN=NS+MS \0
SAKE5030 RETURN SAKE5335
NM=NfI+I\'N SAK[5035 END SAKE5340
Kf(=NN*IIIS SAKE5040
c INCREMENTAL JOINT ROTATICN DR($) SAKE5045
LL=MMM+l SAKE5050
CALL GMPRD (SK(Lll,DD,DR.NN.~S,1)
[10
JOINT R~TATIGN ReT
50 !=l.NN
c., SAI(E50'55
SAKE5060
SAKE5065
50 R~T(NNN+I)=ROT(NNN+I-DRIII S AK E5010
DC 90 IE= l,NM SAKE5015 SUBROUTINE STORE4 IDMX.OMN,FX,fN,DX,DN,WX,WN,NCHANL.NMMB21 SAKE5380
ME=t4E+3 SAKE5080 C SAKE5385
NE=N"'+l SAKE5085 C KEEP TRACK ~f MAXIMUM RESPONSE VALUES. SAKE5390
LE=ISTFOINE SAKE '5090 r. SAKE539S
I(A=KA+2 SAKE5095 CCMMCN ICONST/I~LAG.I\MMBR,NRCWS,NN~DE,NSTORY. SAKE54CO
I<f\=I(A+1 SAKE5100 1 ~BAVS(5I,MSTOR~151 SAKE5405
IF IYF.GT.NBi GO TO 60 SAKE5105 CO~"'CN IELE~NT/EII1001.ESTI3001.FD(2001,FRIZOOI,F2(200I,D2(2001, SAKE5410
c BEAM ME~RER lIE' SAKf5110 1 W22001.ISTFOIIOO"ISTFlI200I,YSTF2IZ001. SAKE5415
IN=YE+(IE-II/MB+NNN SAKE5115 2 LVLIZOOI,KVLI20Cl SAKE5420
OA=(l.O+f\AtLE*SNGLiROTtINIJ+eSILE,$SNGL(ROT(IN+lll SAKE512() COM,..ON IPCNST/ALPH.BETA.CCl,CC2.0DT,DELA,G,GR2 SAKE5425
DB=BAIL1'I*SNGL(ROTIINID+(l.O+BBILEID*SNGl(ROT(IN+ll) SAKE5125 n~UBLE PRFCtSIO~ ALPH,8ETA,CC1,CC2,DOT,DELA SAKE5430
GC TO 70 SAKE'5130 C~MM~N 15TDT/DSP(101,VELIIO),ACCI10I,RSF(101,OPFI101.SKI551. SAKE5435
C C0LUMN MEMBER tIEl SAKE5135 1 ~OT(60).BK(6601.HTI101.SH(101,SMI101.WH(101.WT(101 SAKE 5440
60 IN= H:-N13+NNN SAKE5L40 PRECISION DSP,VFL.ACC,Q~f,DPF.SK.ROT.~K
D~UBLE SAKE5445
IS=II4S-IIE-NP-II/IMB+ll SAK1'5145 IwORK/A,B.X,V,!,K
cnM~CN SAKE5450
IF IIS.LE.ll GO TO 62 SAKF5150 OIMENSICI\ D~X(NCHANL).DMN(NC~Ahll.FXIN~~821.FN(NMMBZI. SA"KE5455
~A=SNGL(CSP(IS-II-D5P(ISII/ALILEI SAKE5155 1 0 X(NII4MB21 ,ON INMMB21, loX! NMMB21 WN( ~M~fl21 SAKE5460
FB=SNGL(qOTIIN+MB+lll SAI<E5160 C SAKE5465
r,o TC 64 SAKE5165 A=GIl2lG SAKE5410
62 FA=-5NGLIDSPIISII/AL(LEI SAKE5110 IF U.GT .O~XH II DpaXC 11=11. SAKE5475
"'f\=0.0 SAKES 115 IF (A.lT.DMNll)) DMNI11=A SAKE54RO
64 ~A=(l.O+RA(LEll*SNGL(RCT(INjl+eB(LEI*FB+FA SAKE51BO X=SNGL(ACC(l')~GRZ SAKE5485
CA=AA If' I*SNGL qOT ( IN) I +( 1. O"t-IHl( LE) I*FA +FA SAI<E'51B5 V=SNGLIOSP(lll SAKE5490
 A=-X*S~(l)+V*WT(ll SAKE549S WRITE (6,6101 UNITI11 SAKE5800
R=-X~S~(ll+Y*WH(l) SAKE5500 WQ IT E (6.612 I IJ.D~XJ+3).C~~(J+31.DMX(J+NSTDRY+31. SAKE5R05
IF (I\JSTORYcEQ.ll GC T'J 12 SAKE5505 1 0~N(J+NSTQRY+31.J=1.NSTORY) SAKF5810
OQ 10 !=2.NSTORY SAKE5510 WRtT~ (6,6141 DMX[21.UNIT{2).D~N(21.UNIT(2). SAKE5815
X=SNGll~(C(II)+GR2 SAKE5515 1 OMX!31, UN IT (2) .UNITIlI,OMND) ,UNIT' 21 ,UNIT! 11 S AK E5820
Y=SNGlIDSP(II-DSP(I-lll SAKE5520 WRITE THE ~AXIMUM MEMBER RESPQNSES SAKE5625
A=A-X~SMiII+Y*WT(I' SAKE5525 J=O SAKF.5630
10 B=B-X*SHII'+V*WH(II SAKE5530 L=O SAKE5835
12 IF (A.GT.OMX(2)' OMX(2)=A SAKE5535 DO 30 IRQW=1.NRO~S SAKE5840
IF U.LT .DMN( 211 OMN(21=A SAKE5540 NE=~ST0RY(IRO~I$(M~AYS(IRDW)+MeAYS(IROWI+1l SAKE5845
IF (8.GT.D~X(311 [)~XI3)=B SAKE5545 WRTTf 16,616 IRCW SAKf5850
IF CfI.lT.DMNnH 01lf-i(3)=!3 SAKE5'55,) WRITE (6,620) UNIT(2P.UNIT(11.UNIT(11 SAKE5855
f'1J 20 !=1.NSTORV SAKE5555 DO Ie i=l,NE SAKE5860
A=(S~GL(AfC(I'D+GR2)/G SAI<.E'5560 J=J+l SAKE5865
1(=3+1 SAKE5565 IST=ISTFHJD SAKf5810
IF A.GT &D~lK DMX(lO=A SAI<f5570 IF IDYIIST)' 5.5,6 SAKf5815
IF U.LT.DMNIKI) DMN( to=A SAI<E5575 5 11.=0.0 SAKE5880
K=I<+I\STOI'Y SAKf'5580 ~=O.O SAKfS885
A=SNGl (DSP II SAKE5585 (;0 TO 7 5AKES890
IF U.GT.OMX(K) I Dp.lX(KI=A SAKE5590 6 A=DXIJI/CY(ISTI SAKf5895
20 IF U .l ,. 0 MN K I ) DMfHIO=A SAKI:'5595 B=DN(J)/DY(ISTI SAKES900
00 30 1(=1.\\f>1~82 SAKE5600 7 WRITE 16,622) I,FXJl.FN(JI.DXIJb.DNIJI,A,B SAKES905
IF (F2(10.GT.FXIKD l FX( Iq=F2 (Kl SAKE5605 J=J+l SAKE5910
IF (F2(K'.LT.FN(l<i} FNOO =F2 (I() SAKE5610 IST=ISTFUJI 5AI<E5915
I'" W200.H.I)X(KI! DX(KI=C2(1<1 SAI<E5615 IF IDV(I5T)1 8,8,9 SAKE5920
IF W2.IKI.lT .. ON(lO) CNIKI:D2(Kl SAKE5620 8 11::0.0 5AKE5925
IF 0.12 (1<1 .. GT .W)( 1<, WXUO=W2 (I<I StlKE56ZS 8=0.0 SAKES930
30 IF iW2(K'.lT.wN(Kll WN(KI"W2(KI SAI<E5630 1';0 TO 10 I-'
SAI(E5935 W
C WRITE (211 GP2.iACC(II.I=1.NSTORYI.(VEL(II.I=I.NSTORY'. SAI<E5635 9 A=OX(JI/CYlIST' SAl<f594C 0
C 1 OSPi II.I=l,I'\SiCRVI. SAI<E5640 R=Ol'JJ 10Y( 1ST! SAI(E594>5
c 2 F2(I).I=1,N~~B21.i02(II.X=1.NMMB2). SAKE5M'5 10 WRITE (6.62~1 FX{JI.fNJi.DXIJI.ON(JI.A.B SAKE5950
c 3 (W2I.I=1.~~~R21 5AI<E5650 WRITE (6,6261 UNIT(21.UNIT(~).UN!T11.UNIT'11 SAKE5955
RETURN SAI<f5655 1')0 20 I=1.NE SAI<E5960
END SAI<E5660 l=l+1 5AI<E5%5
I S1=I STFZ (l SAI(E5970
IF IWY( 1ST 15,15,16 5AI<E5915
15 ,11=0.0 5AI(1E5980
A=O.O 5AI<E5985
GO TO 17 5AKE5990
16 A=WX(ll/~YIISTI SAKE5995
B=WN( ll/~V( I sn SAKE6000
SUBROUTINE MXPSP4 ICMX.OMN.fX.FN.DX.ON.WX.WN.NCHANL,NMMB21 SAKE5100 17 WRITE (6,6221 Y.FXill.FNiLI.WX(l'.WNILI.A,8 5AI<E6005
C SIII(E5705 l=l+l SAl<E6010
C PRINT MA~IMU~ RESPONSE SAKE51l0 IST=ISTFZU SAI(E6015
c SAKE5715 IF (!PlY I ISH) 18,18,19 5AICE6020
CO~MCN ICONST/IFlAG,NM~BR.NROWS.NNODE.NsrORV. 511KI:'5720 18 A=O.O SAKE6025
1 MBAYS5).~STORV(5) SIl.KE512'5 8=0.0 SAI(E6030
CQ~~CN ICURVf/DC20'.DYd201.C~20'.YMI201. SAKE5730 GO TO 20 SAKE6035
1 FOC(20D.FOV{lOI,FDU(20.BTOI201.BT14Z0.BT2120J. SAKE57J5 19 A=WX(L)/~V(ISTI SIIKE6040
2 VV 20 .WV( 20) ,ESP! 20 .PSP 201 ,8T3 (201 SAKE5740 .
B=WNl'/~V!STI SAI<E6045
co~~rN IElE~NT/EI(lOOI.ESr(3001.FDI200.fR2001.F2(2001.02(2001, SAI<E5745 20 WRITE (6.624P FXIll.FNILI.WXIL).WN(ll,A.8 SAI<f'6050
1 W2(200).IS1FOlOOI.ISTF112001.ISTF2(200!. 5 AK E5 750 30 CONTINUE SAKE6055
2 l 'Ill 200 I.KVU 2001 SAI<E5755 C GR2=0.O 5AKf6060
CO~~CN I~TD1/CSP(10).VELIIOI.ACC(10).RSF(101.DPFIIO),SKI551. 5AI(E'5760 C 00 40 I=l,NSTORV SAKF6065
1 ROT{601.~K(660).HT(10.SH(10)tSM(lO.WHIIOD.WTIIO SAKE5765 C Ace I )=('.000 SAKE6070
f10UPlE PRE'CISION DSP.VEL.ACC.RSF.IJPF.SK.ROT.BK SAKf'5770 r. VEl! I 1=0.000 SAKE6015
r~~MCN ITITlf/STNAME(201.EQ~A~E{201.UNIT(3) SAKE51'75 C 40 DSP! I 1=0.000 SAKE6080
CO~~QN IWOR~/I.J.K.l.NE.IROW.IST.A.B 5$1KE5180 C on 50 I=1.NMMB2 S AKE6085
DIMENSIC~ D~X(NCHANlD.DMNCNCHANLI.FXINMM821.FN!NMMB21. SAKE578"i C F2(! 1=0.0 SAKE6090
1 DXN~MB21.DN(~MM821.WX(NMM821.WN(NMMB21 SAKE5HO C 02(Y)=O.C SAI<f6C95
C WRITE T~E MAXIMUM STORY RE'SPONSES SAKl:5195 r 50 W2(!)=0.0 SAKE6100
( r'J 60 J=l,lO SAKE6105 SAKE64H'
60 WRITF (211 GPZ,IACCnl,I=l,IIiSTORVI,IVELIII.I=l,NSTIJRVI. SAKF6110 C(l~~rN ICUQVF./OC(ZOI,DY(ZOI,CMIZOI,YIIli(ZOI. SAKf6415
c 1 (DSP( II,!=1,II.STCRVI. SAKE6115 1 FDC(20),FCV(20),FDU(20),BTO(20),BTl(20),BT2(20), SAI<E6420
c 2 (F2(II,I=1.N~~B2I.(D2(I),I=1,NMMB2D, SAKE6120 Z VY(201.~V(2CI.ESP(ZOI ,PSP(201.BT3(201 SAKE6425
C 3 (W2(!),I=1.NM~B21 SAKE612S SAI(E6430
Cf'1Mr-'CN IEL:'plNT IE H 100 I, EST! 300 ),FOI 7001, FR( ZOO) ,F21 2001,02 (200 I,
C REWI"Jr 21 SAKE6130 1 WZ(ZOO).ISTFOIIOOI.ISTFlI200I,ISTFZ(2001. SAKE6435
P.E TURIll SAKE6135 Z lVLI20~I.KVL(20CI SAKE6440
610 FORMAT (lHl.tI,lH ,lOX, SAI<Eb140 COMMeN IHVSTR/FRO(200I,FRl(200).FR2(2001,FR3(2001. SAKE6445
1 'MAXIMUM ACCELERATION AND DISPLACEMENT RESPONSES',II, SAKE6145 1 DRO(ZOOI.ORlI2001.DR2(200),DR3(200l. SAI<E64S0
2 1HO.T16,'ACCELERATICN, G'.T43,"DISPLACEMENT,',A4,1. SAKE6150 Z XDC(2001.XDl(2CO'.XDZI200I,XD3(200', SAKI:6455
3 1H ,T2.STORV',T14, 'MA)t8,T27,GMIN',T43,'MAX' .T56,'MIN',11 SAKE6155 3 \Ill. OC ZOO I, VRl(200 ),VR2 (2001. IIR 3( 200 I. SAI([6460
61Z FORMAT (lH ,I4.T7,1P2E13.3,T36.2E13.3./J SAKE6160 4 WRO(2001.WR1(200I,WR2(ZOOl.WR3(2001. SAKE6465
614 FORMAT (lHO,lOX.'BASE SI-EAR MAX ... o.lPE13.3.ZX.A4,'. SAKE6165 5 XWOI200I,XWli2CCI.XW2(ZOOI.XW312001. 5AI<E6410
1 1H ,lOX,'BASE SI-EAR IIliIN, :',E1303,2X,A4,1. SAKE6170 6 FDl(200).FD2(200I,FWl(2001. SAKE6475
Z 1HO,10X.'BASE MOMENT MAXe =',E13.3,2X.A4,'-',A4,1. SAKE6175 7 FRM(4001.0RM(4001,VR~(4001,WRMI400) SAI<E6480
3 1H .10X,'RASf ~CMENT ~IN. :",E13.3,2X,A4.'-',A4,/& SAKE6180 SAKE64B5
616 FORMAT (lHl, II.IH ,10X,'MAXIMUM RESPONSE IN FRAME',IZI SAKF6185 HVSTEPESIS RULE OF T~E BIlI~EAR JJINT SPRING (MODIFIED TAKEDA) SAKE64QO
620 FORP-lAT (lH II.IH .20X.'~AXIMUM FLEXURAL RESPJNSE OF MEMBERS',I, SAKF6190 SAKE64Q5
1 1H ,20X,'6S A UNIT-lENGTH SIMPLE MEMBER",I. S AKE61 95 KK=KE+KE-2 SAKE6500
Z IH .20X.'(ElASTIC DEFORMATION INCLUDEDI",II, SAKE6200 "'F=FZ(K:=)+DF SAKE6505
3 1H .T22,MO~ENT',T44.~OTATION""T61.'DUCTILITV'.I. SA~E6Z05 I VL=KVU KE I SAKE6510
4- 1 H T 1. q " , A4. "*", A4. I T 4 3. " 1 0 I' A4. I I , SAKE6210 JIIL= IVL SAI<E6515
5 IH "T2,' MEMBER' T16, ""'IIX' T 28,' MIN' , T'I>O. MAX'. T56, "MIN'. SAKE6215 IF IFRIKEI.LEoO.OI GO TO 11000 SAKE6520
6 T64,MAXD.T16,'~IN',II) SAKE6220 TST=ISTF2(KEI SAK E6525
62Z FORMAT (lH ,Y3,T6.'~-END.lP8EI2e3./1 SAKE6225 IS(;N=2 SAKE6530
f:24 FORMAT (lH ,T6,'B-END',IP8E1203.11 SAKE6230 IF (F2IKEI.LT.0.01 ISG~=l SAI<E6535
6Z6 FORMAT (lH .III,lH .20X.'MAXIMUM JOINT SPRING RESPONSE',II. SAKE6235 rD=FP(KEI*DF SAI<E6540 f-J
1 IH .T22.DMCMENT",T44,"ROTATION",T67,"DUCTILITV,I, SAKE6240 W2IKF)=W2IKEi+DD SAI([6545 W
2 IH TIQ,'" ,A4. ">II' ,A4,')' ,T43.' . ,A4, 'I' .1114.'" .1. SA~E6245 GO Tn 11,Z,3,4.5,6,7,8,9.10,lll.IIIL SAKE6550 f-J
3 IH ,T2, 'i14EMBER ",T16. 'MAX' ,T28, 'MIN" ,T40,'I-1A)(' ,TS2. "MIN'. SAKE6250 c ELASTIC STAGE SI\KE6555
4 T6'1>.MAX.T76,~I~.111 SAKE6255 1 IF (IIY(ISTI-ABSIFFII 100.100.10000 SAKE6560
END SAKE6260 100 VRM(KK+IJ=-VVIISTI SAKE6565
VRMIKK+2'= VYIISTI SI\KE6570
WRM(KK+l'=-wYIYSTI 5AKE6575
WRM(KK+21= ~V(YST) 51\K[6580
110 KVU KE 1=2 SAI(E6585
FR (ICE' )=PSP( ISTI SAI<E6590
FF=C FRC (VV! 1ST) WY I S1I PSP IS TI .W2 KE I I SAI<E6595
Gn TO 10000 SAKE6600
SUElRClUTINE I-lVSTRS (CF,KE.IFlAG i SAKE6300 C lOADI~G ON THE PPIMARY C~RVE AFTER YIELDING SAKE6605
C SAI(r::6305 2 IF (DD*F2!KEIJ 200.200.1COCC SAKE6610
C HYSTERESIS RULES FOR TWO SPRINGS IN SERIES. SAKE6310 200 VRMIKK+ISGNI=FZ!KEI SAKE6615
C (1) FlEXLRAl ROTATIClNAl SPRING (TAKEDA MODEll SAKE631S WRMIK~+ISG~)=COSPIVVIISTI,WY(ISTI.PSPIISTI.F2IKEII SAKE6620
C IZI BILINEAR RQTATIONAl SPRING (MODIFIED TAKEDA MODEll SAKE6320 SN=WRM(I<K+ZI SME6625
C THE TW0 SPRINGS HAVE T~E SAME FORCE LEVEL SAKE6325 IF S~+WRr-'(KK.111.LT.0.01 SN=-WRM(KK+ll SAKE6630
C SAKE6330 FWU KE I=f SJ' II ST 1* (SN/IoV lIST) I ~*BT3( 1ST) SAKE6635
C FOLLOWING PROCEDURES ARE TAKEN IN THIS PROGRAM SAKE6335 KVUKFI=3 SAKE6640
C III GlvE~ INCREME~TAl FORCE rF. CALCULATE INCREMENTAL ROTATION SAKE6340 FRIK'OI=Hl(l<El SAKE6645
C !N T~E AILINEAR JOINT SPRING BASED ON THE PREVIOUS STIFFNESS. SAKE6345 FF=VRM(KK+ISGNI+(W2!KEI-WPM(KK+ISGNII/FRIKEI SAKE6650
C (21 APPL' T~~ HYSTERESIS PULES (MODIFIED TA~EDAI TO FIND NEW SAKE6350 IF IFF~F2(KE" 310,310,10000 SAKE6655
C STIFF~ESS AND INCREMENTAL FCRCE OF JOINT SPRING BASED ON THE SAKE6355 C UNLOAr.ING AFTER YIELDING SAKE6660
C IN(pEME~TAL ROTATION GIVEN IN (lB. SAKE6360 3 IF (FFaF2IKEll 310.310.300 SAKE666'5
C (3) CALCULAT~ INCRE~ENTAL ROTATION IN THE FLEXURAL SPRING FOR THE SAKE6365 300 IF WRpI(KK+ISGNI-W2IKEII*FFI 110,110,10000 SAKE6670
C NEW FORCE LEVEL (21 USING PREVIOUS STIFFNESS, SAKE6370 310 XWO(KEI=wRM(KK+ISGNI-VP~(KK+ISGNI.FR(KEI SAKE6615
c 141 APPLY T~E HYSTERESIS RUl~S (TAKEDA' TO FIND NEW STIFFNESS SAI<E6315 rSGN=3-1 SGN SAKE6680
r AND FQRCE LEVEL OF FLEXUPAL SPRING BASED 0N THE INCREMENTAL SAKF6380 320 n: ((WR/'IIKK+ISGNI-W2(KEII"'FFI 110,110.:no SAKE6685
c PCTATICN (31. SAKE6385 330 KVU KE =4 SAKF6690
C lSI IF T~E FORCE lEVEL IS CHANGED IN STEP (4), APPLV THE SAKE6390 FR(KEI=IWRMIK~+ISr,NI-XWOIKEII/VRM(KK+ISGNI SAKE6695
C HYSTERESIS RULES (MOOIFIED TAKEDAI TO J~INT SPRING TO FIND THE SAKE6395 FF=( W2( KE )-XWO (KEII /FR KE) SAKE6100
C Nfw STIFFNESS AND ROTATICN 8V CHASING THE STRESS HISTORY FROM SAKE6400 (;0 TO HCOO SAKE6105
THE ORIGINAL STRESS LEVEL AT STEP (11. SAKE6405 c LOAO REVERSED FPOM ST~P 3, AND LOADING TOWARD POIIIlT IVRM,WPMI CSAKE61l0
( YH~ p~r~APY (URVE 5AKE6715 FRIKFI=f!,lIKEI SAKF702C
4 IF IDO*~2(1(f) I 400,400.300 SAKE6720 FF=VP3IK~I.(W2IKEI-WR3(KE)I/FR(KEI SAK"'7025
400 VRO(Kfl=F21~F.1 SAKE6725 1010 IFIFF*F2(KEI! 1100,1100.10000 SAKE7030
WR0(KFI=XWOIKEI+f2IKEI*FR(KEl SA'<E6730 ( UNLQADING FRC~ l~VEl 10, peINT IVR3,WR3) SAKE7035
KVlIKFI=5 SAKE6735 11 IF (FF*F2(KEI) 1100.1100.1110 SAI<E1040
FR(I<E'I=Fl1l1ICEI SAKE6740 1100 XW2IKEI=~R3(K~I-VR3IKEI*FR(Kf.1 SAKE7045
FF=VR0IKEI+(~2IKED-WR01KEII/FR(KE) 5AKE6145 G0 rr 710 SIU(E1050
410 IF IFF*F2(KEII 510,510,10000 5AKE6750 1110 IF Wfl3(KEI-W2(KElI*FFI 910.910.10000 SAKE7055
( UNlQACING FRG~ STEP 4, PCI~T (VRO,WROI SIIKE6755 10000 VV=FF SAKE7060
5 IF (FF*f2(KEII 510,510.500 SAKE6760 C SAKE7065
500 IF (WROIKFI-W2(KE II*FF I 320.320.10000 SIII<E6765 C ~YSTEPESIS RULE OF THE FLEXURAL ROTATIONAL SPRING ITAKEDAI SAKE7070
510 ISGN=3-ISGN SAKF6710 C SAKE1C75
XWl(KE)=~RO(KE)-VRO(KE$FR(KEI SAKE6175 11000 IST=ISTFIIKE) SAKE7080
520 IF WRMIKK+ISGNI-W2(1(EII!IkFFI 1l0 . HO.530 SAKE6760 ISG"':2 SAKE708'5
530 KVL( KE I =6 SAKE6785 IF (F21I<EI.LToOoOJ ISGN=1 SAKE1090
FP iKE 1= (WPf'II (I<K+ ISGN 1-)(1141 KE II/VP Me KK+ I SGN I SAKE6790 FDD=FCIKEI.FCCIIST) SAKE7095
FF=(W2!KEI-XWIIKEII/FRIKED SAKF6795 9D=FDO(FF-F2(KEII SAKE7100
r,C T(1 1COOO SAKE6800 D2(KEI=C2(KFI+OD SAKE7105
C LOAD REVERSED FROM STEP 5.AIliD LOADING TOWARD POINT VRM, WRMI SAKE6A05 IVL=LVUKEI SAKE1110
C ON THE PRIMARY CURVE. SAKE6810 GC TO 121,22.23.24.25.26.27,28.29,30,31.32.33,34,35,361.IVL SAKE7115
6 IF (DD*F2(KEII 60C,600,300 SAKf6815 ( ELASTIC TO FORMATION CF (RACKS IN ONE DIRECTION SAKE7120 (
600 VR1IKEI=F2(KEI SAKE6620 21 IF IDC(IST)-ABSID2(KEII) 210(.210(,20000 SAKE1125
WR1(KFI=XWl(KEI+F2(KEI*FRIKEI SIIKE6825 2100 FDl(K~I=FDCIISTI SAKE7l30
KVUKF )=7 SAKE6630 FRM(KK+ll=O.O SAKE7135
FR II<E I=Fld IKE SAKE6835 FP~IKK21=OGO SAKf1140
FF=VRL(Kfl+(W2(KEI-WR1!KEII/FR(KEI SAKE6840 OP.,.. I KK+ 1) =0 eO SAI<1:1145
if IFF*F2IKEI) 100,100.10000 SAKE6845 I)RMI KK+21=0.0 SAKEH50
C UNLOACYNG FRO~ STEP 6. POINT VR1.WRll SAKE6850 2110 IF (DYIISTI-ABSID2IKEIII 2210.2210.2120 SAKE7155
7 IF (FF*F2(KEI 700.70(1.740 SAKE6655 2120 LVU KE I =2 SIIKE7160 f---J
700 XW2(KF)=wR1dKEI-VRl(KEI*FR(KE) LV
SIIKE6860 FCC=FrVI 1ST SAKE7ll.65 [\)
no I SG"=3-I SGN SAKE6865 FF=CFQCICM(ISTI,DCIISTI.FDYIISTl.D2(KE)1 SAKE7110
720 IF WRO(KEI-W2(KEDD*fF) 320.320.730 SAKE6810 r,0 TO leOCa SAKE7175
730 KVUKEI=6 SAKE6875 C LOADING AFTER CRACKING SAKE7lS0
FR(KEI=~PO(KE)-XW2(KE/VRO(I<EI SAI(E6880 22 IF (DD~F2(KEI) 2230.2230,2200 SAKE7185
FF=(W2(KE)-XW2(KEtD/FR(KEI SAKEMl85 220~ IF (DYIISTI-ABS(02(KE'DI 2210.221C.20000 SAKE7190
GO Tr 10000 SAKE6890 2210 lVUKEI=3 SAKE1Jl95
740 IF (iWR1(KEI-W2(KEI'*Ff' 520.520,10000 SAKE6895 s:rO=FOUII5TJ SAKE7200
C LnaD REVERSED FROM STEP 70 AND LOaOYNG TOWARD POINT (VRO.WRO) SAKE6900 FF=CFRC(YMIISrD.DYI!STI.FDUIISTD.D2IKEII SAKE7205
8 IF rO*F2(KEI) aOC.BOO.500 SAI<E6905 GO TO 20000 SAKE1210
80r VR2(KEI=f2(KEI SAKE6911''I 2230 DPMIKK+ISGNI=CDSPICMIISTI,nCIISTI,FDV(ISTI.F2(KEI) SAKE7215
WR2(~FI=XW2(KEI+F2(I<E)$FRKE SAKC::6915 FR~KK+ISGNI=F2(KEl SAKE7220
I!VUKEI=9 SAKF6920 FC1 ( KE I = AP. S ID i< MKK+ IS f N I C( I 1ST I III ABS F2 KE I I +U~ ( 1ST I I SAKE7225
FQ(Kfl=FWUKfl 5111(1"6925 lVUKE)=5 SAKE1?30
FF=VR2(KEI+(W2(KEI-WR2(KEID/FR(KEI SAKE6'BO 2240 Frry=fr.l(~E) SAKE1235
IF (FF*FZ(KEI' 900,900.10000 SAKE6915 FF=FR~(~K+ISGNI+(02(KEi-OR~(KK~ISGNII/FDD 5AKE7240
C UNLOADING FROM STEP 8. POINT (VR2.WR21 5AKI:6940 GO TI'1 20000 5AKE1245
9 IF fF"'F2i Kf I 900,90.0,930 SIII(E6945 C AFTER YIELDING SAI<f72'50
900 X~3(KF'=~R2(KEI-VR2(KEI.FP(I!E SAKE6950 23 IF COD*F2IKEII 2300.2300.200CO SAKF7255
YSGf\J=3-ISGN SAKE6955 2300 FPM(KKvISGN1=F2(KEI SAI<E1260
IF IABS(CFioLToABS(FF-F2IKEIII FF=F2KE+DF SAKE6960 ORM(KKvISGNI=COSP(YMIISTI.DYIISTI.FDU(ISTI.F2(KEII SAKE1265
910 IF (IWR1(KEI-W2(K~)D$FFI 520,520,920 5AKE6965 SN=O~~(1<1<+21 SAKE7270
920 KVUKEI=lO SAKE6970 IF SN+OR~(KI<+llloLToO.OI S",=-DRM(KK+ll SAKE7275
FRKED=~Pl(KEi-XW3IKEII/VR1(KEt 511KE6975 FDUKI5I=(SN/DV(ISTI ."'!IkBTO! ISTI SAKE1280
FF=(~2(KE)-XW3(KE)I/FR(KEi SAK.E6980 1 lIIf)CIISTl+DY(lSTIII(CM(ISTB+VMIISTII SAKE7285
GO TO 10000 SAKE'6985 LVUKEI=4 SAKn290
930 I~ I(WR2(KEI-W2(KEli~FF) 720,720.10000 SAKE6990 GO TC 2240 SAI<E12'95
( LeAD REV~RSEn FROM STFP 9. AND LOADING TOWARD POINT (VR1.WRll SIII<E6995 C UNLOADING AFTER YIELDING FR(M POINT (FRM.DRMI SAKE7300
10 IF I DD*F2 KE I 1000.' 00. HO SAKE7000 24 IF (fF*F2(KEII 2410,2410.2400 SAKE7305
1000 VR3(K~I=F2(KEI SAKE7005 2400 IF DRM(KK+ISGN1-D2(ICEII*FFI 2210.2210.20000 SAKE7310
WP3{KEI=XW3P<EI+F (KEI SAI<=7010 2410 ISGN=3-ISG~ SAKF7315
I<VL (K~ I = 11 SAKE7015 If (CM(IST)-ABS(FR~(KK+ISGNlll 2450,24'50.2420 SAKE7320
 242~ J~ (LVLIKE'-4) 2440.244C,2430 SAKE7325 fF=(D2(KFI-XC2!K!:II/FCD SAKE1630
2430 LVL( Kr: I =14 SAKE'7330 G0 T0 20000 SAI(~7635
GO TO 3410 S tlKE 1335 2940 IF DP1(KEI-D2(KEII*FFI 273C.273C,20000 SAKF7640
2441) l VU KE I =15 SAKE7340 C LeAD PEVERSED FROM STEP 9, aND LOADING TOWARD POINT (FRO,DROI SAKE7645
GO Tf' 3500 SAKEH45 30 IF (DD*F1IK!:)) 3000.30(0.2700 SAK(7651)
2450 XDO(KFI=OP~(KK+3-ISGNI-FP~(KK+3-ISGNI.FOD 5AKQ350 3000 FR2IKfl=F2IKFI SAKE7655
FDX=(ORMIKK+ISGNI-XCO(K~llfFR~(KK+ISGN) SAKE1355 DP2(KEI=XD2(KE)+F2IKEI*FDD SAKE76b0
IF (LVLlKEI.GT.41 GO TO 24eO SAKEi'360 LVL(K~)=11 SAKE7665
IF (I)Y(1STJ.LT.ABStORfJ(KK+ISGNIII GrJ TO 2460 SAKf7365 IF (FDC.GT.FOIIKEI*BT2(ISTII FDD=FDIIKE'*BT2(ISTI SAKE7670
FC'Z= IDY I I ST I +ABS XOO I KE IIIIYM( 1ST I SAKE7370 3010 FF=FP2(KEI+(D2(KEI-DR2(KEII/FDD SAKE7b75
IF (FoX.LT.FDZI GO TO 2460 SAKf1'375 IF (FF~F2(KEII 3100,3100.2CCOO SAKE7680
Ft:12(K!:I=FDI SAKE1380 C UNLOACING FR~M STEP 10. POINT (FR2.DR21 SAKE1685
SN=ISGN+ISGN-3 SAKF.73B5 31 11= (FF*F2IKEII 3100.3100.3130 SAKE1690
FRMIKK+ISGNl=SN*YMIISTI SAK=7390 3100 xn3IKEI=OR2(KEI-FR2IKEI*FDD SAKE7695
ORMIKK+ISGNI=SN*DYIISTI SAKEH95 ISGN=3-ISGN SAKE1700
GO TO 2470 SAKE1400 3110 IF (CDR1IKEI-D2(KEII*FFI 2730,2730.3120 SAKE1705
2460 FC2(K!:I=FDX SAKE7405 3120 LVLIKEI=12 SAKE7710
2470 LVUKEI=6 SAKE7410 FDO=IDRIIKEI-XC3IKEII/FRIIKEI SAKE7715
l=oD:FC2IKEI SAKE7415 FF=ID2IKFI-XD3(KEI"FDD SAKE7720
FF=(02(KEI-XDO(KEII/FOD SAKf7420 GQ ~~ 20000 SAKE7725
GO TO 20000 SAKE'7425 ~130 IF IIOP2(KEI-D2IKEllFFI 2920.2920.20000 SAI<E7731)
C UNLOADING FROM STEP 2. P~I~T (FR~,DRMI SAKE7430 C LrAD REVFRS=O ~RO~ STEP 11, AND LOADING TOWARD POINT IFPl.OR11 SAKE7735
25 IF IFF*F2IKF" 2410,2410,2500 SAKE7435 32 IF (OD*F2IKE)1 3200.3200,2Q40 SAKE1740
2500 IF (IDPMIKK+ISGNI-C2(KEII*FFI 2510,2510,20000 SAKE7440 ~200 FR3(KEI=F2(KEI SAKE7145
251.0 IF IDYIISTI-ABSIDRM(KK+ISGN 2210.2210.2120 SAKE7445 DP3IKFI=XD3(KEI+F2(KEI*FDD SAKE7750
C LeAD REVERSED FROM STEP 4 OR 5, AND LOADING TOWARD POiNT SAKF11450 LVLIKEI=13 SAKE7155
C IFR~.DR~1 CF T~E PRIMARY CURVE. SAKE1455 IF IFDO.GT.FD1(KEI*~T2IISr) FOD=FDIIKE,$BT2IISTI SAKE7160
26 IF IDD$F2(KE 2600.2600.25)0 SAKf74bO 3210 FF=FR3IKEi+(02IKE'-OR3(~EIB/FDC SAKE7765
2600 I=RO(KEI:F2(KEI SAKE11465 IF (FF*F2IKEll 3300.3300,20000 SAKE7770 I--'
DROIKEl=XDO(KEI+F2IKE,$FOD SAKE14l'O C UNLOADING FRaM LEVEL 12, POINT (FR3,VR3' SAKE7775 UJ
LVU KE )=1 SAKEH15 33 IF IFF.F2IKEII 3300.330C,~310 SAKE1780 W
IF (Frr.GT.F01(KE)$BT1IIST)1 FDD=FD1IKE'$BT1IISTI SAKE1480 3300 X02IKEI:DR3IKEI-fR3IKEI*FDO SAKE1785
2610 FF=FROIKEI+ID2IKE'-DROIKED)/FDD SAKEH85 GO TO 2910 SAKE1790
IF (FF*F2(KEJD 2120.2720.2COCO SAKElI490 3310 IF DP3(KEI-02(KEII*FFI 3110,3110,20000 SAKE7795
C UNLOADING FRQM STEP 6. POINT (FRO,DROi SAKE74i95 C LeAD REVERSED FROM STEP 5, AND LOADING TOWARD CRACKING POINT SAKE7800
27 IF (FF*F2IKED) 2720.2720.2700 SAKE1500 34 IF (FF$F2(KEII 3400.3400.3410 SAKE7S05
2700 IF (<<OROIKEI-C2IKE*FFI 211C,271C.20000 SAKE7505 3400 LVL(KEI=5 SAKE1810
2710 IF (DR~(KK+ISGN'-r2KE"$FF 2510.2510.2470 SAKE1510 ISG~=3-ISGN SAKE7815
2720 ISGN=3-ISGN SAKE1515 GO TQ 2500 SAKE78l0
XCl(KEI=CRO(KEI-FROIKEI*FDD SAKE7520 3410 IF (CMIISTI-ABS(FFII 2110,2110.20000 SAKE7825
2730 IF DRMIKK+ISGNI-DZ(KEl'*FF) 2510.2510.2740 SAKE1525 C LOAD REVERSED FROM STEP 4 AFTER YIELDING IN ONLY ONE DIRECTION SAKE7830
2740 LVLU<E)=8 SAKE7530 35 IF (FF*F2(KEII 3520.3~20.350C SAKE7835
FDD= (DIHUKK+ ISGN i-XOi KE ;/FIHI! (KK+l SGNI SAKF7535 3500 IF (CMIIST)-ABSIFFII 3510.3510.20000 SAKE7840
FF=ID2IKEI-XD1IKEDlfFDD SAKE1540 3510 SN=FF/ABSIFFI . SAKE7B45
G8 TIJ 20000 SAKE1545 FDX=DRM(KK+3-ISGN)+(ARS(~PM(KK+3-ISGNII+CMIISTI )*SN*FOD SAKE7850
C LeAD ~EVERSEO FROM STEP 7, AND LOADING TOWARD POINT IFRM,DRM) SAKE7550 LVLIKEI=16 SAKE7855
C eN THE PRIMARY CURVE. SAKE7555 FDo=(rY(ISTI-FOX*SN)/(YMIISTI-C~(ISTII SAKE7860
28 IF IDD$F2CKE'J 2800,280C.250C SAKE7560 I=F=SN*C~(ISTI+(D2(KEI-FDXI/FCC SAKE7865
2800 FPIIKF':F2IKE) SAKE"1565 GOTO 2COOO SAKf7870
rRIIKE):XD1IKEi+F2(KEI*FDD SAKE7570 3520 LVLIKEI=4 SAKE7875
LVUKEI=<J SAKE7575 ISGN=3-ISGN SAKE7880
TF (FOO.GT.FDIIKEI*BT1IIST" FDD=FD1IKEI*BTIIISTI SAKE7580 GO TO 2400 SAKE7885
2810 FF=FRIIKE1+ID2IKI=I-ORIIKEII/FDD SAKE75B5 C LOADING TOWARD YIELDING POINT AFTER STEP 15 SAKE7890
Y~ (FF$f2( I<E II 2900.2900, 200CO SAKE7590 36 IF IOD*F2(KEII 3600.3600,2200 SAKE7895
C UNL'JACING FROid LEVel B, PCINT (FR1.DRlI SAKF7595 3600 SN=ISGN+ISGN-3 SAKE7900
29 11= IFF*F2IKEII 2900,290C.294C SAKE7600 FRMIKK+ISGNI=SN$YM(ISTI SAKE7905
2900 XD2IKEI=CR1(KEI-FPIIKEI*FDD SAKE7605 DRM(KK+ISGNI:SN*DY(ISTI SAKE7910
2910 ISGN=3-ISGN SAKE7610 FC2(KE'=FD~ SAKE7915
2920 IF (I)RC IK EI-C2 (KEII*FF 271 0,271 C.2930 SAKE7615 XOO{KEI=SN*(CYIISTI-YMIISTI$FDDI SAKE7920
2930 LVl(KEI=lO SAKE7620 FROIKFI=F2(KEI SAKE7925
F~D=IDRO(KFI-XD2IK=)I/FRO(KEI SAKE7625 DPO(KF):XCO(KE)+F2(KEl*F02(KEl SAKE7930
 LVl(l<E '=7 SAKf.=7935 4510 Dfl=FF-VRQ (I( E I SAKE8240
FDO=~r1!KEi~BTl[ISTI SAKE7940 r;(l T'J 4411) SAKE8245
TO 2610
r,'] <;AKF7945 4520 XWl(KEI=WROIKEI-VROIKEI*FWIII<EJ SAKE8250
20000 FD(KEI=FCD-FDC(IST) SAK~7950 ISGN=3-ISGN SAKI:8255
IF (IVLG~EaLVlIKEI IFLAG=IFlAG+l SAKE7955 10530 IF VRMCKK+ISGN)-FF'$FFI 41CO,4100,4540 SAKE8260
C SAKt7960 4540 KVUI<EI=6 SAI<E8265
C BACK TO T~E HVSTERESIS RULE FOR THE BILINEAR JOINT SPRING. SAKE7965 FQIKEI=(~RMCKK+ISGN'~XWl(I<EJJ/~RMCKK+ISGNI SAKE8270
C SAl<E7970 W2(KFI=XWl(KEI+FRIKE,*FF SAI<E8Z75
IF (FRI~EI.LF.O.O GO TO 31000 SAKE7975 GO TO 31000 SAKE8280
IF (VV.EQ.FFD GO TO 31~OO SAKE7980 C LeAD IlEVERSED FROM STEP 5 A~D LOADING TOWARD POINT (VRM,WRMJ SAKf8285
IST= ISTF2(1<E SAKE7985 46 IF IOD*F2IKEJI 460C,4600,4530 SAKE8290
! SG\\=2 SAKE7990 4600 VR1(KEI=F2iK~) 5AI<E8295
IF (F2(KEB.LE.0.01 ISGN=l SAKE7'995 WR1(KEI=XW1(KEJ+VRIIKE'*(WRM(KI<+ISGN'-XW1IKEllfVRM(KK+ISGNt SAKE8300
DD=FF-F2(KEJ SAKE8000 IF (1=2(KEI*fFl 4700,4700,4610 SAI<E8305
GO Te 141.42.43,44.45.46.47.48,49.50.51I,JVL SIU<E8005 4610 KVUl<fl=7 SAI<E83LO
C ELASTIC STEP SAKE8010 FRlKEI=FwilKEI SAKE8315
41 IF lVYIISTI-ABSIFF 4100,4100.4110 SAI<E8015 W2IKF)=WRIIKfl+OD*FR(KEI SII.KE8320
4100 I<VUKE 1=2 SAI<E8020 GO TC 31000 SAKE8325
FR(I(EJ=PSPIISn SAKE8025 C UNLOADING FRQM STEP 6, POINT (VPl.WRll SAKE8330
W2(Kfl=CDSP (VYIISTD.WYIISTI.PSPlISTI.FFI SAI<E8030 47 IF (F2(I<EI*FFI 4700.4700,4730 SAKE8335
GO Te 31000 SAKFBC35 4700 XWZ(KEI=~pICKfl-VRl(KEI.FWI(KEI SAKE8340
4110 KVUKEJ=1 SAI(FR040 ISGN=3-YSGf\I SAKE8345
I=RIKEI=Espusn 5AI<E8045 4710 iF (IVROIKEI-FFI*FFI 4330.4330,4720 SAKE8350
W2IKEi=FRIKE'$FF SAKE8C50 4720 KVUKEI=S SAKE8355
GO Tr:." 31000 SAKEB055 FRIKEt= (WROIKE)-XWZIKE II/VRO! I<E I SAKE8360
C LrADI~G CN THE PRIMARY CURVE AFTER YIELDING SAKE8060 W2(KEI=XW2(KEI+FRIKEI*FF SAKE8365
42 IF (OO$F2IKEB' 4200.42fO,4100 SAKEB065 GO TO 31000 SAKE8370
4200 VRM(K~+YSGN'=FZ(KE) SAKE8070 4730 IF VP1(I<EI-FFI*FFI 4530,45~0.4740 SAKE8375 I........
WR~(~K+ISGN)=CDSP(VYIIST).WY(IST'.PSPISTI.F2(KE SAKE8015 4740 OD=FF-VR1(KEI SAKE8380 W
SN=WRM(KK+2) SAKE8080 GO TO 4610 SAKlE8385 ..j::"""
IF SN+WPMIKK+ll!.LTaO.O) SN=-WRMIKK+l1 SAI<E80fll5 C LOAD REVfRSEO FROM STEP 7, aND lOADING TOWARD POINT (VRO,WROI SAKE8390
FWliKED=ESPiIST'$'SN/WV(ISTJ'*ST3IISTI SAKE8090 48 IF (OD$F2(KEI) 4800,4800.4110 SAKE8395
iF IF21~E FFB 43Z0.4~ZO.421C SAKE80'95 4800 V~2KEI=F2CKE) SAKE8400
4210 KVUKEt=3 SAKE8100 WR2(KEI=XW2(KE'+VR2IKEI*WRO(KE'-XW2(KEIDfVROiKEI SAKE8Itt05
FQ(I<EI=FHCKEI SAKE8105 IF IF2(KEI.fFI 4900,4900.4810 SAKE8410
W2(KFI=W~M(KK+ISGN'+DD$FWlI(E) SAKE8HO 4810 KVUKEI=9 SAKE8415
GO TO 31000 SAKi=8U5 FP(KEI=FIIIUKEI SAI<E8420
C U~LOADI~G FROM POINT VR~,WRM) ON THE PRIMARY CURVE SAKERI2f.l W2(Kf)=~P2KED+OD*fR(KE SAKE8425
43 IF (F2IKED*FFD 4320.4~20.~30C SAKE8125 GO TO 31000 SAKE8430
4300 IF VRM(KK+!SGND-fFD$FFD 4100.4100.~310 SAKE8130 C UNLOADING FROM STEP 8, prINT (VR2,WRZI SAKE8435
4310 DD=FF-VRM(KK+YSGNJ SAKE8135 49 IF (1=2(Kf)*FFI 4900,4900,4930 SAKE8440
GO T'J 4210 SAKE8140 4900 XW3IKfl=kR2(KEI-VR2IKEI*FWl(I<EI SAI(E8445
4320 XWO(KEJ=~RM(I(K+ISGND-VR~KK+ISGND*FWICKEJ SAKE8145 ISGI\=3-ISGN SAKE8450
ISGI'.l=3-ISGN SAKEIH50 4910 IF VRIIKE'-FFI$FFB 4530,4530,4920 SAKE8455
4330 IF VR~(KK+ISGN-FF)*FFJ 4100,4100.4340 SAKE8155 4920 KVU KE 1=10 SAKE8460
10340 KVUKED=4 SAKE816/) FR(KEI=(WRlII<EI-XW3(KE)fVRI(I<EI SAKE8465
~R(KE=(WRMKK+ISGND-XWO(KEJ/VR~KK+ISGNJ SAKE8165 W2(KEI=XW3(KEI+fRtKEI*fF SAKE8470
W2CKEI=XWO(KE.+F~(KEI$FF SAKE811'0 GO TG 31000 SAKE8475
GO T~
31000 SAKE8115 4930 IF (VR2(KEI-FF,$FFI 4710,4710,4940 SAKE8480
C LeAD REVERSED FROM STEP 3 AND LOADING TOWARD POINT 'VRM.WR~) SAKE8180 4940 OO=FF-VR2(KEI SAI<E8485
44 IF (DD*F2(I<ED) 4400.4400.4330 5AKE8185 GO Tn 4810 SAKE8490
4400 VRO(I(EJ=F2(I(EJ SAKE81'90 C LeAD REVERSED FROM STEP 9, ~ND LOADI~G TOWARD POINT IVR1,WR11 SAKE8495
WRO(KED=~WO(KED+VROIKEI*(~MCKK+ISGND-XWO(KE'DfVRMIKK+ISGN! SAKEB195 50 IF (DO*F2( I(EI I 5000,5000,4910 SAKE8500
IF (FZ(KED*Ffl 4520,4520.4410 SAKE8200 5000 VR3(KEI=F2(I<EI SAKE8505
44l1'1 I(VUKF 1=5 SAI<E8205 W~31I<EI=XW3(KEJ+VR3(KEI*(WRl(KEi-XW3(KE)lfVR1(KEI SAKE851.0
FR IKE D=FWlC KE SAKE 8210 iF (F2IKEI*FFI 5100,5100.5010 SAK!"8515
W2(KE)=WRO(K~D+DO$I=Wl(KEJ SAKE8215 5010 KVUI<EI=U SAKE8520
GO TO 3100fl SAKE82Z0 FRIKEJ=FWUKEI SAI<E8525
C UNL~ACING FROM POINT(VRO.WROI SAKE8225 WZ(I<EI=WR3(l<fl+DD*FW1(KEI
45 IF IFZIKED*FF, 4520.4520,4500 SII.KE8530
SAKF8230 GO Tf") 31('('0 SAKE8535
4500 IF VRO(~Ei-FF)*FF 4330.4330.4510 SAKE8235 C UNLOADING FROM STEP 10, A~C POINT (VR3,WR31 SAKEB540
 51 IF IF2(I<EI$FFI 5100.51CO.5110 511KE8545
5100 XW2(KFI=WR3(KEI-VR3(KEI$FWlII<EI SAKE8550
tSGN=3-ISGN SAKF8S55
GO T(1 4711) SAKEfl561l
5110 IF (VR3(KEI-FFI$FFI 4910.4910.5120 SAKE8565 SUBROUTINE SPRDG IA.B.C,~,M) SAI<E8B10
5120 nO=FF-VR3(KEI SAKE8510 C (CI=(AI.IBI SAI<E8815
GO TO 5010 SAKE8575 C A(N,N) SY~METRIC MATRI~ (STORAGE MODE=SYMMETRICI SAI<E8880
31000 F2(KfI=FF SAKE8580 C BIN.~) GENERAL MATRIX SI\KE8fl85
IF (JVL.NE.KVlIKEII IFLAG=IFlAG+l 5AKE8585 C (N.M) PESULTANT MATRIx SAKE88'90
RFTURN SAKE8590 DIMENSION A(II,Blll,CI1) SAKE88'95
END 5AKE8595 DOURLE PRECISICN A,B,C SAKEB'900
DO 40 l=l,N SAKE8'905
IJ=I-N SAI<E8910
f(J=O SAKE8915
f)O 40 J=l,fiII SAKE8920
IJ=IJ+N SAKE8925
FUNCTI1N CFRC(CRF~CRO,SLP.DSP) SAKE8625 rtIJI=O.O SAKEB'930
SN=OSP/1I95ICSPD SAKE8630 DO 30 K=l,N SAKE8935
C~RC=(OSP-SN(RD)'5lP+SN.(RF SAKEB6~5 KJ=KJ+l SAKE8'940
RfTURN 511KE8640 IF 11-KI 10.10.20 SAI<E8945
EIllD SAl<EB645 10 IK=f($!K-l1/2+J SAKE8950
r.0 TO 30 SAKE8955
20 IK=I*(1-11/2+K 5A1<E8960
30 CIIJI=C(IJI+A(IKI$SIKJI SAKE8965
40 CONTINUE SAKE8910
RETURN SAKE8915
FUNCTIQN CDSP(CRF,CRD,SLP.FRCI 5AKE8615 END 5AKEB980
SN=FRCI AeS FRC 5AKE8680 I-'
CDSP=(FRC-SN.CRF)$SLP+SN*CRO 5AKE8685 W
RETURN SAK':8690 \Jl
END SAKE8695

SUBROUTINE SYBAN (A,B,NROW,NCCl.NRH51 SAKE9020

C SAKE9025
C SOLUTIJN OF LINEAR EQUATION, IN WHICH T~E COEFFICIENT MATRIX IS SAKE9030
C SYMMETRIC AND BANDED. DEVELOPED BY R. E. MILLER (U. OF 1.1. SAKE9035
SUBROUTINE GMPRD{A.B.R.N.~.l) SAKE8135 C AINROW,NCCLI COEFFICIENT MATRIX SAKE9040
C RC*'=A*'$B($D SAKE8740 C T~E PRINCIPAL DIAGONAL IS STORED IN THE SAKE9045
DIMENSION ACID.fI[l.R(li 5AKE8745 C ~IRST COLUMN AI*.*I. AND THE NONZERO BANDS SAKE9050
nOUBlE PRECISION A,S,R SAKE8150 C ABaVE T~E PRINCIPAL DIAGONAL F1RM THE SAKE9055
c SAKE8155 C BALA.NCE IJF A '.,-). SAKE9060
IR=O SAKE8160 r NROW NUMBER CF EQUATIONS SAKE9065
IK=-M SAKE8765 C NCQL BAND wIDTH OF THE COEFFICIENT ~ATRIX 5AK[9010
Of1 10 K=l.L SAKE8110 c B(NROW,NRH51 PIGHT ~AN9 SIDE C~LUMN VECTORS SAKE9015
IK=Yf<+M SAKE871'5 C UPON RETURN FR~M THIS PROGRAM, THE SAKf9080
01'1 10 J=l,N SAI<E8180 C SIJLUTIC~ ARE STQRED IN B(*.*I 5AKE9085
IP=IIl+l SAI(E8185 C NRHS NUMBER OF RIGHT HAND SIDE COLU~N VECTORS. 5AKE9090
JI=J-N SAKE8190 IMPLICIT REALlIrB(A-H.O-l) SAKE9095
IB=II< 51\KE8195 nI~ENSICN A(NRryW.NCOLI,B(NRCW,NRHS SAKE9100
fl(IRI=O.O SAKFBBOO IF (NROW.LE.l) GO TO 60 SAKE9105
DO 1(1 1=1,14 SI\KE8805 NHALF=NCCL-l SAKE9110
JI=JI+III SAKF8810 NROWM=NROW-l SAKE9115
IR=fA"'l 5AKE8815 NROWH=NROW-NHALF SAKE9120
10 RIIRI=flIIR)+A(JI).BIIBD SAKE8B20 DO 6 I=I,NPCWM SAKE9125
Pf'TUPN SAKF8B25 IF(A(I,111 1.15,1 SAKE9130
E'ND 5AKE8830 ReCIP=I./A! 1.11 5AKE9135
IF! I-NROl>HI 2.2.3 SAIC.E9140
LIMIT=NHAlF 5Af(~9145
GO TO 4 SAKE9150
 LI~IT=NRCW-I SAKE9155 P't =ABS {EPS*SNGU A (KP IV II SAKE9460
4 Cr] /) J=l.LI~IT SAKr:9160 Dr: 11 Y=K,~! SAK!':Q465
JPCW=I+J SAKE9165 IJ'JRK=O.OCO SAKE9470
RATI~=-AII.J+ll*RECIP SAI<=9170 IF (tr:ND) 2,4,2 SAKE9475
1)0 5 /(=j.LI~IT SAKE9115 2 DO 3 L=l.LEND SAH9480
JCCl=K+I-J SAI<E9180 t ANF=KP IV-L SAKE94sr;
5 AIJRQW.JCnl'=AIJqnW,JCJLI+RATIceAll,K+l' SAKF91B5 U NO=l ND-l SAKE9490
[11 /) K=I.IIIRHS SAKE9190 WORK=WORK+AllANF)$A(lINDt SAKE9495
6 PIJPQW,K)=P(JROW,KD+RATIC*SII.K) SAKf9195 4 WORK=A(I~OI-WORK SAKE9500
60 IF(A(NROW.IDJ 7,14,7 SAKE9200 IF IY-/() 10.5.10 SAKE9505
1 RFCIP=l./AINROW.l) SAKE9205 5 IF (SNGlIWORKI-TOll 6,6.9 SAK~9510
f'1J 8 I=1.NflHS SAKE9210 6 IF (WORKI 19,19,1 SAKE9515
8 B(NROW.IJ=BINROW,ID*RECIP SAKE9215 7 IF (IERI 8.8,9 SAKE'9520
IF INROWelE.lI RETURN SAKE9220 8 JER=K-l SAKE'9525
D'C 13 I=l.NRCWM SIl>KE972') 9 DI~=DSQRT(WORK' SAKE9530
IPOW"'NROW-I S~IKE9230 AIKPIVI=DIN SAKE9535
RECIP=l./ACIPOW,I' SAKE9235 D! N= 1.OCO/[l IN SAKE9540
IF(I-NHAlFJ 9,10.10 SAKE9240 GO T(1 11 SAKE9545
q l PHT=I SAKE9245 10 AIIND&=WORK*OIN SAKE'9550
Gn TC 11 SAKE925C 11 I NO= INO+ I SAKE9555
Ie lIMIT=NHAlF SAKF9255 IPIII=N*IN+1)/2 SAKE9560
11 DC 12 J=I.lIMIT SAKF9260 IND=IPIV SAI<E'9565
J~CW=IRO\t+J SAKE9265 DO 16 1=1,111 SAKE9510
JC'll=J+l SAKE9270 DI~:l.DO/A(IPIV' SAKE9575
D'J 12 K=l.NRHS SAKE9215 AIIPIV'=t'yN. SAI<E9580
12 AIPCW.K)=B(IROW.KI-~(IROW.JCCll$B(JROW,KI SAKE92sn IYYN=N SAKE'9585
DO 13 J=l.NRHS SAKE9285 KEND=I-l SAKE9590
13 P(IROW.JI=BIYROW,J'$flECIP SAKE9290 lANF=III-KEND SAKF9595
RETUflN SAKE9295 IF (KENDI 15.15,12 SAKE'9600 f-J
14 I=NPflW SAKE9300 12 J=IND SlKE9605 W
15 WRIT~(6.1011 I SAKE9305 DO 14 1<=1. KENO SAKE9610 0'\
STOP SAKF,9310 W('Rl<=O.OO SAKE9615
101 FORM~T(1164H YOU GOOFEC --- T~ERE IS A lERO ON THE PRINCIPAL DIAGOSAKE9315 MIIII-=I\4IN-l SAKE9620
I~Al IN THE.I4.SH TH ROW.IID SAKE9320 lHQR=IPIV SAKE9625
END SAKF9325 lVER=J SAKE'9630
DO 13 l=lANF,MIN SAKE'9635
lVFR=LVER+l SAKE'9640
Ump=lHCfl+l SAKE9645
13 WQR~=W~RK+AllVER'$A(LHORI SAKE9650
AlJI=-WORI< 2I OII11 SAI(E'9655
14 J:J-f'lIN SAKE'9660
15 IPIV=IPIV-MIN SAKE'9665
SUBPQUTINE CSINVCA.NI SAKE9365 16 II\1D=II11D-l SAKE9670
C SAKF.9370 DO 18 1=1, III SAKE9675
C INVEIlT " GIVEN SYMMETRIC POSITIVE DEFINITE MATRIX SAKF9375 IPIV=IPIV+I 5AKE9680
C IIlICDIFIED FRC~ I~~ SCIENTIFIC SU~ROUTINE PAC~AGE SAKE9380 J=!PIV SAKE9685
C A - OOUBlE PRECISIC~ UPPER TRIANGULAR PART OF GIVEN SAKE9385 D'J 11'1 1<=I.N SAKE9690
C. SYMMETRIC POSITIVE DEFINITE III BY N COEFfICIENT SAKE9390 I>IOPI(=O.DO SAKE9695
C MATR J)(o SAKE9395 lHrp=J SAKE9700
C N - TH~ NUMBER OF ROWS (COLUMNSI IN GIVEN MATRIX. SAKE9400 DO 17 L=K,N SAKE9705
C SAKfC~405 l VEP =LH[}I:!+K- I 5AKE9710
D I ME 1\ S I ON 11. 1 5AKE9410 WOPI<=WORI<+AllHORI$AILVERI SAKE9715
nOUBlE PRECISION A,DIN,WORK SAKE9415 17 lHIlQ=LHnl)+l SAKE9720
!CPS= 1.0E-12 SAKE9420 A(JI=WCRI< SAKE9725
IF IN-I) 19,1,1 SAKE9425 18 J=J+I< SAI<E9130
IER=C SAKE9430 IF (IER, 19.20,19 SAKE9135
KPrv=o 5AKE9435 lq WRITF (6,6001 SAKE9740
Dn 11 K=l.1I! SAKE9440 STOP SAKE9745
KPIV=KP!II+K SIU(E9445 20 PETURN SAKE9150
IN!)=I(PYV
LENn=K-l
SAKE9450
SAKEQ455
600 F~RMIT (IHO.III.IHO.IOX. ERROR .1. SAKE9755
1 IH~,l()X. 'ERROl) IN (,)SINV SU"IRCUTIIIIE') SAKE9760
END <;AKF'ii165
 137

APPENDIX J. EXAMPLE PROBLEM

A three-story one-bay and one-story one-bay structure (Fig. J.l)

was analyzed by the program. Joint springs (mechanical hinge, bilinear

spring or rigid spring), flexural elements, rigid zones were used in

the example problem. The structure had elastic periods of 0.49, 0.22,

and 0.12 sec. for the first three modes. Damping matrix was assumed

to consist of two parts; one part proportional to the mass matrix, and

the other part proportional to a variable stiffness matrix.

The first 5.96 sec. of the 1940 El Centro (NS) record was

amplified by a factor of two, and used as the base motion of the

structure.

Numerical integration was carried out at every O.Onl sec.

Maximum response values were sought at every 0.01 sec. (every 10

numerical integration point,.

The total execution time (GO step) was 108 sec with 13 I/O

requests on the IBM 360/75 computer at the University of Illinois at

Urbana-Champaign.

Moment-rotation relationships of a flexural spring (including

elastic deformation) and a joint spring at an end of the third level

beam are shown in Fig. J.2.

Listings of input data cards and computer output are listed in

the following pages.

LISTING OF INPUT DATA

5. 84 - 3G .441 qz 81 2C
"EXAMPLE PROBLEM FOR SAKE PPOGRAI\4' 5.66-192.2J4174C 5.74 -2.9209540 5.80 -51.5335374 5.96 26.0328941 21
FT KIP SEC 5.88 -26.4503200 S.90 24.t454e33 5.92 17.3412811
3 2 4 '9 5 4SC;600.0 32.2 0.362 0.000945
100 .. 0 70 .. 0 50.0
12 .. 0 10 .. 0 10 .. 0
1 0.5 24.0 0.5
2 0 .. 5 ll9 .. 0 0.5
3 0 .. 5 9.0 0 .. 5
4- 0.5 11 .. 5 0.0
1 0 .. 163 15 .. 0 60.0 66 .. 0 0.00301 0 .. 0301
2 0.661 30 .. !! 120 .. 0 132 .. 0 0 .. 00141 0 .. 0147
3 1 .. 628 45.0 180 .. 0 200.0 0.000903 0.00903
4 0 .. 163 5.0 20.0 22.0 0 .. 001002 0.0100
5 0 .. 0833 20.0 50 .. 0 55.0 0 .. 00490 0.0490
6 0 .. 204 40 .. 0 100.0 HO.O 0.00401 0 .. 0401
7 0.1022 60 .. 0 1'50.0 165 .. 0 0 .. 00290 0 .. 0290
8 0 .. 0633 10 .. 0 25 .. 0 27.5 0.00245 0&0245
'9 0 .. 0833 6.0 15 .. 0 16.5 0 .. 00147 0.0147
1 20 .. 0 60 .. 0 O.N)CC990 0.000891
2 40 .. 0 120.0 0.0000544 0 .. 000490
3 60 .. 0 180 .. 0 0.0000436 0.000392
4 0.,0 0 .. 0 0 .. 0 0.0
5 -1 .. 0 0 .. 0 0 .. 0 0.0
3 1
25 .. 0
1 1 1 1 1 1
2 1 2 2 2 2
3 1 3 3 3 '3 j--I
4 3 5 5 4 4 LA)
en
6
5

7
:3
3
:3
5
6
I/)
5
6
6
..
4

4
4
,.
4

8 4 7 1 4 4
q 4 .., 7 4 4
1 1
20.0
1 2 4 4 4 4
2 4 8 8 -4 4
3 4 '9 4 5
El CENTRO INSJ 1940 '"
H~ST 6 .. 0 SEC ..
4 1 10 0 0 .. 001 0 .. 06" 1.0 0 .. 0 0 .. 0
4(F6.2.Fl20 n
.00 6 .. 1750543 .. 04 .. '9661140 .. 10 10.6543517 .16 -1.2580718 1
.. 22 13.3295023 .26 .. 9212180 .. 28 4.1170753 .. 32 .20878';n 2
.40 1681922156 .44 -23.5160113 .4!! 4.0756390 .50 5.l714943 3
.56 38 .. 2590580 .. 60 9 .. 1'507679 .66 8.7383300 .12 -31 .. 7741089 4
.14 -32 .. 9782553 .. 7f> -56. C824022 .. 80 -32 .. 1'906';176 .88 -27.1072898 5
.96 -80.0238838 1.02 -41e21t3329 1.06 -42.2446327 1.20 86.9263096 6
1.32-168 .. 5985966 1 .. 36 -84.4C6~2 1.42 -92.9193640 1.46 -89 .. 8276711 1
1.50 -95.935~7~ 1.52- 126.4401::31 1 .. 66 233.7307720 1 .. 76 142.7099915 8
1.80 173 .. 7016192 2.02-312.4440460 2.20 292.4184952 2.38-293.4189491 '9
2.42 225 .. 1721142 2 .. 52 -44.04~952 2.60 204.0352402 2.64 167.8269119 10
2.16 -30 .. 0980766 2.82 43 .. 6894221 2.88 102.7769194 2.96 -77.2391470 11
3.02 49.9471'448 3 .. 12-150.1730~16 3.20 7 .. 0102195 3.26-200.7022858 12
3.36 19101768456 3 .. 44 -~1.33'8371 3~52 169.6434650 3.58 -:-7.5690523 13
3.66 64.'U42485 3.72 -&'9.3982135 3.76 -11.1066260 3.82 30.6808448 lit
3.Qo-181.335E574 4.00 22 .. 3432508 4.04 -41.8651042 4 .. 10 20.6223583 15
4.24-223.8068905 ~.40 144.0SC;6504 48~6 -1 .. 3528931 4.56 257.3261871 16
4.66-1~6 .. 4947224 4870 63.,1.969034 4.84-266.3323822 4.92 177.6508579 17
5.00 34.fJ34U25 5.06 216.1215458 5.16 32.0Q060f,5 5.22 131.4880311 18
5.32 -lQ.4329104 5.38 169.3545113 5.46 -99.7622471 5.54 14.6209816 19
LISTING OF PROGI~ OUTPUT
CI~E~SION OF TYPICAL ELE~~NT

NONLINEAR RIC FRAME RESPCNSE PROGRAM (SAKE) NUMflFP (A-END I (MIDDL=) B-~ND I

1 5o00DE-Ol 2.4CDE Cl 5.000E-Ol

2 5.000E-Ol 1.900!" 01 5.000E-Ol
DEVELOP!:O AT UNIVERSITY Of ILLINOIS AT URBANA 3 5.000E-01 '1.000E CO 5o000E-Ol
BY OTANI, S. 4 5.000E-01 1.150e 01 0.0
eN MARCH 26, 1913
SPONSORED BY ~ATIONAL SCIENCE FOUNDATIQN
GRANT NUMBER GI-29934 MOMENT CURVATURE RELATICf\ CF flEXURAL ELEMENTS

MCMENT MCMENT CURVATURE

ELEMENT INERT! A CRACKING YIElr ING UlTIMATE YIELDIIIIG ULTIMATE
PPOfllEM "EXAMPLE PPOBLEM FOR SAKE PR:JGRAM' 1.630E-01 1.500E 01 6.000E 01 6.600E 01 3.010E-03 3.010E-02
1
2 6.670E-Ol 3.rOflE 01 1.200E 02 1.320E 02 1.410E-03 1.470E-C2
3 1.628"" 00 40500E' 01 1.8COE 02 2.000E 02 9.030E-04 9.030E-03
3 4 1.630F-01 5.000E 00 2.000E 01 2.100E 01 1.002E-03 1.000E-01
~UMBER OF STORIES 5 B.330E-02 2.00JE 01 5.CO(lE 01 5.51)OE 01 4.900E-03 4.900E-02
~UMBER CF CONSTITUENT FRAMES 2
NUMBER CF TVPICAL MEMBER DIMENSIONS 4 6 2.040E-Ol 4.000E 01 I.COOE 02 1.100E 02 4.010E-03 4.010E-02
NUMBER OF TYPICAL FLEXURAL ELEMENTS 9 7 4.220E-01 6.000E 01 1.500E 01 1.650E 02 2.900E-03 2.QOOE-01
NUMBER OF TYPICAL ROTATICNAL SPRINGS 5 8 8.33CE-02 1.000E 01 2.500E 01 2.750E 01 2.450E-03 2.450E-02
q B.330E-02 6.000E 00 l.500E 01 1.650f 01 1.41015-03 1.410E-02 I-'
W
Y!:IUNG"S MIJOULUS 4t.B96E 05 \0
E~IRTI- GR/lVny 3 .. 220E 01
FLEXURAL ROTATION AS A UNIT-LENGTH SI~PLE MEMBER
i~LASTIC D~FORMATION INCLUDED'
rnEFFICIENT ALPH ~OR D.MPI~G 3.6200-lll
COEeFIC lENT !'lETA FOR IJllMPYNG '1.4500-04

UNIT OF LENGTH FT ELEMENT CPACKING YIElDIfliG UL TIMATE

U~J1 T OF FORCE KIP
UNIT OF TIME SEC 1 6.265E-05 9.915E-04 2.275 E-03
2 3.062E-05 4.B42E-C4 l.lUE-03
3 1.88217-05 2.974E-04 1.196E-04
4 2.088E-05 3.301E-04 7.564E-C4
') 1.63515-04 i.555E-C3 3. 654E-03
6 1.335E-04 1.2721::-03 2.99015-03
7 9.680E-O,) 9.102E-C4 2.lb3E-03
CATA ALCNG STORY HEIGHT 8 8.113E-05 7.174'.:-C4 lo827E-03
9 4.904E-05 4.664E-04 1.096E-03
STmtY ISEIGHT ... EIGHT
(I( lP D ( FT t

1 i..MOE 02 1.200E 01 CHARACTERISTICS OF TYPICAL JOINT SPqINGS

2 1.000!: 01 1.0eOE 01
3 5aOOOE 01 l.OOOE 01
MC~ENT ROTATI'JN
YIELDING UL TIMATI:: YIelDING UL TIMATE
1 2.000E 01 6.000E 01 9.900E-05 8.910E-04
1 4.CCOE C1 1.200E 02 5.440E-05 4.900 E-04
3 6.000r 01 1.BOO~ 01 4.360'.:-05 3.9201:-04
4 C.O 0.0 0.0 0.0
5 -1.000E 00 0.0 0.0 0.0
 FRAME 1 FR AM!;' 2

NUMBER OF ST~RIES 3 NUMBER OF STORIES

NUM~ER OF BAYS 1 NUMBER 'JF BAYS

STRLCTUR~L DATA ALONG BAY STRUCTURAL DATA ALONG BAY

BAY ~IDTH BAY WIDTH
FT FT
2.500E 01 2.000E 01

MFMBER IDENTIFICATION "'E~BER IDENTIFICATION

FLEXURAL ELEMENT JOINT SPRING JOINT NUMBER STORY NUMBER FLEXURAL ELEMENT JOINT SPRING JOINT NUMBER STORY NUMBER I-'
+:-
DIMENSN (A-ENDBlB-ENOI (A-ENDHR-ENDI (A-ENOI(B-ENOI A-END D B-END I 0
~EMBER MEMRER DIMFNSN (A-ENO)(B-ENDI (A-ENDI(B-ENDI (A-ENOl (B-ENDI (A-ENDHB-ENDI
I 1 1 1 I. 1 1 2 3 3 1 2 4 4 .4 4 1 2 I. I.
2 I. 2 2 2 2 3 4 2 2 2 4 S 8 4 4 1 3 1 0
3 1 ] 3 :3 3 5 6 1 I. 3 4 9 9 4 5 l 4 I. 0
4 3 5 5 4 '+ 1 3 3 7.
5 3 5 5 4 4 2 4 3 2
6 3 6 6 4 I. 3 5 2 1
7 3 6, 6 4 '+ 4 6 2 1
8 4 7 1 4 '+ 5 7 I. 0 MEMBER PROPERTIES
9 4 1 1 4 4 6 8 I. 0
STIFFNESS FLEXIBLE RIG! C ZONES TOTAL LENGTH
MEMBER 0:: I I LENGTH I A-EN!)) ( B-ENDI STR. DMN. ME M. DMN
I. 7.seOE 04 1.900E 01 5.000E-01 5.000E-OI 2.000E 01 2.000E 01
MEMBER PROPERTIES 2 4.078E 04 1.150E 01 5.COO':-01 0.0 1.200E 01 1.200E 01
3 4.018E 04 1.1501" 01 5.000E-01 0.0 i.200E 01 1.200E 01
STIFFNESS FLEXIBLE RIGID IONES TOTAl LENGTH
MEMBER (EIB LENGTH (A-ENOl (B-ENDI STR. DItIN. MEM. OMN
I 7.<;801= 04 2.400E 01 5.000E-01 5.000E-Ol 2.500E 01 2.500E 01
2 3.266E 0'5 2.400E 01 SeOOOE-Ol 5.0aOE-OI 2.500E 01. 2.500E 01
3 1.<;llE 0'5 2.400E 01 5.000E-01 5o000E-Ol 2.'500E 01 l.500E 01
4 4.078E 04 9.000E 00 S.000E-01 S.OOOE-Ol l.OOOE 01 I.OOOE 01
'5 4.071310 04 9.000E 00 S.OOOE-01 5.000E-01 i.OOOE 01 I.OOOE 01
6 9.98BE 04 <;l.OOOE 00 So OC OE-01 S.OOOE-Ol I..OOOE 01 1.000E 01
7 9 988E
0 04 9.000E 00 "i.OOOE-Ol "i.OOOE-Ol l.OOOE 01 I.OOO!: 01
8 2.l'b6E 0'5 i.150E 01 5.~C'OE-Ol o.c 1.200E 01 1.200E 01
9 2.0bbE as i.150E 01 5.00DE-0! 0.0 i.200E 01 1.200F 01
 MaXIMUM A(CFL=RaTI~N AND DISPLACEMENT RESPONSES

ACCELEPAT!'JN. G DISPLACEMENT. FT
"lAS': MOrIeN El CENTRO I~SI 1940 FI~. ST 6.0 SEC. ST'JRY MAX "IN MAX MIN

FORMAT (4(F6.2.Fl2.71) 4.647E-(11 -6.080 E-f'l 8.544E-02 -2.432E-0 1

INPUT CIITA 2 6.779E-01 -5.938E-01 2.014E-Ol -3.804E-01

~UMPER OF rATA POINTS PER CARD 4 3 4.923E-Ol -5.078E-Ol 4.038E-01 -7.483E-01

START POINT OF T~f CALCULATIJN 1
NS~IP FOR RECORDING 10
FLAG F'JR QUAKE [ATA CARDS o BASE SHE~R MAX. 5.481E 01 ~IP
TI~E INTERVAL OF CC~PUT I.GOOE-03 "lASE SHEAR ~IN. -7.018 F 01 KIP
SCALE FACTrR F0R DATA 6.440'=-02
SCALE FACTOR FOR TI~E i.OOOE 00 PASE MCMfNT MAX. I.080E 03 KIP - FT
RASE LINE SHIFT OF Ee. 0.0 PASF MOMENT MIN. -1.244E 03 KIP - FT
TIME INTERVAL OF QUAKE 0.0

EATHQUAKE RECORD
f-J
+=-""
f-J
0.0 6.1750536 0.04 0.9667740 0.10 10.6543512 0.16 -1.2580709
0.22 13.3295021 0.26 0.9212180 0.28 4.11 70750 0.32 0.2087896
0.4(1 1607922(,5R 0.44 -23.5160615 0.48 4.0756388 0.50 5.1714935
0.56 38.2590485 0.60 9.1507617 0.66 8.7383299 0.72 -31.774108Q
0.74 -32o~782410 0.76 -56.0823975 0.80 -32.7906952 O.8A -27.1072845
0.96 -80.0238800 1.02 -41.2363281 1.06 -42.2446289 1.20 86.9263000
1.32-168.5985e70 1.36 -84.40689CC:; 1.42 -92.9193573 1.46 -89.8276672
1.50 -95.~35~741 1.52-126.4401245 1.66 233.7307587 1.76 142.7099915
1.80 173.7016754 2.02-312.4438417 2.20 292.4184570 2.38-293.4189453
2.42 285.1726074 2.52 -44.0480804 2.60 204.0352325 2.64 167.8269043
2.76 -)C.f.980682 2.82 43.6894(113 2.88 102.7769165 2.96 -77.2397461
3.02 49.9477386 3.12-150.1730957 3.20 7.0102186 3.26-200.7022858
3.36 191.1768341 3.44 -91.33982E5 3.52 169.6434631 3.58 -7.5690517
3.66 64.4142456 3.12 -69.39826'>1 3.76 -11.1066256 3.82 30.6808319
3.QO-181.3358459 4.00 22.343241':5 4.04 -41.8650970 4.10 20.6223450
4.24-223.8C6P848 4.40 144.0596466 4.46 -1.3528929 4.56 251.326171Q
4.66-1~6.4q47205 4.70 63.4969025 4.84-266.3322754 4.92 177.6508484
5.00 34.034103-1t 5.06 216.121~3l3 5.16 32.0005951 5.22 131.4880219
5.32 -IG.432GC71 5.38 169.3545074 5.46 -q9.1622375 5.54 14.6209812
5.66-192.2041626 5.74 -2.92C9538 5.RO -51.5335236 5.R4 -30.4419250
5.RR -26.4503174 5.90 24.6454713 5.92 11.3412781 5.96 26.032P827
n.,) 0.(1 0.0 0.0 0.0 n.(\ 1).0 1).0,
 ~I\XY~UM RESPONSE' IN FRAME 1 MI\XI~UM RESPONSE IN FRAME 2

MAXI~UM FLEXURAL RESPONSE OF MEI'II~ERS MAXIMUM FLEXURAL RESPCNSE OF MEMBERS

AS A UNIT-LENGTH SI~PLE MEMBER AS A UNIT-LENGTH SIMPLE MEMBER
ELASTIC nEFOP~ATION INCLUDED I (ELASTIC OEF~R~ATION IIIICLUDED.

MOMENT ROTATIIJN DUCTILITY MOMENT ROTATIIJN DUCTILITY

(KIP'" FT 11.01 FT D KI P *
FT I (l.01 FT I
MAX MIN
MEMSE P III AX MIN MAX MIN MAX MIN "IEMBEP 1<4 AX MIN MAX 141N

A-END 6.1681' 01 -6.327E 01 1.3511'-03 -1.6901'-03 1.3631' 00 -1.704E 00 A-END 2.0641' 01 -2.642E 01 4.6761'-04 -1.698E-IB 1.411E 00 -5 .. 1441' 00
B-EN['I 6.1681' 01 -6.327E 01 1.3511'-03 -1.6901'-03 1.363E 00 -1.704E 00 B-END 1.478E 01 -1.750= 01 2.225E-04 -2.7851'-04 6.7421'-01 -8.4381'-01
2 A-EN'" 1.2201' 07 -1.375E 02 5.. 905E-04 -1.3991'-03 1.219E 00 -2.888E 00 2 A-Et.,lf'! 2.511E 01 -1.960E 01 8.236E-04 -5.269E-04 1.059E 00 -6.778E-Ol
B-END 1.220E 02 -1.31'51' 02 5.9051'-04 -1.399E-03 102191' 00 -2.888= 00 B-E"I'J 3.15!)E Cl -2.602E 01 3.505E-03 -1.2061'-03 4.508E 0(1 -l.551E CO
3 A-END 1.8701: 02 -2.038E 1)2 ".454E-04 -8eOOOE-04 1.497E 00 -2.69(1E OC 3 A-F"ID 1.787F' 01 -1.506E 01 1.6731'-03 -4.902E-04 3.587E 00 -1.0511' 00
B-END I.870E 02 -2.038E 02 1t.454E-04 -8.000E-04 1.497E 00 -2.690': 00 B-Ft.,lO 0.0 0.0 0.0 0.0 0.0 0.0
4 A-H'D 5.916E 01 -5.751E 01 5.404E-03 -4.710E-03 3.475E 00 -3.029E 00
A-END 6.226E Cl -6.373E Cl 607081'-03 -7.322E-1)3 4.315E 00 -4.709E 00
5 II-('I\!D 5.9161: 01 -5.751E 01 5.4041'-03 -4.71OE-03 3.475E 00 -3.029E 00
1\-1'1\11) 6.226E 01 -6.373E 01 6 .. 108E-03 -1.322E-03 4.315E 00 -4.709E 00 ~AXIMU~ JOINT SPRING RESPCNSE
6 A-EIIID 9.831E 01 -1.003E 02 1.240E-03 -1.318E-03 9.748E-Ol -1.036E 00
q-~ND 1.1031' 02 -1.022E 02 3.051E-03 -1.655E-03 2.3981' 00 -i.30tE 00 MOMENT ROTATION DUCTILITY
7 A-END 9.831E 01 -1.003E 02 1.240E-f13 -1.318E-03 9.748E-01 -1.0361:: 00 (KIP'" FT FT I FT
f-J
P-EflJD 1.103E 02 -1.022[ 02 3.0511'-03 -1.655[-03 2.3981' 00 -1.301E 00 IJIE~BER MIN "'AX MIN MAX MIN
A-END 1.626E 02 -1.C89f 02 1.964[-03 -5.458E-04 2.134E 00 -5.931E-01
"'AX +
f\)
B-END 1.820E 02 -1.526E 02 3 .. 513E-03 -1.1321'-03 l.88lE 00 -1.230E 00
9 A-END 1.t:26E 02 -1.08<)E 02 1.964E-03 -5.458E-04 2.D4E 00 -5.931':-01 A-END 2. C64E 01 -2.M2E 01 0.0 D.Q 0.0 000
A-END 1.820E 02 -1.526E 02 3.5731E-03 -1. 13I2E-03 3.883E 00 -1.230E 00 B-END 1.478E 01 -1.150E 01 0.0 0.0 0.0 0.0
2 A-END 2.5111' 01 -1.960E 01 0.0 0.0 0.0 0.0
B-END 3.150E 01 -2.602E 01 0.0 0.0 0.0 0.0
3 A-ENO 1.781E 01 -1.'506E 01 0.0 0.0 0.0 0.0
MAXIMU~ JOINT SPRING RESPONSE B-END 0.0 0.0 0.0 0.0 0.0 0.0

MOMENT ROTATION DUCTILITV

(KIP'" FT FT I FT I
flIIEI'IIBEP MAX MIN MAX MIN MAX MIN

TOTAL NUMBER OF RESFCNSE STEPS 596

A-END 6.168E N -6.327E 01 9.242E-04 -9.551E-04 9.335E 00 -9.653E 00
B-fNO 6.168E 01 -6.3271' 01 9.242E-04 -9.551E-04 9.335E 00 -9.653': 00
2 A-Ef\ID 1.ll0E 02 -1 .. 315E 02 5.0ll.lU:-04 -5.853E-04 9.211E 00 -1.076E 01
S-END 1.220E 02 -1.375E 02 5.011E-0J<. -5.8531:-04 9.211E 00 -1.0161C; 01
3 4-ENI"l 1.870E 02 -1 .. 03S!: 02 4.123[-04 -4.611':-04 9.455E 00 -1.058E 01 ALL DATA PROCESSED
8-END 1.S10E" 02 -2 .. 038E 02 4.123E-04 -4.611E-04 904551' 00 -l.058l: 01
4 A-EN'" 5.9161: 01 -5.1511' 01 0.0 0.0 0.0 0.0
B-END 602261' 01 -6 .. 3130: 01 0.0 0.0 0.0 0.0
5 A-ENfl 5.<)16E 01 -5.751E 01 0.0 GoO 0.0 O.G
B-END 6.226E 01 -6.373E 01 c.o 0.0 0.0 0.0
6 A-BONf) 9.83lE 01 -1.0031' 02 0.0 0.0 0.0 0.0
A-END 1.103E 02 -1.012E 02 0.0 (I.C 0.0 0.0
7 I\-FNf) 9.83lE 01 -1.0031' 02 0.0 0.0 0.0 0.0
8-FI\D 1.103E 02 -1.022E 02 0.0 0.0 0.0 0.0
8 A-E'ND 1.626E 02 -1.C89E 01 0.0 0.0 'J.O (l.0
/'I-ENI:' 1.820E 02 -1.526E 02 0.0 0.0 0.0 0.0
9 A-till':' 1.626E 02 -1.0a9E 02 0.0 0.0 0.0 0.0
A-EN!) 1.820E 02 -1.526F. 02 0 .. 0 0.0 0.0 0.0
 W3 = 50 kip
~

10'

f-
lO'
W = 70 kip
2

I l-.
I""'"
- W = 100 kip
l

i-'
.p-
w
12'
i

,

~ 25' J 1< 20' J
Frame 1 Frame 2

Fig. J.l Example Problem

 144

Moment

Fig. J.2(a) Moment-Rotation Relation of a Flexural Spring

in the rrhird Level Beam
 vu

-wu wu
Moment

-vu

Fig. J.2(b) Moment-Rotation Relation of a Joint Spring

in the Third Level. Be&'ll