model basic -ndm 2 -ndf 2
set E [expr 1.995*1.e6]
set mu 0.22
set Pi [expr 4*atan(1)]
node 1 0. 0.
node 2 20. 0.
node 3 20. 20.
node 4 7. 14.
nDMaterial ElasticIsotropic 1 $E $mu
element quad 1 1 2 3 4 1 PlaneStress 1
fix 1 1 1
fix 2 1 1
fix 3 1 1
recorder Node -file disp.out -node 2 -dof 1 2 disp
pattern Plain 1 Linear {
load 4 [expr -1000000*cos($Pi/6)] [expr 1000000*sin($Pi/6)] 0.
}
constraints Plain
numberer Plain
system BandGeneral
test NormDispIncr 1.e-8 6
algorithm ModifiedNewton
integrator LoadControl 0.1
analysis Static
analyze 10
loadConst -time 0.0
how about this code
model basic -ndm 2 -ndf 2
set E [expr 1.995*1.e6]
set mu 0.22
set Pi [expr 4*atan(1)]
nDMaterial ElasticIsotropic 1 $E $mu
block2D 100 100 1 1 quad "1 PlaneStress2D 1" {
1 0. 0.
2 20. 0.
3 20. 20.
4 7. 14.
}
fix 1 1 1
fix 2 1 1
fix 3 1 1
recorder Node -file disp100.out -node 4 -dof 1 2 disp
recorder Node -file reaction100.out -node 1 2 3 -dof 1 2 reaction
pattern Plain 1 Linear {
load 4 [expr -1000*cos($Pi/6)] [expr 1000*sin($Pi/6)] 0.
}
integrator LoadControl 0.1
test EnergyIncr 1.0e-12 10 0
algorithm Newton
numberer RCM
constraints Plain
system ProfileSPD
analysis Static
In OpenSees
And it is a simple gusset plate that have plane stress conditions and a tensile load is exerted on it's corner and what we want is the displacement of that node.
- Badges
- Science topic