带孔薄板有限元分析

带孔薄板有限元分析
带孔薄板有限元分析

带孔薄板有限元分析

基于ABAQUS/CAE软件对带孔薄板的分析1.模型的建立

1.1模型简为平面应力问题,建立二维平面可变形壳体模型;

1.2模型的草图,单位为 m尺寸如图所示

1.3二维平面模型图

2.赋予材料属性

2.1 sheet 材料为homogeneous isotropic均匀各向同性材料),E=210Gpa, u=0.3,板厚为0.02m;

22上述建立的材料赋给模型,如下图所示

3.约束边界及施加载荷

3.1约束:左端固支;载荷:右端施加均布拉力,大小为60Mpa如下图所示;

4.划分网格

4.1分割模型,孔周围存在应力集中,附近应网格较密,板边界可以网格布置稀疏,有利于电脑求解速度加快,节约时间,而且精度也相应提高;模型分割如下图所示;

4.2 布置边界种子,在圆环边界附近布置种子密集,可以使网格便密,如下图所示;

4.3网格划分单元的选取

CPS8R: An 8-node biquadratic plane stress quadrilateral,

4.4网格划分如下图所示

reduced integration. (8节点四边形二次单元,采用减缩积

分)

7

5.计算求解

5.1建立工作

6.后处理,查看分析结果

6.1 模型Mises应力图,材料力学强度理论中的第四强度理论,机变能密度理论,单位为Pa,如下图所示应力云图,最大值和最小值一再图中标出

最大应力:2.386+eO8Pa,约为238.6Mpa,在孔的最下端处;

最小应力:3.984+eO6Pa,大约为 3.984Mpa;

