Счет рассыпается при различных материальных константах

[AlexLPN 0]

Добрый день. Есть задача о скольжении бруска по бруску и выделении тепла от трения. Вот код, который взят из хэлпа ансиса.

/PREP7

ET,1,PLANE13,4 ! 2-D COUPLED-FIELD SOLID

RECT,,5,,1.25

RECT,,1.25,1.25,2.5

LESIZE,ALL,0.25, , ,1,1

AMESH,ALL

UIMP,1,EX, , ,70000,

UIMP,1,DENS, , ,2.7E-9,

UIMP,1,ALPX, , ,23.86E-6,

UIMP,1,NUXY, , ,0.3,

UIMP,1,MU, , ,0.2,

UIMP,1,KXX, , ,150,

UIMP,1,C, , ,9E8,

TOFFSET,460

TUNIF = 0.00

ET,2,TARGE169 ! 2-D TARGET SEGMENT

ET,3,CONTA171 ! 2-D SURFACE-TO-SURFACE CONTACT

KEYOPT,3,1,1

ASEL,S,,,1

NSLA,S,1

NSEL,R,LOC,Y,1.25

R,1

TYPE,2

ESURF

ALLSEL

ASEL,S,,,2

NSLA,S,1

NSEL,R,LOC,Y,1.25

TYPE,3

ESURF

ALLSEL,ALL

*DIM,PRE,TABLE,2,1,1,TIME

SFL,7,PRES, %PRE%

PRE(1,0,1) = 0

PRE(1,1,1) = 10

PRE(2,0,1) = 10

PRE(2,1,1) = 10

ALLSEL

SAVE

FINISH

/SOLU

NSUB,1

ASEL,S,,,2

NSLA,S,1

NSEL,R,LOC,X,1.25

D,ALL,UX,3.75

ALLSEL

ASEL,S,,,1

NSLA,S,1

D,ALL,UX,0

D,ALL,UY,0

ALLSEL

ANTYPE,TRANS

TIMINT,OFF,STRUC

TINTPR,,,,1.0

NLGEOM,ON

TIME,3.75E-3

AUTO,ON

NSUB,100,10000,100

OUTRES,ALL,-10

NROP,UNSYM

/OUT,SCRATCH

SOLVE

TIME,1

AUTO,ON

TINTPR,,,,1.0

NSUB,100,10000,10

OUTRES,ALL,LAST

SOLVE

Мне необходимо задать другой материал. Тут он задан не в тех размерностях. В общем вот какой мне надо

UIMP,1,EX, , ,2e11,

UIMP,1,DENS, , ,7850,

UIMP,1,ALPX, , ,1.2e-05,

UIMP,1,NUXY, , ,0.3,

UIMP,1,MU, , ,0.2,

UIMP,1,KXX, , ,60.5,

UIMP,1,C, , ,444,

Это для стали. И после такой подстановки уже решение сыпется. Как быть? Как правильно поменять материал?

[Борман 2 376]

Сначала перепроверьте граничные условия. Например это

D,ALL,UX,3.75

И вообще модель... В той ли размерности она набрана.

[AlexLPN 0]

А что перепроверять то?

D,ALL,UX,3.75 - это смещение верхнего бруска.

Размерность так, которая была. Единиц в задаче не было, но уверен, что все правильно т.к. это задача из хэлпа ансиса со ссылкой на автора и публикацию в журнале. Да и к тому же первый код отлично работает без проблем.

Ну если предположить, то модуль Юнга 70000 это ближе к миллиметрам. То есть получается ньютон на мм^2 или ньютон на м^2 * 10^6 и того 7 на 10^11... круто конечно... если сталь 2 * 10^11 ... ну в общем это не принципиально, Ансис то не проверяет размерность, для него это просто числа... так что моя проблема видимо в численном счете... а как тут быть... не разобрался пока

[Борман 2 376]

UIMP,1,DENS, , ,2.7E-9,

Судя по этому - модель набрана в мм, ну т.е. как положено... Ваши новые свойства - в метрах, но АНСИС считает их ""по старому, у вас плотность стала 7850 тонн/мм^3

Обычно такое просто так не проходит...

Совет - приведите в порядок модель.

[AlexLPN 0]

И действительно!!! Поставил материал свой только в мм... все заработало... учиться и учиться))) Спасибо