Двойной маятник Циглера

[Vassa_vl 0]

Здравствуйте!
Я моделирую двойной маятник Циглера под действием мёртвой (постоянного направления) и следящей нагрузках в ANSYS Mechanical APDL. Вопрос именно в следящей нагрузке, я не понимаю, как её приложить. Обе силы должны быть приложены к концу маятника. Я создаю элемент FOLLW201, но как именно через него задать следящую нагрузку? Рассматриваю динамический расчёт 

[Борман 2 375]

6 часов назад, Vassa_vl сказал:

Вопрос именно в следящей нагрузке, я не понимаю, как её приложить.

С этой следящей нагрузкой каждый раз как первый. Заново все вспоминаю.

 

Но она работает это точно.

 

[Vassa_vl 0]

Спасибо, я посмотрела ссылку, там и правда сила следящая, но может быть у Вас сохранился макрос, чтобы понять, как её правильно задать? 
Я попробовала добавить следящую нагрузку, приложу макрос, но у меня она постоянного направления получается... 

 

finish 
/clear 
/BATCH  

/prep7
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
et,1,mass21
keyopt, 1,3,2

et,2,mpc184
keyopt,2,1,1
keyopt,2,2,1

et,3,mpc184
keyopt,3,1,6
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
n,1,0,0
n,2,0.75,0
n,3,1,0
n,4,1,0
n,5,1.75,0
n,6,2,0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!==========
ET,4,FOLLW201
r,4,0,1,0
type,4
real,4
e,6 

esel,s,type,,4
sfe,all,1,pres,,130e3
allsel

/pbc,f,,1
/psf,pres,,2
eplot
!!==========
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
type,1
real,1
m = 390
r,1,m 
E,2
E,5
E,6
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
type,2
real,2
e,1,2
e,2,3
e,4,5
e,5,6
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
local, 11, 0, 0, 0, 0, , , 90
sectype, 3, JOIN, REVO, TESTREVO
secjoin, , 11, 11
TB, JOIN, 3,,, SK4
TBDATA, 3, 3.0E4
type, 3
real, 3
secnum, 3
e,3,4
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
F,6,FX,-100   
eplot
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
mp, ex, 1, 2e11
mp, nuxy, 1, 0.3
mp, density, 1, 7.8e3
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
d, 1, all
ddel, 1, rotz
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
finish

/solu

vel = 2
ic, 6, uy, 0.0, vel
antype, trans
time, 0.1
kbc, 1
nlgeom, on
nsub, 50 ,50, 50
trnopt, full, , , , , HHT
tintp, 0.2
outres, all, all
solve

time, 2
midtol, on, 10
nsub, 200, 1e6, 200
trnopt, full, , , , , HHT
tintp, 0.05
outres, all, all
solve

finish
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/POST1  
/SHOW,WIN32C
SET,FIRST   
/PLOPTS,INFO,3  
/CONTOUR,ALL,18 
/PNUM,MAT,1 
/NUMBER,1   
/REPLOT,RESIZE  
PLDISP,2
SET,,,,,,,26
PLDISP,2
SET,,,,,,,51
PLDISP,2 
/SHOW,WIN32 
/REPLOT,RESIZE  

/pbc,f,,1
/psf,pres,,2
 

[Борман 2 375]

12 часа назад, Vassa_vl сказал:

но может быть у Вас сохранился макрос

Не поверите, сохранился в моей то помойке.

Там фишка буквально в 2-х строчках.

go.txt