Example 1 InputFile: 2D soil column subjected to earthquake base excitation

From OpenSeesWiki
Revision as of 23:45, 11 January 2011 by Foxchameleon (talk | contribs)
Jump to navigation Jump to search
  1. unit. meter, K-N, K-Pa,K-Kg

set h G set g -9.8 model basic -ndm 2 -ndf 2


  1. -----------------------------------------------------
  2. Nodes Definition
  3. -----------------------------------------------------

node 1 272.9 -43.426 -mass 0 0 node 2 272.9 -38.728 -mass 0 0 node 3 272.9 -34.031 -mass 0 0 node 4 272.9 -29.334 -mass 0 0 node 5 272.9 -24.636 -mass 0 0 node 6 272.9 -19.939 -mass 0 0 node 7 272.9 -15.242 -mass 0 0 node 8 272.9 -11.813 -mass 0 0 node 9 272.9 -8.3835 -mass 0 0 node 10 272.9 -4.9543 -mass 0 0 node 11 272.9 -1.5252 -mass 0 0 node 12 272.9 1.904 -mass 0 0 node 13 272.9 5.3331 -mass 0 0 node 14 272.9 8.7623 -mass 0 0 node 15 272.9 12.191 -mass 0 0 node 16 272.9 14.096 -mass 0 0 node 17 272.9 16.001 -mass 0 0 node 18 272.9 17.905 -mass 0 0 node 19 272.9 19.81 -mass 0 0 node 20 272.9 20.877 -mass 0 0 node 21 272.9 21.945 -mass 0 0 node 22 272.9 22.402 -mass 0 0 node 23 272.9 22.859 -mass 0 0 node 24 272.9 24.459 -mass 0 0 node 25 272.9 26.059 -mass 0 0 node 26 272.9 27.659 -mass 0 0 node 27 272.9 29.259 -mass 0 0 node 28 272.9 29.602 -mass 0 0 node 29 272.9 29.945 -mass 0 0 node 30 267.27 -43.426 -mass 0 0 node 31 267.27 -38.728 -mass 0 0 node 32 267.27 -34.031 -mass 0 0 node 33 267.27 -29.334 -mass 0 0 node 34 267.27 -24.636 -mass 0 0 node 35 267.27 -19.939 -mass 0 0 node 36 267.27 -15.242 -mass 0 0 node 37 267.27 -11.813 -mass 0 0 node 38 267.27 -8.3835 -mass 0 0 node 39 267.27 -4.9543 -mass 0 0 node 40 267.27 -1.5252 -mass 0 0 node 41 267.27 1.904 -mass 0 0 node 42 267.27 5.3331 -mass 0 0 node 43 267.27 8.7623 -mass 0 0 node 44 267.27 12.191 -mass 0 0 node 45 267.27 14.096 -mass 0 0 node 46 267.27 16.001 -mass 0 0 node 47 267.27 17.905 -mass 0 0 node 48 267.27 19.81 -mass 0 0 node 49 267.27 20.877 -mass 0 0 node 50 267.27 21.945 -mass 0 0 node 51 267.27 22.402 -mass 0 0 node 52 267.27 22.859 -mass 0 0 node 53 267.27 24.432 -mass 0 0 node 54 267.27 26.012 -mass 0 0 node 55 267.27 27.593 -mass 0 0 node 56 267.27 29.166 -mass 0 0 node 57 267.27 29.497 -mass 0 0 node 58 267.27 29.828 -mass 0 0

  1. massDen G B cohension PeakShearStrain

nDMaterial MultiYieldSurfaceClay 1 2 0.0 54450 1.6e5 33. .1 nDMaterial MultiYieldSurfaceClay 2 2 0.0 33800 1.0e5 26. .1 nDMaterial MultiYieldSurfaceClay 3 2 0.0 96800 2.9e5 44. .1 nDMaterial MultiYieldSurfaceClay 4 2 0.0 61250 1.8e5 35. .1 nDMaterial MultiYieldSurfaceClay 5 2 0.0 180000 5.4e5 60. .1 nDMaterial MultiYieldSurfaceClay 6 2 0.0 369800 1.1e6 86. .1

  1. -----------------------------------------------------
  2. Soil Elememets Definition
  3. -----------------------------------------------------
  4. thick material Mat Press mDensity gravity

