清华大学计算固体力学第四次Lagrangian网格.pptVIP

4 编制程序 Jacobian矩阵 考虑平衡方程 定义残差 通过Newton迭代求解和线性化 内力Jacobian－切线刚度 外力Jacobian－载荷刚度 整体Jacobian为二者之差 ABAQUS的UMAT就是求内力的Jacobian矩阵－切线刚度 Jacobian矩阵是时间的函数， 的逆给出为 4 编制程序 4节点四边形单元采用单元坐标表示的形状函数的梯度为 采用空间坐标表示的形状函数的梯度为 出现在分母中 4 编制程序 速度梯度给出为 内部节点力 上式适用于任意的二维等参单元。因为 是关于单元坐标的有理函数，所以不能对上式进行解析积分，通常使用数值积分。 对于4节点四边形单元，2×2 Gauss积分为完全积分。但是在平面应变问题中，对于不可压缩或几乎不可压缩的材料，完全积分会出现单元自锁，因此必须应用局部减缩积分。 对于4节点四边形单元，沿着每条边的位移是线性的。所以，其外部节点力与3节点三角形单元外部节点力是一致的。 出现在分母中，被积函数 4 编制程序 一致质量矩阵 根据Lobatto积分，使积分点与节点重合，或者将单元的整个质量平均分配到4个节点上可以得到一个集中的对角质量矩阵，分配到4个节点上给出 集中质量矩阵 例题 主－从连接线 4 编制程序 连接线经常用于连接采用不同单元尺度的网格部分，因为它们比使用三角形