Example 2 InputFile: 3D soil block subjected to static push over: Difference between revisions
Jump to navigation
Jump to search
Foxchameleon (talk | contribs) No edit summary |
No edit summary |
||
(One intermediate revision by the same user not shown) | |||
Line 3: | Line 3: | ||
<source lang="tcl"> | <source lang="tcl"> | ||
# | # 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 | #create the ModelBuilder | ||
Line 11: | Line 16: | ||
# massDen G B cohension peakshear | # massDen G B cohension peakshear | ||
nDMaterial MultiYieldSurfaceClay 2 3 1.5 60000.0 240000. 30.0 0.1 | nDMaterial MultiYieldSurfaceClay 2 3 1.5 60000.0 240000. 30.0 0.1 | ||
#updateMaterialStage -material 2 -stage 101 | |||
# create meshes | # create meshes | ||
Line 62: | Line 69: | ||
# -------- define random variable --------------------------------- | # -------- define random variable --------------------------------- | ||
set h K | |||
parameter 1 -element 1 -material $h | parameter 1 -element 1 -material $h | ||
addToParameter 1 -element 2 -material $h | addToParameter 1 -element 2 -material $h | ||
Line 72: | Line 79: | ||
addToParameter 1 -element 7 -material $h | addToParameter 1 -element 7 -material $h | ||
addToParameter 1 -element 8 -material $h | addToParameter 1 -element 8 -material $h | ||
set h G | |||
parameter 2 -element 1 -material $h | parameter 2 -element 1 -material $h | ||
addToParameter 2 -element 2 -material $h | addToParameter 2 -element 2 -material $h | ||
Line 82: | Line 90: | ||
addToParameter 2 -element 7 -material $h | addToParameter 2 -element 7 -material $h | ||
addToParameter 2 -element 8 -material $h | addToParameter 2 -element 8 -material $h | ||
set h cohesion | |||
parameter 3 -element 1 -material $h | parameter 3 -element 1 -material $h | ||
addToParameter 3 -element 2 -material $h | addToParameter 3 -element 2 -material $h | ||
Line 92: | Line 103: | ||
addToParameter 3 -element 7 -material $h | addToParameter 3 -element 7 -material $h | ||
addToParameter 3 -element 8 -material $h | addToParameter 3 -element 8 -material $h | ||
# -------------- recorder ------------------------------------------------ | # -------------- recorder ------------------------------------------------ | ||
recorder Node -file | recorder Node -file ddm3k.out -time -node 3 -dof 1 2 3 -precision 16 "sensitivity 1" | ||
recorder Node -file | recorder Node -file ddm3G.out -time -node 3 -dof 1 2 3 -precision 16 "sensitivity 2" | ||
recorder Node -file | 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 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 stress.out material 2 stress | ||
recorder Element -ele 1 -time -file strain.out material 2 strain | recorder Element -ele 1 -time -file strain.out material 2 strain | ||
Line 112: | Line 126: | ||
load 21 0.0 2.5 0 | load 21 0.0 2.5 0 | ||
} | } | ||
system BandGeneral | system BandGeneral | ||
Line 119: | Line 135: | ||
algorithm Newton | algorithm Newton | ||
numberer RCM | numberer RCM | ||
analysis Static | |||
sensitivityAlgorithm -computeAtEachStep | sensitivityAlgorithm -computeAtEachStep | ||
set startT [clock seconds] | set startT [clock seconds] | ||
Line 129: | Line 144: | ||
set endT [clock seconds] | set endT [clock seconds] | ||
puts "Execution time: [expr $endT-$startT] seconds." | puts "Execution time: [expr $endT-$startT] seconds." | ||
</source> | </source> |
Latest revision as of 17:49, 17 October 2018
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."