element quadWithSensitivity 1 57 28 29 58 6.0957 "PlaneStrain" 1 0 2.0 0 $g element quadWithSensitivity 2 56 27 28 57 6.0957 "PlaneStrain" 1 0 2.0 0 $g element quadWithSensitivity 3 55 26 27 56 6.0957 "PlaneStrain" 2 0 2.0 0 $g element quadWithSensitivity 4 54 25 26 55 6.0957 "PlaneStrain" 2 0 2.0 0 $g element quadWithSensitivity 5 53 24 25 54 6.0957 "PlaneStrain" 2 0 2.0 0 $g element quadWithSensitivity 6 52 23 24 53 6.0957 "PlaneStrain" 2 0 2.0 0 $g element quadWithSensitivity 7 51 22 23 52 6.0957 "PlaneStrain" 3 0 2.0 0 $g element quadWithSensitivity 8 50 21 22 51 6.0957 "PlaneStrain" 3 0 2.0 0 $g element quadWithSensitivity 14 15 16 45 44 6.0957 "PlaneStrain" 4 0 2.0 0 $g element quadWithSensitivity 13 16 17 46 45 6.0957 "PlaneStrain" 4 0 2.0 0 $g element quadWithSensitivity 12 17 18 47 46 6.0957 "PlaneStrain" 4 0 2.0 0 $g element quadWithSensitivity 11 18 19 48 47 6.0957 "PlaneStrain" 4 0 2.0 0 $g element quadWithSensitivity 10 19 20 49 48 6.0957 "PlaneStrain" 4 0 2.0 0 $g element quadWithSensitivity 9 20 21 50 49 6.0957 "PlaneStrain" 4 0 2.0 0 $g element quadWithSensitivity 15 44 43 14 15 6.0957 "PlaneStrain" 5 0 2.0 0 $g element quadWithSensitivity 16 43 42 13 14 6.0957 "PlaneStrain" 5 0 2.0 0 $g element quadWithSensitivity 17 42 41 12 13 6.0957 "PlaneStrain" 5 0 2.0 0 $g element quadWithSensitivity 18 41 40 11 12 6.0957 "PlaneStrain" 5 0 2.0 0 $g element quadWithSensitivity 19 40 39 10 11 6.0957 "PlaneStrain" 5 0 2.0 0 $g element quadWithSensitivity 20 39 38 9 10 6.0957 "PlaneStrain" 5 0 2.0 0 $g element quadWithSensitivity 21 38 37 8 9 6.0957 "PlaneStrain" 5 0 2.0 0 $g element quadWithSensitivity 22 37 36 7 8 6.0957 "PlaneStrain" 5 0 2.0 0 $g element quadWithSensitivity 23 36 35 6 7 6.0957 "PlaneStrain" 6 0 2.0 0 $g element quadWithSensitivity 24 35 34 5 6 6.0957 "PlaneStrain" 6 0 2.0 0 $g element quadWithSensitivity 25 34 33 4 5 6.0957 "PlaneStrain" 6 0 2.0 0 $g element quadWithSensitivity 26 33 32 3 4 6.0957 "PlaneStrain" 6 0 2.0 0 $g element quadWithSensitivity 27 32 31 2 3 6.0957 "PlaneStrain" 6 0 2.0 0 $g element quadWithSensitivity 28 31 30 1 2 6.0957 "PlaneStrain" 6 0 2.0 0 $g

fix 1 1 1 fix 30 1 1

equalDOF 29 58 1 2 equalDOF 28 57 1 2 equalDOF 27 56 1 2 equalDOF 26 55 1 2 equalDOF 25 54 1 2 equalDOF 24 53 1 2 equalDOF 23 52 1 2 equalDOF 22 51 1 2 equalDOF 21 50 1 2 equalDOF 20 49 1 2 equalDOF 19 48 1 2 equalDOF 18 47 1 2 equalDOF 17 46 1 2 equalDOF 16 45 1 2 equalDOF 15 44 1 2 equalDOF 14 43 1 2 equalDOF 13 42 1 2 equalDOF 12 41 1 2 equalDOF 11 40 1 2 equalDOF 10 39 1 2 equalDOF 9 38 1 2 equalDOF 8 37 1 2 equalDOF 7 36 1 2 equalDOF 6 35 1 2 equalDOF 5 34 1 2 equalDOF 4 33 1 2 equalDOF 3 32 1 2 equalDOF 2 31 1 2

  1. ------------------------- add sensitivity ------------------------

reliability

set h G

  1. --- G1 ---

parameter 1 -element 1 -material $h addToParameter 1 -element 2 -material $h

  1. --- G2 ---

parameter 2 -element 3 -material $h addToParameter 2 -element 4 -material $h addToParameter 2 -element 5 -material $h addToParameter 2 -element 6 -material $h

  1. --- G3 ---

