Elastic-Perfectly Plastic Material

[kyungtae79]

Hi, all

I am just checking the elastoplastic material as a hook for a expansion joint by using the zero-length element.

A model is very simple following :

wipe
model BasicBuilder -ndm 2 -ndf 3

# Column
set Ac 3.3428e+1
set Ic 2.195e+1
set Ec 20.0e+10

# Shear keys
uniaxialMaterial ElasticPP 13 1.08e+6 3.4259e-3

node 1 0. 0.
node 2 0. 10.
node 3 0. 10.

fix 1 1 1 1

set transfColumn 2
geomTransf Linear $transfColumn

element elasticBeamColumn 1 1 2 $Ac $Ec $Ic $transfColumn
element zeroLength 2 2 3 -mat 13 -dir 1
equalDOF 2 3 2 3

pattern Plain 100 Linear {
load 3 100 0. 0.
}

recorder Node -file disp.txt -node 1 2 3 -dof 1 2 3 disp
recorder Element -file force.txt -ele 2 force

constraints Penalty 1.0e+18 1.0e+18
numberer Plain
system BandGeneral
integrator DisplacementControl 3 1 0.0001
#integrator LoadControl 0000.1

test EnergyIncr 1.e-6 10 0
algorithm Newton
analysis Static

set ok [analyze 1000]

------------------------------------------
The output just comes out only with a elastic range.
To look at a plastic part, how can I do?

I will appreciate any comments on that.

[vesna]

I changed a little bit your file and it works. Look at differences.

wipe
model BasicBuilder -ndm 2 -ndf 3

# Column
set Ac 3.3428e+1
set Ic 2.195e+1
set Ec 20.0e+10

# Shear keys
uniaxialMaterial ElasticPP 13 1.08e+6 3.4259e-3

node 1 0. 0.
node 2 0. 10.
node 3 0. 10.

fix 1 1 1 1

set transfColumn 2
geomTransf Linear $transfColumn

element elasticBeamColumn 1 1 2 $Ac $Ec $Ic $transfColumn
element zeroLength 2 2 3 -mat 13 -dir 1
equalDOF 2 3 2 3

pattern Plain 100 Linear {
#load 3 100 0. 0.
sp 3 1 1.0
}

recorder Node -file disp.txt -node 1 2 3 -dof 1 2 3 disp
recorder Element -file force2.txt -ele 2 force
recorder Element -file force1.txt -ele 1 force

#constraints Penalty 1.0e+18 1.0e+18
constraints Transformation
numberer Plain
system BandGeneral
#integrator DisplacementControl 3 1 0.0001
integrator LoadControl 0.00005

test EnergyIncr 1.e-6 100 0
algorithm Newton
analysis Static

set ok [analyze 1000]