CADCAM第6章总结.docxVIP

CAD/CAM第六章学习总结计算机辅助工程分析是涉及有限元分析技术、数值计算技术、产品优化设计方法和工程分析与仿真等在内的一个综合性系统。计算机辅助工程分析的核心技术是工程问题的模型化和数值实现方法 。计算机辅助工程分析的理论基础是有限元法、边界元法及结构优化设计技术等计算力学方法 。一．有限元分析技术1.1 有限元分析方法概述 Finite Element Analysis 是一种用于求解各类工程实际问题的数值计算方法 ，广泛应用于解决力学问题和场问题 。简称为FEA。有限元分析的思路和方法：物体离散化。单元特性分析。单元组集。求解未知节点位移。1. 物体离散化将某个工程结构离散为由各种单元（每种单元可以是一维、二维或三维的情况）组成的计算模型，这一步称作单元剖分。离散后单元与单元之间利用单元的节点相互连结起来。单元节点的设置、性质、数目等应根据问题的性质、描述变形形态的需要和计算精度而定(一般情况，单元划分越细则描述变形情况越精确，即越接近实际变形，但计算量越大)。 2.单元特性分析 选择位移模式：可选择节点位移作为基本未知量；可选择节点力作为基本未知量；亦可取一部分节点力和一部分节点位移作为基本未知量。 分析单元的力学性质：根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式。 计算等效节点力：用等效的节点力来替代所有作用在单元上的力。3.单元组集 利用结构力的平衡条件和边界条件把各个单元按原来的结构重新联接起来，形成整体的有限元方程K是整体结构的刚度矩阵；q是节点位移列阵；F是载荷列阵。4.求解未知节点位移 解有限元方程式得出位移。这里，可以根据方程组的具体特点来选择合适的计算方法。有限元法的基本思想是“一分一合”，“分”是为了进行单元分析，“合”则是为了对整体结构进行综合分析。1.2 有限元分析方法中单元特性的导出方法直接方法。虚功原理法。能量变分原理方法1.直接方法 图中所示是xoy平面中的简支梁弯曲简图，EI为梁的抗弯刚度。现在，以它为例用直接方法建立单元的刚度矩阵。梁在横向外载荷(可以是集中力或力矩或分布载荷等)作用下产生弯曲变形，对于平面弯曲问题，每个点(包括支承点)处的位移有两个，即挠度和倾角；相应地也有两个节点力，即与挠度对应的剪力和与倾角对应的弯矩。我们规定挠度和剪力向上为正，倾角和弯矩逆时针方向为正。为使问题简化，把图示的梁看成是一个单元。当令左支承点为节点i，右支承点为节点j时，则节点位移和节点力可以分别写成Vi、θzi、Vj、θzj和Fyi、Mzi、Fyj、Mzj。也可写成矩阵形式 梁的节点力和节点位移是有联系的。在弹性小位移范围内，这种联系是线性的。它代表了单元的载荷与位移之间(或力与变形之间)的联系，称为单元的有限元方程。式中 [K]称为单元刚阵，它是单元的特性矩阵。单元刚度矩阵中任一元素kij表示j号节点的单位位移对i号节点力的贡献。如 [K]中第1列各元素就分别代表当在i节点处挠度方向产生单位位移i=1时，它们对其他各位移(包括vi)方向上引起的节点力Fyi，Mzi，Fyj，Mzj的贡献。由功的互等定理有kij = kji，所以单元刚度矩阵是对称的。对于梁单元平面弯曲问题，可以计算出各系数kij的数值。假设则梁变形情况如图所示：这时，平面弯曲梁单元的刚度矩阵或单元特性矩阵为： 2.虚功原理法 解题基本步骤：设定位移函数由位移函数求应变根据虎克定律，通过应变求应力由虚功原理求单元的刚度矩阵以平面问题中的三角形单元为例，说明其方法步骤。设三节点三角形单元内的位移函数为：当单元很小时，单元内一点的位移可以通过节点的位移数值来表示。假设单元内位移为x, y的线性函数，即：写成矩阵形式：同理，单元的三个节点i，j，k上的位移也可用它来表示：即其中式中A是三角形面积上式相乘后，可写成也可简写成此式即为单元内某点的位移用节点位移插值表示的多项式。称[N]为形状函数。由弹性力学知根据虎克定律，对于平面问题，有其中对于平面问题为根据虚功原理，当结构受载荷作用处于平衡状态时，在任意给出的节点虚位移下，外力(节点力){F}及内力{σ}所做的虚功之和应等于零，即经推导得平面应力问题三角形单元刚度矩阵为：1.3 有限元法的解题步骤1 单元剖分和插值函数的确定根据构件的几何特性、载荷情况及所要求的变形点，建立由各种单元所组成的计算模型。再按单元的性质和精度要求，写出表示单元内任意点的位移函数u(x,y,z)，v(x,y,z)， w(x,y,z)或{d }=[S(x,y,z)] {a}。上式是用节点位移表示单元体内任意点位移的插值函数式。2 单元特性分析根据位移插值函数，由弹性力学中给出应变和位移关系，计算出应变。由物理关系，得应变与应力的关系式。由虚位移原理，可得单元的有限元方