parameter 3 -element 7 -material $h addToParameter 3 -element 8 -material $h

  1. --- G4 ---

parameter 4 -element 9 -material $h addToParameter 4 -element 10 -material $h addToParameter 4 -element 11 -material $h addToParameter 4 -element 12 -material $h addToParameter 4 -element 13 -material $h addToParameter 4 -element 14 -material $h

  1. --- G5 -----

parameter 5 -element 15 -material $h addToParameter 5 -element 16 -material $h addToParameter 5 -element 17 -material $h addToParameter 5 -element 18 -material $h addToParameter 5 -element 19 -material $h addToParameter 5 -element 20 -material $h addToParameter 5 -element 21 -material $h addToParameter 5 -element 22 -material $h

  1. ---- G6 -----

parameter 6 -element 23 -material $h addToParameter 6 -element 24 -material $h addToParameter 6 -element 25 -material $h addToParameter 6 -element 26 -material $h addToParameter 6 -element 27 -material $h addToParameter 6 -element 28 -material $h

  1. -----

set h cohesion

  1. ---- h7 --------

parameter 7 -element 1 -material $h addToParameter 7 -element 2 -material $h

  1. ---- h8 --------

parameter 8 -element 3 -material $h addToParameter 8 -element 4 -material $h addToParameter 8 -element 5 -material $h addToParameter 8 -element 6 -material $h

  1. ---- h9 --------

parameter 9 -element 7 -material $h addToParameter 9 -element 8 -material $h

  1. ---- h10 --------

parameter 10 -element 9 -material $h addToParameter 10 -element 10 -material $h addToParameter 10 -element 11 -material $h addToParameter 10 -element 12 -material $h addToParameter 10 -element 13 -material $h addToParameter 10 -element 14 -material $h

  1. ---- h11 --------

parameter 11 -element 15 -material $h addToParameter 11 -element 16 -material $h addToParameter 11 -element 17 -material $h addToParameter 11 -element 18 -material $h addToParameter 11 -element 19 -material $h addToParameter 11 -element 20 -material $h addToParameter 11 -element 21 -material $h addToParameter 11 -element 22 -material $h

  1. ---- h12 --------

parameter 12 -element 23 -material $h addToParameter 12 -element 24 -material $h addToParameter 12 -element 25 -material $h addToParameter 12 -element 26 -material $h addToParameter 12 -element 27 -material $h addToParameter 12 -element 28 -material $h


  1. -----

set h K

  1. ------ k13 --------

parameter 13 -element 1 -material $h addToParameter 13 -element 2 -material $h

  1. ------ k14 --------

parameter 14 -element 3 -material $h addToParameter 14 -element 4 -material $h addToParameter 14 -element 5 -material $h addToParameter 14 -element 6 -material $h

  1. ------ k15 --------

parameter 15 -element 7 -material $h addToParameter 15 -element 8 -material $h

  1. ------ k16 --------

parameter 16 -element 9 -material $h addToParameter 16 -element 10 -material $h addToParameter 16 -element 11 -material $h addToParameter 16 -element 12 -material $h addToParameter 16 -element 13 -material $h addToParameter 16 -element 14 -material $h

  1. ------ k17 --------

parameter 17 -element 15 -material $h addToParameter 17 -element 16 -material $h addToParameter 17 -element 17 -material $h addToParameter 17 -element 18 -material $h addToParameter 17 -element 19 -material $h addToParameter 17 -element 20 -material $h addToParameter 17 -element 21 -material $h addToParameter 17 -element 22 -material $h

  1. ------ k18 --------

parameter 18 -element 23 -material $h addToParameter 18 -element 24 -material $h addToParameter 18 -element 25 -material $h addToParameter 18 -element 26 -material $h addToParameter 18 -element 27 -material $h addToParameter 18 -element 28 -material $h


  1. ------------------------ recorder ----------------------------------

recorder Node -file disp29.out -time -node 29 -precision 16 -dof 1 disp recorder Node -file disp27.out -time -node 27 -dof 1 disp recorder Node -file disp23.out -time -node 23 -dof 1 disp recorder Node -file disp21.out -time -node 21 -dof 1 disp recorder Node -file disp15.out -time -node 15 -dof 1 disp recorder Node -file disp7.out -time -node 7 -dof 1 disp