6.2通过后处理,拉伸平面0.02m 后模型的三维应力云图;Sj Mimes (Avg: 75%) +?.3?^+D5 +?.l?Qfl+Ot : +l,?95fl+08 +1.7^9fl+0ti +1.64Mfl+D0 +l.-q06fl+fl6 +1怎1站+骼 +l 1(Sl?e i +£l6 +6■总畐站+D A +41S0B?+DT +i :s M3c+D7 1—1- +J.WIe+De Mtxs 斗孔3駅 Elim! SHEETS-Ut htodff ; 64 Min* +3:.理?十 D 百 日OH%; SHEET-1 ^3 陶则:150 Maw: +ZJOGE+DOD Anflyfif th* *h ?*t with 4 心阳? hoi* in th* rniddi* OCB: shee^odt! Abdqus/Stjndard 6.9-3 Thu Jun 2B 20:30:21 GMT-it-aBiao M1Z 5t?P; S^p-1 Encrement 1: Step Tirrie Primary Var: £ Mtsof Defcrrned V”; U 。讲射怕試审“ ScA 刊曲匕 *1-,95S?+Q2 6.2 模型位移云图(位移单位:m ; 最大位移:5.106e-04m ,约为 51.06mm 最小位移:0; g t :K :s :K ^^ le-L e-i l s - :r "*-?-.t ■■*.■: ?=?1: 4-- -■ ■- "- ' -T a L/Ainr L r d l r a u r ;" 一 Ln*

1 J h f l * f i d - n k -1 I s f l * F l

2 y -

7总结

通过ABAQUS/CA 有限元软件简单的对一个各向同性均匀的板, 进行了应力 分析,通过分析我们可以知道,板的最大应力及变形值,这个值可以供我们参 考,可以采取一些措施来使模型的变形及应力在材料所允许的安全范围之类, 这种分析方法在工程实际问题中具有一定的实际意义! Mat; 4 J U --fl-ii 4 MfideE

l .1?J MH!中口4QRWQ rp-dt- ^HErr-iz

r-

-z'

? .■」_■: c

r 5 - L r 吕 r E H 7

1 J h f l * f i

J - n k -1 I s f l * h l -2 u ■L-L ‘

ANSYS 有限元分析 平面薄板

《有限元基础教程》作业二:平面薄板的有限元分析 班级:机自101202班 姓名:韩晓峰 学号:201012030210 一.问题描述: P P h1mm R1mm 10m m 10mm 条件:上图所示为一个承受拉伸的正方形板,长度和宽度均为10mm ,厚度为h 为1mm ,中心圆的半径R 为1mm 。已知材料属性为弹性模量E=1MPa ,泊松比为0.3,拉伸的均布载荷 q =1N/mm 2。根据平板结构的对称性,只需分析其中的二分之一即可,简化模型如上右图所 示。 二.求解过程: 1 进入ANSYS 程序 →ANSYS 10.0→ANSYS Product Launcher →File management →input job name: ZY2→Run 2设置计算类型 ANSYS Main Menu: Preferences →select Structural → OK 3选择单元类型 ANSYS Main Menu: Preprocessor →Element Type →Add/Edit/Delete →Add →select Solid Quad 4node 42 →OK → Options… →select K3: Plane Strs w/thk →OK →Close 4定义材料参数 ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX: 1e6, PRXY:0.3 → OK 5定义实常数以及确定平面问题的厚度 A NSYS Main Menu: Preprocessor →Real Constants …→Add/Edit/Delete →Add →Type 1→OK →Real Constant Set No.1,THK:1→OK →Close 6生成几何模型 a 生成平面方板 ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Rectangle →By 2 Corners →WP X:0,WP Y:0,Width:5,Height:5→OK b 生成圆孔平面 ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Circle →Solid Circle →WPX=0,WPY=0,RADIUS=1→OK b 生成带孔板 ANSYS Main Menu: Preprocessor →Modeling →Operate →Booleans → Subtract →Areas →点击area1→OK →点击area2→OK 7 网格划分 A NSYS Main Menu: Preprocessor →Meshing →Mesh Tool →(Size Controls) Global: Set →SIZE: 0.5 →OK →iMesh →Pick All → Close

有限元分析薄板挠度(附C程序)

1问题描述 某周边简支非均匀的矩形(或圆形)板在均布载荷作用下挠度过大。结合实际,提出集中改进设计方案,并进行对比分析。 2.问题分析 不均匀板有两种主要的情况,结构不均匀和材料不均匀,结构不均匀是指板的厚度不是常量,材料不均匀体现在板的弹性模量和泊松比是变化的。另外,有的板可以是以上两种情况的混合情形。 不均匀板与均匀板的有限元问题有哪些差别呢?下面从均匀板问题推导出非均匀板有限元问题的解决方法。 2.1应力应变 先以结构不均匀板为例来讨论。假设一矩形板长为2,宽为2,厚度沿x ,y 不均匀,由一函数()h ,h x y =描述,但仍然符合薄板假设。对于均匀板,显然h 是一个常数。设挠度为()=x,y ωω,则板内应变向量可以表示为 {}2222211==z 1 2x x y y xy xy x z y x y ρεεεω εγγ?????????????????????????? ?=-???????????????????????? ?????????? 应力应变关系为 {}1p z D σρ????=? ????? 弯矩扭矩矩阵 {}{}()() h ,2h ,2 x y x y M zdz σ-=? 这里就体现出不均匀板和均匀板的区别了。积分完毕后,可以得到 {}[]1M D ρ?? =????

其中薄板的弯曲系数矩阵 []()()()3 21 ,101210 1/2Eh x y D μ μμμ?? ??=??-??-?? 是关于薄板总体坐标的函数,所以对各个分单元都是不同的。 各单元的弯曲系数矩阵可以采用单元中心处的代替。那么就可以得出一系列的弯曲系数矩阵[]D e i 。如果单元划分得足够细,是可以代替真实解的。 2.2单元分析 可以将板分为边长为0.25的矩形小单元,每一个单元都是一样的。对于任何一个单元的节点,都有3项独立的位移 {}i i i xi i yi i w w w y w x δθθ???? ? ???????????? ==???? ??????????? ??????- ???????? 位移模式 ()223123456722333 89101112,w x y x y x xy y x x y xy y x y xy αααααααααααα=+++++++ ++++ 形状函数矩阵是一个112?的行向量 ()[],k l m n N x y N N N N =???? 其中 222222222 2 22222211128111111i i i i i i i i i i i i i x x y y x x y y x y N a b a b a b x x y y y y x x y y x x y x a b b a b a ? ??????=++++--?? ? ????????????? ? ????????????++--++-? ??? ? ? ????????????????? (),,,i k l m n = 单元刚度矩阵 [][][][]1212e e T S k B D B dxdy ?=? 很明显,积分式中包含了弹性系数矩阵,而不同单元的弹性系数矩阵是不同的,所以, 即便单元划分相同,得到的单元刚度矩阵也不同。对于均匀板,相同形式的单元,刚度矩阵

有限元模态分析报告实例

ANSYS模态分析实例 5.2ANSYS建模 该课题研究的弹性联轴器造型如下图5.2: 在ANSYS中建立模型,先通过建立如5.2所式二分之一的剖面图,通过绕中轴线旋转建立模拟模型如下图5.3

5.3单元选择和网格划分 由于模型是三给实体模型,故考虑选择三维单元,模型中没有圆弧结构,用六面体单元划分网格不会产生不规则或者畸变的单元,使分析不能进行下去,所以采用六面体单元。经比较分析,决定采用六面体八结点单元SOLID185,用自由划分的方式划分模型实体。课题主要研究对象是联轴器中橡胶元件,在自由划分的时候,中间件2网格选择最小的网格,smart size设置为1,两端铁圈的smart size设置为6,网格划分后模型如图5.4。 5.4边界约束 建立柱坐标系R-θ-Z,如5-5所示,R为径间,Z为轴向

选择联轴器两个铁圈的端面,对其面上的节点进行坐标变换,变换到如图5.5所示的柱坐标系,约束节点R,Z方向的自由度,即节点只能绕Z轴线转 5.5联轴器模态分析 模态分析用于确定设计中的结构或者机器部件振动特性(固有频率和振型),也是瞬态变动力学分析和谐响应分析和谱分析的起点。 在模态分析中要注意:ANSYS模态分析是线性分析,任何非线性因素都会被忽略。因此在设置中间件2的材料属性时,选用elastic材料。 5.5.1联轴器材料的设置 材料参数设置如下表5-1: 表5.1材料参数设置 表5.1材料参数设置 铁圈1 中间件2 铁圈3 泊松比0.3 0.4997 0.3 弹性模量Mpa 2E5 1.274E3 2E5 密度kg/m 7900 1000 7900 5.5.2联轴器振动特性的有限元计算结果及说明 求解方法选择Damped方法,频率计算结果如表5-2,振型结果为图5.6: 表5.2固有频率 SET TEME/FREQ LOAO STEP SUBSTEP CUMULATIVE 1 40.199 1 1 1 1 73.63 2 1 2 2 3 132.42 1 3 3 4 197.34 1 4 4

有限元分析报告样本

《有限元分析》报告基本要求: 1. 以个人为单位完成有限元分析计算,并将计算结果上交;(不允许出现相同的分析模型,如相 同两人均为不及格) 2. 以个人为单位撰写计算分析报告; 3. 按下列模板格式完成分析报告; 4. 计算结果要求提交电子版,报告要求提交电子版和纸质版。(以上文字在报告中可删除) 《有限元分析》报告 一、问题描述 (要求:应结合图对问题进行详细描述,同时应清楚阐述所研究问题的受力状况和约束情况。图应清楚、明晰,且有必要的尺寸数据。) 一个平面刚架右端固定,在左端施加一个y 方向的-3000N 的力P1,中间施加一个Y 方向的-1000N 的力P2,试以静力来分析,求解各接点的位移。已知组成刚架的各梁除梁长外,其余的几何特性相同。 横截面积:A=0.0072 m2 横截高度:H=0.42m 惯性矩:I=0.0021028m4x 弹性模量: E=2.06x10n/ m2/ 泊松比:u=0.3 二、数学模型 (要求:针对问题描述给出相应的数学模型,应包含示意图,示意图中应有必要的尺寸数据;如进行了简化等处理,此处还应给出文字说明。) (此图仅为例题)

三、有限元建模(具体步骤以自己实际分析过程为主,需截图操作过程) 用ANSYS 分析平面刚架 1.设定分析模块 选择菜单路径:MainMenu—preference 弹出“PRreferences for GUI Filtering”对话框,如图示,在对话框中选取:Structural”,单击[OK]按钮,完成选择。 2.选择单元类型并定义单元的实常数 (1)新建单元类型并定 (2)定义单元的实常数在”Real Constants for BEAM3”对话框的AREA中输入“0。0072”在IZZ 中输入“0。0002108”,在HEIGHT中输入“0.42”。其他的3个常数不定义。单击[OK]按 钮,完成选择 3.定义材料属性 在”Define Material Model Behavier”对话框的”Material Models Available”中,依次双击“Structural→Linear→Elastic→Isotropic”如图

带孔平板拉伸作业

带孔平板有限元分析 本文采用有限元法,对带圆孔的矩形平板进行了弹塑性受力分析,分析了圆孔处的应力集中现象,为其设计和应用提供了参考依据。 1. 研究问题概述 本文研究带圆孔矩形平板在轴对称拉力作用下的平面应力问题。平板开孔的应力问题是弹塑性力学平面中的一个经典的问题,也是实际工程中常见的问题。平板长200mm ,宽50mm ,厚8mm ,具体几何参数及受力见图1。 图1 平板几何参数及受力 2.弹性力学方法解答 由弹性力学知识知,在距圆孔圆心()r ρρ>处的径向正应力、环向正应力、切应力分别为: 222222 1c o s 211322p r p r r ρσψρρρ?????? =-+-- ? ????????? 22221cos 21322p r p r ?σψρρ????=+-+ ? ???? ? 2222sin 21132p r r ρψψρ ττψρρ???? ==--+ ?????? ? 沿着y 轴,90ψ=。,环向正应力为: 242413122r r p ?σρρ?? =++ ???

max 3q ?σ=由上表可知: ()max = 3K q ψ σ=故应力集中因子: 可见孔边最大应力比无孔时提高了3倍,应力集中系数k=3,如图2所示。 图2 孔边应力集中 3.有限元分析 3.1模型建立 图3 有限元模型 3.2边界条件和载荷 为避免在计算时平板产生移动引发计算问题,必须对试件的外部边界条件进行限定。对平板左侧进行铰接约束,示意图如下

图4 平板约束示意图 由于我们只关注孔附近的应力分布情况,根据圣维南原理,载荷的具体分布只影响载荷作用区附近的应力分布。故我们用均布力代替集中力施加在平板右侧的作用面上,其大小为225P MPa ,为负值。 图5 平板载荷示意图 3.3材料 平板的弹性模量为200GPa ,泊松比为0.3。其塑性的应力应变参数见下图 图6 塑性应力应变参数 3.4有限元网格划分 网格划分是非常重要的过程,它会对计算速度、精度、可靠性产生重要影响。网格划分主要包括两方面:尺寸、单元类型。

薄板圆孔的ANSYS分析

板中圆孔的应力集中 问题:如图所示为一个承受单向拉伸的无限大板,在其中心位置有一个小圆孔。材料属性为弹性模量E=211Pa,泊松比为0.3,拉伸载荷q=1000Pa,平板厚度t=0.1. 1、定义工作名和工作标题 (1)定义工作文件名:在弹出的Change Jobname对话框中输入Plate。选择New log and error files复选框,单击OK按钮。 (2)定义工作标题:在弹出的的Change Title对话框中输入The analysis of plate stress with small circle,单击OK按钮。 (3)重新显示:执行replot命令。 2、定义单元类型和材料属性 (1)选择单元类型:在弹出的Element Type中,单击Add按钮,弹出所示对话框,选择Structural Solid和Quad 8node 82选项,单击OK,然后 单击close。 (2)设置材料属性:在弹出的define material models behavior窗口中,双击structural/linear/elastic/isotropic选项,弹出linear isotropic material properties for material number 1对话框,EX和PRXY分别输入2e11和 0.3,单击OK,执行exit命令。 (3)保存数据:单击SAVE_DB按钮。 3、创建几何模型 (1)生成一个矩形面:执行相应操作弹出create rectangle by dimensions对话框,输入数据,单击OK,显示一个矩形。 (2)生成一个小圆孔:执行创建圆的操作弹出对话框,输入数据,单击OK,生成一个圆。 (3)执行面相减操作:执行Booleans/Subtract/Areas命令,生成结果如图示。(4)保存几何模型:单击SAVE_DB按钮。 4、生成有限元网格(自由网格划分) (1)设置网格的尺寸大小:执行size cntrlsl-global-size命令,弹出对话框,在element edge lenge文本框中输入0.5,单击OK. (2)采用自由网格划分:执行mesh/areas/free命令,生成网格模型如图示。

学习模态分析要掌握的的知识

模态分析中的几个基本概念 一、模态定义:物体按照某一阶固有频率振动时,物体上各个点偏离平衡位置的位移是满足一定的比例关系的,可以用一个向量表示。 模态分析一般是在振动领域应用,每个物体都具有自己的固有频率,在外力的激励作用下,物体会表现出不同的振动特性: 一阶模态是外力的激励频率与物体固有频率相等的时候出现的,此时物体的振动形态叫做一阶振型或主振型; 二阶模态是外力的激励频率是物体固有频率的两倍时候出现,此时的振动外形叫做二阶振型,以依次类推。 一般来讲,外界激励的频率非常复杂,物体在这种复杂的外界激励下的振动反应是各阶振型的复合。 二、模态分析:模态是结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得,这样一个计算或试验分析过程称为模态分析。 有限元中模态分析的本质是求矩阵的特征值问题,所以“阶数”就是指特征值的个数。将特征值从小到大排列就是阶次。 实际的分析对象是无限维的,所以其模态具有无穷阶。但是对于运动起主导作用的只是前面的几阶模态,所以计算时根据需要计算前几阶的。 一个物体有很多个固有振动频率(理论上无穷多个),按照从小到大顺序,第一个就叫第一阶固有频率,依次类推。所以模态的阶数就是对应的固有频率的阶数。 三、振型是指体系的一种固有的特性。它与固有频率相对应,即为对应固有频率体系自身振动的形态。每一阶固有频率都对应一种振型。振型与体系实际的振动形态不一定相同。振型对应于频率而言,一个固有频率对应于一个振型。按照频率从低到高的排列,来说第一振型,第二振型等等。此处的振型就是指在该固有频率下结构的振动形态,频率越高则振动周期越小。在实验中,我们就是通过用一定的频率对结构进行激振,观测相应点的位移状况,当观测点的位移达到最大时,此时频率即为固有频率。实际结构的振动形态并不是一个规则的形状,而是各阶振型相叠加的结果。 四、模态扩展是为了是结果在后处理器中观察而设置的,原因如下: 求解器的输出内容主要是固有频率,固有频率被写到输出文件Jobname.OUT及振型文件Jobnmae.MODE中,输出内容中也可以包含缩减

带孔平板的线性静力分析

带孔平板的线性静力分析 本示例将对一个给定的带孔平板几何模型创建有限元模型、施加边界条件、进行有限元分析并在HyperView中观察受载平板的变形和应力结果。 本示例包括以下步骤: ?在HyperMesh中建立有限元模型 ?施加载荷和边界条件 ?求解 ?观察结果 1.在HyperMesh中建立有限元模型 (1)载入OptiStruct用户界面并打开模型文件 1)启动HyperMesh。 2)在User Profile对话框中选择OptiStruct,点击OK。 这就加载了OptiStruct用户界面,它包括OptiStruct模板、宏菜单等。简化了与OptiStruct 使用相关的HyperMesh功能。 User Profiles…可以从下拉式菜单中的Preferences中进入。 3)在工具条选择按钮。 弹出Open file…窗口。 4)选择plate_hole.hm文件,模型位于/tutorials/os/。 5)点击Open。 plate_hole.hm的数据被载入当前的HyperMesh中,替代了原有的数据。数据仅包含几何。 注意此时plate_hole.hm的路径显示在file:文本框中。 6)点击Return。 (2)定义材料属性、单元属性卡片及component 1)点击定义材料。 2)在面板左边选择create子面板。 3)点击name =并输入steel。 4)点击card image =并从弹出菜单中选择MAT1 5)点击create/edit。 弹出MAT1 的卡片信息。 如果括号中的量下面没有值,表示其处于关闭状态。要改变该状态,点击括号中的量,

