24小时热门版块排行榜    

查看: 1865  |  回复: 1

Mr_Husky

新虫 (初入文坛)

[求助] ansys提取刚度问题

我用简支梁进行ansys整体刚度矩阵K和整体质量矩阵M的提取
ansys提取的模态是按照质量归一话的振型phi
按理说phi'*M*phi = I(单位矩阵)
           phi'*K*phi = Λ(圆频率平方的矩阵)

但是我提取出来的,matlab算只有phi'*M*phi = I(单位矩阵) 这个满足
phi'*K*phi 算出来的平动位移上差别很大,然而扭转却是正确的,感觉很纳闷,希望大牛解惑,不甚感激

简支梁长L = 300m ,宽B=40m,竖向抗弯刚度EI_z = 2.1E12 Pa*m4,横向抗弯刚度EI_y = 1.8E13 Pa*m4 扭转刚度GI_t = 4.1E11 Pa*m4
桥面质量 m = 2000kg/m3 ,质量惯性矩 I_m = 4.5E6 kg*m2/m
命令流如下


FINI
/UIS,MSGPOP,3
/CLEAR
/PREP7
ET,1,BEAM4
ET,2,MASS21
N,1,0,0
N,31,300,0,0
FILL,1,31
R,1,10,10,85.714,40,0.25,,,     !面积10,IZZ10,IYY85.714, TKZ40,TKY0.25
RMORE,0,5.076,,,,,,             !IXX5.076
R,2,,,,45000000,,,,                 !TKZ4500000
MP,EX,1,2.1E11                   !弹性模量  
MP,DENS,1,2000 !质量密度
MP,GXY,1,8.077E10  !剪切模量
TYPE,1
REAL,1
*DO,I,1,30
E,I,I+1
*ENDDO
TYPE,2          !集中质量
REAL,2
*DO,I,2,30
E,I
*ENDDO
/SOLU
D,1,UX,0,0,,,UY,UZ,ROTX
D,31,UY,0,0,,,UZ,ROTX
ANTYPE,MODAL       !Specifies the analysis type and restart status.
MODOPT,LANB,10     !Specifies modal analysis options.
MXPAND,10       !Specifies the number of modes to expand
LUMPM,ON
SOLVE !无阻尼特征值分析,获取前10阶模态频率初值


/post1
set,list
回复此楼

» 收录本帖的淘帖专辑推荐

ANSYS命令流

» 猜你喜欢

O Ever Youthful O Ever Weeping
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xing_yan

铁虫 (初入文坛)

楼主搞定没,我想问问您那个ANSYS中模型建好后怎么样计算横向抗弯刚度和纵向抗弯刚度以及扭转刚度呀
2楼2018-07-11 15:48:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Mr_Husky 的主题更新
信息提示
请填处理意见