我这有个和你类似的例子,希望对你有帮助:
梁的正弦加载的瞬态分析
模型描述:
1)悬臂梁顶端施加X向正弦变化的载荷,进行瞬态动力求解;
2)单元为beam188;
3)/post26中能提取末端的位移结果。
模型如下图,已经通过模态、谐响应计算(如下图),但瞬态计算总是无法得到理想的结果。
问题:
我认为,梁顶端的位移输出为一正弦曲线,但ANSYS计算结果的时间-位移历程总是无法得到正弦曲线,不知是何原因?我觉得是力载荷的加载有问题,主要出现在激励频率的确定上,但总是调试不出来。请高手指点。
命令流:
/FILNAME,BTSanlysis
/UNITS,SI
/TITLE,Bond Tool计算
D=3e-3!直径
L=15e-3 !长度
S=3.141593*(D/2)*(D/2)!截面积
I=3.141593*D*D*D*D/64
/PREP7!建模
ET,1,BEAM188!梁单元 3,188,189
MP,EX,1,620.53E9!弹性模量
MP,DENS,1,14880
MP,PRXY,1,0.29
!MP,DAMP,1,0.0085
secnum,1
sectype,1,beam,CSOLID
secdata,D/2,10,3!
N,1,0,0
N,21,0,-L
FILL,1,21
E,1,2
EGEN,20,1,1,1,1
/PNUM,NODE,1!图示节点及编号
/PNUM,ELEM,1!图示单元及编号
EPLOT!显示单元
FINISH
/SOLU
ANTYPE,2!模态分析
EQSLV,SPAR
MXPAND,100, , ,1
LUMPM,0
PSTRES,0
MODOPT,LANB,100,5000,100000, ,OFF
OUTPR,BASIC,ALL
d,1,uy,,,,,uz,rotx,roty,rotz
SOLVE
FINISH
/POST1
!做谐响应分析
/SOLU
ANTYPE,HARMIC! HARMONIC RESPONSE ANALYSIS
HROUT,OFF! Print results as amplitudes and phase angles
HROPT,full
OUTPR,BASIC,1
harfrq,5000,100000! FREQUENCY RANGE FROM 0 TO 7.5 HZ
NSUBST,200! 30 INTERVALS WITHIN FREQ. RANGE
KBC,1! STEP BOUNDARY CONDITION
d,1,uy,,,,,uz,rotx,roty,rotz!very important!!
!D,1,UX,10
F,1,FX,200
ALLS
SOLVE
FINISH
/POST26
NSOL,2,21,U,X,UX! STORE UX DISPLACEMENT OF NODE 21作为变!量2
!ESOL,3,21,,SMISC,1,F1! STORE 摩擦力FORCE F1 OF ELEMENT 1 AS VARIABLE 3
!由K1引起的弹性力F1.
!PRVAR,2,3! PRINT VARIABLES 2 AND 3
/GRID,1! TURN GRID ON
/AXLAB,Y,DISP! Y AXIS LABEL AS DISP
/GTHK,CURVE,2! CURVE LINES THICKNESS RATIO OF 2
PLVAR,2
!做瞬态分析
/SOLu
antype,trans
trnopt,full
outpr,all,all
outres,all,all
!alphad,0.014
!betad,0.000092
autots,on
kbc,1
!tintp,,0.25,0.5,0.5
d,1,uy,,,,,uz,rotx,roty,rotz
*do,i,1,2000,1
t=i/100000
time,t
nsubst,50
!施加约束条件
Pressure=10*sin(2*3.14159*91743*t) !如何设置激励力的频率?f为其中的FX方 向振动的固有频率?
F,1,FX,Pressure
solve
lswr,i
*enddo
fini
/POST26
NSOL,2,21,U,X,UX! STORE UX DISPLACEMENT OF NODE 21作为变!量2
!ESOL,3,21,,SMISC,1,F1! STORE 摩擦力FORCE F1 OF ELEMENT 1 AS VARIABLE 3
!由K1引起的弹性力F1.
!PRVAR,2,3! PRINT VARIABLES 2 AND 3
/GRID,1! TURN GRID ON
/AXLAB,Y,DISP! Y AXIS LABEL AS DISP
/GTHK,CURVE,2! CURVE LINES THICKNESS RATIO OF 2
PLVAR,2 |