模态分析有限元仿真分析学习心得

有限元仿真分析学习心得 1 有限元分析方法原理 有限元分析(FEA,Finite Element Analysis)利用数学近似的方法对真实物理系统(几何和载荷工况)进行模拟。还利用简单而又相互作用的元素,即单元,就可以用有限数量的未知量去逼近无限未知量的真实系统。 有限元法是随着电子计算机发展而迅速发展起来的一种工程力学问题的数值求解方法。20世纪50年代初,它首先应用于连续体力学领域—飞机结构静、动态特性分析之中,用以求得结构的变形、应力、固有频率以及阵型。由于其方法的有效性,迅速被推广应用于机械结构分析中。随着电子计算机的发展,有限元法从固体力学领域扩展到流体力学、传热学、电磁学、生物工程学、声学等。 随着计算机科学与应用技术的发展,有限元理论日益完善,随之涌现了一大批通用和专业的有限元计算软件。其中,通用有限元软件以ANSYS,MSC公司旗下系列软件为杰出代表,专业软件以ABAQUS、LS-DYNA、Fluent、ADAMS 为代表。 ANSYS作为最著名通用和有效的商用有限元软件之一,集机构、传热、流体、电磁、碰撞爆破分析于一体,具有强大的前后处理及计算分析能力,能够进行多场耦合,结构-热、流体-结构、电-磁场的耦合处理求解等。 有限元分析一般由以下基本步骤组成: ①建立求解域,并将之离散化成有限个单元,即将问题分解成单元和节点; ②假定描述单元物理属性的形(shape)函数,即用一个近似的连续函数描述每个单元的解; ③建立单元刚度方程; ④组装单元,构造总刚度矩阵; ⑤应用边界条件和初值条件,施加载荷; ⑥求解线性或者非线性微分方程组得到节点值,如不同节点的位移; ⑦通过后处理获得最大应力、应变等信息。 结构的离散化是有限元的基础。所谓离散化就是将分析的结构分割成为有限

