1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184 | print("==========================")
from openseespy.opensees import *
print("Starting RCFrameGravity example")
# Create ModelBuilder (with two-dimensions and 3 DOF/node)
model('basic', '-ndm', 2, '-ndf', 3)
# Create nodes
# ------------
# Set parameters for overall model geometry
width = 360.0
height = 144.0
# Create nodes
# tag, X, Y
node(1, 0.0, 0.0)
node(2, width, 0.0)
node(3, 0.0, height)
node(4, width, height)
# Fix supports at base of columns
# tag, DX, DY, RZ
fix(1, 1, 1, 1)
fix(2, 1, 1, 1)
# Define materials for nonlinear columns
# ------------------------------------------
# CONCRETE tag f'c ec0 f'cu ecu
# Core concrete (confined)
uniaxialMaterial('Concrete01', 1, -6.0, -0.004, -5.0, -0.014)
# Cover concrete (unconfined)
uniaxialMaterial('Concrete01', 2, -5.0, -0.002, 0.0, -0.006)
# STEEL
# Reinforcing steel
fy = 60.0; # Yield stress
E = 30000.0; # Young's modulus
# tag fy E0 b
uniaxialMaterial('Steel01', 3, fy, E, 0.01)
# Define cross-section for nonlinear columns
# ------------------------------------------
# some parameters
colWidth = 15
colDepth = 24
cover = 1.5
As = 0.60 # area of no. 7 bars
# some variables derived from the parameters
y1 = colDepth / 2.0
z1 = colWidth / 2.0
section('Fiber', 1)
# Create the concrete core fibers
patch('rect', 1, 10, 1, cover - y1, cover - z1, y1 - cover, z1 - cover)
# Create the concrete cover fibers (top, bottom, left, right)
patch('rect', 2, 10, 1, -y1, z1 - cover, y1, z1)
patch('rect', 2, 10, 1, -y1, -z1, y1, cover - z1)
patch('rect', 2, 2, 1, -y1, cover - z1, cover - y1, z1 - cover)
patch('rect', 2, 2, 1, y1 - cover, cover - z1, y1, z1 - cover)
# Create the reinforcing fibers (left, middle, right)
layer('straight', 3, 3, As, y1 - cover, z1 - cover, y1 - cover, cover - z1)
layer('straight', 3, 2, As, 0.0, z1 - cover, 0.0, cover - z1)
layer('straight', 3, 3, As, cover - y1, z1 - cover, cover - y1, cover - z1)
# Define column elements
# ----------------------
# Geometry of column elements
# tag
geomTransf('PDelta', 1)
# Number of integration points along length of element
np = 5
# Lobatto integratoin
beamIntegration('Lobatto', 1, 1, np)
# Create the coulumns using Beam-column elements
# e tag ndI ndJ transfTag integrationTag
eleType = 'forceBeamColumn'
element(eleType, 1, 1, 3, 1, 1)
element(eleType, 2, 2, 4, 1, 1)
# Define beam elment
# -----------------------------
# Geometry of column elements
# tag
geomTransf('Linear', 2)
# Create the beam element
# tag, ndI, ndJ, A, E, Iz, transfTag
element('elasticBeamColumn', 3, 3, 4, 360.0, 4030.0, 8640.0, 2)
# Define gravity loads
# --------------------
# a parameter for the axial load
P = 180.0; # 10% of axial capacity of columns
# Create a Plain load pattern with a Linear TimeSeries
timeSeries('Linear', 1)
pattern('Plain', 1, 1)
# Create nodal loads at nodes 3 & 4
# nd FX, FY, MZ
load(3, 0.0, -P, 0.0)
load(4, 0.0, -P, 0.0)
# ------------------------------
# End of model generation
# ------------------------------
# ------------------------------
# Start of analysis generation
# ------------------------------
# Create the system of equation, a sparse solver with partial pivoting
system('BandGeneral')
# Create the constraint handler, the transformation method
constraints('Transformation')
# Create the DOF numberer, the reverse Cuthill-McKee algorithm
numberer('RCM')
# Create the convergence test, the norm of the residual with a tolerance of
# 1e-12 and a max number of iterations of 10
test('NormDispIncr', 1.0e-12, 10, 3)
# Create the solution algorithm, a Newton-Raphson algorithm
algorithm('Newton')
# Create the integration scheme, the LoadControl scheme using steps of 0.1
integrator('LoadControl', 0.1)
# Create the analysis object
analysis('Static')
# ------------------------------
# End of analysis generation
# ------------------------------
# ------------------------------
# Finally perform the analysis
# ------------------------------
# perform the gravity load analysis, requires 10 steps to reach the load level
analyze(10)
# Print out the state of nodes 3 and 4
# print node 3 4
# Print out the state of element 1
# print ele 1
u3 = nodeDisp(3, 2)
u4 = nodeDisp(4, 2)
results = open('results.out', 'a+')
if abs(u3 + 0.0183736) < 1e-6 and abs(u4 + 0.0183736) < 1e-6:
results.write('PASSED : RCFrameGravity.py\n')
print("Passed!")
else:
results.write('FAILED : RCFrameGravity.py\n')
print("Failed!")
results.close()
print("==========================")
|