[AlexBud01]

Учусь моделировать нелинейный бетон в Ансис. Решил сделать расчет, аналогичный представленному в данном видео:

Нажмите для просмотра

https://www.youtube.com/watch?v=Vb-xVXdY3Go&index=27&list=PLRyAsVf7brpS2Y5sNGde-7utlM8TZ3BcB . Разобрался в командах apdl; однако при расчете выдает одни и те же ошибки. Перепроверил все много раз, не могу понять в чем проблема. Помогите пожалуйста найти решение (прикладываю проект). Версия 16.2

Команды apdl:

ДЛЯ ГЕОМЕТРИИ:

Temperature

TOFFST,273.15
TREF,22.0

Material model (concrete):

*SET,MAT_ID,101

DENS_=2400.0
EX_=3.6E+10
NUXY_=0.18

MP,DENS,MAT_ID,DENS_,
MP,EX,MAT_ID,EX_,
MP,NUXY,MAT_ID,NUXY_,

Rb_c=42E+6
Rb_t=2.5E+6


C1=0.3
C2=0.7
C3=Rb_t
C4=-1
C9=0.6

TB,concr,MAT_ID,1,,1

TBTEMP,22.0

tbdata,1,C1
tbdata,2,C2
tbdata,3,C3
tbdata,4,C4
tbdata,9,C9

Element type (Solid65)

*SET,ELEM_TYPE_ID,101,
ET,ELEM_TYPE_ID,SOLID65
KEYOPT,ELEM_TYPE_ID,1,0
KEYOPT,ELEM_TYPE_ID,3,2
KEYOPT,ELEM_TYPE_ID,5,0
KEYOPT,ELEM_TYPE_ID,6,0
KEYOPT,ELEM_TYPE_ID,7,1
KEYOPT,ELEM_TYPE_ID,8,0

Material model (Steel)

*SET,MAT_ID,102

DENS_=7850.0
EX_=2.1E+11
NUXY_=0.3

MP,DENS,MAT_ID,DENS_,
MP,EX,MAT_ID,EX_,
MP,NUXY,MAT_ID,NUXY_,

Yield_strenght=2.5E+08
Tangent_modulus=1.4E+09

TB,BKIN,MAT_ID,1,,0

TBTEMP,22.0

TBDATA,1,Yield_strenght
TBDATA,2,Tangent_modulus

Element type (Link180)

*SET,ELEM_TYPE_ID,102

ET,ELEM_TYPE_ID,LINK180

KEYOPT,ELEM_TYPE_ID,2,0
KEYOPT,ELEM_TYPE_ID,3,0
KEYOPT,ELEM_TYPE_ID,12,0

Sections (Link180)

*SET,SEC_ID,301

SECTYPE,SEC_ID,LINK,CSOLID,,

R_=0.005
S_=3.141592*R_*R_
SECDATA,S_

*SET,SEC_ID,302

SECTYPE,SEC_ID,LINK,CSOLID,,

R_=0.007
S_=3.141592*R_*R_
SECDATA,S_

ДЛЯ STATIC STRUCTURAL:

Commands_PREP7

fini
/prep7


*SET,ELEM_TYPE_ID,101

CMSEL,S,NS_concr, ELEM

EMODIF,ALL,TYPE,ELEM_TYPE_ID

*SET,MAT_ID,101

EMODIF,ALL,MAT,MAT_ID
*set,matid,MAT_ID

ALLSEL,ALL



*SET,ELEM,TYPE_ID,102

CMSEL,S,NS_reinf_all,ELEM

EMODIF,ALL,TYPE,ELEM_TYPE_ID

ALLSEL,ALL

*SET,MAT_ID,102

CMSEL,S,NS_reinf_D10,ELEM
CMSEL,S,NS_reinf_D14,ELEM

EMODIF,ALL,MAT,MAT_ID
*set,matid,MAT_ID

ALLSEL,ALL


*SET,SEC_ID,301

CMSEL,S,NS_reinf_D10,ELEM

EMODIF,ALL,SECNUM,SEC_ID

ALLSEL,ALL

*SET,SEC_ID,302
CMSEL,S,NS_reinf_D14,ELEM

EMODIF,ALL,SECNUM,SEC_ID

ALLSEL,ALL



CMSEL,S,NS_concr,ELEM
CMSEL,A,NS_reinf_all,ELEM
NSLE,S,ALL
cpintf,all,0.0004

ALLSEL,ALL

fini

Nonlinear solution control

NLGEOM,OFF

SOLCONTROL,ON,OFF,NOPL

nropt,unsym

PRED,OFF,,OFF

LNSRCH,OFF

neqit,100

KBS,0

cutcontrol,plslimit,0.5


cutcontrol,cutbackfactor,0.25
cutcontrol,noshape,1

SHPP,WARN

pivchek,off

Equation Solver control

EQSLV,SPRASE,1.0E-8,3.0,,
BCSOPTION,,INCORE,16000,,,OFF

Step control [1]

AUTOTS,ON

NSUBST,3,100,3,OFF

time,25.0

outres,all,all

STABILIZE,OFF,,,,

CNVTOL,F,,0.5,2,0.01

Step control [2]

AUTOTS,ON

NSUBST,20,1000,20,OFF

time,35.0

outres,all,all

STABILIZE,OFF,,,,

CNVTOL,F,,0.5,2,1.0E-3

ДЛЯ SOLUTION:

PLCRACK

/GRAPHICS,FULL

/SHOW,PNG

/gfile,600

/device,victor,1

/PLOPTS,INFO,3

/GROPT,CURL,ON

!/GST,ON

!/CPLANE,1

!WPSTYLE,,,,,,,,0

!/auto,1

/FOCUS,1,0.0,0.2,1.5

/VIEW,1,1.0,0.0,0.0

/ANGLE,1,1.0,0.0,0.0

/DIST,1,0.8,0

*SET,step_number,2

*SET,k_,1

*SET,substep_number,0

esel,ALL

*IF,substep_number,EQ,0,THEN

SET,step_number

*get,n_substep,active,,solu,NCMSS

*do,i,1,n)substep,k_

SET,step_number,i,,,,,

PLCRACK,0,0

*enddo

*ENDIF

*IF,substep_number,GT,0,THEN

SET,step_number,substep_number,,,,,

PLCRACK,0,0

[MrWhite]

На какой стадии ошибка то? Просто тут вот такие вещи есть например: *do,i,1,n)substep,k_. Но вряд ли до них доходит.
И вообще зачем Вам это (решил сделать расчет...)?