理论教育 ANSYS13.0入门实战:有限元分析成果

ANSYS13.0入门实战:有限元分析成果

时间:2023-10-25 理论教育 版权反馈
【摘要】:5.网格划分打开网格划分工具面板GUI:Preferences|Meshing|MeshTool设置单元尺寸选择网格划分工具面板Size Controls中的Global|Set,在弹出的对话框中输入SIZE=0.003,单击OK。

ANSYS13.0入门实战:有限元分析成果

1.定义文件名

GUI:Utility Menu|File|Change Jobname

操作后弹出一个对话框,在输入栏中输入“Optimal Design of Disk”,单击OK。

2.定义几何参数

GUI:Utility Menu|Parameters|Scalar Parameters

在弹出对话框中的Selection中输入pi=3.1415,单击Accept完成对圆周率参数的定义,按照此方法继续定义转盘内半径R1=0.12,转盘轮毂厚度W1=0.015,转盘轮毂边缘到轴心的距离R2=R1+W1,转盘轮毂高度H1=0.05,转盘边缘到轴心的距离R4=0.25,转盘轮缘厚度W2=0.01,转盘轮缘到轴心的距离R3=R4-W2,转盘轮缘高度H2=0.02,转盘轮廓线切角SEITA1=90,转盘轮廓线切角SEITA2=90,转盘中间X坐标Xmid=(R3-R2)/2,转盘中间Y坐标Ymid=(H1+H2)/2,定义完毕后单击Close。

3.定义单元和材料常数

(1)定义单元

GUI:Main Menu|Preprocessor|Element Type|Add/Edit/Delete

弹出对话框,单击Add,在弹出对话框的左边选择Solid,然后在右边选择Quad4node182单元,单击OK。

(2)定义单元属性

单击单元类型面板Element Types的Options按钮,在弹出的对话框中设置单元关键字K3为轴对称Axisymmetric,单击OK,单击Close。

(3)定义单元材料常数

GUI:Main Menu|Preprocessor|Material Props|Material Models

弹出对话框,单击Structural|Liner|Elastic|Isotropic,又弹出一个输入材料属性的对话框,输入EX=3E11,PRXY=0.3,单击OK,单击Structural|Density设置材料密度,在弹出的对话框中设置密度(DENS)为7200,选择Material|Exit,退出面板。

4.建立模型

(1)定义矩形

GUI:Main Menu|Preprocessor|Modeling|Create|Areas|Rectangle|By 2 Corners

在弹出的对话框中输入X=R1,Y=0,Width=W1,Height=H1,单击Apply;继续输入X=R3,Y=0,Width=W4,Height=H2,单击OK。

(2)定义关键点

GUI:Main Menu|Preprocessor|Modeling|Create|Keypoints|In Active CS

在弹出的对话框中输入NPT=10,X=R2+Xmid,Y=Ymid,Z=0,单击OK。

(3)激活总体柱坐标系

GUI:Utility Menu|WorkPlane|Change Active CS to|Global Cylindrical

(4)定义样条曲线

GUI:Main Menu|Preprocessor|Modeling|Create|Lines|Splines|With Options|Spline thru KPs

弹出拾取对话框,用鼠标拾取关键点K3、K10和K8,即轮毂右顶点关键点、转盘中间位置关键点和轮缘左顶点关键点,单击OK,在弹出的对话框中输入XV1=1,YV1=90+SEITA1,ZV1=0,XV6=0,YV6=90-SEITA2,ZV6=0,单击OK。

(5)激活总体笛卡儿坐标系

GUI:Utility Menu|WorkPlane|Change Active CS to|Global Cartesian

(6)定义线

GUI:Main Menu|Preprocessor|Modeling|Create|Lines|Lines|In Active Coord

弹出拾取对话框,拾取关键点K2和K5,即轮毂右下点和轮缘左下点,单击OK。

(7)通过线定义面

GUI:Main Menu|Preprocessor|Modeling|Create|Areas|Arbitrary|By Lines

弹出拾取对话框,拾取线L9,L2,L10和L8,单击OK。

5.网格划分

(1)打开网格划分工具面板

GUI:Preferences|Meshing|MeshTool

(2)设置单元尺寸

选择网格划分工具面板Size Controls中的Global|Set,在弹出的对话框中输入SIZE=0.003,单击OK。

(3)进行网格划分

选择网格划分工具面板中的Mesh按钮,弹出拾取对话框,单击Pick All。

6.静力学分析

(1)设置分析类型

1)定义静力学分析

GUI:Main Menu|Solution|Analysis Type|New Analysis

在弹出的对话框中选中Static,含义为进行静力学问题求解。

2)设置分析选项

GUI:Main Menu|Solution|Analysis Type|Analysis Options

在弹出的对话框中设置PETRES为Prestress ON,含义为应力分析,单击OK。

(2)定义边界条件

1)定义对称边界条件

GUI:Main Menu|Solution|Define Loads|Apply|Structural|Displacement|Symmetry B.C.|On Lines

弹出拾取对话框,拾取线L1,L5和L10,即转盘的对称线,单击OK。

2)定义转速

GUI:Main Menu|Solution|Define Loads|Apply|Structural|Inertia|Angular Veloc|Global

