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 | print("==========================")
print("Start RCFrameEarthquake Example")
# Units: kips, in, sec
#
# Written: Minjie
from openseespy.opensees import *
import ReadRecord
import numpy as np
import matplotlib.pyplot as plt
wipe()
# ----------------------------------------------------
# Start of Model Generation & Initial Gravity Analysis
# ----------------------------------------------------
# Do operations of Example3.1 by sourcing in the tcl file
import RCFrameGravity
print("Gravity Analysis Completed")
# Set the gravity loads to be constant & reset the time in the domain
loadConst('-time', 0.0)
# ----------------------------------------------------
# End of Model Generation & Initial Gravity Analysis
# ----------------------------------------------------
# Define nodal mass in terms of axial load on columns
g = 386.4
m = RCFrameGravity.P/g
mass(3, m, m, 0.0)
mass(4, m, m, 0.0)
# Set some parameters
record = 'elCentro'
# Permform the conversion from SMD record to OpenSees record
dt, nPts = ReadRecord.ReadRecord(record+'.at2', record+'.dat')
# Set time series to be passed to uniform excitation
timeSeries('Path', 2, '-filePath', record+'.dat', '-dt', dt, '-factor', g)
# Create UniformExcitation load pattern
# tag dir
pattern('UniformExcitation', 2, 1, '-accel', 2)
# set the rayleigh damping factors for nodes & elements
rayleigh(0.0, 0.0, 0.0, 0.000625)
# Delete the old analysis and all it's component objects
wipeAnalysis()
# Create the system of equation, a banded general storage scheme
system('BandGeneral')
# Create the constraint handler, a plain handler as homogeneous boundary
constraints('Plain')
# 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 )
# Create the solution algorithm, a Newton-Raphson algorithm
algorithm('Newton')
# Create the DOF numberer, the reverse Cuthill-McKee algorithm
numberer('RCM')
# Create the integration scheme, the Newmark with alpha =0.5 and beta =.25
integrator('Newmark', 0.5, 0.25 )
# Create the analysis object
analysis('Transient')
# Perform an eigenvalue analysis
numEigen = 2
eigenValues = eigen(numEigen)
print("eigen values at start of transient:",eigenValues)
# set some variables
tFinal = nPts*dt
tCurrent = getTime()
ok = 0
time = [tCurrent]
u3 = [0.0]
# Perform the transient analysis
while ok == 0 and tCurrent < tFinal:
ok = analyze(1, .01)
# if the analysis fails try initial tangent iteration
if ok != 0:
print("regular newton failed .. lets try an initail stiffness for this step")
test('NormDispIncr', 1.0e-12, 100, 0)
algorithm('ModifiedNewton', '-initial')
ok =analyze( 1, .01)
if ok == 0:
print("that worked .. back to regular newton")
test('NormDispIncr', 1.0e-12, 10 )
algorithm('Newton')
tCurrent = getTime()
time.append(tCurrent)
u3.append(nodeDisp(3,1))
# Perform an eigenvalue analysis
eigenValues = eigen(numEigen)
print("eigen values at end of transient:",eigenValues)
results = open('results.out','a+')
if ok == 0:
results.write('PASSED : RCFrameEarthquake.py\n');
print("Passed!")
else:
results.write('FAILED : RCFrameEarthquake.py\n');
print("Failed!")
results.close()
plt.plot(time, u3)
plt.ylabel('Horizontal Displacement of node 3 (in)')
plt.xlabel('Time (s)')
plt.show()
print("==========================")
|
