车辆结构强度基本理论4-6章.ppt

图4.36 计算成果的整理 图4.37 主应力方向 图4.38 结构内部的应力处理 图4.39 两个相邻单元的应力处理 图4.40 边界应力推算 4.12 平面矩形单元 单元节点位移有8个分量， 其列阵为： 单元节点力也有8个分量，其列阵为： 图4.41 平面矩形单元 (1)单元位移模式 图4.42 单元边界受表面力 5 线性代数方程组的解法 由前面的讲述可知，有限元法分析到最后总是归结为解线性代数方程组Kδ=P，即使是求解振动问题和场问题也都经常 用到。 5.1 高斯消元法 (1)向前消元 (2)向后回代 5.2 三角分解法(乔列斯基分解) (1)基本思想 (2)求解过程 (3)例题 (4)一般公式 5.3 波前法 针对这一情况，人们想出了两种方法：一种是分块消元法，另一种是波前法。 1)波前法的思路 图5.1所示为一平面问题，现以该题为例对高斯消元法进行再分析。 图5.1 受面内集中 力三角形平板 (2)第二轮消元(把k32,k42,k52变为0) 2)波前法的步骤 现以图5.1为例介绍波前法的步骤。 第1步：选择波源点。 第2步：选择消元路径。 第3步：消元运算。 第4步：回代求解 图5.2 波前法消 元过程图 6 等参数单元 6.1 平面4节点四边形等参数单元 图6.1 任意四边形单元 6.1.1 单元位移模式和坐标变换 相应的形函数为： 若用ξi,ηi(i=1,2,3,4)表示单元节点i的局部坐标，则式(6.2)可缩写为： 容易看出，式(6.2)或式(6.2a)满足形函数的基本要求： 仿照式(4.65)可写出插值函数形式的单元位移模式： 或写成： 式中的形函数矩阵为： 单元节点位移列阵为： 用同样的形函数式(6.2)，写出单元内任一点坐标(x，y)的单元几何模式： 以节点1为例，在节点1处，ξ=η=-1，只有N1=1，其余N2=N3=N4=0，故有： 式中，a,b,c,d是2,3两节点的整体坐标x2,y2, (a)单元几何形状 (b)斜边的参数坐标 (c)形函数N1的图形表达 图6.2 四节点双线性单元 x3,y3所确定的常数。消去以上两式中的η，可得23边在整体坐标系中的方程，为如下任意斜直线： 由于点的位移是点的坐标的函数，故可假设任意节点i处的位移ui中，包含了与其节点坐标xi ,yi的如下关系： 代入位移模式(6.3)中，单元内任一点的位移有： 因为等参单元有： 6.1.2 单元的应变、应力及雅可比矩阵 单元应变由几何方程可表示成两个矩阵相乘的形式： 由复合函数的求导法则，有： 将以上两式写成矩阵形式： ②在单元内任一点的三个形函数之和等于1，即： 图4.17 三角形单元的形函数 4.4.3 解答的收敛性 (1)必须包含单元的常量应变 在一般情况下，单元中的应变包含常量应变和变量应变。 证明：把式(4.13)代入几何方程(4.5)中，得： (2)必须包含单元的刚体位移 图4.18 悬臂梁的变形 图4.19 刚体转动 图4.20 单元在公共边界不协调 4.5 单元内的应变和应力 4.5.1 单元应变与单元节点位移的关系 4.5.2 单元应力与单元节点位移的关系 4.6 单元刚度矩阵 4.6.1 单元节点力分量 4.6.2 单元平衡方程 单元平衡方程即单元节点力和单元节点位移间的关系，可利用虚功方程来推导。 图4.21 三角形单元节点力分量 图4.22 三角形单元节点位移分量 则式(4.35)变为： 由于虚位移是任意的，因此，上式两边与{δ*}eT相乘的矩阵应当相等，故有： 4.6.3 单元刚度矩阵的性质 ①单元刚度矩阵中任一个元素都是一个刚度系数。 ②单元刚度矩阵是对称矩阵。 ③单元刚度矩阵是奇异矩阵。 ④单元刚度矩阵只与单元的几何形状、大小、方位及材料的性质有关，而与单元的位置无关，不随单元或坐标轴的平移而 改变。 4.7 载荷向节点的移置 图4.23 三角形单 元载荷移置 图4.24 图4.25 4.8 总刚度矩阵的形成 4.8.1 对总刚度矩阵的初步认识 (1)写出单元刚度矩阵的一般形式 图4.26 平面离散结构 图4.27 单元编号与节点编号 图4.28 三角形单元的局部 (2)将每个单元刚度矩阵下标轮换成整体编号，并写出其平衡方程 (3)写出每个节点的受力平衡方程 4.8.2 结构的总刚度方程一般论述 (1)总体思