[vsae]

Добрый день, сразу оговорюсь что в ANSYS я недавно. Работаю исключительно с твердотелкой в APDL mechanical.
Была поставлена задача рассчета электрического напряжения пьезокерамики сложной модели. Параллельно шёл её расчет на бумаге. Результаты сильно рознились, поэтому я взялся за более простые тестовые задачи.
Суть задачи в том, чтобы к массе и скрепленному с ней пьезоэлементу приложить ускорение, получив эл. напряжение. Построил простую модель, задал постоянное ускорение системы в целом. Тут возникла проблема следующего характера, при заданном ускорении на массу с определенным сечением должна действовать одна сила, в постпроцессинге значится совсем другая, на порядок меньше.
Далее вместо ускорения была задана непосредственно сила, масса была удалена, остался лишь пьезоэлемент. Результаты постпроцессинга совпали с бумажными.
Вопрос: Как нужно задавать ускорение чтобы при массе m груза на пьезоэлемент создавалась сила F=ma, равная той силе что была задана мануально без ускорения и массы.

Пример прилагается

FINISH ! Make sure we are at BEGIN level
/CLEAR,NOSTART ! Clear model since no SAVE found

/FILNAME,MULTILAYER1,0
/TITLE,MULTILAYER PIEZOELECTRIC DEVICE
/PREP7
*SET,L,40e-3
*SET,D1,30e-3
D2=0.027
S=L*0.1
*SET,m,1e-3
CYLIND,D1/2,0,0,L-S,0,360
CYLIND,D1/2,D2/2+m,L-S,L-m,0,360
CYLIND,0,D2/2,L-S,L,0,360


! Материальные константы пьезокерамики ЦТС-19
RO1=7.51e3
C11E=12.2e10
C12E=7.4e10
C13E=7.1e10
C33E=11.0e10
C66E=(C11E-C12E)/2.
C44E=C33E/5.6
E31=-1.7
E33=16.4
E15=10.3
EPS11=913
EPS33=873

! Материал 1 - поляризованная керамика
MP,DENS,1,RO1
TB,ANEL,1
TBDATA,1,C11E,C12E,C13E
TBDATA,7,C11E,C13E
TBDATA,12,C33E
TBDATA,16,C66E
TBDATA,19,C44E
TBDATA,21,C44E
TB,PIEZ,1
TBDATA,3,E31
TBDATA,6,E31
TBDATA,9,E33
TBDATA,14,E15
TBDATA,16,E15
MP,PERX,1,EPS11
MP,PERZ,1,EPS33

!Материал 2 - Абстрактная инерционная Масса
RO2=10e3
E2=1e10
NU2=0.3
Q2=10
MP,DENS,2,RO2
MP,EX,2,E2
MP,NUXY,2,NU2
MP,DAMP,2,omdi/Q2

! Реальные константы массы
R,2,Mass

MP,DAMP,3,0

ET,1,SOLID5,3
ET,2,SOLID45
ET,3,MASS21,,,4

VSEL,S,VOLU,,1 $ VATT,2,,2 $ VSEL,ALL
VSEL,S,VOLU,,2 $ VATT,2,,2 $ VSEL,ALL
VSEL,S,VOLU,,3 $ VATT,1,,1 $ VSEL,ALL

! склеивание объектов, отличающихся на величину, меньшую 1e-8
NUMMRG,ALL,1e-8
NUMCMP,ALL


LESIZE,ALL,,,10
VSWEEP,ALL


!*Нанесение электродов
CSYS,0
ESEL,S,MAT,,1
NSLE,S
NSEL,R,LOC,Z,L
CP,1,VOLT,ALL
CP,3,UZ,ALL
CSYS,0
ESEL,S,MAT,,1
NSLE,S
NSEL,R,LOC,Z,L-S
CP,2,VOLT,ALL
CP,4,UZ,ALL
NSEL,S,CP,,1
*GET,N_VOLT1,NODE,,NUM,MIN
NSEL,S,CP,,2
*GET,N_VOLT2,NODE,,NUM,MIN
NSEL,ALL

NSEL,S,LOC,Z,L
D,ALL,UX,0,,,,UY,UZ
NSEL,ALL
ALLSEL


/solu
ANTYPE,TRANS
OUTRES,BASIC,ALL
TINTP,,0.25, 0.5, 0.5
KBC,1

ACEL,0,0,0
NSUBST,10
TIME,0.01
SOLVE

ACEL,0,0,-2e+5
NSUBST,200
TIME,0.02
SOLVE

ACEL,0,0,0
NSUBST,100
TIME,0.05
SOLVE

В примере сила каким то невероятным образом исчезает... Она равна ~ 3500Н, когда расчетная должна быть 52000Н
В задаче без массы и ускорения но с приложенной силой все ок.
Задачи абсолютно эквивалентные, и расчетные силы в них практически одинаковы в теории

Offtop: PS: Если кто укажет на косяки не относящиеся к сабжу буду только рад

[СергейД]

вы задавали ускорение несколькоми способами, но ни разу правильно
надо (система СИ)
mp,dens,1,10000
acel,0,0,10

[vsae]

Мне необходимо задать ускорение 200000 м/с2, а не 10 м/с2
Не понимаю зачем менять плотность материала
В чем фишка поясните пожалуйста