recorder Node -file ddm29G1.out -time -node 29 -dof 1 -precision 16 "sensitivity 1" recorder Node -file ddm27G1.out -time -node 27 -dof 1 "sensitivity 1" recorder Node -file ddm23G1.out -time -node 23 -dof 1 "sensitivity 1" recorder Node -file ddm21G1.out -time -node 21 -dof 1 "sensitivity 1" recorder Node -file ddm15G1.out -time -node 15 -dof 1 "sensitivity 1" recorder Node -file ddm7G1.out -time -node 7 -dof 1 "sensitivity 1"

recorder Node -file ddm29G2.out -time -node 29 -dof 1 "sensitivity 2" recorder Node -file ddm27G2.out -time -node 27 -dof 1 "sensitivity 2" recorder Node -file ddm23G2.out -time -node 23 -dof 1 "sensitivity 2" recorder Node -file ddm21G2.out -time -node 21 -dof 1 "sensitivity 2" recorder Node -file ddm15G2.out -time -node 15 -dof 1 "sensitivity 2" recorder Node -file ddm7G2.out -time -node 7 -dof 1 "sensitivity 2"

recorder Node -file ddm29G3.out -time -node 29 -dof 1 "sensitivity 3" recorder Node -file ddm27G3.out -time -node 27 -dof 1 "sensitivity 3" recorder Node -file ddm23G3.out -time -node 23 -dof 1 "sensitivity 3" recorder Node -file ddm21G3.out -time -node 21 -dof 1 "sensitivity 3" recorder Node -file ddm15G3.out -time -node 15 -dof 1 "sensitivity 3" recorder Node -file ddm7G3.out -time -node 7 -dof 1 "sensitivity 3"

recorder Node -file ddm29G4.out -time -node 29 -dof 1 "sensitivity 4" recorder Node -file ddm27G4.out -time -node 27 -dof 1 "sensitivity 4" recorder Node -file ddm23G4.out -time -node 23 -dof 1 "sensitivity 4" recorder Node -file ddm21G4.out -time -node 21 -dof 1 "sensitivity 4" recorder Node -file ddm15G4.out -time -node 15 -dof 1 "sensitivity 4" recorder Node -file ddm7G4.out -time -node 7 -dof 1 "sensitivity 4"

recorder Node -file ddm29G5.out -time -node 29 -dof 1 "sensitivity 5" recorder Node -file ddm27G5.out -time -node 27 -dof 1 "sensitivity 5" recorder Node -file ddm23G5.out -time -node 23 -dof 1 "sensitivity 5" recorder Node -file ddm21G5.out -time -node 21 -dof 1 "sensitivity 5" recorder Node -file ddm15G5.out -time -node 15 -dof 1 "sensitivity 5" recorder Node -file ddm7G5.out -time -node 7 -dof 1 "sensitivity 5"

recorder Node -file ddm29G6.out -time -node 29 -dof 1 "sensitivity 6" recorder Node -file ddm27G6.out -time -node 27 -dof 1 "sensitivity 6" recorder Node -file ddm23G6.out -time -node 23 -dof 1 "sensitivity 6" recorder Node -file ddm21G6.out -time -node 21 -dof 1 "sensitivity 6" recorder Node -file ddm15G6.out -time -node 15 -dof 1 "sensitivity 6" recorder Node -file ddm7G6.out -time -node 7 -dof 1 "sensitivity 6"

recorder Node -file ddm29t1.out -time -node 29 -dof 1 "sensitivity 7" recorder Node -file ddm27t1.out -time -node 27 -dof 1 "sensitivity 7" recorder Node -file ddm23t1.out -time -node 23 -dof 1 "sensitivity 7" recorder Node -file ddm21t1.out -time -node 21 -dof 1 "sensitivity 7" recorder Node -file ddm15t1.out -time -node 15 -dof 1 "sensitivity 7" recorder Node -file ddm7t1.out -time -node 7 -dof 1 "sensitivity 7"

recorder Node -file ddm29t2.out -time -node 29 -dof 1 "sensitivity 8" recorder Node -file ddm27t2.out -time -node 27 -dof 1 "sensitivity 8" recorder Node -file ddm23t2.out -time -node 23 -dof 1 "sensitivity 8" recorder Node -file ddm21t2.out -time -node 21 -dof 1 "sensitivity 8" recorder Node -file ddm15t2.out -time -node 15 -dof 1 "sensitivity 8" recorder Node -file ddm7t2.out -time -node 7 -dof 1 "sensitivity 8"

recorder Node -file ddm29t3.out -time -node 29 -dof 1 "sensitivity 9" recorder Node -file ddm27t3.out -time -node 27 -dof 1 "sensitivity 9" recorder Node -file ddm23t3.out -time -node 23 -dof 1 "sensitivity 9" recorder Node -file ddm21t3.out -time -node 21 -dof 1 "sensitivity 9" recorder Node -file ddm15t3.out -time -node 15 -dof 1 "sensitivity 9" recorder Node -file ddm7t3.out -time -node 7 -dof 1 "sensitivity 9"

recorder Node -file ddm29t4.out -time -node 29 -dof 1 "sensitivity 10" recorder Node -file ddm27t4.out -time -node 27 -dof 1 "sensitivity 10" recorder Node -file ddm23t4.out -time -node 23 -dof 1 "sensitivity 10" recorder Node -file ddm21t4.out -time -node 21 -dof 1 "sensitivity 10" recorder Node -file ddm15t4.out -time -node 15 -dof 1 "sensitivity 10" recorder Node -file ddm7t4.out -time -node 7 -dof 1 "sensitivity 10"

recorder Node -file ddm29t5.out -time -node 29 -dof 1 "sensitivity 11" recorder Node -file ddm27t5.out -time -node 27 -dof 1 "sensitivity 11" recorder Node -file ddm23t5.out -time -node 23 -dof 1 "sensitivity 11" recorder Node -file ddm21t5.out -time -node 21 -dof 1 "sensitivity 11" recorder Node -file ddm15t5.out -time -node 15 -dof 1 "sensitivity 11" recorder Node -file ddm7t5.out -time -node 7 -dof 1 "sensitivity 11"

recorder Node -file ddm29t6.out -time -node 29 -dof 1 "sensitivity 12" recorder Node -file ddm27t6.out -time -node 27 -dof 1 "sensitivity 12" recorder Node -file ddm23t6.out -time -node 23 -dof 1 "sensitivity 12" recorder Node -file ddm21t6.out -time -node 21 -dof 1 "sensitivity 12" recorder Node -file ddm15t6.out -time -node 15 -dof 1 "sensitivity 12" recorder Node -file ddm7t6.out -time -node 7 -dof 1 "sensitivity 12"

  1. ----

recorder Node -file ddm29k1.out -time -node 29 -dof 1 "sensitivity 13" recorder Node -file ddm29k2.out -time -node 29 -dof 1 "sensitivity 14" recorder Node -file ddm29k3.out -time -node 29 -dof 1 "sensitivity 15" recorder Node -file ddm29k4.out -time -node 29 -dof 1 "sensitivity 16" recorder Node -file ddm29k5.out -time -node 29 -dof 1 "sensitivity 17" recorder Node -file ddm29k6.out -time -node 29 -dof 1 "sensitivity 18"

  1. -----

recorder Element -element 1 -time -file stress1.out material 1 stress -dT 0.01 recorder Element -element 1 -time -file strain1.out material 1 strain -dT 0.01 recorder Element -element 15 -time -file stress15.out material 1 stress -dT 0.01 recorder Element -element 15 -time -file strain15.out material 1 strain -dT 0.01 recorder Element -element 22 -time -file stress22.out material 1 stress -dT 0.01 recorder Element -element 22 -time -file strain22.out material 1 strain -dT 0.01 recorder Element -element 28 -time -file stress28.out material 1 stress -dT 0.01 recorder Element -element 28 -time -file strain28.out material 1 strain -dT 0.01


  1. --------------------------------------------------------
  2. Soil Gravity Analysis Model Definition
  3. --------------------------------------------------------

system BandGeneral test NormDispIncr 1.E-12 50 2 constraints Transformation integrator LoadControl 1 1 1 1 algorithm Newton numberer RCM

  1. ---------------- add sensitivity ---------------------

sensitivityIntegrator -static sensitivityAlgorithm -computeAtEachStep

analysis Static analyze 1

puts "soil gravity nonlinear analysis completed ..."

wipeAnalysis

constraints Transformation test NormDispIncr 1.E-12 50 2 algorithm Newton numberer RCM system BandGeneral

integrator NewmarkWithSensitivity 0.55 0.2756 sensitivityIntegrator -definedAbove sensitivityAlgorithm -computeAtEachStep analysis Transient

pattern UniformExcitation 1 1 -accel "Series -factor 1 -filePath acce.txt -dt 0.01"

set startT [clock seconds] analyze 1000 0.01 set endT [clock seconds] puts "Execution time: [expr $endT-$startT] seconds."

Retrieved from "https://opensees.ist.berkeley.edu/wiki/index.php?title=Example_1_InputFile:_2D_soil_column_subjected_to_earthquake_base_excitation&oldid=7430"