有道理,我采用solid45单元可以看到有频率的改变了.放上来给有需要的人参考.
:lol
/COM, Structural
!*
/PREP7
BLOCK,0,100,0,10,0,10,
!------------
ET,1,SOLID45
!------------
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7980
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,1.90e11
MPDATA,PRXY,1,,0.30
!------------------------
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
MSHAPE,0,3d
MSHKEY,1
VMESH,_Y1
MSHKEY,0
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!-------------------------*
ANTYPE,0
FLST,2,1,5,ORDE,1
FITEM,2,5
!*
/GO
DA,P51X,ALL,0
FLST,2,2,1,ORDE,2
FITEM,2,75
FITEM,2,-76
!*
/GO
F,P51X,FY,4000
!FLST,2,1,5,ORDE,1
!FITEM,2,6
!/GO
!*
!SFA,P51X,1,PRES,400000
GPLOT
FINISH
!------------------------------------------------------------------------------
/SOL !静力大变形求解
ANTYPE,0
NLGEOM,1
PSTRES,ON !计及预应力效果
LUMPM,ON
TIME,1
AUTOTS,1
NSUBST,20, , ,1
KBC,0
SOLVE
FINISH
/SOLUTION
ANTYPE,2 !进行模态求解
MSAVE,0
MODOPT,LANB,10
MXPAND,10, , ,0 !取前十阶模态
PSTRES,1 !打开预应力效应
LUMPM,ON
MODOPT,LANB,10,0,0, ,OFF
UPCOORD,1,ON !修正坐标以得到正确的应力
PSOLVE,TRIANG !三角化矩阵
PSOLVE,EIGLANB !提取特征值和特征向量
FINISH
/SOLU
EXPASS,1 !扩展模态解
PSOLVE,EIGEXP
FINISH
/POST1
SET,LIST !观察结果
FINISH |