有限元理论基础

《有限元理论基础》由会员分享，可在线阅读，更多相关《有限元理论基础（7页珍藏版）》请在装配图网上搜索。

1、有限元理论基本21数值模拟技术21.1数值模拟技术简介在工程技术领域中许多力学问题和场问题，实质上就是在一定的边界条件下求解某些微分方程。对于少数简朴问题,人们可以通过建立它们的微分方程与边界约束求出该问题的解析解。但是对于比较复杂的数学方程问题以及不规则的边界条件通过激吻戏法往往难以求解，而需要借助多种数值模拟措施活的相应的工程数值解，这就是所谓的数值模拟技术。在实际工程领域中,用数值模拟技术可以对复杂的工程构造进行受力和响应分析,这样可以在设计或者加工前预知实体构造工作状态下的大概状况。目前在工程实际应用中,常用的数值求解措施有:有限单元法、有限差分法、边界元等但从实用性和使用范畴来说,有

2、限单元法则是随着计算机技术的发展而被广泛应用的一种行之有效的数值计算措施。2.22 有限元法有限元法是一种基于能量原理的数值计算措施,是解决工程实际问题的一种有效的数值计算工具。它是里茨法的另一种表达形式，它可应用里茨法分析的所有弹性理论。限元法是解决持续的构造体离散或有限个单元集合,也就是将持续的求解域离散为一定数量的单元集合体。且每个单元都具有一定的节点,相邻单元通过节点互相持续，同步使用等效节点力替代作用于单元上的力和选定场函数的节点值作为基本未知量。并在每一单元中假设一种近似插值函数以表达单元中场函数的分布规律:进而运用力学中的某些变分原理去建立用以求解节点未知量的有限元法方程，从而将

3、一种持续域中的无限自由度问题化为离散域中的有限自由度问题。求解后，可运用解出的节点值和设定的插值函数拟定整个单元集体上的场函数。有限元求解问题中的单元分析：式中:单元节点作用力。 :单元刚度矩阵。 ：单元节点位移。通过单元分析拟定单元刚度矩阵,建立单元节点作用力和单元为伊关系。有限元求解问题时建立的构造整体平衡方程：式中：P构造整体等效点力载荷 K构造总体刚度矩阵 U构造节点位移阵列单元内力的计算：式中:弹性矩阵 P应变矩阵整个构造的有限元分析就是一句上述方程而进行的具体的有限元求解过程如图结构离散生成有限元网格求解线性代数方程构造总体刚度矩阵构造节点位移阵列引入约束条件计算并输出单元应输出节

4、点位移 2.2有限元法的基本理论2.2有限元法理论在有限元法中,单元的应变位移关系可表达为：式中:t应变向量 u位移向量 应变位移变换矩阵单元的应力应变关系表达为：式中:应力向量 材料有关系数在线弹性材料条件下,D矩阵是一种常量：在非线性弹性材料中,矩阵上市应变t的函数。有限元刚度方程为：式中：P构造总体刚度矩阵 K单元刚度矩阵其中单元刚度矩阵为:式中:V积分域 对于非线性弹性材料而言,D矩阵和单元刚度矩阵K均是应变t和位移u的函数。在小变形问题中，矩阵与位移u没有有关关系。而在大变形问题时,矩阵B和单元刚刚度 矩阵K则均是位移的函数。在轮胎分析中，轮胎由于充气和垂直载荷等作用，轮胎构造会产生

5、较大的变形,轮胎集合构造的这种变形属于集合非线性问题。轮胎构造 自身又是多种材料构成的复合体,其材料属性既有各向同性又有各项异性,这属于材料非线性问题。在轮胎的静态解触、自由滚动和动态接触状态下,轮胎与路面之间的接触，又波及到接触非线性问题提。2.22应力应变理论有限元是里茨法的另一种表达形式，它可应用里茨法分析的所有弹性理论,而应力应变理论则是里茨分析法的弹性理论的基本。因此在有限元分析中，一般使用弹性理论俺就载荷作用下物体中的内力状态和变形规律。1. 应力 物体收到外力的作用时发生变形，这种变形变化了物体内各分子间的间距，在物体内形成了一种内立场。当内力和外力互相平衡时，变形不再继续，物体

6、达到稳定平衡状态。 这种由于物体受外力的作用因其物体的变形，而导致内部各部分之间因相对位置变化而因其的互相作用，这种互相作用成为内力。所谓应力，就是指分布内力系在物体内某一点处的强弱的互相作用。为了研究物体内某一点处的内力，假设用以通过点C的截面mm将物体分开，在这选用涉及点C的一种部分进行研究。如图2.2所示。环绕点取微笑面积,上存在着分布内力系的合力,如图2.2b所示。的大小和方向与点C的位置和的大小密切有关。的比值成为平均应力。 为一矢量,表达在范畴内，单位面积上内力的平均集度,称作平均应力。随着的逐渐变小,的大小和方向都将逐渐变化。当趋近于零时,的大小和方向都将趋近于一定极限P,即C点

