ansys作业.docx
《ansys作业.docx》由会员分享,可在线阅读,更多相关《ansys作业.docx(29页珍藏版)》请在冰豆网上搜索。
ansys作业
目录
问题描述………………………………………………………….2
ANSYS分析步骤………………………………………………….3
加载并求解…………………………………………………….14
结果分析…………………………………………………….18
建议和体会………………………………………………………24
一、问题描述
在旋转机械中通常会遇到轴与轴承、轴与齿轮、轴与盘连接的问题,根据各自的不同情况可能有不同的连接形式。
但大多数连接形式中存在过盈配合,也就是涉及到接触问题的分析。
这里以某转子中轴和盘的连接为例,分析轴和盘的配合应力以及将轴从盘中拔处时盘轴连接处的应力情况。
本实例的轴为一等直径空心轴,盘为等厚度圆盘,其结构及尺寸如图1所示。
由于模型和载荷都是轴对称的,可以用轴对称方法进行分析。
这里为了后处理时观察结果更直观,我们采用整个模型的四分之一进行建模分析,最后将其进行扩展,来观察整个结构的变形及应力分布、变化情况。
盘和轴用同一种材料,其性质如下:
弹性模量
EX=2.1e5
泊松比
NUXY=0.3
接触摩擦系数
MU=0.2
图1.盘轴结构图
二、ANSYS分析步骤
2.1设置分析标题
首先我们通过完成如下工作来建立本实例的有限元模型,需要完成的工作有:
指定分析标题,定义单元类型,定义材料性能,建立结构几何模型、进行网格划分等。
根据本实例的结构特点,我们将首先建立代表盘和轴的两个1/4圆环面,然后对其进行网格划分,得到有限元模型。
本实例为进行如图1所示的盘轴结构的接触分析,属于非线性结构分析范畴。
跟前面实例一样,为了在后面进行菜单方式操作时的方便,需要在开始分析时就指定本实例分析范畴为“Structural”。
本实例的标题可以命名为:
“AnalysisofaAxisContactingaholeinaDisc”,具体的操作过程如下:
1.选取菜单路径UtilityMenu|File|ChangeJobname,将弹出ChangeJobname(修改文件名)对话框,如图2所示。
在Enternewjobname(输入新文件名)文本框中输入文字“CH20”,为本分析实例的数据库文件名。
并单击Newloganderrorfiles(新的日志和错误文件)单选框,使其变为“Yes”,为本实例的分析过程创建新的日志。
单击按钮关闭对话框,完成文件名的修改。
图2.修改文件名对话框
2.选取菜单路径UtilityMenu|File|ChangeTitle,将弹出ChangeTitle(修改标题)对话框,如图3所示。
在Enternewtitle(输入新标题)文本框中输入文字“AnalysisofaAxisContactingaholeinaDisc”,为本分析实例的标题名。
单击
按钮,完成对标题名的指定。
图3.修改标题对话框
3.选取菜单路径UtilityMenu|Plot|Replot,指定的标题“AnalysisofaAxisContactingaholeinaDisc”将显示在图形窗口的左下角(图略)。
4.选取菜单路径MainMenu|Preference,将弹出PreferenceofGUIFiltering(菜单过滤参数选择)对话框。
单击Structual(结构)选项使之被选中,以将菜单设置为与结构分析相关的选项。
单击
按钮关闭,完成分析范畴的指定。
2.2定义单元类型
本实例分析的问题中涉及到大变形,故选用Solid185单元类型来建立本实例的模型。
本接触问题属于面面接触,目标面和接触面都是柔性的,将使用接触单元TARGE170和CONTA174来模拟接触面。
接触单元在分析过程中使用接触向导时可以自动添加,这里就不再添加。
下面为定义单元类型的具体操作过程。
1.选取菜单路径MainMenu|Preprocessor|ElementType|Add/Edit/Delete,将弹出ElementTypes(单元类型)对话框。
单击对话框中的
按钮,将弹出LibraryofElementTypes(单元类型库)对话框,如图4所示。
图4.单元类型库对话框
2.在单元类型库对话框中,靠近左边的列表中,单击“StructuralSolid”一次,使其高亮度显示,指定添加的单元类型为结构实体单元。
然后,在靠近右边的列表中,单击“Brick8node185”一次,选定单元类型Solid185为第一类单元。
单击对话框中的
按钮,关闭LibraryofElementTypes对话框。
然后,ElementTypes(单元类型)对话框会重新显示,且在对话框中的单元列表框中列出了定义的单元类型1:
Solid185,如图5所示。
图5.定义的单元类型
3.单击ElementTypes(单元类型定义)对话框中的
按钮,关闭对话框,完成单元类型的定义。
2.3定义材料特性
本例中盘轴使用同一种材料类型,所以只需定义一种材料就可以了。
下面是具体的操作过程。
1.选取菜单路径MainMenu|Preprocessor|MaterialProps|MaterialModels,将弹出DefineMaterialModelBehavior(材料模型定义)对话框,如图6所示。
图6.材料模型定义对话框
2.在对话框的右边选项框中,依次双击Structural||Linear||Elastic||Isotropic,将弹出1号材料的弹性模量EX和泊松比PRXY的定义对话框,如图7所示。
图7.线性各向同性材料定义对话框
3.在线性各向同性材料属性对话框中的EX(弹性模量)文本框中输入“2.1E5”,PRXY(泊松比)文本框中输入0.3。
单击对话框中的按钮关闭对话框。
4.在DefineMaterialModelBehavior(材料模型定义)对话框的左边列表框中将列出定义的材料1的属性。
在对话框选取路径Material|Exit关闭对话框,完成对材料模型的定义。
2.4建立几何模型并分网
下面来建立本实例的轴对称几何模型,并进行合理的分网。
本实例的轴对称模型比较简单,可以用多种途径很方便地建立。
前面我们讲过对于这种旋转体几何模型可以先建立一个形面并对其进行网格划分,然后将这个形面绕其对称轴旋转需要的角度而得到;也可以直接建立圆环,再对其进行实体网格划分而得到整个模型的网格。
这里我们采用后面一种方法,具体操作过程如下。
1.创建四分之一圆盘。
选择菜单路径MainMenu|Preprocessor|Modeling|Create|Volumes|Cylinder|PartialCylinder命令,将打开PartialCyclinde(r创建部分圆环)对话框,如图8所示。
图8。
创建部分圆环对话框
2.在创建部分圆环对话框中的输入圆心坐标为:
WPX=0、WPY=0,内径Rad-1=34,起始角度为Theta-1=0,外径为Rad-2=100,结束角度为Theta-2=90,圆盘厚度为Depth=25,如图20.8所示。
单击
按照设置的值建立圆盘模型,在ANSYS图形显示对话框中将会显示刚创建圆盘(图略)。
3.创建代表轴的圆环。
接着,在创建部分圆环对话框中输入相应的值来建立代表四分之一轴的圆环。
输入的值为:
WPX=0、WPY=0,内径Rad-1=25,起始角度为Theta-1=0,外径为Rad-2=35,结束角度为Theta-2=90,圆盘厚度为Depth=150。
然后单击对话中的
按钮关闭对话框。
在ANSYS图形显示对话框中将会显示刚创建圆盘。
4.调整所建实体的视图。
选择菜单路径UtilityMenu|PlotCtrls|PanZoomRotate,打开Pan-Zoom-Rotate(平移-缩放-旋转)对话框,并单击对话框右上角的
按钮,将建立的模型调到比较合适的位置。
调整后的两个圆环几何体如图9所示。
读者可以单击对话框下部的动态显示选择框,然后通过鼠标右键来调整视图位置。
图9.建立的盘轴几何模型
5.将圆环2(轴)移动到合适的位置。
选择菜单路径MainMenu|Preprocessor|Modeling|Move/Modify|Volumes,将弹出MoveVolume(移动体)拾取对话框,在图形输出窗口中单击代表轴的圆环将其选中,然后单击对话框中的
按钮关闭对话框。
将弹出MoveVolume(移动体)对话框,如图10所示。
图10.移动体对话框
6.在对话框中的Z-offsetinActiveCS(Z轴方向平移量)文本框中输入“-10”,将该体沿Z轴负方向平移10mm,然后单击
按钮关闭此对话框。
7.对创建的体进行网格划分。
选取菜单路径MainMenu|Preprocessor|Meshing|MeshTool打开分网工具对话框,如图11所示。
8.对端面上的线进行分网控制。
在网格划分工具对话框(MeshTool)中的尺寸控制(SizeControls)区中,单击Lines(线单元)的
按钮,将弹出ElementSizesonPickedLines(选定线的单元尺寸定义)拾取对话框。
在图形输出窗口中单击轴某个端面上的两条圆弧线,然后单击拾取对话框中的
按钮,将弹出ElementSizesonPickedLines(选定线的单元尺寸)定义对话框,如图12所示。
图11.分网格工具对话框
图12.对线进行网格控制
9.在线单元尺寸控制对话框中的No.ofelementdivisions(单元划分个数)文本框中输入“15”,指定轴沿周向划分15个单元。
单击对话框中的
按钮对设置进行确认。
10.对轴的网格进行控制。
重复上面的步骤8、9的过程,在图形输出窗口中选则同一端面上的径向边,将它们划分为2份。
再选择轴的高度上的某条线,将其划分20份。
11.完成对轴的网格划分。
单击分网工具对话框中的Mesh下拉框中的“Volume”,指定分网对象为体。
再单击Shape(分网形状)控制区的“Hex/Wedge”单选按钮,指定形状为六面体。
单击其下面的“Sweep”单选按钮,指定分网方式为扫掠。
再单击对话框中的按钮,将弹出VolumeSweeping(扫掠体)拾取对话框,单击图形显示窗口中的轴,将其选中,单击拾取对话框中的OK按钮,对轴进行网格划分。
12.对盘进行网格划分。
重复步骤8~11,将盘周向划分10份,径向划分8份,轴向划分3份,同样用扫掠的方式对其进行网格划分,最后单击分网工具对话框中的Close按钮关闭对话框。
至此,完成了盘轴结构的有限元建模的全部工作,通过上面的工作建立的有限元模型如图13所示。
图13.盘轴结构的有限元模型
2.5创建接触对
由于轴和盘在连接时是过盈配合,轴的外表面和盘心的表面之间将构成面面接触对。
ANSYS6.1的接触对生成向导可以使用户非常方便地生成分析需要的接触对。
下面我们将利用接触对生成向导来生成本实例需要的接触对。
在生成接触对的同时,ANSYS程序将自动给接触对分配实常数号。
1.打开接触管理器。
选取菜单路径MainMenu|Preprocessor|Modeling|Create|ContactPair,打开ANSYS6.1中的ContactManager(接触管理器),如图14所示。
图14.接触管理器
1.单击接触管理器中的工具条上的最左边按钮,将弹出AddContactPair(添加接触对)对话框,如图15所示。
图15.添加接触对对话框
3.单击对话框中的“Areas”单选按钮,指定接触目标表面为面,然后单击按钮来选择具体的目标面。
将弹出SelectAreaforTarget(选择目标面)拾取对话框。
在图形输出窗口中单击圆盘的盘心面将其选定,然后单击拾取对话框中的按钮将其关闭。
这时,AddContactPair(添加接触对)对话框中的按钮将被激活,单击按钮进入下一步,将弹出选中接触面的对话框。
4.单击对话框中的“Areas”单选按钮,指定接触表面为面,然后单击按钮,来选择具体的接触面。
将弹出SelectAreaforContact(选择目标面)拾取对话框。
在图形输出窗口中单击轴的外环面将其选定,然后单击拾取对话框中的按钮将其关闭。
这时,AddContactPair(添加接触对)对话框中的按钮将被激活,单击按钮进入下一步,对接触对属性进行设置。
其对话框如图16所示。
图16.设置接触对属性
5.在对话框中,单击“Includeinitialpenetration”选择框将其选中,使分析中包括初始渗透。
单击MaterialID(材料代号)下拉框中的“1”,指定接触材料属性为定义的一号材
料。
并在CoefficientofFriction(摩擦系数)文本框中输入“0.2”,指定摩擦系数为0.2。
单击按钮,来对接触问题的其它选项进行设置,弹出的对话框如图17所示。
图17.基本设置选项
6.在对话框中的NormalPenaltyStiffness(正则处罚刚度)文本框中输入“0.1”,指定接触刚度的处罚系数为0.1。
然后单击对话框上部的Friction(摩擦)标签,打开对摩擦选项设置的选项卡,如图18所示。
图18.设置摩擦选项
7.单击选项卡中的Stiffnessmatrix(刚度矩阵)下拉框中的“Unsymmetric”选项,将其选中,指定本实例的接触刚度为非对称矩阵。
其余的设置保持缺省,单击按钮关闭对话框,完成对接触选项的设置。
8.单击图16所示的按钮,ANSYS程序将根据前面的设置来创建接触对。
然后弹出如图19所示的对话框。
图19.完成接触对的创建
9.查看图19所示对话框中的信息,然后单击按钮关闭对话框。
在ANSYS的接触管理器的接触对列表框中,将列出刚定义的接触对,其实常数为3。
关闭接触管理器。
在图形输出窗口中显示的接触对,如图20所示。
图20.定义的接触对
10.单击ANSYSToolbar(工具条)上的按钮,保存数据库文件。
11.选取菜单路径UtilityMenu|Plot|Areas,对建立的模型进行从新显示,在ANSYS图形输出窗口中对所建模型的面进行显示。
至此,完成了本实例有限元模型的全部工作,下面将进行加载求解工作。
三、加载并求解
因为本实例是对盘轴结构中,将轴从盘心拔出的过程进行非线性分析。
根据条件知道在这个过程中盘的外缘节点的自由度该全部约束。
又由于建模时为了节约费用,根据模型结构的特点,只建立了四分之一的模型,所以需要在分析时定义对称边条,来模拟真实的情况。
本实例的分析过程由两个载荷步组成,第一个载荷步为过盈分析,求解盘轴过盈安装时的应力情况。
第二个载荷步为将轴从盘心拔出时的接触分析,分析在这个过程中盘心面和轴的外表面之间的接触应力。
它们都属于大变形问题,属于非线性问题。
在分析时需要定义一些非线性选项来帮助问题的收敛。
下面进行本实例的加载和求解操作。
3.1定义对称边条和位移约束
前面我们讲过,为了节约建模费用,我们根据模型的结果特点只建立了四分之一模型,这样在分析时就需要定义轴对称边条来模拟实际的情况。
另外根据问题的描述知道,圆盘外缘的节点的所有自由度应该被约束。
下面为具体的操作过程。
1.定义轴对称边条。
选取菜单路径MainMenu|Solution|DefineLoads|Apply|Structural|Displacement|SymmetryB.C.|OnAreas,将弹出ApplySYMMonAreas(在面上施加轴对称边条)拾取对话框。
在图形输出窗口中单击选取盘和轴的四个径向截面,如图21所示。
然后,单击拾取对话框中的OK按钮关闭对话框,对它们施加轴对称边条。
2.定义盘外缘的位移约束。
选取菜单路径MainMenu|Solution|DefineLoads|Apply|Structural|Displacement|OnAreas,将弹出ApplyU,ROTonAreas(在面上施加位移约束)拾取对话框。
在图形输出窗口中单击盘的外缘面,然后单击拾取对话框中的按钮关闭拾取对话框,将弹出ApplyU,ROTonAreas(在面上施加位移约束)对话框,如图22所示。
图22.定义位移约束
2.在对话框中,单击DOFstobeconstrainted(约束自由度)列表框中的“AllDOF”,使其高亮度显示,将其选中。
其余设置保持缺省值(缺省的位移值为0),单击按钮关闭拾取对话框,完成对位移约束的定义。
施加的轴对称约束和位移约束如图23所示。
图23.定义的轴对称约束和位移约束
3.2定义并求解第一个载荷步
对于本实例,第一个载荷步是盘轴连接时的过盈配合分析,它属于结构静力分析的大变形分析。
这里需要进行的工作是指定分析类型,载荷步选项,以及输出文件控制。
1.指定分析类型。
选择菜单路径MainMenu|AnalysisType|Solution|NewAnalysis,将会弹出NewAnalysis(新分析)对话框(图略),单击对话框中的“Static”单选按钮,指定分析类型为静力分析。
2.设定分析选项。
选择菜单路径MainMenu|Solution|AnalysisType|Sol'nControls,将弹出求解控制对话框。
3.单击Basic(基本)选项卡左边的AnalysisOptions(分析选项)下拉框中的“LargeDisplacementStatic”将其选中,使分析中考虑大变形影响。
然后在TimeControl(时间控制区)中设定载荷步结束时间Timeatendofloadstep=100,并关掉自动时间步(Automatictimestepping为“Off”)。
其余设置保持缺省,然后单击对话框中的按钮关闭对话框。
方法对比:
通过下面方法同样可以对时间步进行设置。
单击求解控制(SolutionControls)对话框中的AdvancedNL(高级非线性卡)右下部的按钮,打开隐藏的其余求解选项菜单。
然后通过菜单路径:
MainMenu|Solution|LoadStepOpts|Time/Frequenc|Time-TimeStep来设置。
另外,在进行结构静力分析时ANSYS缺省的加载方式为斜坡(Ramped)加载方式。
3.进行第一载荷步求解。
选择菜单路径MainMenu|Solution|Solve|CurrentLS,将弹出STATUSCommand(求解设置状态)窗口(见图25),和SolveCurrentLoadStep(求解当前载荷步)对话框(见图26)。
图25.求解设置状态对话框
4.仔细阅读求解状态窗口中的信息,检查所有的选项和设置是否都正确。
如果正确,单击图26所示对话框框中的按钮,对当前载荷步进行求解。
如果有不对的地方,单击对话框中的按钮,取消这次求解,对正确的设置进行修改并从新求解。
5.ANSYS程序会在求解时弹出不同的状态对话框显示正在进行的求解内容。
经过一定的运算求解之后,会弹出求解完成提示对话框,单击按钮,完成第一载荷步的求解。
6.求解完成之后ANSYS图形显示窗口中显示的是求解过程的迭代曲线,选择菜单路径UtilityMenu|Replot,可以对窗口中的内容重新显示成盘轴结果的有限元模型。
3.3定义并求解第二载荷步
本实例中,第二载荷步是求解将轴从盘心拔出过程中轴和盘的接触应力情况。
在这个载荷步中需要定义轴的位移值(沿轴向移动的距离)。
同时,需要定义多个载荷子步来进行迭代求解。
下面是定义并求解第二载荷步的具体操作过程。
1.设定分析选项和输出控制选项。
选择菜单路径MainMenu|Solution|AnalysisType|Sol'nControls,将弹出SolutionControls(求解控制)对话框,如图27所示。
图27.求解控制对话框
2.在求解控制对话框中的Basic(基本)选项卡中,将AnalysisOptions(分析选项)设置为“LargeDisplacementStatic”,指定为大变形分析。
将TimeControl(时间控制)区的Timeatendofloadstep(载荷步结束时间)设置为“250”,并打开Automatictimestepping(自动时间步选项),设置为“On”)。
并且将Numberofsubsteps(载荷子步数)文本框设置为“150”,Maxno.ofsubsteps(最大子步数)为“10000”,Minno.ofsubsteps(最小载荷子步数)为“10”。
3.单击对话框右边WriteItemstoResultsFile(结果输出项)设置区下面的Frequency(输出频率)下拉框中的“Writeeverysubstep”,将每个载荷子步结果都输出到结果文件中。
然后单击按钮关闭对话框。
4.施加位移载荷(将轴沿轴向平移40mm,拔出盘孔)。
选取菜单路径UtilityMenu|Select|Entities,将弹出SelectEntities(选择实体)对话框。
单击对话框中最上面下拉框中的“Node”,指定选择对象为节点。
接着在下面的下拉框中单击“ByLocation”,指定选择方式为根据坐标值来选取。
单击“Zcoordinates”单选按钮,在下面的文本框中输入“140”,指定选取Z坐标值为140的所有节点。
单击按钮,然后单击按钮关闭对话框,完成选取。
5.选取菜单路径MainMenu|Solution|DefineLoads|Apply|Structural|Displacement|OnNodes,将弹出施加节点位移载荷拾取对话框,单击对话框中的按钮,将弹出ApplyU,ROTonNOdes(施加节点位移载荷)对话框,如图20.29所示。
单击对话框中约束自由度列表框中的“UZ”,使其高亮度显示。
然后在Displacementvalue(位移值)文本框中输入“40”,其余设置保持缺省,单击按钮关闭对话框,完成位移载荷的施加。
6.选取菜单路径UtilityMenu|Select|Everythings,选取所有的有限元元素。
7.由于大变形影响和加载方式在第一载荷步中都已经设置,这里不需要再重新定义。
下面直接求解第二载荷步。
选择菜单路径MainMenu|Solution|Solve|CurrentLS,将弹出STATUSCommand(求解设置状态)窗口和求解SolveCurrentLoadStep(当前载荷步)对话框(图略)。
8.仔细阅读求解状态窗口中的信息,检查所有的选项和设置是否都正确,然后关闭它。
如果信息正确,单击求解当前载荷步对话框框中的按钮,对当前载荷步进行求解。
如果有不对的地方,单