Example 2 InputFile: 3D soil block subjected to static push over
Jump to navigation
Jump to search
Return to Tcl input files for the provided application examples:
# unit. meter, K-N, K-Pa,K-Kg
# refer to "Quan Gu, Joel P. Conte, Ahmed Elgamal and Zhaohui Yang. Response Sensitivity Analysis of a Multi-Yield-Surface J2 Plasticity Model by
# Direct Differentiation Method.June, 2009. Vol 198(30-32), Computer Methods in Applied Mechanics and Engineering, Pages 2272-2285.
# test by Lei Wang, July 2018
#create the ModelBuilder
model BasicBuilder -ndm 3 -ndf 3
reliability
# massDen G B cohension peakshear
nDMaterial MultiYieldSurfaceClay 2 3 1.5 60000.0 240000. 30.0 0.1
#updateMaterialStage -material 2 -stage 101
# 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
# 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
# -------- 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
# -------------- 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."