ansys导线舞动APDL

SECOFFSET, CENT
SECDATA,R,0,0,0,0,0,0,0,0,0,0,0
K,1,,,,
K,2,L1,,,
K,3,L1+L2,,,
K,4,L1+L2+L3,,,
K,5,L1+L2+L3+L4,,,
K,6,,,JG,
K,7,L1,,JG,
*SET,%_FNCNAME%(0,2,1), 0.0, -1, 9, 1, -1, 0, 0
*SET,%_FNCNAME%(0,3,1), 0, -2, 0, 1, 17, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, 3, 4, 18
*SET,%_FNCNAME%(0,5,1), 0.0, -1, 10, 1, -1, 0, 0
*DIM,_FNC_C3,,1
*SET,_FNC_C1(1),-fl
*SET,_FNC_C2(1),v
*SET,_FNC_C3(1),fd
*SET,_FNCCSYS,0
! /INPUT,fh.func,,,1
*DIM,%_FNCNAME%,TABLE,6,8,1,,,,%_FNCCSYS%
*DIM,_FNC_C3,,1
*SET,_FNC_C1(1),fl
*SET,_FNC_C2(1),v
*SET,_FNC_C3(1),fd
*SET,_FNCCSYS,0
! /INPUT,fl.func,,,1
*DIM,%_FNCNAME%,TABLE,6,8,1,,,,%_FNCCSYS%
K,8,L1+L2,,JG,
K,9,L1+L2+L3,,JG,
K,10,L1+L2+L3+L4,,JG,
LSTR, 1, 2
LSTR,2,3
LSTR,3,4
LSTR,4,5
LSTR,6,7
LSTR,7,8
LSTR,8,9
LSTR,9,10
FH=1.458
FV=1.39
MT=2.28
FD=MD*xc*V*V*CD/2/1000000
Fl=MD*xc*V*V*Cl/2/1000000
Fm=MD*xc*V*V*Cm/2/1000000
/PREP7
ET,1,BEAM188
SECTYPE, 1, BEAM, CSOLID, , 0
LDQD=LDL/S !导线拉断强度
MD=1.349/S !导线密度
TXML=4.5E9 !导线弹性模量
BSB=0.35 !导线泊松比
JG=1/MD/S !中间棒的长度
!函数计算参数
cd=2.887
cl=1.575
cm=-5.55e-1
fini
save
LMESH, all
!边界条件
DK,1, , , ,0,ALL, , , , , ,
DK,6, , , ,0,ALL, , , , , ,
DK,5, , , ,0,ALL, , , , , ,
DK,10, , , ,0,ALL, , , , , ,
ACEL,0,9.8/100000,0,
*DEL,fh
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNCCSYS
*SET,_FNCNAME,'fh'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
!
! Begin of equation: FL*cos({Y}/V)-FD*sin({Y}/V)
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
FINISH
/CLEAR
R=0.02682/2 !导线的半径
LDL=103900 !导线的拉断力
PI=3.1415926
XC=204.3 !导线总长
V=10 !风速
L1=50
L2=45
L3=55
L4=50
!换算
S=R*R*PI !导线截面积
!-->
*DEL,fv
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNCCSYS
*SET,_FNCNAME,'fv'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*SET,%_FNCNAME%(0,6,1), 0.0, -3, 0, 1, 19, 3, -1
*SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 1, -2, 2, -3
*SET,%_FNCNAME%(0,8,1), 0.0, 99, 0, 1, -1, 0, 0
! End of equation: FL*cos({Y}/V)-FD*sin({Y}/V)
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 1, 3, 4, 18
LSTR,2,7
LSTR,3,
LSTR,4,9
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,txml
MPDATA,PRXY,1,,bsb
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,md
ESIZE,XC/30,0,
/REPLOT
!****************************************舞动过程中的位移场结果
/post1
*DO,J,1,100
SET,,, ,,, ,j
PLNSOL, U,SUM, 0,1.0
*ENDDO
/seg,off,,EQV,0.4
anim,1,1 !*********************************生成动画完毕
OUTRES,ERASE
OUTRES,ALL,ALL
TIME,1
F,all,FX, %FH%
F,all,FY, %FV%
alls
solve
!*********************舞动动画 *********************
/post1
/seg,dele !允许图形数据存储到局部终端存储器
!-->
/solu
ANTYPE,0
TIME,1
alls
solve
/PREP7
UPGEOM,1,LAST,LAST,'file','rst',' '
/solu
ANTYPE,4
TRNOPT,FULL
LUMPM,0
NSUBST,100,200,100
*SET,%_FNCNAME%(0,2,1), 0.0, -1, 10, 1, -1, 0, 0
*SET,%_FNCNAME%(0,3,1), 0, -2, 0, 1, 17, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1, 3, 4, 18
*SET,%_FNCNAME%(0,5,1), 0.0, -1, 9, 1, -1, 0, 0
/seg,multi,wudong,0.4 !存动画,分节存储后来的显示,动画名字,每个画面间延迟系数
!*************************************************位移场结果 !从结果文件读,载荷步,子步,比例因子实部,比例因子虚部,时间
/dscale,1,1 !显示位移时,设置放大系数
avprin,0,0 !规定主项量和总向量如何被计算
avres,1 !规定当powergraphies可用时,结果数据如何被平均
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 1, 3, 4, 18
*SET,%_FNCNAME%(0,6,1), 0.0, -3, 0, 1, 19, 3, -1
*SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 1, -2, 2, -3
*SET,%_FNCNAME%(0,8,1), 0.0, 99, 0, 1, -1, 0, 0
! End of equation: -FL*SIN({Y}/V)-FD*COS({Y}/V)
!
! Begin of equation: -FL*SIN({Y}/V)-FD*COS({Y}/V)
*SET,%_FNCNAME%(0,0,1), 0.0, -999