基础激励的谐响应分析——详解

╄→尐↘猪︶ㄣ 提交于 2020-03-08 04:39:16

=================命令==============================

1.转速与角速度
w = 2*pi*f = 2*pi*n
w角速度 rad/s
f 频率 1/s
n转速 r/s

2.定义单元截面
SECTYPE,截面ID,截面类型,截面子类型
sectype命令的格式由单元决定。
当截面类型为beam时,可以定义截面子类型。其中截面子类型为csolid表示实心圆面。

3.定义截面的几何形状
SECDATA,半径,周向划分的数量,径向划分的数量
secdata命令的格式由单元决定。
当截面为beam,且子类型为csolid时,命令格式如上所示。

4.指定单元
TYPE,单元ID

5.指定单元的截面
SECNUM,截面ID

6.定义线的划分的数量,三种命令

定义某条线被划分的数量
命令:LESIZE,线ID,,,每条线上的单元数量

定义关键点最近的单元的边长
命令:KESIZE,关键点ID,单元边长
只影响关键点附近线的长度,且只划分没有被LESIZE定义的线。
如果使用了SMRTSIZE智能划分,那么之前定义的KESIZE可能失效。

定义每条线被划分的数量
ESIZE,单元边长
ESIZE,,每条线上的单元数量
只划分除KESIZE LESIZE命令定义之外的线。
如果使用了智能划分,那么定义的ESIZE失效。

所以优先级LESIZE>KESIZE>ESIZE。 (我粗略地认为)
先定义好如何被划分,然后再使用划分命令,进行划分网格。

7.划分网格

划分线,生成节点与线单元
LMESH,线ID1,线ID2,增量INC
从线1到线2,增量为INC,这些线都将被划分

划分面,生成节点与面单元
AMESH,面ID1,面ID2,增量INC 

8.定义实常数
R,实常数ID,实常数1,实常数2,实常数3,实常数4,实常数5,实常数6

R命令只能定义最多6个实常数,当还有实常数需要定义时,可在R命令后面继续用RMORE命令。
RMORE,实常数7,实常数8,实常数9,实常数10,实常数11,实常数12
注意RMORE必须紧跟在R后面。

9.指定单元实常数
REAL,实常数ID

10.定义一个刚性区域
刚性区域可以是一个刚性连接(两个节点),刚性面,刚体。
此例中,先在原点建立了一个质量点单元,此质量点代表刚性夹具,质量点有惯性属性:即质量与转动惯量。
然后选择转轴端面圆心处的节点,再将该节点与原点处的节点建立刚性连接。
选择盘的内径处的所有结点,再将这些节点与原点处的节点建立一个刚性面。
(注:我并未区分结点与节点)

命令:
cerig,主节点,从节点,all   在主节点与从节点之间建立刚性区域,建立所有自由度关联的等式,即建立6个等式。

11.node函数
表示距离原点最近的结点编号
ncent = node(0,0,0)  将原点处的节点编号赋给ncent

12.谐响应分析选项
HROPT,谐响应分析方法,谐响应计算的最大模态数
谐响应分析方法包括:完全法FULL 振型叠加法MSUP 变分法(VT 、VTPA、VTRU)

13.控制求解结果写入数据库
outres,all,none 除了积分点位置LOCI与状态变量SVAR外,其余的求解结果都不保存到数据库。
outres,nsol,all 将节点自由度解保存到数据库,每个子步保存一次。
保存频率:ALL 这是谐响应分析和扩展通道默认的    LAST这是静力学与瞬态分析默认的。
ALL指每个子步都进行求解结果的保存。
LAST指只在每个载荷步的最后一个子步保存求解结果。

14.定义加速度
ACEL,X,Y,Z     定义x y z方向的加速度

15.定义载荷步内,子步施加的形式
分为阶跃与斜坡
kbc,0 斜坡
kbc,1 阶跃

16.定义谐响应分析的频率范围
HARFRQ,freqb,freqe
即从freqb开始分析,每间隔(FREQE-FREQB) / NSBSTP分析一次,直到freqe。

HARFRQ,freqb,freqe,--,对数选项
则表示从freqb开始分析,每间隔(log(FREQE)-log(FREQB)) / (NSBSTP-1)分析一次,直到freqe。

其中nsbstp是由NSUBST命令定义的,表示每个载荷步内的子步数量。
--表示保留的参数位置,以后会开发出来的。

此例中
harfrq,,fexcit 
nsubst,500 
表示 谐响应分析的频率范围为:0,1,2,3,4,....,直到500Hz

17.阻尼矩阵
betad, 刚度矩阵乘数   即阻尼矩阵C=系数*K
alphad,质量矩阵乘数  即C=系数*M

18.对于振型叠加法进行载荷向量换算
在振型叠加法谐响应分析与振型叠加法瞬态分析中,所有的载荷向量需要在第一个载荷步之前进行换算。
LVSCALE,1

19.扩展通道

打开扩展通道
EXPASS,on
扩展振型叠加法和子结构分析的求解结果
NUMEXP,all   扩展所有载荷步的子步求解结果

20.模态扩展

定义模态分析与屈曲分析中模态被扩展,同时被写入结果文件的数量。
MXPAND,30 扩展前30阶模态,同时写入结果文件。

21.后处理
nsol,2,node(0,0,0),U,X,uxTip    2号变量=原点处的节点的X方向的位移
/gropt,logy,on     y轴方向的间隔是对数
/axlab,x,FREQUENCIES  x轴的标题
/axlab,y,DISPLACEMENTS (m)   y轴的标题
plvar,2,3,4  显示2 3 4 号变量的曲线
extreme,2,4,1 列出2 3 4号变量画出曲线的极值   

EXTREME,变量1,变量n,增量inc  列出变量1,变量1+inc,变量1+2*inc......,变量n的极值。


=====================感想=============================

1.(在此例中)模态分析中,考虑转速,阻尼,并在y方向加一个加速度 -1 m/s2,然后进行模态分析,分析得到30阶模态下的结果。
接着进行谐响应分析,将模态分析中的计算出的载荷向量原封不动的引入到谐响应中(使用命令LVSCALE,1),同时取消加速度,因为已经把模态分析中的y轴负方向的加速度引入了。然后进行一个频率范围内的激励(方向为:沿着y轴负方向)。最后得到谐响应结果。
2.只有运行solve命令时,才开始执行载荷步。
 

参考
http://tieba.baidu.com/p/5208651536 因为我自己翻译出来后比较困惑,然后查看了一下别人的说法,验证一下。
http://blog.sina.com.cn/s/blog_700544f10100wh7n.html

标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!