Learning Forum

[Mauricio Casarotto]

I have the following script to model and study an isotropic hardening model:

---------------------------------------------------------------------------------------------------

FINISH

/CLEAR

 

!Asignando nombres al analisis

 

/FILNAME,viga-ep-script,1

 

/Title,Analisis-ElastoPlastico

 

!Definiendo unidades (m,kg,s,K)

/UNITS,SI

 

!Entrando al preprocesador

/PREP7

 

!Definiendo tipo de análisis

ANTYPE,STATIC,NEW

 

 

!Definiendo propiedades del material

MAT,1

MP,EX,1,2E11

MP,PRXY,1,0.3

TB,BISO,1

TBDATA,1,4E8,1E9

 

!Definiendo keypoints

K,1,0,0,0

K,2,0.003,0,0

K,3,0.003,0.001,0

K,4,0,0.001,0

 

 

!Definiendo lineas

LSTR,1,2

LSTR,2,3

LSTR,3,4

LSTR,4,1

 

!Definiendo area

A,1,2,3,4

 

!Definiendo nodos

NKPT,1,1

NKPT,13,2

NKPT,65,3

NKPT,53,4

 

FILL,1,13

FILL,13,65,3,26,13

FILL,1,53,3,14,13

FILL,14,26

FILL,27,39

FILL,40,52

FILL,53,65

 

!Definiendo elemento

ET,1,PLANE182,0,,2,,,1

 

!Asignando malla al area

AATT,1,,,,1

 

!Estableciendo la malla

MSHAPE,0,2D

MSHKEY,0

LESIZE,1,,,12

LESIZE,2,,,4

LESIZE,3,,,12

LESIZE,4,,,4

AMESH,ALL

NUMMRG,NODE

 

!Estableciendo las cargas y restricciones

D,1,ALL,0

D,14,UX,0,,53,13

D,65,UX,0

D,65,UY,-0.001

 

 

!Iniciando el solucionador

/SOLU

 

!Incluyendo no linealidades

NLGEOM,ON

AUTOTS,OFF

DELTIM,1

NSUBST,80

 

KBC,0

 

!Solucionando

OUTRES,ALL,ALL,,,,

SOLVE

FINISH

 

/POST1

 

 

SET,LAST,80

 

PLDISP,1

 

PLNSOL,S,EQV,0

 

PLNSOL,EPPL,EQV,0,,,,AUTO

 

PLNSOL,EPTO,EQV,0,,,,AUTO

 

/POST26

 

!Curva U-L

NUMVAR,200

TIMERANGE,0,1

NSOL,2,65,U,Y,DESPLAZAMIENTO

XVAR,2

ABS,2,2,,,DESPLA,,,1

PLTIME,0,1

ESOL,3,48,65,F,Y,CARGA

/AXLAB,X,DISPLACEMENT

/AXLAB,Y,LOAD

PLVAR,3

---------------------------------------------------------------------------------------------------

 

While the deformed shape is right and compares to the ones in previous studies, when I graph the load vs displacement I get the right numbers but a thousand times bigger than the previous studies. The values of forces should be  in the range of 60 to 65.

 

I dont know what am I doing wrong.

 

[Chandra Sekaran]

My guess is that it must be something inconsistent with the units. Can you add units as comments to the input deck? What is the reference you are comparing to?

[Mauricio Casarotto]

I tried commenting the /UNITS command, but the output was the same. Here you can find the reference A comparison of triangle and quadrangle finite elements for large strain elastoplastic formulations