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 | print("==========================")
print("Start RCFramePushover Example")
# Units: kips, in, sec
#
# Written: GLF/MHS/fmk
# Date: January 2001
from openseespy.opensees import *
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
# ----------------------------------------------------
# ----------------------------------------------------
# Start of additional modelling for lateral loads
# ----------------------------------------------------
# Define lateral loads
# --------------------
# Set some parameters
H = 10.0 # Reference lateral load
# Set lateral load pattern with a Linear TimeSeries
pattern('Plain', 2, 1)
# Create nodal loads at nodes 3 & 4
# nd FX FY MZ
load(3, H, 0.0, 0.0)
load(4, H, 0.0, 0.0)
# ----------------------------------------------------
# End of additional modelling for lateral loads
# ----------------------------------------------------
# ----------------------------------------------------
# Start of modifications to analysis for push over
# ----------------------------------------------------
# Set some parameters
dU = 0.1 # Displacement increment
# Change the integration scheme to be displacement control
# node dof init Jd min max
integrator('DisplacementControl', 3, 1, dU, 1, dU, dU)
# ----------------------------------------------------
# End of modifications to analysis for push over
# ----------------------------------------------------
# ------------------------------
# Start of recorder generation
# ------------------------------
# Stop the old recorders by destroying them
# remove recorders
# Create a recorder to monitor nodal displacements
# recorder Node -file node32.out -time -node 3 4 -dof 1 2 3 disp
# Create a recorder to monitor element forces in columns
# recorder EnvelopeElement -file ele32.out -time -ele 1 2 forces
# --------------------------------
# End of recorder generation
# ---------------------------------
# ------------------------------
# Finally perform the analysis
# ------------------------------
# Set some parameters
maxU = 15.0 # Max displacement
currentDisp = 0.0
ok = 0
test('NormDispIncr', 1.0e-12, 1000)
algorithm('ModifiedNewton', '-initial')
while ok == 0 and currentDisp < maxU:
ok = analyze(1)
# if the analysis fails try initial tangent iteration
if ok != 0:
print("modified newton failed")
break
# print "regular newton failed .. lets try an initail stiffness for this step"
# test('NormDispIncr', 1.0e-12, 1000)
# # algorithm('ModifiedNewton', '-initial')
# ok = analyze(1)
# if ok == 0:
# print "that worked .. back to regular newton"
# test('NormDispIncr', 1.0e-12, 10)
# algorithm('Newton')
currentDisp = nodeDisp(3, 1)
results = open('results.out', 'a+')
if ok == 0:
results.write('PASSED : RCFramePushover.py\n')
print("Passed!")
else:
results.write('FAILED : RCFramePushover.py\n')
print("Failed!")
results.close()
# Print the state at node 3
# print node 3
print("==========================")
|