带孔板的建模及有限元分析Word版

基于SolidWorks带孔板的建模及有限元分析 李军 摘要:利用SolidWorks对带孔矩形板进行虚拟建模,通过赋予板材材质、载荷后进行网格划分,进而进行有限元分析,得出其应力、应变和位移的分布图,并对结果进行分析研究对板材安全性的影响。 关键词:SolidWorks;带孔板;建模;有限元分析 0 SolidWorks简介 Solidworks是一款优秀的三维设计软件,具有十分强大的零件设计功能及装配模块,同时也拥有丰富的后置处理模块。由于其功能强大,新手上手快,应用领域广,所以成为了主流的三维造型软件。经过17年的发展,在全球已经拥有30多万的客户,最新版本为SolidWorks 2011版。在中国SolidWorks在计算机辅助设计、计算机辅助工程、计算机辅助制造、计算机辅助工艺、数据管理等方面为企业提供了强大的动力,使企业在管理、设计和制造方面有了很大的提升。 1 带孔板的模型建立 矩形板材的尺寸为300*180*10mm,孔位于中心,直径为50mm,模型如图1。 图1 带孔矩形板模型 2前置处理 2.1在Command Manager中点击SIMULATION选项,建立新算例,名称默认,确认。 2.2赋予板材材料属性 材料为AISI304,材料属性如表1

