《弹性力学与有限元法分析及实例讲解》由会员分享，可在线阅读，更多相关《弹性力学与有限元法分析及实例讲解（15页珍藏版）》请在金锄头文库上搜索。
1、弹性力学与有限元法分析弹性力学是固体力学的一个重要分支，是研究弹性固体在受外力作用、温度改变、边界约束或其他外界因素作用下而发生的应力、形变和位移状态的科学。有限单元法是力学、数学、物理学、计算方法、计算机技术等多种学科综合发展和结合的产物，是随着计算机技术的广泛应用而迅速发展起来的一种数值分析方法。有限元法的基本思想就是化整为零，分散分析，再集零为整。即用结构力学方法求解弹性力学问题，实质是将复杂的连续体划分为有限多个简单的单元体，单元体之间仅仅通过结点相连，实现化无限自由度问题为有限稀有度问题，将连续场函数的（偏）微分方程的求解问题转化为有限个参数的代数方程组的求解问题。有限元方法经过近半
2、个世纪的发展，目前已经成为各种工程问题特别是结构分析问题的标准分析方法，而有限元软件也已成为现代结构设计中不可缺少的工具。有限元软件是有限元理论通向实际工程应用的桥梁，它的应用极大地提高了力学学科解决自然科学和工程实际问题的能力，进一步促进了有限元方法的发展。ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件，广泛用于机械制造、石油化工、航空航天、汽车交通、土木工程、造船、水利等一般工业及科学研究。ANSYS软件的组成：（一）前处理模块该模块为用户提供了一个强大的实体建模及网格划分工具，可以方便的构造有限元模型，软件提高了100种以上的单元类型，用来模拟工程中的各
3、种结构和材料。包括：1. 实体建模：参数化建模，布尔运算及体素库，拖拉、旋转、拷贝、蒙皮、倒角等。2. 自动网格划分，自动进行单元形态、求解精度检查及修正。3. 在集合模型上加载：点加载、分布载荷、体载荷、函数载荷。4. 可扩展的标准梁截面形状库。（二）分析计算模块该模块包括结构分析（可进行线性分析、非线性分析和高度非线性分析）、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析，可模拟多种物理介质的相互作用，具有灵敏度分析及优化分析能力。（三）后处理模块将计算结果以彩色等值线、梯度、矢量、粒子流、立体切片、透明及半透明等图形方式显示出来，也可以用图表、曲线形式显示或输出。由
4、于现在只是对ANSYS工程软件有初步的了解和掌握，所以本次作业仅以（1）结构静力学分析为例，运用ANSYS软件对汽车连杆进行受力分析；（2）三杆桁架的优化设计为例。举例1对连杆进行受力分析结构静力学分析是ANSYS工程软件中应用最广泛的一种分析方法，“结构”不仅包含像桥梁、建筑物等建筑工程结构，也包括像活塞、机械零件和工具等机械零部件一样的船只、航空和机械结构，如船的外壳、航空器、机器的机架等。结构静力学分析用来求解稳态外载荷引起的系统或部件的位移、应变、应力和力,稳态载荷包括外部施加的力和压力、稳态的惯性力(重力和旋转速度、施加位移、温度等)。静力分析很适合求解惯性和阻尼对结构的影响并不显著
5、的问题，如确定结构中的应力集中现象，可以进行线性分析，也可以进行非线性分析，如塑性、蠕变、膨胀、大应变及接触分析。静力分析过程分为三个步骤：图1.静力分析过程三步骤如图所示为汽车连杆的几何模型，连杆的厚度为0.5in，图中所注尺寸均为英寸，在小头孔的内侧90度范围内承受p=1000psi的面载荷作用，试利用有限元分析该连杆的受力状态。连杆的材料属性为：杨氏模量E=30x106psi，泊松比为0.3。连杆的几何模型1.建立计算模型本例是对一个弹性力学平面问题进行结构分析，所以只选择结构分析模块。由于连杆的结构、载荷均对称，因此在分析时只要采用一半进行分析即可。采用由底向上的建模方式，具体的操作过
6、程如下：(1) 生成圆环面和两个矩形，包括创建关键点、直线、样条曲线、倒圆角等；通过对面进行布尔操作产生真实模型。JUIiS烹D0920:42:31AREASftPUArawThestresscacmlatingofLiangan图1.1生成连杆大头圆环面、矩形ANThe9tre9scaculatingofLianganAREASJULS2D09P.PHAm_w2D:44:11图1.2生成连杆小头圆环面，进行布尔运算A?AREASJULS2D0920:Thestresscaculatingofliangan图1.3创建样条曲线ANLIWESLINEm_HJUL52D0920:Thegtress
7、caculatingofliangan图1.4生成连杆轮廓线ANJULS200?Tn*srrsazacmlariniofliann图1.5由边界线住成面(2)划分网格，包括定义单元类型：选择单元为solid,quad8-node,设置单元尺寸为0.2,采取自由网格划分用20节点的S0LID95单元划分网格。ANELBHENTSThmaculatingofLiananJUL5200920:51:22图1.6采用自由划分网格，生成2D网格ELEHENTiiANJULS20092CH舅注nThestresscaculatingofliangan图1.7生成3D网格(3)定义材料属性对平面连杆进行线弹
8、性分析计算，所定义的材料属性有弹性模量E=30X106psi,泊松比为03。2.施加约束和载荷并求解(1)施加约束ELEKENTiiANJULS2D092Q:59:20ThestcesamaculatingofHangan图2.1大孔内表面施加对称约束图22Y=0的所有面上施加对称约束(2)施加载荷ELEEr4TiANJUL52D0921:06:14ThestresscaculatingofLiangan2.3小孔内表面施加面载荷(3)选择PCG求解器，求解运算JULS2D0921:06:42NFORRFORPRESLanaThestceaamaculatingofHangan图3求解运算结果
9、显示图3.查看显示结果ANJUL5200921!人11!57WODRL9OLUTICWSTEP1SUE-1T=lS&JW(AVG1DMX.L33E-n33MH-&21066SMS=1097.SZLOBbAZL.93A643-2人3lZfii1636图4连杆体的受力分布显示由受力图可知，当连杆小头内侧受到面分布载荷时，小头孔内侧靠近大头孔部分受力最大，在一定的条件下，可能最先出现裂纹，也最容易断裂。举例2:三杆桁架的优化设计优化设计是一种寻找确定最优设计方案的技术，人们总希望在一切可能的方案中选择一个最好的，设计方案的任何方面都可以优化，进行优化设计首先要把实际的实际问题用数学表达式加以描述，转
10、化成数学模型，然后根据数学模型的特性，选择某种适当的优化计算方法及其程序，通过计算机求得最优解。优化设计的步骤：(2) (1)生成循环所用的分析文件，建立优化过程中的参数；进入OPT处理器，指定分析文件，指定优化变量；选择优化工具或优化方法；指定优化循环控制方式；进行优化分析；查看设计序列结果。如图所示为一个由三根杆组成的桁架结构，它承受纵向和横向载荷，求该桁架的最小质量25图1桁架结构模型已知桁架特性如下表所示:表1材料特性几何属性弹性模量E=2X50Mpa横截面变化范围0.6x10-30.645m3泊松比0.3基本尺寸B变化范围1025m比重7.8X10kg/m集中载荷F=9X10N许用应
11、力2.76Mpa横截面积A15A2,A3假设结构的初始重量为500t，缺省允差由计算机自动选择，为了便于收敛，一阶方法的优化分析中将目标函数的允差定为1。根据分析问题的性质，选择设计变量，状态变量，目标函数如下：表2设计变量A,人2小3B状态变量6目标函数MinWt该问题的优化数学模型为:Minf(E)X=Xi,X2,X3,X4】=Ai,A2,A3汉10Ai兰0.64516i=1,2,3葩兰B兰250_max(j_2.76Mpaj=1,2,31.1参数化建立几何模型定义参数和材料属性参数初始值：B=25,A仁A2=A3=0.645;材料属性:E=2e11,PRXY=0.3,DENS=7800定
12、义单元类型及属性A1,A2,A3。单元类型：structurallink-2Dspar1;定义实常数：(1) 建立有限元模型hl-ut-EBNTiDKMIIPIOPTofhangjid图1.1生成单元JULb2(口932(4)施加约束和载荷lU-OC-EBANJUL(口20=4S：32OPTofhangjid图1.2施加边界约束N-uC-EBNTiDEJUEL石E(J口92D:45:32OPTofhangjid图1.34节点施加集中载荷1.2求解运算2.1提取并指定状态变量和目标函数WT=481520.422计算单元体积的总和VTOT=61.7333874;计算初始重量显示单元形状和大小：2.
13、2进入优化处理器，指定分析文件指定分析文件为：hangjia_optgw，指定设计变量:A1、A2、A3、B；设置状态变量:sigl、sig2、sig3；设置目标函数：wt,允差为1；指定一阶优化方法：first-order。3.1查看优化结果信息窗口中将显示了迭代序列的结果，包括每个设计变量、状态变量、目标函数的值。所有序列的结果见程序文件。图3.1迭代序列优化运行结果3.2显示目标函数的变化规律ANJU1D21109Zl：57：05LE33S34-UIrerationMuffiDecUD13JilA1工0113:cmMnp0aupisuFTafhaaqiia1ZCQaizu30X1图3.2目标函数WT变化规律=3.3显示基本尺寸的变化规律图3.3B的变化规律3.4显示杆横截面的变化规律ANCTTTMIZATICWJU162皿皿;叩.尢rRIteration为自tOFTaf图3.4Ai的变化规律3.5显示杆中应力的变化规律QPTIMIZATICHSIGI3IG2SIG3SIGZ5IG3(HWIIJZOT1LCUJ2093.3.S34-ftIterationMumDer0PTafhaagjia图3.5勺的变化规律ANJUXo21109Z2!D1：39
