Seismic analysis platform for structures
##########################################################################
Author: Junjun Guo(HomePage)
E-mail: jjguo2@bjtu.edu.cn/guojj_ce@163.com
##########################################################################
- Tutorial-1:Install SAPStructure and view OpenSeesPy model
- Tutorial-2:Auxiliary modules
- Tutorial-3:Quickly generate model nodes and elements with AutoCAD dxf file
- Tutorial-4:2D soil profile meshing with AutoCAD dxf file
- Tutorial-5:Near fault ground motion simulation
- Download the zip file
- Run the example model (eg. exmple 1)
- Download SAPStructure from (https://www.yunpan.com/surl_yGNgf2BXaD9 （code：b00f）) or (https://drive.google.com/file/d/1qRovfA5kCcRkKkifJXYBlciovm06GX9c/view?usp=drive_link)
- When encounter error, just close the window, and reopen it.
- Select SAPStructure and right click the mouse, then click the properties and choose running the program as an administrator.
- Click loadResultDB button, and load the result database
- Then display the model and conduct post process.
Prepare your own openseespy model by referring the examples
- a class for calculate ground motion intensity measure. please use the command print(help(CalculateGroundMotionIMs)) to check the structure and the usage of the class
- class for ground motion baseline correction ,fltering, and conversion among acceleration,velocity and displacement, please use the command print(help(GroundMotionProcess)) to check the structure and the usage of the class
- a class for the visualization and quick construction of OpenSeesPy model. please use the command print(help(OpenSeesPyX)) to check the structure and the usage of the class
- a class for calculating the section properties. please use the command print(help(SectionPropertyCalculate)) to check the structure and the usage of the class
- a class for section moment curvature analysis. please use the command print(help(SectMCAnalysis)) to check the structure and the usage of the class
- a class for horizontally rotate FE model,it is convinient to get the rotated node coordinates use this class. please use the command print(help(ExciteAnyDirectionOpenSees)) to check the structure and the usage of the class
- a class for python interacting with the SAP2000 program.
- a class for processing shake table test data, such as calculating the periods and damping ratios of a structure
- Clicking the "ModelPath" button to select an folder to save the model database and the generated model nodes and elements.(There generates a template dxf file named "cadModelTemplate.dxf" with several specific layers,such as girder,pier,etc., plot each model part in the corresponding layer.), see Figure 2.1.
- Setting the number of segments for arc and spline, the elements length of girder, pier and cap beam.
- Clicking the "DXFModelLoad" button to load your model dxf file. (You can reference an example girder model dxf file constructed with the generated template file, called "girderBridgeExample.dxf", see Figure 2.2)
- Clicking the "generateModel" button to automatically generate the model nodes and elements, and save them into the model database and txt files in the selected folder.
- Go to the main windown of SAPStructure, clicking the button "loadModelDB" to load the generated model database, and display the model, see Figures 2.3 and 2.4. In addition, the model can be visualized with node and element tags, see Figures 2.5 and 2.6.
Figure 3.1
Figure 3.2
Figure 3.3
Figure 3.4
Figure 3.5
Figure 3.6
- Open the 2D soil profile divide window, see Figure 3.1.
- Click "resultsSavePath" button to select an folder to save the generated nodes and elements files.
- Click "dxfModelLoad" button to load the soil profile dxf file (see the example dxf file-soilProfile.dxf).
- Enter the layer name of the rectangular border of the soil profile (Make sure the soil profile is within the rectangular border).
- Specify the width and height of the divided elements.
- Click the "meshSoilProfile" button to divide the rectangular region (see Figure 3.2).
- Then, entering the layer name of the closed boundary of a local soil region, and enter the number of the soil region.
- Click "soilIdentify" button to identify the generated elements that within the specified soil region.
- Repeat steps 7 and 8 until meshing all the soil regions (see Figure 3.3).
- Finally, click "saveNodesEles" button to save the nodes and the elements of each soil region to the selected folder.
Figure 4.1
Figure 4.2
Figure 4.3
- The near_fault ground motions are simulated based on the method proposed by Dabaghi and Der Kiureghian((2014 PEER report) "Stochastic Modeling and Simulation of Near-Fault Ground Motions for Performance-Based Earthquake Engineering)
- The input parameters:
faultType: type of faulting,including ["strikeSlip","reverseAndReverseOblique"]
simType: type of simulatated grond motions, including ["pulseAndNoPulse","onlyPulse","onlyNoPulse"]
momentMag: moment magnitude of an Earthquake
zTOR(km): depth to the top of the rupture plane in kilometers
Rrup(km): the closest distance from the site to the fault rupture in kilometers
Vs30(m/s): site soil shear wave average velocity over the top 30 meters in meters per second
s_or_d(km): directivity parameters s or d (input the large of s and d)
theta_or_phi(degrees): directivity angle theta or phi (input corresponding to s_or_d)
numSimMotions(int): the number of simulated ground motion time histories
Figure 5.1
Figure 5.2
Figure 5.3