Модель материала Drucker-Prager

/PREP7

!Тип КЭ
ET,1,182,,,1

!Модель материала
MP,EX,1,12E6
MP,NUXY,1,0.3
MP,DENS,1,1700
MP,ALPX,1,11E-6

*AFUN,DEG

C=1500
FI=25
DEL=0
ALPHA1=(6*SIN(FI))/(3-SIN(FI))
ALPHA2=(6*SIN(DEL))/(3-SIN(DEL))
SIGMA=(6*C*COS(FI))/(3-SIN(FI))

TB,EDP,1,1,2,LYFUN          
TBDATA,1,ALPHA1,SIGMA

TB,EDP,1,1,2,LFPOT          
TBDATA,1,ALPHA2

UIMP,1,REFT,,,0

!Геометрия
K,1,0,0,0
K,2,2,0,0
K,3,1.5,1,0
K,4,0,1,0

A,1,2,3,4

!КЭ сетка
ESIZE,0.05
MSHKEY,2
MSHAPE,0,2D

TYPE,1
MAT,1
AMESH,ALL

ALLS

FINISH

 

/SOL

!Стадия №1

!Закрепление
NSEL,S,LOC,Y,-0.01,0.01
D,ALL,UY,0
NSEL,S,LOC,X,-0.01,0.01
D,ALL,UX,0
ALLS

!Нагрузка
ACEL,0,10,0
BFE,ALL,TEMP,1,1000

ALLS
SBCTRAN

!Настройка

ANTYPE,0
NLGEOM,1

OUTRES,ALL,ALL
NSUBST,100,1000,50
NEQIT,10
AUTOTS,ON
NROPT,UNSYM,,ON    

SOLVE

!Стадия №2

!Нагрузка

ALLS
BFE,ALL,TEMP,1,0,,,

ALLS
SOLVE

!Стадия №3

!Нагрузка

ALLS
BFE,ALL,TEMP,1,1000,,,

ALLS
SOLVE

 

/POST1  
PLNSOL,U,SUM,0,1.0
PLNSOL,U,Y,0,1.0
PLNSOL,U,X,0,1.0