轴对称问题的有限元分析.docx
《轴对称问题的有限元分析.docx》由会员分享,可在线阅读,更多相关《轴对称问题的有限元分析.docx(52页珍藏版)》请在冰点文库上搜索。
![轴对称问题的有限元分析.docx](https://file1.bingdoc.com/fileroot1/2023-5/1/e379d133-e2cd-4a61-93ad-54f47d3344c9/e379d133-e2cd-4a61-93ad-54f47d3344c91.gif)
轴对称问题的有限元分析
第11章轴对称问题的有限元分析
第1节基本知识
本节的有限元对象为轴对称问题,目的是学习将3D问题转化为2D问题分析的轴对称方
法,涉及如何选取轴对称单元、建模规律、载荷的施加方法和后处理技术。
一、轴对称问题的定义
轴对称问题是指受力体的几何形状、约束状态,以及其它外在因素都对称于某一根轴(过该轴的任一平面都是对称面)。
轴对称受力体的所有应力、应变和位移均对称于这根轴。
二、用ANSYS解决2D轴对称问题的规定
用ANSYS解决2D轴对称问题时,轴对称模型必须在总体坐标系XOY平面的第一象限中
创建,并且Y轴为轴旋转的对称轴。
求解时,施加自由约束、压力载荷、温度载荷和Y方向的加速度可以像其它非轴对称模
型一样进行施加,但集中载荷有特殊的含义,它表示的是力或力矩在360°范围内的合力,
即输入的是整个圆周上的总的载荷大小。
同理,在求解完毕后进行后处理时,轴对称模型输出的反作用力结果也是整个圆周上的合力输出,即力和力矩按总载荷大小输出。
在ANSYS中,X方向是径向,Z方向是环向,受力体承载后的环向位移为零,环向应力和应变不为零。
常用的2D轴对称单元类型和用途见表11-1。
表11-12D轴对称常用结构单元列表
板单元名称
说明
Shell51
圆锥壳单元,可以生成圆柱壳或环行圆盘单元,具有非线形材料性质,用于解决轴对称问题
Shell61
圆锥壳单元,可以生成圆柱壳或环行圆盘单元,没有非线形材料性质,用于解决轴对称问题
Plane42
二维结构实体单元,用于解决平面应力问题、平面应变问题和轴对称问题
Plane82
Plane42的高阶单元
Plane182
二维结构实体单元,用于解决平面应力问题、平面应变问题和轴对称问题
Plane183
Plane182的高阶单元
在利用ANSYS进行有限元分析时,将这些单元定义为新的单元后,设置单元配置项
KEYOPT3)为Axisymmetric(Shell51和Shell61单元本身就是轴对称单元,不用设置该项),
单元将被指定按轴对称模型进行计算。
后处理时,可观察径向和环向应力,它对应的是SX与SZ应力分量,并且在直角坐标系
下观察即可。
可以通过轴对称扩展设置将截面结果扩展成任意扇型区域大小的模型,以便更加真实地
观察总体模型的各项结果。
第2节2D轴对称问题有限元分析实例
、案例1――圆柱筒的静力分析
如图11-1所示,圆柱筒材质为A3钢,受1000N/m的压力作用,其厚度为0.1m,直
径12m高度为16m,并且圆柱筒壳的下部轴线方向固定,其它方向自由,试计算其变形、径向应力和轴向应力。
弹性模量为x1011N/m2,泊松比为。
以圆柱筒底部中心为坐标原点,建立直角坐标系如图11-1所示,标出主要点(1点和2
点)的坐标,为实体造型做好准备。
制定分析方案。
分析类型为线弹性性材料,结构静力分析,轴对称问题,由于受力题为圆柱壳,选用Shell51单元,筒的厚度为0.1m为单元的实常数;边界条件为圆柱筒下部施加轴线方向固定支撑,2点的受力为1000*12*n等于37699N。
1.ANSYS^析开始准备工作
(1)清空数据库并开始一个新的分析选取Utility>Menu>File>Clear&StartNew,
弹出CleardatabaseandStartNew对话框,单击OK按钮,弹出Verify对话框,单击OK按钮完成清空数据库。
(2)指定新的工作文件名指定工作文件名。
选取Utility>Menu>File>Change
Jobname,弹出ChangeJobname对话框,在EnterNewJobname项输入工作文件名,本例中输入的工作文件名为"yuanzhutong”,单击OK按钮完成工作文件名的定义。
(3)指定新的标题指定分析标题。
选取Utility>Menu>File>ChangeTitle,弹出ChangeTitle对话框,在EnterNewTitle项输入标题名,本例中输入"shellproblem”为标题名,然后单击OK按钮完成分析标题的定义。
(4)重新刷新图形窗口选取Utility>Menu>Plot>Replot,定义的信息显示在图形窗
口中。
(5)定义结构分析运行主菜单MainMenu>Preferences,出现偏好设置对话框,赋值分析模块为Structure结构分析,单击OK按钮完成分析类型的定义。
2.定义单元
运行主菜单MainMenu>Preprocessor>ElementType>Add/Edit/Delete命令,弹出
ElementTypes对话框,单击"Add”按钮新建单元类型,弹出LibraryofElementTypes
对话框,先选择单元大类为Shell,接着选择2DaxisSymmetr51(Shell51),单击OK按钮,
完成单元类型选择,单击Close按钮完成设置,如图11-2所示。
dnlputriKturia]*刑右曲are
LKhrAjrofRlle*^alTypes
图11-2定义单元类型Shell
3•定义实常数
运行主菜单MainMenu>Preprocessor>RealConstantsAdd/Edit/Delete命令,弹出实
常数定义对话框,单击Add按钮进入实例常量输入对话框,在TK(I)项输入板厚(实例常
数)0.1米,单击OK按钮完成实例常量输入。
回到实例常量对话框,此时显示出新建编号为1的实例常量,单击Close按钮完成输入,如图11-3所示。
4.定义材料属性
运行主菜单MainMenu>Preprocessor>MaterialProps>MaterialModels命令,系统显
示材料属性设置对话框,在材料属性对话框中依次选择Structure、Linear、Elastic、
Isotropic,如图11-4所示。
完成选择后,弹出材料属性输入对话框,分别输入弹性模量2e11,泊松比,如图11-5
所示,单击0K按钮完成材料属性输入并返回图11-4。
完成材料属性设置后,单击对话框右上方“X”按钮离开材料属性设置。
命令弹出创建关键点对话
Menu>Preprocessor>Modeling>Create>Keypoints>InactiveCS
框,在对话框中分别输入关键点编号1,X=6,Y=0,Z=0,单击Apply按钮定义第一点,如
图11-6所示。
在对话框中分别输入关键点编号2,X=6,Y=16,Z=0,单击OK按钮定义第二点,如图
11-7所示。
(2)创建圆柱筒母线运行主菜单MainMenu>Preprocessor>Modeling>Create>lines>
Lines>StraightLine命令,在弹出的拾取对话框中分别拾取1点和2点,单击OK按钮完
成圆柱筒母线绘制,图11-8所示。
图11-6定义关键点1
图11-7定义关键点2
6.划分网格
运行主菜单MainMenu>Preprocessor>Meshing>MeshTool(网格划分工具)命令,出现
MeshTool菜单,单击Lines中的Set按钮,出现选择拾取对话框,单击PickAll按钮选择
所有线,在单元尺寸控制对话框中的No.ofelementdivisions项中输入线的等分数为8,
其它项默认,单击OK按钮确定。
单击Mesh按钮划分网格,在出现的MeshLines对话框中单击PickAll按钮,系统将自动完成网格划分。
运行菜单UtilityMenu>PlotCtrls>Numbering弹出PlotNumberingControls选择
对话框,显示单元编号,设置如图11-9所示,单击OK按钮完成设置,划分网格结果如图
11-10所示。
图11-9设置显示单元编号
图11-10划分网格结果
7.加载
(1)施加约束运行主菜单MainMenu>Solution>DefineLoads>Apply>Structural>
Displacement>OnKeypoints命令,出现拾取菜单,选择关键点1,单击OK按钮出现约束定
义对话框,如图11-11所示,选择UY约束Y方向自由度,在DisplacementValue选项输入0,其它项默认,再单击0K按钮,完成约束定义。
(2)施加载荷运行主菜单MainMenu>Solution>DefineLoads>Apply>Structural>
Force/Moment>OnKeypoints命令,出现拾取菜单,拾取关键点2,单击OK按钮出现载荷定
义对话框,如图11-12所示,载荷类型为FY,载荷数值为-37699N,再单击OK按钮完成载荷的施加,结果如图11-13所示。
图11-11施加约束
[FK1ApplyPorc^/loB<@ntoriKe’poirit険
LM0irection□f
Ifvaluerthan:
VALUEForce/ioi^rit
图11-12施加载荷
&求解
运行主菜单MainMenu>Solution>CurrentLS命令,出现SolveCurrentLoadStep对
话框,单击/STATCommand窗口菜单/STATCommand>File>Close关闭/STATCommanc窗口,然后单击SolveCurrentLoadStep菜单中OK按钮确定,计算机开始进行求解,求解完成后出现"Solutionisdone”提示表示求解完成,单击Close按钮完成求解。
选择菜单路径MainMenu>Finish退出求解器。
8查看分析结果
(1)扩展获得360度模型运行下拉菜单UtilityMenu〉PlotCrls>Style>Symmetry
Expansion>2DAxis-Symmetric弹出轴对称扩展设置对话框,选择Fullexpansion,单击
OK按钮完成操作,图形窗口将显示扩展后的图形。
(2)打开单元形状显示开关运行下拉菜单UtilityMenu〉PlotCrls>Style>Sizeand
Shape弹出SizeandShape对话框,将Displayofelementshapesbasedonrealconstantdescriptions设置为On即打开状态,将ReplotuponOK/Apply设置为Replot,单击OK按
钮完成操作,图形窗口将显示图形的截面形状。
(3)显示节点(单元)位移云图运行主菜单MainMenu>GeneralPostproc>Plot
Results>ContourPlot>NodalSolu(orElementSolu)命令,选择DOF
Solution>DisplacementVectorsum合位移,单击OK按钮,节点位移云图如图11-14所示。
(4)显示径向应力(SX)运行主菜单MainMenu>GeneralPostproc>PlotResults>
ContourPlot>NodalSolu(orElementSolu)命令,选择Stress>X-Componentstress径
向应力,单击OK按钮,节点应力云图如图11-15所示。
(5)显示环向应力(SZ)运行主菜单MainMenu>GeneralPostproc>PlotResults>
ContourPlot>NodalSolu(orElementSolu)命令,选择Stress>Z-Componentstress
环向应力,单击OK按钮,节点应力云图如图11-16所示。
5TEP»1
-I
心*JUtDKS
U^UHfAWl
B£X5-0
Dffi-s^SE-0«
Sbfit*.EBiBE-O7
sm.S05E-O«|i
AN
UH2€SMft;歸
.3B6E-O7
1HE-04
.24«E-€4
.127E-&6
.40TE-O6.5C>gE-06.?
^£-4S
.407£->01«_«46E-46,IO5EM9ti
图11-14显示变形图
flaaALSOLUTI-CIII
STEP«l
SUB-1
tz»o«l
sxCiWj
R£¥S«4>
DMX■4师EC-0首smbw-sna
SMI七Id
-1434
■Wi
-1071
121«S
图11-15径向应力(SX)等值云图
STW-1
S'UBwl
T£l«»l
£Z1AVQ>
RST£-^
DHX-a7l«£-0«却■■一豹酊
5MX*13-.2f
-9T85
-esv«
-7506
-4517
-5<2B-3.249--1070
-4331-215ft1«.2B
图11-16环向应力(SZ)等值云图
、案例2——飞轮的静力分析
如图11-17所示是一飞轮的截面图。
飞轮角速度为rad/s,飞轮边缘受压力作用,压力p为1MPa飞轮轴孔固定,飞轮尺寸如图(单位为mr),试进行静力分析。
弹性模量为210GPa泊松比为,密度为7800kg/m3。
2
R5
3V
6
5■
50
图11-17飞轮截面图
军题过程
建立直角坐标系如图11-17所示,算出1—7点的坐标为实体造型做好准备。
制定分析方案。
分析类型为线弹性性材料,结构静力分析,轴对称问题,采用板单元,选用Plane82单元;边界条件飞轮的角速度引起的惯性和飞轮边缘受的压力。
1.ANSY盼析开始准备工作
(1)清空数据库并开始一个新的分析选取Utility>Menu>File>Clear&StartNew,
弹出CleardatabaseandStartNew对话框,单击OK按钮,弹出Verify对话框,单击OK按钮完成清空数据库。
(2)指定新的工作文件名指定工作文件名。
选取Utility>Menu>File>Change
Jobname,弹出ChangeJobname对话框,在EnterNewJobname项输入工作文件名,本例中输入的工作文件名为“feilun”单击OK按钮完成工作文件名的定义。
(3)指定新的标题指定分析标题。
选取Utility>Menu>File>ChangeTitle,弹出
ChangeTitle对话框,在EnterNewTitle项输入标题名,本例中输入“axisymmetric
problem"为标题名,然后单击OK按钮完成分析标题的定义。
(4)重新刷新图形窗口选取Utility>Menu>Plot>Replot,定义的信息显示在图形窗
口中。
2.确定分析类型
运行主菜单MainMen>preferences弹出分析类型设定对话框,选择分析模块为
Structural结构分析,然后单击OK按钮完成分析类型定义。
3.定义单元及材料属性
(1)新建单元类型运行主菜单MainMenu>Preprocessor>ElementType>Add/Edit/
Delete命令,弹出ElementTypes对话框,单击Add按钮新建单元类型,弹出Libraryof
ElementTypes对话框,先选择单元大类为Solid,接着选择Quad8node82(plane82),
单击OK按钮,完成单元类型选择,如图11-18所示。
接着单击Option按钮进入单元设置选项,在Elementbehavior(K3)栏中更改选项为
Axisymmetric(轴对称),在单击OK按钮返回ElementType对话框,单击Close按钮完成设置,如图11-19所示。
(2)定义材料属性运行主菜单MainMenu>Preprocessor>MaterialProps>Material
Models命令,系统显示材料属性设置对话框,在材料属性对话框中依次选择Structure、
Linear、Elastic、Isotropic,如图11-20所示。
完成选择后,弹出材料属性对话框图11-21,分别输入弹性模量210e9,泊松比,单击
OK返回图11-20,选择Densiy,弹出密度定义对话框图11-22,输入密度7800,单击OK按
钮返回。
完成材料属性设置后,单击对话框右上方“X”按钮离开材料属性设置。
OplionvforPI.ANEA2^E)easntTyp此Rrf.No.9
图11-19设置单元配置项
图11-20进入材料属性设置
LinearI&otroolcHaterLalPropertiesforlaterIalMuab&rI
T1
T&bjj1
PRtY
AddT£'BperH.turt-IhjLuttTi'iipeiH.turc
图11-21设置材料属性
Densit$rforKaleriaJNunberL
11
AddTe»perature|BeletcTeiiteralure|Graj)h
VOKJJHelp
图11-22定义材料密度
4.建立几何图形
(1)建立关键点关键点1—7各点的坐标见表11-2。
(2)绘制矩形运行主菜单MainMenu>Preprocessor>Modeling>Create>Areas>Rec-
tangle>ByDimensions命令,在对话框中分别输入1点坐标X仁50,丫仁0,2点坐标X2=55,Y2=50,,单击Apply按钮完成第一个矩形绘制;同理,输入3点坐标X1=55,丫1=24,6点坐标X2=75,Y2=16,,单击Apply按钮完成第二个矩形绘制;输入4点坐标X1=75,Y1=40,5点坐标X2=80,Y2=5,,单击OK按钮完成第三个矩形绘制,如图11-23所示。
表11-2几何模型关键点坐标
关键点号
在活动坐标系的X,Y,Z坐标
1
(50,0,0)
2
(55,50,0)
3
(55,24,0)
4
(75,40,0)
5
(80,5,0)
6
(75,16,0)
图11-23绘制矩形
(3)布尔操作合并图形运行主菜单MainMenu>Preprocessor>Modeling>Operate>
Booleans>Add出现选择拾取对话框,单击PickAll按钮选择所有面,完成布尔加操作。
(4)设置显示方式为了下一步倒圆角方便,设置显示图中的各直线编号。
运行菜单
UtilityMenu〉PlotCtrls>Numbering弹出PlotNumberingControls选择对话框,设置
LineNumbers为On,其它项默认,设置如图11-23所示,单击OK按钮完成设置。
运行菜单UtilityMenu>Plot>Line,设置显示方式为直线,以便下一步为直线倒圆角,
为如图11-24所示。
(5)倒圆角运行主菜单MainMenu>Preprocessor>Modeling>Create>Lines>Line
Fillet弹出拾取对话框,操作如下:
拾取L7、L18,单击OK按钮,出现圆角半径设置对话框,FilletRadius项输入5,其它项默认,单击Apply按钮完成第一个圆角绘制,如图11-25所示;同理,拾取L7、L20,
单击OK按钮,出现圆角半径设置对话框,FilletRadius项输入5,其它项默认,单击Apply
按钮完成第二个圆角绘制;拾取L5、L17,单击OK按钮,出现圆角半径设置对话框,FilletRadius项输入5,其它项默认,单击Apply按钮完成第三个圆角绘制;拾取L5、L19,单击OK按钮,出现圆角半径设置对话框,FilletRadius项输入5,其它项默认,单击OK按钮
完成圆角绘制,如图11-26所示。
尢P」胞1也C.V'-UiZtXfjl
1/PHUiJPlatNumberingtkmtrols
KFKevpoLntnuibers
roff
LINELinenhibers
FOn
AREAAreanuabers
Off
VOLUVolujmenuiibers
Off
MODENodentubers
roff
Elco/Attribruabcri
MoniiHhering审
TARNT7J雌N)u■住幕
Off
SVAl.Hurler1ccodttruf
Off
If/MUM;Hiinberinftahornwith
GoLors*nij-Bbcn▼
r/kEPLOT]KcploluphiOft/AppLrl
Rep1(31皆
OK:
ApplyCancel|Help
图11-23设置显示方式
图11-25设置圆角半径
LII&
LliK1TLH
L3
laX
AN
2120恥01:
21:
10
L4
■-L2
1^2
LIO
jF
L1.7
图11-26绘制圆角
(6)生成面运行菜单UtilityMenu>Plot>Line设置直线显示方式,运行主菜单Main
Menu>Preprocessor>Modeling>Create>Areas>Arbitrary>Bylines,出现选择拾取对话框,
拾取L2、L8、L6单击Apply按钮完成第一个面的生成,如图11-27所示;同理,完成其它
三个面的生成,单击0K按钮结束。
AN
LIVES
Af战21
LIVEHW
Ki4BAl
LI
MF
Li4i
LI4
Lie
L21
■2