[原创]材料力学与有限元入门之ANSYS分析
2017-02-17 by:CAE仿真在线 来源:互联网
导读:材料力学刘鸿文版本下册,讲到矩阵位移法,可做为了解有限元法的入门。结构分析中,以位移为基本未知量,借助于平衡方程求解的方法称为位移法,该方法又可以称为杆件结构的有限元法,ANSYS软件的理论基础就是有限元法。本次ANSYS分析要结合矩阵位移法了解有限元法的思想,熟练定义杆单元的单元属性(单元类型、实常数和材料),熟练进行后处理(约束反力,节点位移,单元的轴力和应力)。
一、问题描述
一杆两端固定,长度a =0.1 m,横截面积均为A= 100 mm2,F2= 3 kN,F3=6 kN。材料的弹性模量E= 200 GPa,泊松比为0.3。求杆两端的约束反力,各节点的位移,各杆的轴力和应力。
二、理论计算
位移法的求解过程是:
首先将结构分散成若干个杆件,以后把分散的每一根杆件称为一个单元。把单元的内力用位移表示。
其次,当再把各单元组成原来的结构时,结构各节点的平衡方程便可由位移表示出。从平衡方程中解出位移,于是又可由位移求出各单元的内力。
位移法中如以矩阵为运算工具,就是矩阵位移法。因它易于编制程序,便于利用计算机进行计算,所以受到广泛重视。
本算例采用矩阵位移法的解答过程如下
计算结果如下
节点位移
节点 |
u /mm |
1 |
0 |
2 |
0.02 |
3 |
0.025 |
4 |
0 |
各杆应力
单元号 |
应力/MPa |
1 |
40 |
2 |
10 |
3 |
50 |
约束反力
节点 |
Fx/kN |
1 |
-4 |
4 |
-5 |
参考教材
刘鸿文. 材料力学 II(第5版) [M]. 北京: 高等教育出版社, 2011: 180-184.
三、GUI操作步骤
1.进入ANSYS
程序→ ANSYS → ANSYSProduct Launcher → 改变working directory到指定文件夹→ 在job name输入:Link → Run。job name可默认为file。
2.定义工作文件名及工作标题
(1)定义工作文件名:UtilityMenu > File > Change Jobname → 输入文件名Link→ OK。可不用输入,默认为file。
(2)定义工作标题:UtilityMenu > File > Change Title → 输入Link→ OK。可不用输入。
3.定义单元属性
(1)定义单元类型:MainMenu > Preprocessor > Element Type > Add/Edit/Delete → Add → 在左列表框中选择Link,在右列表框中选择3Dfinit stn 180 → OK → Close。
(2)定义实常数:MainMenu > Preprocessor > Real Constants > Add → Type 1 → OK → RealConstant Set No.:1,AREA:100 → OK → Close。ANSYS无单位,需自己统一,本算例采用N、mm和MPa单位制。
(3)设置材料属性:MainMenu > Preprocessor > Material Props > Material Models → Structural→ Linear → Elastic → Isotropic → EX:2.0E5,PRXY:0.3 → OK。
4.建立模型
(1)定义节点:MainMenu > Preprocessor > Create > Nodes > In Active CS → 依次输入4个节点坐标1 (0,0,0),2 (0.1e3,0,0),3 (0.2e3,0,0),4 (0.3e3,0,0) → OK。
(2)定义单元:
①定义单元:MainMenu > Preprocessor > Create > Elements > Auto Numbered > ThruNodes 按照杆的编号,依次拾取节点。
拾取节点1和2(生成单元1)→ Apply →
拾取节点2和3→ Apply →
拾取节点3和4→ OK。
各杆单元编号,先拾取的节点为I节点,后拾取的节点为J节点。
②打开节点和单元编号:Utility Menu > PlotCtrls > Numbering→ NODE numbers:ON,Elem/Attribnumbering: Element numbers。→ OK。
③显示单元和节点:Utility Menu > Plot > Element,查看节点和单元编号,见下图。
④打开杆单元的单元形状:Utility Menu > PlotCtrls> Style > Size and Shape → [/ESHAPE]:On。
5.施加边界条件并求解
(1)施加约束:MainMenu > Solution > Define loads > Apply > Structural > Displacement> On Nodes →
拾取节点1→ OK → Lab2:AllDOF → Apply →
拾取节点4→ Lab2:AllDOF → OK。
(2)施加载荷:MainMenu > Solution > Define loads > Apply > Structural > Force/Moment> On Nodes →
拾取节点2→ OK → Lab:FX,VALUE:3E3→ Apply →
拾取节点3→ OK → Lab:FX,VALUE:6E3→ OK。
(3)求解前保存模型:UtilityMenu > Files > Save as → 输入Link_Load.db→ OK。
(4)求解:MainMenu > Solution > Solve > Current LS → File > Close → Solve CurrentLoad Step → OK → Solution is done → Close。
(5)保存结果文件:UtilityMenu > Files > Save as → 输入Link_Solve.db→ OK。
6.后处理
(1)查看各节点位移:
①云图显示:MainMenu > General Postproc > Plot Results > Contour Plot > Nodal Solu→ 选择DOFsolution。 → 选择X-componentof displacement(X方向的位移)→OK。云图中的UX表示x方向位移。
②列表显示:MainMenu > General Postproc > List Results > Nodal solu → 选择DOFsolution 下的X-component of displacement → OK。
(2)查看各单元轴力:
①查定义LINK180轴力单元表的方法:打开ANSYSHelp,搜索LINK180。
找到LINK180的单元介绍页面后,找Table180.1: LINK180 Element OutputDefinitions,在Name栏中FORCE便是轴力,Sxx是轴向应力。
找Table180.2: LINK180 Item and Sequence Numbers,表中OutputQuantity Name列找到FORCE,Item为SMISC,编号为1。
②定义轴力单元表:MainMenu > General Postproc > Element Table > Define Table,→ Lab:输入FN→ Item:选择Bysequence num → Comb:选择SMISC,在SMISC后面输入“1” → OK → Close。
③轴力列表显示:Main Menu > General Postproc > Element Table > List Element Table→ 选择FN → OK → 记录各个单元的轴力 → File → Close。
④画轴力图:MainMenu > General Postproc > Plot Results > Contour Plot > Line ElemRes → LabI选择FN,LabJ选择FN → OK。
(3) 查看各单元应力:
①画轴向应力图:MainMenu > General Postproc > Plot Results > Contour Plot > Nodal Solu→ Nodal Solution → Stress → X – Component of stress → OK。
打开杆单元的单元形状:Utility Menu > PlotCtrls> Style > Size and Shape → [/ESHAPE]: On。显示节点应力云图之前需要打开单元形状,否则不能显示该结果。
②轴向应力列表显示:MainMenu>General Postproc>List Results>Element Solution → Element Solution→ Stress → X – Component of stress → OK。
单元应力列表结果
说明:杆单元的一个单元只有一个轴力,云图显示节点应力与列表显示单元应力的结果一致。杆单元的应力是在单元坐标系下显示的,因此应选择X – Component of stress。
(5)列表查看约束反力:MainMenu > General Postproc > List Results > Reaction Solu → All strucforc F → OK。
节点1和节点3的约束反力结果如下:
7.退出ANSYS软件
Utility Menu> File > Exit → Quit-No Save → OK。
相关标签搜索:[原创]材料力学与有限元入门之ANSYS分析 Ansys有限元培训 Ansys workbench培训 ansys视频教程 ansys workbench教程 ansys APDL经典教程 ansys资料下载 ansys技术咨询 ansys基础知识 ansys代做 Fluent、CFX流体分析 HFSS电磁分析 Abaqus培训