|
第一步：首先建立模型，如图，一个很简单的角形结构梁
其中上面的号码是节点号，连接1，10的就是调频器，节点2与1均有质量。
finish
/clear,start
/prep7
k,1
k,2,10000
k,3,5000,10000
et,1,beam44
R,1,186160,178980000000,56238000000,1278,680,128980000000,
RMODIF,1,21,1278,680,
RMODIF,1,25,76160,110000,0,0,7.0678e-009,0,
UIMP,1,EX, , ,210000,
UIMP,1,DENS, , ,1.1261e-5,
l,2,3,4
l,3,1,4
lsel,all
lmesh,all
et,2,mass21
r,2,1000,1000,1000
type,2
real,2
e,1
e,2
dk,1,all
finish
第二步：进行模态分析
finish
/solu
antype,modal
modopt,subsp,10
mxpand,10
solve
finish
自己进入后处理器，察看结果
TIME 1 TIME
TIME
0.13996 0.139960
0.23500 0.235003
0.25406 0.254063
0.43206 0.432056
0.95366 0.953658
1.6187 1.61867
1.8375 1.83753
2.0020 2.00205
2.5345 2.53452
3.0464 3.04638
再进入扩展模态，找出对应想抑制震动的主自由度
在上面就以二阶和三阶为例吧，他们的振型是在y，z平面，这也好理解，这个方向刚度很低，所以低频率就出现在这个振型上。所以上面加的tmd的方向也很讲究，也就是节点10的位置是垂直于y，z平面，不然，在其他方向就没有好的效果了。
加tmd：
finish
/prep7
et,3,combin14
r,3,100,10
r,4,100,100,100
csys,4
kwpave,2
n,10,0,0,-1000
type,3
real,3
e,1,10
type,2
real,4
e,10
finish
然后进入阻尼模态分析，由于关心的频率范围是前两阶，因此，设置了频率范围在0.1到0.2，
在此强调一下，在damp中频率范围一定要设置，不然感兴趣的结果有很多遗漏。
finish
/solu
antype,modal
modopt,damp,10
MXPAND,10, , ,0
MODOPT,DAMP,10,0.1,0.2, ,OFF
solve
finish
然后进入后处理器
0.13772
0.16125
有兴趣者可以对tmd命令流中的r,3,100,10中的阻尼进行变化，
就可发现阻尼的作用了
对于低频作用不大，但对高频有很好的抑制效果
弹簧刚度的作用相反。
|