表1 材料的属性 模型参考属性零部件 名称:AISI 304 模型类型:线性弹性同向性 默认失败准则:最大von Mises 应力屈服强度: 2.06807e+008 N/m^2 张力强度: 5.17017e+008 N/m^2 弹性模量: 1.9e+011 N/m^2 泊松比:0.29 质量密度:8000 kg/m^3 抗剪模量:7.5e+010 N/m^2 热扩张系数: 1.8e-005 /Kelvin SolidBody 1(凸台-拉伸1)(aisi304带孔矩形钢板静力分析) 曲线数据:N/A 2.3网格生成 在SIMULATION选项中选择“运行”中的“生成网格”,使用默认网格划分。网格 信息如表2,网格信息细节如表3,网格划分后的模型如图2。 表2 网格信息 网格类型实体网格 所用网格器: 基于曲率的网格 雅可比点 4 点 最大单元大小7.44196 mm 最小单元大小7.44196 mm 网格品质高 表3 网格信息细节 节点总数23523 单元总数13612 最大高宽比例 3.9347 单元(%),其高宽比例< 3 99.7 单元(%),其高宽比例> 10 0 扭曲单元(雅可比)的% 0 完成网格的时间(时;分;秒): 00:00:03 计算机名: PC-201009062016

带孔平板的应力集中分析

有限元方法 Finite Element Method ——基于ANSYS的有限元建模与分析 姓名吴威 学号20100142 班级10级土木茅以升班2班 西南交通大学 2014年4月

综合练习——带孔平板的应力分布及应力集中系数的计算一、问题重述 计算带孔平板的应力分布及应力集中系数。 二、模型的建立与计算 在ANSYS中建立模型,材料的设置属性如下 分析类型为结构(structural),材料为线弹性(Linear Elastic),各向同性(Isotropic)。弹性模量、泊松比的设定均按照题目要求设定,以N、cm为标准单位,实常数设置中设板厚为1。

采用solid 4 node 42板单元,Element Behavior设置为Plane strs w/thk。 建立模型时先建立完整模型,分别用单元尺度为5cm左右的粗网格和单元尺度为2cm左右的细网格计算。 然后取四分之一模型计算比较精度,为了使粗细网格单元数与完整模型接近,四分之一模型分别用单元尺度为2.5cm左右的粗网格和单元尺度为1cm左右的细网格计算。 (1) 完整模型的计算 ①粗网格

