CODE:
model BasicBuilder -ndm 2 -ndf 3
set P 100
set l 121.9;
set h 91.44;
set E 20600;
set a1 1.33;
set a2 1.49;
set k 1;
set I 1 ;
node 1 0.0 0.
node 2 $l 0.0
node 3 $l $h
node 4 0 $h
node 5 0 $h
fix 1 1 1 0
fix 2 0 1 0
equalDOF 3 5 1 3
uniaxialMaterial Elastic 1 $E
puts materijal
geomTransf Linear 1
element elasticBeamColumn 1 1 2 $a1 $E $I 1
element elasticBeamColumn 2 2 3 $a2 $E $I 1
element elasticBeamColumn 3 3 4 $a2 $E $I 1
element elasticBeamColumn 4 5 1 $a1 $E $I 1
element elasticBeamColumn 5 4 2 $a1 $E $I 1
element elasticBeamColumn 6 1 3 $a1 $E $I 1
element zeroLength 7 4 5 -mat 1 -dir 2
pattern Plain 1 Linear {
load 4 $P 0.0 0.0
}
recorder Element -file 7.txt -time -ele 7 localForce
recorder Element -file 71.txt -load -ele 7 deformation
recorder Node -file D.txt -time -node 3 5 -dof 1 2 3 disp
constraints Plain
numberer RCM
test NormUnbalance 1.0e-2 250 0
integrator LoadControl 0.025 1 0.025 0.025
algorithm Newton
system ProfileSPD
analysis Static
integrator LoadControl 0.05 3 0.05 0.05
analyze 16
integrator LoadControl 0.02 3 0.02 0.02
analyze 10
print ele 1 7