在弹出的对话框中输入OMEGY=2*pi*8500/60,因为ANSYS转速的单位是Hz,所以需进行转换。

(3)求解

GUI:Main Menu|Solution|Current LS

弹出求解对话框,单击OK。

(4)进入后处理提取参数

提取参数使用命令流比较方便,如下:

NSORI,S,EQV

*GET,smax,SORT,,MAX !提取模型中最大等效应力并赋予参数SMAX

*GET,NELEM,ELEM,,COUNT,, !提取模型单元的数量并赋予参数NELEM

ETABLE,eseqv,S,EQV !使用单元表定义单元等效应力

*DIM,sarray,ARRAY,NELEM,1,1… !定义维数为NELEM,1,1数组

!求模型中单元等效应力之和(www.daowen.com)

SUM=O

*DO,I,1,NELEM,1

*GET,sarray(I),ELEM,I,ETAB,E SEQV

SUM=sarray(I)+SUM

*ENDDO

!求模型中单元平均应力

MS=SUM/NELEM

!求模型等效应力的标准差并赋予参数SD

SD=O

*DO,I,1,NELEM,1

SD-(sarray(I)-MS)**2+SD

*ENDDO

SD=SQRT(SD/(NELEM-1))

7.模态求解

(1)设置分析类型

1)定义模态分析分析

GUI:Main Menu|Solution|Analysis Type|New Analysis

在弹出的对话框中选中Modal,含义为进行模态问题求解。

2)设置分析选项

GUI:Main Menu|Solution|Analysis Type|Analysis Options

在弹出的对话框中勾选模态求解方法为Block Lanczos,提取阶数(No.of modes to extract)为3,扩展模态阶数(NMODE No.of modes to expand)为3,并勾选预应力选项(PSTRES),单击OK,保持弹出对话框的默认值,单击OK。

(2)定义边界条件

1)删除对称位移约束

GUI:Main Menu|Solution|Define Loads|Delete|All Load Data|All Loads&Opts

2)定义转盘对称线为反对称约束

GUI:Main Menu|Solution|Define Loads|Apply|Structural|Displacement|AntisymmB.C.|OnLines

弹出拾取对话框,拾取线L1,L5和L10,即转盘的对称线,单击OK。

(3)求解

GUI:Main Menu|Solution|Current LS

弹出求解对话框,单击OK。

(4)提取固有频率

提取参数使用命令流比较方便,如下:

*GET,FREQl,MODE,1,FREQ

*GET,FREQ2,MODE,2,FREQ

*GET,FREQ3,MODE,3,FREQ

8.优化设计

(1)整理出命令流并复制到ANSYS的工作目录

本实例设置命令流文件名为OPT.txt。

(2)读入OPT.txt

GUI:Utility Menu|File|Read Input from,在工作目录中选中OPT.txt。

(3)指定分析文件

GUI:Main Menu|Design Opt|Analysis File|Assign

弹出选择文件对话框,单击Browse,选择文件OPT.MAC,单击OK。

(4)定义设计变量

GUI:Main Menu|Design Opt|Design Variables

弹出定义设计变量对话框,单击Add,在弹出的对话框中选择NAME中的SEITA1,同时设置MIN=70,MAX=90,单击Apply;继续选择NAME中的SEITA2,同时设置MIN=70,MAX=90,单击Apply;继续选择NAME中的XMID,同时设置MIN=0.03,MAX=0.08,单击Apply;继续选择NAME中的YMID,同时设置MIN=0.03,MAX=0.08,单击OK,单击Close。

(5)定义状态变量

GUI:Main Menu|Design Opt|State Variables

弹出定义状态变量对话框,单击Add,在弹出的对话框中,选择NAME中的FREQ2,同时设置MIN=2500,单击Apply;继续选择NAME中的SMAX,同时设置MAX=230E6,单击OK,单击Close。

(6)定义目标函数

GUI:Main Menu|Design Opt|Objective

弹出定义目标函数对话框,单击Add,在弹出的对话框中,选择NAME中的SD,单击OK,单击Close。

(7)定义优化方法

GUI:Main Menu|Design Opt|Method/Tool

弹出指定优化方法(Specify Optimization Method)面板,选择随机优化方法(Random Designs),单击OK,在弹出的对话框中设置最大迭代次数(Maximum iterations)为50,单击OK。

(8)查看优化结果

1)查看最优设计。

GUI:Main Menu|Design Opt|Design Sets|List

弹出列举设计序列对话框,选择最优序列(BESTSet),单击OK。

2)查看设计变量与目标函数的曲线

GUI:Main Menu|Design Opt|Design Sets|Graphs/Tables

在弹出的对话框中选择XVAROPT中的XMID,选择Y-variable params中的SD,单击Apply,如图26-3所示;继续选择XVAROPT中的SEITA1,选择Y-variable params中的SD,单击OK,如图26-4所示。

978-7-111-35546-5-Chapter26-5.jpg

图26-3 目标函数与XIMD关系曲线

978-7-111-35546-5-Chapter26-6.jpg

图26-4 目标函数与SEITA1关系曲线

免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。

我要反馈