使用轴对称方法求解模态

njweiwei2

该事例对一个圆柱体进行模态求解，分别以建立整个模型和只建立1/4模型进行求解，然后进行比较。
/PREP7  
ET,1,SOLID95
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,2.5E5  
MPDATA,PRXY,1,,0.3  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,7.85E-9  
CYLIND,1,0.5,0,0.5,0,360,   
FLST,5,12,4,ORDE,4  
FITEM,5,1   
FITEM,5,-4  
FITEM,5,9   
FITEM,5,-16
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
LESIZE,_Y1, , ,12, , , , ,1
FLST,5,4,4,ORDE,2   
FITEM,5,17  
FITEM,5,-20
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
LESIZE,_Y1, , ,6, , , , ,1  
CM,_Y,VOLU  
VSEL, , , ,       1
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
VSWEEP,_Y1  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
ALLSEL  
CSYS,1  
NROTAT,ALL  
FINISH  
/SOL
ANTYPE,2
MSAVE,0
MODOPT,LANB,12  
EQSLV,SPAR  
MXPAND,12, , ,0
LUMPM,0
PSTRES,0
MODOPT,LANB,12,0,0, ,OFF
FLST,2,2,5,ORDE,2   
FITEM,2,5   
FITEM,2,-6  
/GO
DA,P51X,UX,
FLST,2,1,5,ORDE,1   
FITEM,2,1   
/GO
DA,P51X,UZ,
ALLSEL  
SOLVE   
FINISH  
/POST1  
SET,LIST
FINISH

!!!!!!!!以下为建立1/4模型，用轴对称方法求解

/PREP7  
VCLEAR,       1
wpro,,-90.000000,   
VSBW,       1   
wpro,,,-90.000000   
FLST,2,2,6,ORDE,2   
FITEM,2,2   
FITEM,2,-3  
VSBW,P51X   
FLST,2,3,6,ORDE,3   
FITEM,2,1   
FITEM,2,4   
FITEM,2,-5  
VDELE,P51X, , ,1
ET,2,MESH200
KEYOPT,2,1,7
KEYOPT,2,2,0
MSHAPE,0,2D
MSHKEY,0
CM,_Y,AREA  
ASEL, , , ,       8
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
AMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
FLST,5,4,4,ORDE,4   
FITEM,5,18  
FITEM,5,20  
FITEM,5,23  
FITEM,5,-24
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
LESIZE,_Y1, , ,6, , , , ,1  
CM,_Y,AREA  
ASEL, , , ,       8
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
ACLEAR,_Y1  
AMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
FLST,3,1,5,ORDE,1   
FITEM,3,8   
AGEN,2,P51X, , , ,90, , ,0  
NUMMRG,ALL, , , ,LOW
CM,_Y,VOLU  
VSEL, , , ,       6
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
VSWEEP,_Y1  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
NROTAT,ALL  
CYCLIC  
FINISH  
/SOL
ANTYPE,2
MSAVE,0
MODOPT,LANB,8   
EQSLV,SPAR  
MXPAND,8, , ,0  
LUMPM,0
PSTRES,0
MODOPT,LANB,8,0,0, ,OFF
FLST,2,1,5,ORDE,1   
FITEM,2,27  
/GO
DA,P51X,UX,
FLST,2,1,5,ORDE,1   
FITEM,2,25  
/GO
DA,P51X,UZ,
ALLSEL  
SOLVE   
FINISH  
/PREP7  
FINISH  
/SOL
FINISH  
/POST1  
SET,LIST, , , , , , ,ORDE

njweiwei2

两者结果对比：
建立整个模型的结果：
  *****  INDEX OF DATA SETS ON RESULTS FILE  *****

   SET   TIME/FREQ    LOAD STEP   SUBSTEP  CUMULATIVE
     1 0.46004E-01         1         1         1
     2 0.94863E+06         1         2         2
     3 0.94863E+06         1         3         3
     4 0.17707E+07         1         4         4
     5 0.17707E+07         1         5         5
     6 0.23377E+07         1         6         6
     7 0.23377E+07         1         7         7
     8 0.24158E+07         1         8         8
     9 0.24158E+07         1         9         9
    10 0.24175E+07         1        10        10
    11 0.24313E+07         1        11        11
    12 0.24313E+07         1        12        12
建立1/4模型的结果：
  *****  INDEX OF DATA SETS ON RESULTS FILE  *****

   SET   TIME/FREQ    LOAD STEP   SUBSTEP  CUMULATIVE  HRM-INDEX
     1 0.71644E-01         1         1         1         0
     9 0.94863E+06         2         1         1         1
    10 0.94863E+06         2         2         2         1
    17 0.17707E+07         3         1         1         2
    18 0.17707E+07         3         2         2         2
    11 0.23376E+07         2         3         3         1
    12 0.23376E+07         2         4         4         1
    13 0.24158E+07         2         5         5         1
    14 0.24158E+07         2         6         6         1
     2 0.24174E+07         1         2         2         0
    19 0.24313E+07         3         3         3         2
    20 0.24313E+07         3         4         4         2

njweiwei2

可以看到，两者结果几乎没有任何差别。对于大型轴对称结构，使用其中一部分模型进行求解可以节省大量资源，大大提高效率

inndoor

不错，用命令流的都是高手，学习了

YW3624021

我觉的楼主的结果是错误的，不可能没有区别，1/4模型的质量比整体的小，质量的变化必定会引起频率的变化
看看你下面的数据，set下面的，好像是作假了，呵呵！
SET   TIME/FREQ    LOAD STEP   SUBSTEP  CUMULATIVE  HRM-INDEX
     1 0.71644E-01         1         1         1         0
     9 0.94863E+06         2         1         1         1
    10 0.94863E+06         2         2         2         1
    17 0.17707E+07         3         1         1         2
    18 0.17707E+07         3         2         2         2
    11 0.23376E+07         2         3         3         1
    12 0.23376E+07         2         4         4         1
    13 0.24158E+07         2         5         5         1
    14 0.24158E+07         2         6         6         1
     2 0.24174E+07         1         2         2         0
    19 0.24313E+07         3         3         3         2
    20 0.24313E+07         3         4         4         2

htbbzzg

1  1/4 模型不是轴对称模型，而是一般的面对称模型，概念上有问题；
2  5 楼质疑者的概念也是不清楚，简单的说： 1/4 模型的质量和刚度同样缩小，因此频率应该不变。更深入的说，频率结果还与约束条件等有关。采用一般面对称模型计算频率时，对称面需要分别考虑对称约束和反对称约束的各种组合，才能和整体的结果进行比较；如果确实是轴对称模型，需要考虑轴对称振动以及各阶谐波的振动才能与整体结果比较。

Wallace998

谢谢分享！！

victao

那我可不可以先做一个矩形面的模态分析，然后利用三维应力扩展到一周呢？