单元网格的划分及约束荷载的施加如图(单元尺度为5cm) 约束施加时在模型左侧边界所有节点上只施加x方向的约束,即令U X=0,在左下角节点上施加x、y两个方向的约束,即U X=0、U Y=0。荷载施加在右侧边界上,大小为100。 对模型进行分析求解得到: 节点应力云图(最大值222.112)

单元应力云图(最大值256.408) 可看出在孔周围有应力集中现象,其余地方应力分布较为均匀,孔上部出现最大应力。 ②细网格 单元网格的划分及约束荷载的施加如图(单元尺度为2cm)

ANSYS有限元分析与实体建模

第五章实体建模 5.1实体建模操作概述 用直接生成的方法构造复杂的有限元模型费时费力,使用实体建模的方法就是要减轻这部分工作量。我们先简要地讨论一下使用实体建模和网格划分操作的功能是怎样加速有限元分析的建模过 程。 自下向上地模造有限元模型:定义有限元模型顶点的关键点是实体模型中最低级的图元。在构造实体模型时,首先定义关键点,再利用这些关键点定义较高级的实体图元(即线、面和体)。这就是所谓的自下向上的建模方法。一定要牢记的是自下向上构造的有限元模型是在当前激活的坐标系内 定义的。 图5-1自下向上构造模型 自上向下构造有限元模型:ANSYS程序允许通过汇集线、面、体等几何体素的方法构造模型。当生成一种体素时,ANSYS程序会自动生成所有从属于该体素的较低级图元。这种一开始就从较高级的实体图元构造模型的方法就是所谓的自上向下的建模方法。用户可以根据需要自由地组合自下向上和自上向下的建模技术。注意几何体素是在工作平面内创建的,而自下向上的建模技术是在激活的坐标系上定义的。如果用户混合使用这两种技术,那么应该考虑使用CSYS,WP或CSYS,4命令强迫坐标 系跟随工作平面变化。 图5-2自上向下构造模型(几何体素) 注意:建议不要在环坐标系中进行实体建模操作,因为会生成用户不想要的面或体。

运用布尔运算:可以使用求交、相减或其它的布尔运算雕塑实体模型。通过布尔运算用户可直接用较高级的图元生成复杂的形体。布尔运算对于通过自下向上或自上向下方法生成的图元均有效。 图5-3使用布尔运算生成复杂形体。 拖拉或旋转:布尔运算尽管很方便,但一般需耗费较多的计算时间。故在构造模型时,如果用拖拉或旋转的方法建模,往往可以节省计算时间,提高效率。 图5-4拖拉一个面生成一个体〔VDRAG〕 移动和拷贝实体模型图元:一个复杂的面或体在模型中重复出现时仅需要构造一次。之后可以移动、旋转或拷贝到所需的地方。用户会发现在方便之处生成几何体素再将其移动到所需之处,这样 往往比直接改变工作平面生成所需体素更方便。 图5-5拷贝一个面 网格划分:实体建模的最终目的是为了划分网格以生成节点和单元。在完成了实体建模和建立了单元属性,网格划分控制之后,ANSYS程序可以轻松地生成有限元网格。考虑到要满足特定的要求,用户可以请求映射网格划分生成全部都是四边形、三角形或块单元。

ANSYS_有限元分析_平面薄板

: P P h 1mm R1mm 10m m 10mm 条件:上图所示为一个承受拉伸的正方形板,长度和宽度均为10mm ,厚度为h 为1mm ,中心圆的半径R 为1mm 。已知材料属性为弹性模量E=1MPa ,泊松比为0.3,拉伸的均布载荷q = 1N/mm 2。根据平板结构的对称性,只需分析其中的二分之一即可,简化模型如上右图所示。 二.求解过程: 1 进入ANSYS 程序 →ANSYS12.0→ANSYS Product Launcher →File management →input job name: ZY2→Run 2设置计算类型 ANSYS Main Menu: Preferences →select Structural → OK 3选择单元类型 ANSYS Main Menu: Preprocessor →Element Type →Add/Edit/Delete →Add →select Solid Quad 4node 42 →OK → Options… →select K3: Plane Strs w/thk →OK →Close 4定义材料参数 ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX: 1e6, PRXY:0.3 → OK 5定义实常数以及确定平面问题的厚度 ANSYS Main Menu: Preprocessor →Real Constants …→Add/Edit/Delete →Add →Type 1→OK →Real Constant Set No.1,THK:1→OK →Close 6生成几何模型 a 生成平面方板 ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Rectangle →By 2 Corners →WP X:0,WP Y:0,Width:5,Height:5→OK b 生成圆孔平面 ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Circle →Solid Circle →WPX=0,WPY=0,RADIUS=1→OK b 生成带孔板 ANSYS Main Menu: Preprocessor →Modeling →Operate →Booleans → Subtract →Areas →点击area1→OK →点击area2→OK 7 网格划分 ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool →(Size Controls) Global: Set →SIZE: 0.5 →OK →iMesh →Pick All → Close 8 模型施加约束 a 分别给左边施加x 和y 方向的约束 ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →

开孔薄板有限元分析报告

