首页/文章/ 详情

有限元基础知识:嵌入式区域(Embedded Region)

2月前浏览944

在Abaqus中Embedded region通常用来模拟不同维度单元之间的相互作用,比如钢筋混凝土,其钢筋用1D单元(桁架、梁)进行模拟,而混凝土则用3D的实体单元进行模拟,而真实中钢筋混凝土一起承受重量,这时如果采用兼容网格的方式进行网格划分,往往非常难以实现,则用嵌入式区域的方式,将3D的混凝土设为被嵌入的区域,而1D的钢筋模型则嵌入其中。


钢筋
混凝土

Embedded Region的基本思想

其实这个算法的原理也很简单,如上图所示的在母材料区域中嵌入一个钢筋,那么最为简单的考虑两者的共同作用的方式就是在交接面上划分兼容网格,或者不兼容网格然后绑定,这两种方法在钢筋和母材料(混凝土)都为实体建模的时候比较有效,而工程中更为高效的方法则是将钢筋用1D(truss或者beam)的方式建模出来,那么上述的方式就不太好用了。整体来说我们更期望下面的这种方式:

而这其中又有两种方式可以采用:

  1. 考虑钢筋与母材之间的滑移 (ERS method)
  2. 认为钢筋和母材之间完全没有相对滑移 (ER method)

Abaqus中默认采用的使第二种方法,这种方法其实简述起来很是简单,就分三步:

  1. 通过几何搜索确认钢筋积分点在母材料中的位置 (逆映射)
  2. 得到该点母材料的位移场插值结果
  3. 计算钢筋的贡献

由于第一条主要就是个几何搜索与逆映射的过程,第二条与第三条其实就可以用两条公式表达

 

其中     为钢筋1D单元上的位移场,而     则为母材的节点位移,     则为母材形函数在钢筋积分点处的插值。那么既然已经知道钢筋积分点处位移的表达方式,则可以通过位移-应变的表达,与应变-应力之间的表达写出其整体的贡献, 以线弹性小变形的truss-in-solid为例:

 

Embedded region支持的单元组合

观察上述图片,虽然Abaqus也存在solid-in-solid的这种组合方式,而在真实应用中并不常用(起码我自己不常用到),我个人来说更多的是用Beam-in-Solid ,Truss-in-Solid 或者 Shell-in-Solid。这里需要注意的几点是:

  • 对于转动自由度,采用Embedded region的方式并不能将其约束住
  • 母单元中对应位置的材料并没有被替换为被嵌入的材料,也就是说那里的刚度其实被计算了两次(钢筋一次,混凝土一次,而现实中有钢筋的地方其实没有混凝土)
  • 在计算嵌入单元的位移等变量时,对于母单元为二阶单元的,其插值可能只选用角节点,而不考虑中间节点。

笔者博士期间写了一篇论文,解决了上述第一、二点问题,并考虑了两种不同材料之间材料属性的不同引入了弱不连续性(    连续但    不连续),非常完美的解决了很多土木、FRP(Fiber reinforced polymer)与FRC (Fiber reinforced composite)中用这种方式计算出来应变、应力不太准的问题。


其实Abaqus中还有另外一种方式,就是可以通过单元技术在实体、壳中按照一定间隔、角度引入Rebar(钢筋层),以此来直接影响刚度矩阵的计算,这个我们后续也来说说。


来源:大狗子说数值模拟
Abaqus材料单元技术
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-12-13
最近编辑:2月前
大狗子说数值模拟
博士 传播国际一流的数值模拟算法
获赞 6粉丝 6文章 45课程 0
点赞
收藏
作者推荐

有限元基础知识:偏导、方向导数、微分与变分

偏导数与方向导数假设有如下的函数那么其两个方向的偏导数则定义为与一般我们也就表示为与单变量函数的求导的物理意义类似,偏导求得的是对某一个变量的变化率,对于多个变量的偏导数可以写为而与偏导数不同,方向导数则是关注函数沿着特定方向的变化率,这里一个非常常见的问题,那比如上边的公式,不就各自是个轴的方向吗?这两者有啥不同?其实这个疑问就特别合理,那我们来看方向导数的定义其中是该方向导数选取的方向,根据这个可以看到,其实偏导数是方向导数的一种特殊形式,诸如上述的公式当的时候就是,而当的时候就是故此我们可以看到方向导数是一种更为通用的形式,在有限元分析中,比如求形函数在参数空间上对某个参数的变化率就是在求偏导,而优化分析中对于多个可能的方向进行求导则是在求方向导数。微分与变分微分大家都很熟悉了,对于一个变量,其微分为,同时我们可以通过偏导数将函数的微分与变量的微分联系起来:而变分则是泛函分析的概念,现在我们有的不是一个变量,而是一个变化的函数,那么我们对其求变分,得到我们现在要求的可能是这样一个泛函那么现在我们将整个函数视为一个变化的量,关注函数在所定义的函数空间上的变化引起的变化,那么上述积分的变分就可以写为那么我们就可以来看一下有限元的概念,对于我们有限元来说,每个单元都有多个自由度,都是独立变量,我们以最为简单的线弹性材料为例,我们看每个单元的应变就是一个对节点位移的函数那么根据最简单的弹性理论,我们的应变能可以表示为那么我们同样的根据能量最小原则,考虑这个函数的变化,对能量这个标量的影响,求极值,就得到了我们常见的有限元变分公式而这也就是变分理论在有限元推到中最常用的地方,这里面大家可以看到我们并不是对某一个节点位移求导,而是对“应变函数”这个函数在泛函空间上求变分。变分和微分的运算法则基本相同,我觉得我不说大家很多也都会那么用的。最后强调一下这几个概念的不同符号:变分(\delta),微分,偏导(\partial),方向导数概念是死的,人是活的,活学活用,懂了就行,好用就行。来源:大狗子说数值模拟

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