自己编译最新版的OpenSees (地震工程模拟开放系统)

 source truss.tcl

# ------------------------------# Start of model generation# ------------------------------# Remove existing modelwipe# Create ModelBuilder (with two-dimensions and 2 DOF/node)model BasicBuilder -ndm 2 -ndf 2# Create nodes# ------------# Create nodes & add to Domain - command: node nodeId xCrd yCrdnode 1   0.0  0.0node 2 144.0  0.0node 3 168.0  0.0node 4  72.0 96.0# Set the boundary conditions - command: fix nodeID xResrnt? yRestrnt?fix 1 1 1 fix 2 1 1fix 3 1 1# Define materials for truss elements# -----------------------------------# Create Elastic material prototype - command: uniaxialMaterial Elastic matID EuniaxialMaterial Elastic 1 3000# Define elements# ---------------# Create truss elements - command: element truss trussID node1 node2 A matIDelement Truss 1 1 4 10.0 1element Truss 2 2 4  5.0 1element Truss 3 3 4  5.0 1# Define loads# ------------#create a Linear TimeSeries (load factor varies linearly with time): command timeSeries Linear $tagtimeSeries Linear 1# Create a Plain load pattern with a linear TimeSeries: command pattern Plain $tag $timeSeriesTag { $loads }pattern Plain 1 1 {        # Create the nodal load - command: load nodeID xForce yForce    load 4 100 -50}# ------------------------------# Start of analysis generation# ------------------------------# Create the system of equation, a SPD using a band storage schemesystem BandSPD# Create the DOF numberer, the reverse Cuthill-McKee algorithmnumberer RCM# Create the constraint handler, a Plain handler is used as homo constraintsconstraints Plain# Create the integration scheme, the LoadControl scheme using steps of 1.0integrator LoadControl 1.0# Create the solution algorithm, a Linear algorithm is createdalgorithm Linear# create the analysis object analysis Static # ------------------------------# Start of recorder generation# ------------------------------# create a Recorder object for the nodal displacements at node 4recorder Node -file example.out -time -node 4 -dof 1 2 disp# Create a recorder for element forces, one in global and the other local systemrecorder Element -file eleGlobal.out -time -ele 1 2 3 forcesrecorder Element -file eleLocal.out -time -ele 1 2 3  basicForces# ------------------------------# Finally perform the analysis# ------------------------------# Perform the analysisanalyze 1# ------------------------------# Print Stuff to Screen# ------------------------------# Print the current state at node 4 and at all elementsputs "node 4 displacement: [nodeDisp 4]"print node 4print ele