几何大变形下的两种朗格朗日格式

一、理论推导

在这里插入图片描述在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

二、TL和UL的高性能代码实现

TL和UL虽然是两种基于不同时刻构型的求解方法,但它们理论以及计算结果到最后都是完全一致的,上面的每一项都是相互对应且相等的。以上形式是在学术程序之中经常用到的,但是如果编写成代码需要考虑以下两点:
1.减少大规模变量的存储。TL格式下每个积分点上形函数偏导虽然只需计算一次,但是存储起来进行每步的计算花费实际上代价更高。(我个人更推荐UL格式)
2.尽量减少变量的计算次数。在计算两种格式下整体刚度矩阵的时候需要分别计算材料刚度矩阵和几何刚度矩阵,而形函数偏导有可能计算了两次,这对代码的运行效率也有着一定不利。
这说明以上的理论知识,编写出的代码是有可能效率不高的。在查阅很多相关文献,在Nonlinear solid mechanics—A Continuum Approach for Engineering(chapter 8, pg 396-399)这本书找到了答案,将几何刚度项纳入到材料刚度项中,使得以前的两个刚度项相加,变为只需一个刚度项相加,形函数偏导只需计算一次
在这里插入图片描述
在这里插入图片描述

三、注意事项(06-24)

UL格式下,每一迭代步都进行坐标更新的话,会计算精度更高点?如果是每个载荷步更新位移,计算有以下好处

  • 计算内存更小,因为为了每迭代步更新位移而且为了不能收敛返回到原来值,必须存储以前的坐标值,必定会消耗内存,如果节点数据多的话。
  • 4
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值