在线做家装设计的网站/网络营销seo培训
旋转软化与应力刚化
平板旋转
两种单元:shell181,solid185。分析5种情况:静止,绕y=-0.1处的x轴转3200rpm,绕y=-0.1处的x轴转9600rpm,绕y=-0.1处的z轴转3200rpm,绕y=-0.1处的z轴转9600rpm。
只提取一阶模态频率,因此得到10个结果。
验证手册指出:平板绕x轴旋转时,应力刚化与旋转软化都影响弯曲模态;绕z轴旋转时旋转软化对弯曲模态的影响可以忽略。
这里的应力刚化与旋转软化都是旋转时的离心力造成的,先进行 nonlinear static analysis ,再进行 linear perturbation modal analysis,计算特征值用的是Block-Lanczos eigensolver。
举个例子:
一根杆有弹性无质量,一端有一个质量球,绕着一根轴旋转,离心力会把杆拉长,离心力就会变大,这就是旋转软化;但是杆轴向受力使得横向刚度增大,横向弯曲模态频率也就提高,这就是应力刚化。
SHELL181适合分析溥的及中等厚度的板壳结构零件。它一个单元有四个结点,每个结点有六个自由度,即沿X,Y,Z方向的移动自由度和转动自由度。(如果选了“薄膜”选项,则单元只能用移动自由度)。退化“三角形”选项,只能过滤单元化分网中使用。
SHELL181 支持线性,大扭转和大应变,变厚度非线性分析。它既能用完全法也可用缩减法,可用于分布载荷及应力刚化。
SHELL181可用于多层材料模型。板壳建模的精度由剪切变形理论(Mindlin-Reissner shell theory)决定。
在遇到有收敛困难的时候,用SHELL181代替SHELL43单元,该单元详细介绍,请参看ANSYS, In
solid185单元有8节点,每个节点有三个移动自由度。因此对比上下两个图,y=0处的节点自由度都被约束了,但上图多了黄色的部分是旋转自由度。
代码
由于有*VWRITE命令,只能批处理运行
/clear
/COM,ANSYS MEDIA REL. 150 (11/8/2013) REF. VERIF. MANUAL: REL. 150
/VERIFY,VM281
/TITLE,VM281,EFFECT OF STRESS STIFFENING AND SPIN SOFTENING ON A ROTATING PLATE
/COM,
/COM, REFERENCE: NASA TECHNICAL MEMORANDUM 89861,PAGE 14,FIGURE 7
/COM,
/COM, =====================================================================
/COM, MODELING THE PLATE USING 4-NODE STRUCTURAL SHELL ELEMENT (SHELL181)
/COM, ======================================================================
/PREP7
lg = 152.4e-3 ! LENGTH IN METERS
wd = 50.8e-3 ! WIDTH IN METERS
tk = 2.54e-3 ! THICKNESS IN METERS
K,1,-wd/2,0,0 ! Defines a keypoint.
K,2,+wd/2,0,0
K,3,+wd/2,lg,0
K,4,-wd/2,lg,0
A,1,2,3,4 ! Defines an area by connecting keypoints.
/COM, MATERIAL PROPERTIES OF STRUCTURAL STEEL PLATE
MP,EX,1,2.0e+11 ! 1号材料的弹性模量
MP,DENS,1,7850 ! 1号材料的密度
MP,NUXY,1,0.3 ! 1号材料的泊松比
ET,1,SHELL181 ! 4-NODE STRUCTURAL SHELL ELEMENT
SECTYPE,1,SHELL ! 定义壳单元截面
SECDATA,tk,1,10,3 ! THICKNESS 厚度 层坐标系相对单元坐标系的角度 积分点
TYPE,1 !设置当前单元
MAT,1 !设置当前材料
SECNUM,1 !设置当前截面
ESIZE,5e-3 !划分单元的长度
AMESH,1 !划分面单元
NSEL,S,LOC,Y,0.0 !选中y=0的所有节点
D,ALL,ALL ! CONSTRAINING ALL DOF AT BASE
ALLSEL,ALL
CM,PLATE,ELEM ! FORMING A COMPONENT WITH ALL ELEMENTS
SAVE,SHELL,DB ! SAVE THE DATABASE
FINISH/COM, ----------------------
/COM, NO ROTATION 静止
/COM, ----------------------
/OUT,SCRATCH
/SOLUTION
ANTYPE,MODAL ! 模态分析
MODOPT,LANB,4 ! EXTRACT 4 MODES USING BLOCK-LANCZOS
MXPAND,4 !http://blog.sina.com.cn/s/blog_9b0ec2070102z2mu.html
SOLVE
/OUT,
*GET,FREQ1_SH,MODE,1,FREQ ! *GET Solution Items, Entity = MODE 提取一阶模态频率赋给FREQ1_SH
FINISH
/DELETE,,rst ! 删除rst文件/COM, --------------------------------------
/COM, PLATE ROTATED ALONG X AXIS AT 3200 RPM 绕x轴3200rpm 模态分析加入了离心力产生的预应力
/COM, --------------------------------------
/OUT,SCRATCH ! /OUTPUT的Notes
/SOLUTION
ANTYPE,STATIC ! 静力分析
OFFSET= -0.1
pi = acos(-1)
CMOMEGA,PLATE,3200*pi/30,,,0,OFFSET,0,1,OFFSET,0 ! ROTATING THE PLATE ALONG X AXIS
NLGEOM,ON ! NON-LINEAR ANALYSIS
TIME,1.0 ! 载荷步时间为1
AUTOTS,ON ! 自动载荷步
NSUBS,10,100,1 ! 子步:10步 最少1步 最大100步
OUTRES,ALL,ALL
SOLVE
FINISH
/SOLUTION
ANTYPE,STATIC,RESTART,,,PERTURB ! LINEAR PERTURBATION ANALYSIS 预应力模态分析
PERTURB,MODAL, ! MODAL ANALYSIS,RETAIN CMOMEGA
SOLVE,ELFORM ! REGENERATE STIFFNESS MATRIX 重新生成单元矩阵
MODOPT,LANB,4
MXPAND,4
SOLVE
/OUT,
*GET,FREQ2_SH,MODE,1,FREQ ! *GET Solution Items, Entity = MODE 提取一阶模态频率赋给FREQ2_SH
FINISH
/DELETE,,rst
/DELETE,,rstp/COM, --------------------------------------
/COM, PLATE ROTATED ALONG X AXIS AT 9600 RPM
/COM, --------------------------------------
/OUT,SCRATCH
/SOLUTION
ANTYPE,STATIC
OFFSET= -0.1
pi = acos(-1)
CMOMEGA,PLATE,9600*pi/30,,,0,OFFSET,0,1,OFFSET,0 ! ROTATING THE PLATE ALONG X AXIS
NLGEOM,ON ! NON-LINEAR ANALYSIS
TIME,1.0
AUTOTS,ON
NSUBS,10,100,1
OUTRES,ALL,ALL
SOLVE
FINISH/SOLUTION
ANTYPE,STATIC,RESTART,,,PERTURB ! LINEAR PERTURBATION ANALYSIS
PERTURB,MODAL, ! MODAL ANALYSIS,RETAIN CMOMEGA
SOLVE,ELFORM ! REGENERATE STIFFNESS MATRIX
MODOPT,LANB,4
MXPAND,4
SOLVE
/OUT,
*GET,FREQ3_SH,MODE,1,FREQ ! *GET Solution Items, Entity = MODE 提取一阶模态频率赋给FREQ3_SH
FINISH
/DELETE,,rst
/DELETE,,rstp/COM, --------------------------------------
/COM, PLATE ROTATED ALONG Z AXIS AT 3200 RPM
/COM, --------------------------------------
/OUT,SCRATCH
/SOLUTION
ANTYPE,STATIC
OFFSET= -0.1
pi = acos(-1)
CMOMEGA,PLATE,3200*pi/30,,,0,OFFSET,0,0,OFFSET,1 ! ROTATING THE PLATE ALONG Z AXIS
NLGEOM,ON ! NON-LINEAR ANALYSIS
TIME,1.0
AUTOTS,ON
NSUBS,10,100,1
OUTRES,ALL,ALL
SOLVE
FINISH
/SOLUTION
ANTYPE,STATIC,RESTART,,,PERTURB ! LINEAR PERTURBATION ANALYSIS
PERTURB,MODAL, ! MODAL ANALYSIS,RETAIN CMOMEGA
SOLVE,ELFORM ! REGENERATE STIFFNESS MATRIX
MODOPT,LANB,4
MXPAND,4
SOLVE
/OUT,
*GET,FREQ4_SH,MODE,1,FREQ
FINISH
/DELETE,,rst
/DELETE,,rstp/COM, --------------------------------------
/COM, PLATE ROTATED ALONG Z AXIS AT 9600 RPM
/COM, --------------------------------------
/OUT,SCRATCH
/SOLUTION
ANTYPE,STATIC
OFFSET= -0.1
pi = acos(-1)
CMOMEGA,PLATE,9600*pi/30,,,0,OFFSET,0,0,OFFSET,1 ! ROTATING THE PLATE ALONG Z AXIS
NLGEOM,ON ! NON-LINEAR ANALYSIS
TIME,1.0
AUTOTS,ON
NSUBS,10,100,1
OUTRES,ALL,ALL
SOLVE
FINISH
/SOLUTION
ANTYPE,STATIC,RESTART,,,PERTURB ! LINEAR PERTURBATION ANALYSIS
PERTURB,MODAL, ! MODAL ANALYSIS,RETAIN CMOMEGA
SOLVE,ELFORM ! REGENERATE STIFFNESS MATRIX
MODOPT,LANB,4
MXPAND,4
SOLVE
/OUT,
*GET,FREQ5_SH,MODE,1,FREQ
FINISH/OUT,SCRATCH
*DIM,LABEL,CHAR,5,2
*DIM,VALUE,,5,3
LABEL(1,1) = 'None'
LABEL(1,2) = '0rpm'
LABEL(2,1) = 'X axis'
LABEL(2,2) = '3200rpm'
LABEL(3,1) = 'X axis'
LABEL(3,2) = '9600rpm'
LABEL(4,1) = 'Z axis'
LABEL(4,2) = '3200rpm'
LABEL(5,1) = 'Z axis'
LABEL(5,2) = '9600rpm'
*VFILL,VALUE(1,1),DATA,90
*VFILL,VALUE(1,2),DATA,FREQ1_SH
*VFILL,VALUE(1,3),DATA,(90/FREQ1_SH)
*VFILL,VALUE(2,1),DATA,108
*VFILL,VALUE(2,2),DATA,FREQ2_SH
*VFILL,VALUE(2,3),DATA,(108/FREQ2_SH)
*VFILL,VALUE(3,1),DATA,195
*VFILL,VALUE(3,2),DATA,FREQ3_SH
*VFILL,VALUE(3,3),DATA,(195/FREQ3_SH)
*VFILL,VALUE(4,1),DATA,121
*VFILL,VALUE(4,2),DATA,FREQ4_SH
*VFILL,VALUE(4,3),DATA,(121/FREQ4_SH)
*VFILL,VALUE(5,1),DATA,250
*VFILL,VALUE(5,2),DATA,FREQ5_SH
*VFILL,VALUE(5,3),DATA,(250/FREQ5_SH)
SAVE,TABLE_1
FINISH/CLEAR,NOSTART
/OUT,
/COM, =====================================================================
/COM, MODELING THE PLATE USING 8-NODE STRUCTURAL SOLID ELEMENT (SOLID185)
/COM, ======================================================================
/PREP7
lg = 152.4e-3 ! LENGTH IN METERS
wd = 50.8e-3 ! WIDTH IN METERS
tk = 2.54e-3 ! THICKNESS IN METERS
K,1,-wd/2,0,0
K,2,+wd/2,0,0
K,3,+wd/2,lg,0
K,4,-wd/2,lg,0
K,5,-wd/2,0,tk
K,6,+wd/2,0,tk
K,7,+wd/2,lg,tk
K,8,-wd/2,lg,tk
V,1,2,3,4,5,6,7,8
/COM, MATERIAL PROPERTIES OF STRUCTURAL STEEL PLATE
MP,EX,1,2.0e+11 MP,DENS,1,7850
MP,NUXY,1,0.3
ET,1,SOLID185 ! 3-D 8 NODE STRUCTURAL SOLID ELEMENT
KEYOPT,1,2,3 ! SIMPLIFIED ENHANCED STRAIN FORMULATION
TYPE,1
MAT,1
ESIZE,5e-3
VMESH,1
NSEL,S,LOC,Y,0.0
D,ALL,ALL ! CONSTRAINING ALL DOF AT BASE
ALLSEL,ALL
CM,PLATE,ELEM ! FORMING A COMPONENT WITH ALL ELEMENTS
SAVE,SOLID,DB ! SAVE THE DATABASE
FINISH
/COM, ----------------------
/COM, NO ROTATION
/COM, ----------------------
/OUT,SCRATCH
/SOLUTION
ANTYPE,MODAL
MODOPT,LANB,4 ! EXTRACT 4 MODES USING BLOCK-LANCZOS
MXPAND,4
SOLVE
/OUT,
*GET,FREQ1_SO,MODE,1,FREQ
FINISH
/DELETE,,rst
/COM, --------------------------------------
/COM, PLATE ROTATED ALONG X AXIS AT 3200RPM
/COM, --------------------------------------
/OUT,SCRATCH
/SOLUTION
ANTYPE,STATIC
OFFSET= -0.1
pi = acos(-1)
CMOMEGA,PLATE,3200*pi/30,,,0,OFFSET,0,1,OFFSET,0 ! ROTATING THE PLATE ALONG X AXIS
NLGEOM,ON ! NON-LINEAR ANALYSIS
TIME,1.0
AUTOTS,ON
NSUBS,10,100,1
OUTRES,ALL,ALL
SOLVE
FINISH
/SOLUTION
ANTYPE,STATIC,RESTART,,,PERTURB ! LINEAR PERTURBATION ANALYSIS
PERTURB,MODAL, ! MODAL ANALYSIS,RETAIN CMOMEGA
SOLVE,ELFORM ! REGENERATE STIFFNESS MATRIX
MODOPT,LANB,4
MXPAND,4
SOLVE
/OUT,
*GET,FREQ2_SO,MODE,1,FREQ
FINISH
/DELETE,,rst
/DELETE,,rstp
/COM, --------------------------------------
/COM, PLATE ROTATED ALONG X AXIS AT 9600RPM
/COM, --------------------------------------/OUT,SCRATCH
/SOLUTION
ANTYPE,STATIC
OFFSET= -0.1
pi = acos(-1)
CMOMEGA,PLATE,9600*pi/30,,,0,OFFSET,0,1,OFFSET,0 ! ROTATING THE PLATE ALONG X AXIS
NLGEOM,ON ! NON-LINEAR ANALYSIS
TIME,1.0
AUTOTS,ON
NSUBS,10,100,1
OUTRES,ALL,ALL
SOLVE
FINISH
/SOLUTION
ANTYPE,STATIC,RESTART,,,PERTURB ! LINEAR PERTURBATION ANALYSIS
PERTURB,MODAL, ! MODAL ANALYSIS,RETAIN CMOMEGA
SOLVE,ELFORM ! REGENERATE STIFFNESS MATRIX
MODOPT,LANB,4
MXPAND,4
SOLVE
/OUT,
*GET,FREQ3_SO,MODE,1,FREQ
FINISH
/DELETE,,rst
/DELETE,,rstp
/COM, --------------------------------------
/COM, PLATE ROTATED ALONG Z AXIS AT 3200RPM
/COM, --------------------------------------
/OUT,SCRATCH
/SOLUTION
ANTYPE,STATIC
OFFSET= -0.1
pi = acos(-1)
CMOMEGA,PLATE,3200*pi/30,,,0,OFFSET,0,0,OFFSET,1 ! ROTATING THE PLATE ALONG Z AXIS
NLGEOM,ON ! NON-LINEAR ANALYSIS
TIME,1.0
AUTOTS,ON
NSUBS,10,100,1
OUTRES,ALL,ALL
SOLVE
FINISH
/SOLUTION
ANTYPE,STATIC,RESTART,,,PERTURB ! LINEAR PERTURBATION ANALYSIS
PERTURB,MODAL, ! MODAL ANALYSIS,RETAIN CMOMEGA
SOLVE,ELFORM ! REGENERATE STIFFNESS MATRIX
MODOPT,LANB,4
MXPAND,4
SOLVE
/OUT,
*GET,FREQ4_SO,MODE,1,FREQ
FINISH
/DELETE,,rst
/DELETE,,rstp
/COM, --------------------------------------
/COM, PLATE ROTATED ALONG Z AXIS AT 9600RPM
/COM, --------------------------------------
/OUT,SCRATCH
/SOLUTION
ANTYPE,STATIC
OFFSET= -0.1
pi = acos(-1)
CMOMEGA,PLATE,9600*pi/30,,,0,OFFSET,0,0,OFFSET,1 ! ROTATING THE PLATE ALONG Z AXISNLGEOM,ON ! NON-LINEAR ANALYSIS
TIME,1.0
AUTOTS,ON
NSUBS,10,100,1
OUTRES,ALL,ALL
SOLVE
FINISH
/SOLUTION
ANTYPE,STATIC,RESTART,,,PERTURB ! LINEAR PERTURBATION ANALYSIS
PERTURB,MODAL, ! MODAL ANALYSIS,RETAIN CMOMEGA
SOLVE,ELFORM ! REGENERATE STIFFNESS MATRIX
MODOPT,LANB,4
MXPAND,4
SOLVE
/OUT,
*GET,FREQ5_SO,MODE,1,FREQ
FINISH
/OUT,SCRATCH
*DIM,LABEL,CHAR,5,2
*DIM,VALUE,,5,3
LABEL(1,1) = 'None'
LABEL(1,2) = '0rpm'
LABEL(2,1) = 'X axis'
LABEL(2,2) = '3200rpm'
LABEL(3,1) = 'X axis'
LABEL(3,2) = '9600rpm'
LABEL(4,1) = 'Z axis'
LABEL(4,2) = '3200rpm'
LABEL(5,1) = 'Z axis'
LABEL(5,2) = '9600rpm'
*VFILL,VALUE(1,1),DATA,90
*VFILL,VALUE(1,2),DATA,FREQ1_SO
*VFILL,VALUE(1,3),DATA,(90/FREQ1_SO)
*VFILL,VALUE(2,1),DATA,108
*VFILL,VALUE(2,2),DATA,FREQ2_SO
*VFILL,VALUE(2,3),DATA,(108/FREQ2_SO)
*VFILL,VALUE(3,1),DATA,195
*VFILL,VALUE(3,2),DATA,FREQ3_SO
*VFILL,VALUE(3,3),DATA,(195/FREQ3_SO)
*VFILL,VALUE(4,1),DATA,121
*VFILL,VALUE(4,2),DATA,FREQ4_SO
*VFILL,VALUE(4,3),DATA,(121/FREQ4_SO)
*VFILL,VALUE(5,1),DATA,250
*VFILL,VALUE(5,2),DATA,FREQ5_SO
*VFILL,VALUE(5,3),DATA,(250/FREQ5_SO)
SAVE,TABLE_2
FINISH
/NOPR
RESUME,TABLE_1 !
/GOPR
/OUT,vm281,vrt
/COM,
/COM, -----------------------------VM281 RESULTS COMPARISON-----------------------
/COM,
/COM, | TARGET | MECHANICAL APDL | RATIO
/COM,
/COM,
/COM, USING SHELL181 ELEMENTS
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3) ! https://www.cnblogs.com/gaozihan/p/11200897.html
(1X,A8,' ',A8,' ',F12.3,' ',F12.3,' ',F8.3)
/COM,
/COM, USING SOLID185 ELEMENTS
/COM,
/NOPR,
RESUME,TABLE_2
/GOPR,
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,' ',A8,' ',F12.3,' ',F12.3,' ',F8.3)
/COM,
/COM,---------------------------------------------------------------------------------
/OUT,
*LIST,vm281,vrt ! 将结果输出到vm281.vert
FINISH
http://blog.sina.com.cn/s/blog_153e0668a0102x1zg.html 大变形分析中(NLGEOM,ON)包含应力刚化效应(SSTIF,ON)
https://www.bridgehead.com.cn/thread-3367-1-1.html
输出结果
-----------------------------VM281 RESULTS COMPARISON-----------------------| TARGET | MECHANICAL APDL | RATIOUSING SHELL181 ELEMENTSNone 0rpm 90.000 90.831 0.991X axis 3200rpm 108.000 108.273 0.997X axis 9600rpm 195.000 195.139 0.999Z axis 3200rpm 121.000 120.686 1.003Z axis 9600rpm 250.000 252.359 0.991USING SOLID185 ELEMENTSPRINTOUT RESUMED BY /GOPNone 0rpm 90.000 90.875 0.990X axis 3200rpm 108.000 108.324 0.997X axis 9600rpm 195.000 195.237 0.999Z axis 3200rpm 121.000 120.725 1.002Z axis 9600rpm 250.000 252.419 0.990-------------------------------------------------------------------------------
再补充一个例子
一个小插曲:命令d,1,all,10,9 误输入为 d,1,all,10.9
D, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6
d,1,all,10.9 表示将1-11节点的自由度都约束
d,1,all,10.4 表示将1-10节点的自由度都约束,因此nend这个参数是四舍五入的,ninc步长默认为1
d,1,all,10.9,3 表示将1,4,7,10节点的自由度都约束
d,1,all,10,9 表示将节点1与10约束
fini
/clear/prep7
et,1,shell63,,,,,,,1
r,1,3e-3
mp,ex,1,217e9
mp,nuxy,1,0.3
mp,dens,1,7850
n,1,-0.014,,0.15
n,9,-0.014,,0.478
fill
ngen,2,9,1,9,1,0.028
e,1,2,11,10
egen,8,1,-1
finish! no stress stiff and no spin soft
/solu
d,1,all,,,10,9
antype,modal
modopt,subsp,5
solve
!/post1
!set,list
*GET,FREQ1,MODE,1,FREQ
fini! with stress stiff and no spin soft
/solu
antype,static
pstres,on
omega,314.159265
outpr,,1
solve
fini/solu
antype,modal
modopt,subsp,5
pstres,on
omega,314.159265,,,0
solve
!/post1
!set,list
*GET,FREQ2,MODE,1,FREQ
fini! with stress stiff and with spin soft
/solu
antype,modal
modopt,subsp,5
pstres,on ! 包含应力刚化
omega,314.159265,,,1 ! 包含旋转软化
solve
!/post1
!set,list
*GET,FREQ3,MODE,1,FREQ
fini
一阶固有频率
no stress stiff and no spin soft:23.96
with stress stiff and no spin soft:72.21
with stress stiff and with spin soft:52.01
对比
两段都是考虑应力刚化与旋转软化,但使用的代码不同
shell63不支持预应力模态分析,改用shell181进行对比
方1
fini
/clear/prep7
et,1,shell181,,,,,,,1
r,1,3e-3
mp,ex,1,217e9
mp,nuxy,1,0.3
mp,dens,1,7850
n,1,-0.014,,0.15
n,9,-0.014,,0.478
fill
ngen,2,9,1,9,1,0.028
e,1,2,11,10
egen,8,1,-1
finish! 方1
/solu
d,1,all,,,10,9
antype,static
omega,314.159265
NLGEOM,ON
TIME,1.0 ! 载荷步时间为1
AUTOTS,ON ! 自动载荷步
NSUBS,10,100,1 ! 子步:10步 最少1步 最大100步
outpr,,1
solve
fini/SOLUTION
ANTYPE,STATIC,RESTART,,,PERTURB
PERTURB,MODAL
SOLVE,ELFORM
modopt,subsp,5
solve
*GET,method1_freq1,MODE,1,FREQ
*GET,method1_freq2,MODE,2,FREQ
*GET,method1_freq3,MODE,3,FREQ
*GET,method1_freq4,MODE,4,FREQ
*GET,method1_freq5,MODE,5,FREQ
fini
方2
fini
/clear/prep7
et,1,shell181,,,,,,,1
r,1,3e-3
mp,ex,1,217e9
mp,nuxy,1,0.3
mp,dens,1,7850
n,1,-0.014,,0.15
n,9,-0.014,,0.478
fill
ngen,2,9,1,9,1,0.028
e,1,2,11,10
egen,8,1,-1
finish! 方2
/solu
d,1,all,,,10,9
antype,static
pstres,on
omega,314.159265
outpr,,1
solve
fini/solu
antype,modal
modopt,subsp,5
pstres,on ! 包含应力刚化
omega,314.159265,,,1 ! 包含旋转软化
solve
*GET,method2_freq1,MODE,1,FREQ
*GET,method2_freq2,MODE,2,FREQ
*GET,method2_freq3,MODE,3,FREQ
*GET,method2_freq4,MODE,4,FREQ
*GET,method2_freq5,MODE,5,FREQ
fini
NLGEOM只能用于静力与全瞬态分析。
NLGEOM,ON
Large-deflection effects有两种:大挠度(大旋转)效应,大应变效应,取决于单元类型。打开NLGEOM,也会自动包含应力刚化效应
对比两种方法前五阶模态频率
略微有区别
结果为何不同
两段代码都是验证手册的,前者是15版的,后者没在15版中找到,我猜是旧的,(查看2021版命令)现在OMEGA命令只有三个参数,没有第四个参数设置旋转软化,而在 CORIOLIS提到旋转软化ansys默认考虑。
热流固间的耦合
流固耦合
热固耦合
热流耦合
热流固耦合