地球讲义13,14课－计算土力学.ppt

计算土力学 主讲教师：张爱军 §4.4 常应变三角形单元（补充） 单元刚度矩阵形成 有了刚度值，就可以由位移求出力（非应力） 有了弹性模量值，就可以由应变求出应力 整体各个单元刚度矩阵分析 总体刚度矩阵的叠加 整体刚度矩阵的性质 整刚表示整体位移与力的关系 总体刚度矩阵的行数和列数＝总结点数×单元结点自由度数的方阵，其排列的顺序按照结点标号依次从小到大排列，每个结点的按照自由度的顺序排列。上例中总刚度矩阵排列见后，这样排列是为了与单元结点位移向量对应。 结点的多少决定着整刚的大小。 约束处理 就是将整刚中相应的约束结点的约束自由度上的行、列划去，将荷载向量中相应项划去即可。 其实质是：强制使得u＝0，或v＝0；或u＝0，v=0 对于给定位移约束条件，即：u＝c1或v＝c2 是将整体刚度矩阵中相应的主元乘以一个大值A，等效结点向量相应的值设为A×c1或A×c2，利用“大数吃小数”的原理，得到u＝c1或v＝c2。在上例中，若整体结点号为2的水平位移u2＝c2，则有： §4.5 等参元 以上三角形单元分析中可以看出，从位移模式到形函数的构造上面，需要解三元一次方程，若位移模式中取得的项数进一步增加到4～20项，其形函数的构造非常麻烦，而且不能够统一编程，因此需要找更方便的办法 在三角形单元中，对于等效结点荷载向量的求解等就很复杂了，对于其他复杂单元其求解将更加困难。 这就引出来另外一种思路。即：先将实际单元通过坐标转换函数转化成一个母元（即：形状规则、简单的单元），所有位移模式、形函数、单元分析，等效结点荷载向量建立等均在该母元上进行，使得求解简单化，并便于变成。然后将求出的量再转化成实际单元的量，求出最后的结果。这种思路概念清楚、便于理解，并且便于标准化。这就是等参元的思路。 其中：坐标转换函数＝母元的形函数时称为等参元。 等参元的形函数 母元 对于平面问题 母元位移模式与形函数的构造 位移模式的待定参数的个数应该与结点个数一致，才能保证待定参数可以解出来。对于4边形单元，位移模式取为（在局部坐标中分析）： 将母元各个结点的局部坐标值代入位移模式中，得到： 以上公式就是平面4结点等参元的形函数 对于空间8结点 坐标转换函数的建立 坐标转换函数的作用在于将几何形状不规则的真实单元转换成为标准单元（数学上称为映射），等参元是将坐标转换函数与形函数一致，即： 我们分析一下形函数的性质： 坐标转换函数形式： 将η=1（3-4边）代入到坐标转换方程： 化解得到： 将联立方程中ζ消去得到： 雅克比矩阵(Jacobi) [J] 就是真实坐标对母元局部坐标的导数矩阵 对于平面问题： 对于空间问题： Jacobi矩阵地意义 得到： 那么： [J]阵通过以下公式求得： 微元的面积dS和微元的体积dV 在有限元计算中，对于三维问题，等效荷载向量的计算中我们要用到面积积分，而形成单元刚度矩阵时要用到体积积分，这里讲授微元面积与微元体积的计算问题。为形成{R}、{k}准备。 预备知识 设空间上矢量a，b，其表示成为分量的形式为： 该标量乘积为一个标量，表示矢量a在b上的投影长度（也叫模）与矢量b长度的乘积。 a,b两个矢量的向量乘积，即所谓叉乘： 行列式的值： 设η，ζ为平面上的曲线坐标，dη是与曲线η＝1相切的矢量，dζ是与曲线ζ＝1相切的矢量，如图示： 则有： 得到平面问题的微元面积计算公式： 得到： 得到对于空间曲面微元的秘面积为： 微元的体积计算 同样设在空间上的曲线坐标为：ξηζ， ξηζ=1为三个曲面，沿着η＝1 和ζ=1 曲面的交线作dξ，沿ξ＝1 和ζ=1 曲面交线作矢量dη，沿ξ＝1 和 η＝1交线设矢量dζ，如图所示： 微元体的体积就是： Any questions？ 大学课件出品 版权归原作者所有 联系QQ ：910670854 如侵权，请告知，吾即删 更多精品文档请访问我的个人主页 /611696569 附赠人生心语 人生太短，聪明太晚 人生太短，聪明太晚(1) 我们都老得太快?却聪明得太迟 把钱省下来，等待退休后再去享受 结果退休后，因为年纪大，身体差，行动不方便，哪里也去不成。钱存下来等养老，结果孩子长大了，要出国留学，要创业做生意，要花钱娶老婆，自己的退休金都被拗走了。 人生太短，聪明太晚(2) 当自己有足够的能力善待自己时，就立刻去做，老年人有时候是无法做中年人或是青少年人可以做的事，年纪和健康就是一大因素。小孩子从小就告诉他，养你到高中，大学以后就要自立更生，要留学，创业，娶老婆，自己想办法，自己要留多一点钱，不要为了小