开孔薄板有限元分析报告 一、有限元分析的目的 通过对两种模型(一个上边开口的和另一个上下两边开口的模型)的静力分析,比较与其对应的理论解的不同,了解有限元仿真软件与理论计算存在的,进一步熟悉workbench求解有限元问题的一般步骤。 二、实体建模(两个模型) 建立如下所示的模型,其中,边长300mm,宽80mm,厚5mm,边缘为半径是10mm的半孔。 上边开口的实体模型(模型A) 上下两边开口的模型(模型B)

模型A 模型B

模型采用的单元类型 模型A: 1 386 2 LID186 (20 Node Quadratic Hexahedron) 2 3858 SOLID186 (20 Node Quadratic Wedge) 3 112 CONTA17 4 (Quadratic Quadrilateral Contact) 4 112 TARGE170 (Quadratic Quadrilateral Target) 5 64 SURF154 (3D Quadratic Quadrilateral) 模型B: 1 3856 SOLID186 (20 Node Quadratic Hexahedron) 2 3866 SOLID186 (20 Node Quadratic Wedge) 3 100 CONTA17 4 (Quadratic Quadrilateral Contact) 4 100 TARGE170 (Quadratic Quadrilateral Target)

5 64 SURF154 (3D Quadratic Quadrilateral) 2.载荷与约束的施加方法(绘图表示并说明); 两模型施加的载荷与约束相同 约束:单击static structural,选择长方体的左侧面,鼠标右键选择“insert>fixed support” 载荷:选择长方体的左侧面,鼠标右键选择“insert>force”,大小为50N。 四、计算结果;(变形图,应力等色线图,约束反力列表等) (1)模型A的求解结果: 总的变形位移图

根据MARC的含圆孔正方形薄板四周受力性能的有限元分析

基于MARC的含圆孔正方形薄板四周受

学院: 班级: 学号: 姓名: 标题:针对含圆孔的正方形板四周受力性能的有限元分析 摘要:采用通用的有限元程序MARC研究含圆孔的正方形板四周受力问题。 在工件工作时,小孔的边缘会产生应力集中的现象,极端情况下甚至 会发生破坏,导致失效。通过对该模型的分析,计算出其最大应力、 最大位移及所发生的位置,得出其承载能力和变形特征,使该力学模 型更好服务于建造等工程方面。 关键词:圆孔、正方形板、受均布力、最大应力、最大位移、位置、四分之一

Title: hole for a square plate with four weeks of the force Finite Element Analysis Abstract: In view of daily life, building structure, mechanical steel structure of the existence of multi-shaped plate with a circular hole is the mechanical model, its bearing capacity and design studies and calculations of concern. In this paper, general finite element program MARC square hole of the plate four weeks with the force the issue. Through analysis of the model to calculate the maximum stress, maximum displacement and the location of occurrence, reached its carrying capacity and deformation characteristics. So that the mechanical model to better serve the construction and other projects. Keywords: round hole, square plate, force, maximum stress, maximum displacement, position, deformation characteristics,horizontal direction, vertical direction, a quarter 正文 1.引言:

有限元分析

隔板对悬臂梁力学性能影响的静力学分析 (byTYH 机自) 摘要:本文基于现代设计技术课程,结合课上所学到的有限元分析技术及理论,运用ansys workbench软件对模型进行静力分析,获得采用不同类型隔板的空心悬臂梁受力后的变形情况,分析其力学性能,验证以前学到的理论知识。 正文: 一.模型 悬臂梁模型一。如图1所示,其基本尺寸为:400mm×100mm×100mm,壁厚为10mm,其中一端固定,另一端为自由状态。为了便于在自由端施加作用力,在自由端增加一个尺寸为:100mm×20mm×5mm的凸台。 图1.悬臂梁模型一 悬臂梁模型二在模型一的基础上添加纵向隔板,如图2所示。 图2.悬臂梁模型二 悬臂梁模型三在模型一的基础上添加斜向隔板隔板,如图3所示。 图3.悬臂梁模型三 悬臂梁模型四在模型一的基础上添加横向隔板隔板,如图4所示。 图4.悬臂梁模型四 为了更易于分析,以上四个模型先在3维绘图软件solidworks中绘制出来,在分析时依次导入使用。 二.有限元分析

启动Ansys Workbench进入工作界面,要做的分析类型为静态结构分析,因此双击toolbox中的在工具箱中的Analysis System→Static StStatic新建一个项目。 项目建好后,首先需要编辑材料参数。所用材料为45号钢,查相关资料可知45号钢的密度为7890 kg/m^-3,杨氏模量为2.09E+11,泊松比为0.269。 双击项目框中的Engineering Data项,进入材料参数设置界面,新建材料并命名45,选中Density和IsotropicElastidty选项,然后输入相应参数,如图5所示。材料设置好后退回workbench主界面。 图5.编辑材料参数 导入模型,双击项目框中的Geometry,进入建模界面。由于模型已经提前建好,因此这里只需导入即可,如图6所示。完成之后退回workbench主界面。 图6.导入模型 分析预处理。双击项目框中的Model,进入操作界面。由于软件默认材料为结构钢,首先需要定义模型材料,将材料选为45号钢,如图7。 图7.定义材料 划分网格,这里我将使用智能网格划来划分网格。选中project中的mesh,在details of mesh中设置网格参数,右键选择“Generate Mesh”即可完成网格划分。网格划分完成后如图8所示。

