有限元的数学求解原理三大步骤预案.ppt

⑷最小势能原理法 * 试函数只要求满足位移边界条件，对函数连续性要求相对较低。 三大步骤： ⑴.结构离散 ⑵.单元分析 ⑶.整体分析 三、有限元法的基本流程 结构离散:就是用假想的线或面将连续物体分割成有限个单元组成的集合体且单元之间仅在节点处连接，单元之间的作用仅由节点传递。（基本要求） 注意的问题 单元：满足一定几何特性和物理特性的最小结构域 节点：单元与单元间的连接点 节点力：单元与单元间通过节点的相互作用力 节点载荷：作用于节点上的外载 ⑴.结构离散 1).选择插值（位移）函数。 插值函数：用以表示单元内物理量变化（如位移或位移场）的近似函数。由于该近似函数常由单元节点物理量值插值构成，故称为插值函数，如单元内物理量为位移，则该函数称为位移函数。 选择位移函数的一般原则 位移函数在单元节点的值应等于节点位移（即单元内部是连续的）； 所选位移函数必须保证有限元的解收敛于真实解。 位移函数一般采用多项式形式，在单元内选适当阶次的多项式可得到与真实解接近的近似解 ⑵.单元分析 构造位移函数：如平面问题位移函数的一般形式为 ⑴.多项式项数越多，则逼近真实位移的精度越高，项数的多少由单元的自由度数决定。 ⑵.多项式选取应由低阶到高阶，尽量选择完全多项式以提高单元精度。 ⑶.选取多项式时，还应使所选取的多项式具有坐标的对称性，即按Pascal（帕斯卡）三角形来选择 2).构造位移函数 位移函数构造方法： ⑴.广义坐标法： ⑵.插值函数法：即将位移函数表示为各个节点位移与已知插值基函数积的和 单元特性分析的基本任务就是建立单元的平衡方程，也称为刚度方程。在选择了单元类型和相应的位移函数后，即可按弹性力学的几何方程、物理方程导出单元应变与应力的表达式，最后利用虚位移原理或最小势能原理或直接法或加权残值法建立单元的平衡方程，即单元节点力与节点位移间的关系。 3).单元特性分析 整体分析的基本任务包括建立整体平衡方程，引入边界条件，完成整体方程求解。 整体平衡方程的建立有多种方法，可基于能量原理（势能变分或虚位移原理）推导，也可基于节点力平衡得到。 在引入边界条件之前，整体平衡方程是奇异的，这意味着整体方程是不可解的。 方程求解包括边界条件引入和数值计算，一旦利用适当的数值方法求出未知的节点位移，则可按前述的应力应变公式计算出各个单元的应变、应力等物理量 ⑶.整体分析 刚度由使其产生单位变形所需的外力值来量度，刚度是指零件在载荷作用下抵抗弹性变形的能力。 单元的刚度矩阵：单元刚度矩阵反应的是单元节点力与单元节点位移的关系； 总刚度矩阵反应的是整体的节点力与节点位移的关系；刚度矩阵将总体坐标下的节点位移与整个结构的总体力联系在一起。 四、刚度矩阵 单元刚度矩阵行数等于位移向量的分量个数，列数等于为位移的列向量的分量个数，由于两者相等所以单刚是个方阵。 结构的总体刚度矩阵即结构的原始刚度矩阵，每1个元素的物理意义就是当其所在列对应的节点位移分量等于单位位移（其余结点位移分量为0）时，其所在行对应的节点力的数值。 表示由于第j个自由度的单位位移dj在第i个自由度需要的力 a. 单元刚度矩阵仅与单元的几何特征和材料性质有关。仅与单元的横截面积A、惯性矩I、单元长度l、单元的弹性模量E有关。 b. 单元刚度矩阵是一个对称阵。在单元刚度矩阵对角线两侧对称位置上的两个元素数值相等，即，根据是反力互等定理。 c. 单元刚度矩阵是一个奇异阵。 d. 单元刚度矩阵可以分块矩阵的形式表示。具有确定的物理意义。 单元的刚度矩阵的性质 整体刚度矩阵 中位于主对角线上的子块 ，称为主子块，其余 为副子块。 a. 中主子块 由结点i的各相关单元的主子块扩展之后叠加求得，即 b. 当结点i、 j为单元e的相关结点时， 中副子块 为该单元e相应的副子块，即 。 c. 当结点i、 j为非相关结点时， 中副子块 为零子块，即 。 d. 仅与各单元的几何特性、材料特性，即A、I、l、E等因素有关。 e. 为对称方阵， f. 为奇异矩阵，其逆矩阵不存在，因为建立整体刚度矩阵时没有考虑结构的边界约束条件。 * 整体刚度矩阵的性质 g. 为稀疏矩阵，整体刚度矩阵中的非零元素分布区域的宽度与结点编号有关，非零元素分布在以对角线为中心的带状区域内，称为带状分布规律，见图a。在包括对角线元素在内的区域中，每行所具有的元素个数叫做把半带宽，以d表示。 最