Мультилинейная модель материала

[AlexeyT 0]

Помогите пожалуйста, где-то я видно запутался.

Имею материала модуль юнга E=2e11 Па, коэффициент пуассона mu=0.3, условный предел текучести Sig02=5e8 Па, коэффициент упрочнения m=0.1

Пытаюсь расчитать точки для задания кривой деформирования (в ANSYS)

1. определяю предел текучести

SigT =(Sig02/(0,002*E+Sig02))^(1/(1-m))

2. определяю деформацию для предела текучести EpsT = SigT/E

3. пересчитываю на истиные деформации и напряжения

EpsTi = ln(1 + EpsT)

SigTi = SigT(1 + EpsT)

4. задаюсь произвольным шагом по Sig и нахожу Eps по формуле

Eps = EpsTi(Sig/SigTi)^(1/m)

Где-то получаю ошибку. Не могу найти где, просто не сходится реальная диаграмма с моей аппроксимацией.

Посмотрите пожалуйста где я заблудился. Или если у кого-то есть кода на APDL решающий данную задачу, поделитесь пожалуста.

Спасибо.

[Артем Кулаченко 0]

Алексей, у Вас что действительно такие огромные деформации что надо переводить. Кстати напряжение так пересчитывать можно не всегда.

Я как то вот так задавал зависимость, правда другую

! Multilinear Isotropic Hardening

NNN=20

TB,MISO,1,,NNN

*DO,i,1,NNN

a1=0.05/NNN*i

a2=A*a1+B*tanh(C*a1)*1E6

*IF,i,EQ,1,THEN

MP,EX,1,a2/a1 ! Material constants

*ENDIF

TBPT,,a1,a2

*ENDDO

[Islander 2]

Для одновременного учета физической и геометрической нелинейности.

*SET,s02,345000 !Условный предел текучести

*SET,fs,490000 !Предел прочности

*SET,psi,63 !Относительное сужение при разрыве

*SET,E,2.1e8 !Модуль Юнга

*SET,m,0.75*log(fs*(1+1.4*(psi/100))/s02)/log(log(100/(100-psi))/((s02/E)+0.002))

*SET,ys,(s02/(E*0.002+s02)**m)**(1/(1-m))

*SET,sm,ys*(log(100/(100-psi))/(ys/E))**m !Напряжение, соответствующее предельной деформации

/SOLU !Вход в решатель

ANTYPE,0

NLGEOM,1 !Учет больших перемещений

NSUBST,20,0,0 !Число шагов нагружения

OUTRES,ALL,ALL !Сохранение промежуточных результатов для всех этапов нагружения

TIME,1 !Время конца нагружения

FINISH

/PREP7 !Вход в препроцессор

TB,MISO,1,1,20 !Закон упрочнения материала (мультилинейное изотропное упрочнение)

TBTEMP,0 !Температура

TBPT,,(ys/E),ys !Начальная точка графика σ =f(e)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(1/19))**(1/m),ys+(sm-ys)*(1/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(2/19))**(1/m),ys+(sm-ys)*(2/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(3/19))**(1/m),ys+(sm-ys)*(3/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(4/19))**(1/m),ys+(sm-ys)*(4/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(5/19))**(1/m),ys+(sm-ys)*(5/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(6/19))**(1/m),ys+(sm-ys)*(6/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(7/19))**(1/m),ys+(sm-ys)*(7/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(8/19))**(1/m),ys+(sm-ys)*(8/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(9/19))**(1/m),ys+(sm-ys)*(9/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(10/19))**(1/m),ys+(sm-ys)*(10/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(11/19))**(1/m),ys+(sm-ys)*(11/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(12/19))**(1/m),ys+(sm-ys)*(12/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(13/19))**(1/m),ys+(sm-ys)*(13/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(14/19))**(1/m),ys+(sm-ys)*(14/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(15/19))**(1/m),ys+(sm-ys)*(15/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(16/19))**(1/m),ys+(sm-ys)*(16/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(17/19))**(1/m),ys+(sm-ys)*(17/19)

TBPT,,(ys/E)*(1+((sm/ys)-1)*(18/19))**(1/m),ys+(sm-ys)*(18/19)

TBPT,, log(100/(100-psi)),sm !Конечная точка графика σ =f(e)

***

[Gravyx 0]

asds,

Это конечно не принципиально, но можно было бы последний кусок циклом сделать:

TBPT,,(ys/E),ys !Начальная точка графика σ =f(e)

*do,i,1,18,1

TBPT,,(ys/E)*(1+((sm/ys)-1)*(i/19))**(1/m),ys+(sm-ys)*(i/19)

*enddo

TBPT,, log(100/(100-psi)),sm !Конечная точка графика σ =f(e)

[AlexeyT 0]

Огромное спасибо ВСЕМ !

Вы мне очень помогли

2Артем Кулаченко:

дело в том, что в дальнейшем планируется расчитывать слоистый композит, и заранее неизвестно как поведут себя некоторые слои.