有限元薄板弯曲问题分析

第4章 弹性薄板弯曲问题的有限元法 薄板弯曲问题在理论上和应用上都具有重要意义,并有专门著作加以论述(如杨耀乾《平板理论》)。象其它弹性力学问题一样,用微分方程、差分法等经典方法所能求解的薄板问题很有限,一般只能解决等厚、小孔口、支承情况较简单的单跨板。故工程设计中以往多采用简化、近似、图表等方法来解决板的设计问题。 在板的分析中,常取板的中面为xoy 平面(如图)。平板结构按其厚度t 与短边a 的比值大小而分为: 厚板(Thick plate )和 薄板(Thin plate)两种。 当1<t w 为绝对柔性板。) 4.1 基本理论 一、基本假定 1、略去垂直于中面的法向应力。(0=z σ),即以中面上沿Z 方向的挠度W 代表板的挠度) 2、变形前垂直中面的任意直线,变形后仍保持为垂直中面的直线。(─法向假定 0=zx τ,0=zy τ) 3、板弯曲时,中面不产生应力。(─中面中性层假定) 上述假定常称为薄板小挠度问题假定(or 柯克霍夫假定)。符合上述假定的平板即为刚性板。 二、基本方法

模态分析中的几个基本概念 模态分析中的几个基本概念 分析

模态分析中的几个基本概念 物体按照某一阶固有频率振动时,物体上各个点偏离平衡位置的位移是满足一定的比例关系的,可以用一个向量表示,这个就称之为模态。模态这个概念一般是在振动领域所用,你可以初步的理解为振动状态,我们都知道每个物体都具有自己的固有频率,在外力的激励作用下,物体会表现出不同的振动特性。一阶模态是外力的激励频率与物体固有频率相等的时候出现的,此时物体的振动形态叫做一阶振型或主振型;二阶模态是外力的激励频率是物体固有频率的两倍时候出现,此时的振动外形叫做二阶振型,以依次类推。一般来讲,外界激励的频率非常复杂,物体在这种复杂的外界激励下的振动反应是各阶振型的复合。模态是结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得,这样一个计算或试验分析过程称为模态分析。有限元中模态分析的本质是求矩阵的特征值问题,所以“阶数”就是指特征值的个数。将特征值从小到大排列就是阶次。实际的分析对象是无限维的,所以其模态具有无穷阶。但是对于运动起主导作用的只是前面的几阶模态,所以计算时根据需要计算前几阶的。一个物体有很多个固有振动频率(理论上无穷多个),按照从小到大顺序,第一个就叫第一阶固有频率,依次类推。所以模态的阶数就是对应的固有频率的阶数。振型是指体系的一种固有的特性。它与固有频率相对应,即为对应固有频率体系自身振动的形态。每一阶固有频率都对应一种振型。振型与体系实际的振动形态不一定相同。振型对应于频率而言,一个固有频率对应于一个振型。按照频率从低到高的排列,来说第一振型,第二振型等等。此处的振型就是指在该固有频率下结构的振动形态,频率越高则振动周期越小。在实验中,我们就是通过用一定的频率对结构进行激振,观测相应点的位移状况,当观测点的位移达到最大时,此时频率即为固有频率。实际结构的振动形态并不是一个规则的形状,而是各阶振型相叠加的结果。 固有频率也称为自然频率( natural frequency)。物体做自由振动时,其位移随时间按正弦或余弦规律变化,振动的频率与初始条件无关,而仅与系统的固有特性有关(如质量、形状、材质等),称为固有频率,其对应周期称为固有周期。 物体做自由振动时,其位移随时间按正弦规律变化,又称为简谐振动。简谐振动的振幅及初相位与振动的初始条件有关,振动的周期或频率与初始条件无关,而与系统的固有特性有关,称为固有频率或者固有周期。 物体的频率与它的硬度、质量、外形尺寸有关,当其发生形变时,弹力使其恢复。弹力主要与尺寸和硬度有关,质量影响其加速度。同样外形时,硬度高的频率高,质量大的频率低。一个系统的质量分布,内部的弹性以及其他的力学性质决定 模态扩展是为了是结果在后处理器中观察而设置的,原因如下: 求解器的输出内容主要是固有频率,固有频率被写到输出文件Jobname.OUT 及振型文件Jobnmae.MODE 中,输出内容中也可以包含缩减的振型和参与因子表,这取决于对分析选项和输出控制的设置,由于振型现在还没有被写到数据库或结果文件中,因此不能对结果进行后处理,要进行后处理,必须对模态进行扩展。在模态分析中,我们用“扩展”这个词指将振型写入结果文件。也就是说,扩展模态不仅适用于Reduced 模态提取方法得到的缩减振型,而且也适用与其他模态提取方法得到的完整振型。因此,如果想在后处理器中观察振型,必须先扩展模态。谱分析中的模态合并是因为激励谱是其实是由一系列的激励组合成的一个谱,里面的频率不会是只有一个,而不同的激励频率对于结构产生的结果是不一样的,对于结果的贡献也是不一样的,所以要选择模态组合法对模态进行组合,得到最终的响应结果。

相关文档
最新文档