7、应力P为应力P是分布内力系在点C处的集度，反映了分布内力系在点C处的强弱成都。对于应力P，威力表征其与无建议的形变或者材料的有关性,一般将应力P分解成垂直于截面的分量正应力和切与截面的分量切应力。2. 应变 应变表达物体收到外力的作用时发生变形的强弱限度。在图2.a中，物体中的点因变形位移到点,为无体统变形时M点的位移。这里假设物体受到约束,没有刚性位移,M点的位移全是由变形引起的,假想在点附近取平面与坐标平面平行的正六面体(当正六面体的边长趋于无限小时为单元体)，设该六面体的棱边边长分别为。变形后其边长和棱边夹角都发生微小变化将单元体投影到xy平面,如图2.3b所示。变形前单元体平行于x轴的

8、边长长度为,变形后，点M和点N分别位移到点和点。的长度为，且比值：表达单元体边长的平均长度变化,成为平均应变。当长度趋近于0时，则的极限为：2.23 大变形的有限元数学描述在轮胎构造分析中,轮胎变形问题属于集合非线性范畴，此时小变形状况下的几何方程不再适合轮胎的构造分析,为此必须重新定义新的平衡方程。这重要是方程和平衡方程不能消除刚性运动的影响,从而无法度量大变形物体的形态。在度量物体的变形时,需要选用一种特定的构造为基准在变形问题中，一般有两种参照构型。设初始时刻=0,时质点坐标为，表达任意时刻质点的位置。当选择初始构型为参照构型时则为自变量。应变体现式为:式中：为格林应变 为克朗内记号。若

9、选择现时构型，即作为自变量，其应变体现式为:式中：为阿尔曼西应变。2.2.4大变形状态下的平衡方程描述在变形状况下，应力是和应变有关的，在外部载荷的作用下，物体发生形变，并在变形后重新达到平衡。因此，用变形后状态即现时构型构成的auhy应力平衡方程更是合适。设物体在现时构型中占据的区域为,它的边界为,其中A由两部分构成,同事设现时构型的单位体积的体积载荷和表面载荷分别为，如图所示。这时在区域的Cauhy应力平衡方程为:应力边界条件：式中：为现时构型边界的外法线方向余弦。结合Cuchy应力张量和Kirchohf应力张量可以将大变形问题的平衡方程表述为:.5大变形增量问题 轮胎的充气和与路面接触分

10、析都与其变形过程有着密切关系,对于这些非线性问题要想获得一种合适解,就必须采用增量法。所谓增量法，即将时间历程离散或时间序列,即在的增量阶段，需要选用一种参照构型，一般使用两种方式：一种是完全拉格朗日方程，即选用某一时刻的构型为参照构型。一般在轮胎此类大变形问题中,一般选用拉格朗日方程求解比较合适。.26接触问题在轮胎的静态接触、自由滚动和动态接触状态下，轮胎与路面之间发生接触行为。在NSYS中，接触属于一种高度的状态非线性行为，它是状态非线性类型中一种特殊而又比较重要的子集部分。接触问题重要体目前（1)接触区域的大小、互相接触物体的相对位置以及接触的具体状态事先都是都是未知的;(2)接触无建

11、议之间不可互相侵入，接触界面间的法向作用只存在压力，切向可以滑移,这些条件共同构成高度非线性的单边性不等式约束。在状态非线性分析中,刚度矩阵K相有关位移矩阵x。 式中:力与位移的关系是非线性的。其关系如图所示接触一般分为刚柔接触和柔柔接触,接触方式一般分为三种形式：面面接触、点面接触和点点接触。 轮胎路面接触属于典型的刚柔面面接触。在刚柔面面接触中,需要把一种个面设立为“目的”面,而把另一种面设立为“接触面”，一般选择刚性路面为“目的面”,而接触面选择柔性的胎面花纹表面。这两个面共同构成“接触对”,使用相似的实常数号，有限元程序则是通过这组相似的实常数号来辨认该接触对。接触分析中，接触区域需要同步满足运动和力学两方面约束条件。运动约束条件是指接触的两个物体在接触边界上法向无相对运动，即无建议之间不存在互相渗入。力学约束是指接触面法向方向上只能存在接触压力，而在接触面的切向方向上会产生大小相等方向相反的摩擦力。两个接触物体在接触过程中也许会浮现三种情形：粘式、滑移和分离。粘式指接触的两个无建议在接触点粘在一处，不存在相对运动。滑移是指接触的两个物体在接触边界上发生切向运动。分离是指两个接触的物体在变形前存在接触，但在变形后发生分离。23本章小结: