24小时热门版块排行榜

返回列表

Mr_Husky

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 17
帖子: 32
在线: 7.1小时
虫号: 5303599
注册: 2016-12-01
性别: GG
专业: 结构工程

[求助] 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

回复此楼