《有限元理论基础(共5页).doc》由会员分享，可在线阅读，更多相关《有限元理论基础(共5页).doc（5页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上有限元理论基础2.1 数值模拟技术2.1.1数值模拟技术简介在工程技术领域中许多力学问题和场问题，实质上就是在一定的边界条件下求解一些微分方程。对于少数简单问题，人们可以通过建立它们的微分方程与边界约束求出该问题的解析解。但是对于比较复杂的数学方程问题以及不规则的边界条件通过激吻戏法往往难以求解，而需要借助各种数值模拟方法活的相应的工程数值解，这就是所谓的数值模拟技术。在实际工程领域中，用数值模拟技术可以对复杂的工程结构进行受力和响应分析，这样可以在设计或者加工前预知实体结构工作状态下的大概情况。目前在工程实际应用中，常用的数值求解方法有：有限单元法、有限差分法、边界
2、元等但从实用性和使用范围来说，有限单元法则是随着计算机技术的发展而被广泛应用的一种行之有效的数值计算方法。2.2.2 有限元法有限元法是一种基于能量原理的数值计算方法，是解决工程实际问题的一种有效的数值计算工具。它是里茨法的另一种表示形式，它可应用里茨法分析的所有弹性理论。限元法是处理连续的结构体离散或有限个单元集合，也就是将连续的求解域离散为一定数量的单元集合体。且每个单元都具有一定的节点，相邻单元通过节点相互连续，同时使用等效节点力代替作用于单元上的力和选定场函数的节点值作为基本未知量。并在每一单元中假设一个近似插值函数以表示单元中场函数的分布规律：进而利用力学中的某些变分原理去建立用以求
3、解节点未知量的有限元法方程，从而将一个连续域中的无限自由度问题化为离散域中的有限自由度问题。求解后，可利用解出的节点值和设定的插值函数确定整个单元集体上的场函数。有限元求解问题中的单元分析：式中：单元节点作用力。 :单元刚度矩阵。 ：单元节点位移。通过单元分析确定单元刚度矩阵，建立单元节点作用力和单元为伊关系。有限元求解问题时建立的结构整体平衡方程：式中：P结构整体等效点力载荷 K结构总体刚度矩阵 U结构节点位移阵列单元内力的计算：式中：D弹性矩阵 P应变矩阵整个结构的有限元分析就是一句上述方程而进行的具体的有限元求解过程如图结构离散生成有限元网格求解线性代数方程结构总体刚度矩阵结构节点位移阵
4、列引入约束条件计算并输出单元应输出节点位移 2.2有限元法的基础理论2.2.1 有限元法理论在有限元法中，单元的应变位移关系可表示为：式中：t应变向量 u位移向量 B应变位移变换矩阵单元的应力应变关系表示为：式中：应力向量 D材料相关系数在线弹性材料条件下，D矩阵是一个常量：在非线性弹性材料中，D矩阵上市应变t的函数。有限元刚度方程为：式中：P结构总体刚度矩阵 K单元刚度矩阵其中单元刚度矩阵K为：式中：V积分域 对于非线性弹性材料而言，D矩阵和单元刚度矩阵K均是应变t和位移u的函数。在小变形问题中，矩阵B与位移u没有相关关系。而在大变形问题时，矩阵B和单元刚刚度 矩阵K则均是位移u的函数。在轮
5、胎分析中，轮胎由于充气和垂直载荷等作用，轮胎结构会产生较大的变形，轮胎集合结构的这种变形属于集合非线性问题。轮胎结构 本身又是多种材料构成的复合体，其材料属性既有各向同性又有各项异性，这属于材料非线性问题。在轮胎的静态解触、自由滚动和动态接触状态下，轮胎与路面之间的接触，又涉及到接触非线性问题提。2.2.2应力应变理论有限元是里茨法的另一种表示形式，它可应用里茨法分析的所有弹性理论，而应力应变理论则是里茨分析法的弹性理论的基础。因此在有限元分析中，一般使用弹性理论俺就载荷作用下物体中的内力状态和变形规律。1. 应力 物体收到外力的作用时发生变形，这种变形改变了物体内各分子间的间距，在物体内形成
6、了一个内立场。当内力和外力相互平衡时，变形不再继续，物体达到稳定平衡状态。 这种由于物体受外力的作用因其物体的变形，而导致内部各部分之间因相对位置改变而因其的相互作用，这种相互作用成为内力。所谓应力，就是指分布内力系在物体内某一点处的强弱的相互作用。为了研究物体内某一点C处的内力，假设用以经过点C的截面mm将物体分开，在这选取包含点C的一个部分进行研究。如图2.2所示。围绕点C取微笑面积，上存在着分布内力系的合力，如图2.2b所示。的大小和方向与点C的位置和的大小密切相关。的比值成为平均应力。 为一矢量，表示在范围内，单位面积上内力的平均集度，称作平均应力。随着的逐渐变小，的大小和方向都将逐渐
7、变化。当趋近于零时，的大小和方向都将趋近于一定极限P，即C点应力P为应力P是分布内力系在点C处的集度，反应了分布内力系在点C处的强弱成都。对于应力P，威力表征其与无提议的形变或者材料的相关性，通常将应力P分解成垂直于截面的分量正应力和切与截面的分量切应力。2. 应变 应变表示物体收到外力的作用时发生变形的强弱程度。在图2.2a中，物体中的M点因变形位移到点，为无体统变形时M点的位移。这里假设物体受到约束，没有刚性位移，M点的位移全是由变形引起的，假想在M点附近取平面与坐标平面平行的正六面体（当正六面体的边长趋于无限小时为单元体），设该六面体的棱边边长分别为。变形后其边长和棱边夹角都发生微小变化
8、将单元体投影到xy平面，如图2.3b所示。变形前单元体平行于x轴的边长长度为，变形后，点M和点N分别位移到点和点。的长度为，且比值：表示单元体边长的平均长度变化，成为平均应变。当长度趋近于0时，则的极限为：2.2.3 大变形的有限元数学描述在轮胎结构分析中，轮胎变形问题属于集合非线性范畴，此时小变形情况下的几何方程不再适合轮胎的结构分析，为此必须重新定义新的平衡方程。这主要是方程和平衡方程不能消除刚性运动的影响，从而无法度量大变形物体的形态。在度量物体的变形时，需要选取一个特定的结构为基准在变形问题中，一般有两种参考构型。设初始时刻t=0,时质点坐标为，表示任意时刻t质点的位置。当选择初始构型
9、为参考构型时则为自变量。应变表达式为：式中：为格林应变 为克朗内记号。若选择现时构型，即作为自变量，其应变表达式为：式中：为阿尔曼西应变。2.2.4 大变形状态下的平衡方程描述在变形情况下，应力是和应变相关的，在外部载荷的作用下，物体发生形变，并在变形后重新达到平衡。因此，用变形后状态即现时构型构成的Cauchy应力平衡方程更是合适。设物体在现时构型中占据的区域为V，它的边界为A，其中A由两部分构成，同事设现时构型的单位体积的体积载荷和表面载荷分别为，如图所示。这时在区域V的Cauchy应力平衡方程为：应力边界条件：式中：为现时构型边界的外法线方向余弦。结合Cauchy应力张量和Kirchoh
10、off应力张量可以将大变形问题的平衡方程表述为：2.2.5 大变形增量问题 轮胎的充气和与路面接触分析都与其变形过程有着密切关系，对于这些非线性问题要想获得一个合适解，就必须采用增量法。所谓增量法，即将时间历程离散或时间序列，即在的增量阶段，需要选用一个参考构型，一般使用两种方式：一种是完全拉格朗日方程，即选取某一时刻t的构型为参考构型。通常在轮胎这类大变形问题中，一般选用拉格朗日方程求解比较合适。2.2.6接触问题在轮胎的静态接触、自由滚动和动态接触状态下，轮胎与路面之间发生接触行为。在ANSYS中，接触属于一种高度的状态非线性行为，它是状态非线性类型中一个特殊而又比较重要的子集部分。接触问
11、题主要体现在（1）接触区域的大小、相互接触物体的相对位置以及接触的具体状态事先都是都是未知的；（2）接触无提议之间不可相互侵入，接触界面间的法向作用只存在压力，切向可以滑移，这些条件共同构成高度非线性的单边性不等式约束。在状态非线性分析中，刚度矩阵K相关于位移矩阵x。 式中：力与位移的关系是非线性的。其关系如图所示接触一般分为刚柔接触和柔柔接触，接触方式一般分为三种形式：面面接触、点面接触和点点接触。 轮胎路面接触属于典型的刚柔面面接触。在刚柔面面接触中，需要把一个个面设置为“目标”面，而把另一个面设置为“接触面”，通常选择刚性路面为“目标面”，而接触面选择柔性的胎面花纹表面。这两个面共同组成“接触对”，使用相同的实常数号，有限元程序则是通过这组相同的实常数号来识别该接触对。接触分析中，接触区域需要同时满足运动和力学两方面约束条件。运动约束条件是指接触的两个物体在接触边界上法向无相对运动，即无提议之间不存在相互渗透。力学约束是指接触面法向方向上只能存在接触压力，而在接触面的切向方向上会产生大小相等方向相反的摩擦力。两个接触物体在接触过程中可能会出现三种情形：粘式、滑移和分离。粘式指接触的两个无提议在接触点粘在一处，不存在相对运动。滑移是指接触的两个物体在接触边界上发生切向运动。分离是指两个接触的物体在变形前存在接触，但在 变形后发生分离。23本章小结：专心-专注-专业