【请教】两个构件模态分析频率结果

ljq2008

最近想用ANSYS做下两个构件的振动分析。其实就是两个简单的构件的不同方向的振动模态叠加问题。
（1）第一个是一个杆件的扭转振动形态，单独情况下，用PIPE16单元 分析结果一阶扭转频率为6.98
以下是命令流
/PREP7
ET,1,PIPE16
ET,2,MASS21
R,1,19E-3,0.9E-3
R,2,,,,,,0.32
MP,EX,1,199.16E9
MP,DENS,1,7750
MP,NUXY,1,0.3
/VIEW,1,1,1,1
N,1
N,NODENUM,,,520E-3
FILL
E,1,2
EGEN,(NODENUM-1),1,1
TYPE,2
REAL,2
E,NODENUM
FINISH
!进入求解模块
/SOLU
ANTYPE,MODAL
MODOPT,SUBSP,FREQY
D,1,UX,,,NODENUM,,UY,UZ,ROTX,ROTY
D,1,ROTZ
!设置输入信息 并开始求解
OUTPR,,1
SOLVE
FINISH
!进入后处理模块 提前前两阶频率的计算结果
/POST1
*GET,FREQ1,MODE,1,FREQ
*GET,FREQ2,MODE,2,FREQ
*GET,FREQ3,MODE,3,FREQ
（2）第二个杆件是一个简单的悬臂梁，单独情况下用beam3单元分析得到前两阶的模态频率为0.5579 13.42
以下是命令流
*SET,NODENUM,50
*SET,FREQY,5
/PREP7
ET,1,BEAM3
ET,2,MASS21,,,4
R,1,15E-6,1.25E-12,15E-3
R,2,0.42
MP,EX,1,200E9
MP,DENS,1,7750
MP,NUXY,1,0.3
!定义模型节点和单元连接关系
N,1
N,NODENUM,520E-3
FILL
E,1,2
EGEN,(NODENUM-1),1,1
TYPE,2
REAL,2
E,NODENUM
FINISH
!进入求解模块
/SOLU
ANTYPE,MODAL
MODOPT,SUBSP,FREQY
D,1,ALL
!设置输入信息 并开始求解
OUTPR,,1
SOLVE
FINISH
（3）将这两个构件放在一起，由于都是一维单元，并且PIPE16的自由度可以兼容beam3的自由度，所以对公共节点没有进行任何处理
结果得到的 两个构件的前三阶的频率分别为 0.56 7.0 13.43
以下是命令流
FINISH
/CLE
*SET,NODENUM,6
*SET,FREQY,6
/PREP7
ET,1,PIPE16
ET,2,MASS21
ET,3,BEAM3
ET,4,MASS21,,,4
R,1,19E-3,0.9E-3
R,2,,,,,,0.32
R,3,15E-6,1.25E-12,15E-3
R,4,0.42
MP,EX,1,199.16E9
MP,DENS,1,7750
MP,NUXY,1,0.3
MP,EX,2,200E9
MP,DENS,2,7750
MP,NUXY,2,0.3
!定义模型节点和单元连接关系

N,1
N, NODENUM,,,520E-3
FILL
TYPE,1
REAL,1
MAT,1
E,1,2
EGEN,(NODENUM-1),1,1
TYPE,2
REAL,2
E,NODENUM

N,2*NODENUM-1,520E-3,0,520E-3
FILL,NODENUM,2*NODENUM-1
TYPE,3
REAL,3
MAT,2
E,NODENUM,NODENUM+1
EGEN,NODENUM-1,1,NODENUM+1
TYPE,4
REAL,4
E,2*NODENUM-1
FINISH
!进入求解模块
/SOLU
ANTYPE,MODAL
MODOPT,SUBSP,FREQY
D,1,UX,,,NODENUM,,UY,UZ,ROTX,ROTY
D,1,ROTZ
OUTPR,,1
SOLVE
FINISH

问题提出
明显 7.0就是杆的扭转频率，但是如果用理论分析的话，由于梁本身的质量和梁末端质量的影响，实际导致杆的转动惯量变大，而杆的扭转频率应该降低？这是什么原因？希望高手给指点下！

ljq2008

有没有高手能解答下！