# Example 2 InputFile: 3D soil block subjected to static push over

Return to Tcl input files for the provided application examples:


```tcl

  1. unit. meter, K-N, K-Pa,K-Kg
  2. refer to "Quan Gu, Joel P. Conte, Ahmed Elgamal and Zhaohui Yang. Response Sensitivity Analysis of a Multi-Yield-Surface J2 Plasticity Model by
  3. Direct Differentiation Method.June, 2009. Vol 198(30-32), Computer Methods in Applied Mechanics and Engineering, Pages 2272-2285.
  4. test by Lei Wang, July 2018
  1. create the ModelBuilder

model BasicBuilder -ndm 3 -ndf 3 reliability

  1. massDen G B cohension peakshear

nDMaterial MultiYieldSurfaceClay 2 3 1.5 60000.0 240000. 30.0 0.1

  1. updateMaterialStage -material 2 -stage 101
  1. create meshes

node 1 0.00000 0.0000 0.00000 node 2 0.00000 0.0000 0.5000 node 3 0.00000 0.0000 1.000000 node 4 0.00000 0.5000 0.00000 node 5 0.00000 0.5000 0.5000 node 6 0.00000 0.5000 1.000000 node 7 0.00000 1.00000 0.00000 node 8 0.00000 1.00000 0.5000 node 9 0.00000 1.00000 1.000000 node 10 0.5000 0.0000 0.00000 node 11 0.5000 0.0000 0.5000 node 12 0.5000 0.0000 1.000000 node 13 0.5000 0.5000 0.00000 node 14 0.5000 0.5000 0.5000 node 15 0.5000 0.5000 1.000000 node 16 0.5000 1.00000 0.00000 node 17 0.5000 1.00000 0.5000 node 18 0.5000 1.00000 1.000000 node 19 1.000000 0.0000 0.00000 node 20 1.000000 0.0000 0.5000 node 21 1.000000 0.0000 1.000000 node 22 1.000000 0.5000 0.00000 node 23 1.000000 0.5000 0.5000 node 24 1.000000 0.5000 1.000000 node 25 1.000000 1.00000 0.00000 node 26 1.000000 1.00000 0.50000 node 27 1.000000 1.00000 1.000000

  1. ele Node mat

element bbarBrickWithSensitivity 1 1 10 13 4 2 11 14 5 2 element bbarBrickWithSensitivity 2 2 11 14 5 3 12 15 6 2 element bbarBrickWithSensitivity 3 10 19 22 13 11 20 23 14 2 element bbarBrickWithSensitivity 4 11 20 23 14 12 21 24 15 2 element bbarBrickWithSensitivity 5 4 13 16 7 5 14 17 8 2 element bbarBrickWithSensitivity 6 13 22 25 16 14 23 26 17 2 element bbarBrickWithSensitivity 7 14 23 26 17 15 24 27 18 2 element bbarBrickWithSensitivity 8 5 14 17 8 6 15 18 9 2

fix 1 1 1 1 0 0 0 fix 4 1 1 1 0 0 0 fix 7 1 1 1 0 0 0 fix 10 1 1 1 0 0 0 fix 13 1 1 1 0 0 0 fix 16 1 1 1 0 0 0 fix 19 1 1 1 0 0 0 fix 22 1 1 1 0 0 0 fix 25 1 1 1 0 0 0

  1. -------- define random variable ---------------------------------

set h K

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

set h G

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

set h cohesion

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

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

recorder Node -file ddm3k.out -time -node 3 -dof 1 2 3 -precision 16 "sensitivity 1" recorder Node -file ddm3G.out -time -node 3 -dof 1 2 3 -precision 16 "sensitivity 2" recorder Node -file ddm3cohesion.out -time -node 3 -dof 1 2 3 -precision 16 "sensitivity 3"

recorder Node -file disp3.out -time -node 3 -dof 1 2 3 -precision 16 disp

recorder Element -ele 1 -time -file stress.out material 2 stress recorder Element -ele 1 -time -file strain.out material 2 strain

pattern Plain 1 "Sine 0.0 1000.0 10.0 -factor 0.8 " { load 3 2.5 0.0 0 load 6 2.5 0.0 0 load 9 2.5 0.0 0 }

pattern Plain 2 "Sine 0.0 1000.0 10.0 -shift 1.5708 -factor 0.8" { load 3 0.0 2.5 0 load 12 0.0 2.5 0 load 21 0.0 2.5 0 }

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

sensitivityAlgorithm -computeAtEachStep

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

```