Skip to content

Commit

Permalink
Update TDmodel.py
Browse files Browse the repository at this point in the history
  • Loading branch information
BijanSeif committed Dec 2, 2022
1 parent dd84a9d commit 6d7e6b6
Showing 1 changed file with 6 additions and 31 deletions.
37 changes: 6 additions & 31 deletions Examples/Python/Models/TDmodel.py
Original file line number Diff line number Diff line change
Expand Up @@ -7,7 +7,6 @@
## Updated - 09/10/2020
##################################################################

import opsvis as opsplt
import openseespy.opensees as ops
import numpy as np
import topspy.dynamic as bjd
Expand Down Expand Up @@ -111,24 +110,13 @@
nodeTag1 += 1
nodeTag1 += 1

# calculate eigenvalues & print results
# numEigen = 7
# eigenValues = ops.eigen(numEigen)
# PI = 2 * asin(1.0)

for r in [186, 192, 198, 204, 210, 181, 187, 193, 199, 205, 182, 183, 184, 185, 206, 207, 208, 209]:
ops.equalDOF(196,r,1,2)

for r in [156,162,168,174,151,157,163,169,175,176,177,178,179,180,152,153,154,155]:
ops.equalDOF(165,r,1,2)

###################################
#### Display the active model with node tags only
# opsplt.plot_model("nodes")

#### Display specific mode shape with scale factor of 300 using the active model
# opsplt.plot_modeshape(5, 300)

###################################
# To save the analysis output for deformed shape, use createODB command before running the analysis
# The following command saves the model data, and output for gravity analysis and the first 3 modes
# in a folder "3DFrame_ODB"

# opsplt.createODB("3DFrame", "Gravity", Nmodes=3)

# Define Static Analysis
# create SOE
Expand All @@ -149,19 +137,6 @@

print(ops.nodeReaction(6))

# IMPORTANT: Make sure to issue a wipe() command to close all the recorders. Not issuing a wipe() command
# ... can cause errors in the plot_deformedshape() command.



####################################
### Now plot mode shape 2 with scale factor of 300 and the deformed shape using the recorded output data

# opsplt.plot_modeshape(2, 300, Model="3DFrame")

# opsplt.plot_defo()



G = 3860
[dt,TH]=bjd.LAT2('models/1')
Expand Down

0 comments on commit 6d7e6b6

Please sign in to comment.