以连续体模型为基准,研究接触非线性模型的计算精度,包括压力和位移加载方式的影响、渗透容差的影响、非匹配网格的影响、接触面和目标面选取的影响。
目标面和接触面的选取原则:
- 如果一个表面比另一个表面硬,则硬面为目标面。
- 如果一个表面大于另一个表面,则大的表面为目标面。
- 如果一个表面为高阶另一个表面为低阶,则低阶表面为目标面。
- 如果一个凸表面与一个平面或凹面接触,应该选择平面或凹面为目标面。
- 如果一个表面有粗糙的网格而另一个表面网格细密,则选择粗糙的网格表面为目标面。
![657b5c8c7d37f45ee5155838f314ae76.png](https://img-blog.csdnimg.cn/img_convert/657b5c8c7d37f45ee5155838f314ae76.png)
![55c9bc47e28e388de75f2e847b452cbc.png](https://img-blog.csdnimg.cn/img_convert/55c9bc47e28e388de75f2e847b452cbc.png)
一、问题描述
为分析接触非线性问题,设计了2个模型。第1个模型为1块板;第2个模型是将第1个板切成2块板,切割成的2块板之间创建接触关系。第1个为连续体模型,以此为基准,研究第2个接触非线性模型的计算精度。
![3732c6cc93f23e80b3dcb4f704b64fac.png](https://img-blog.csdnimg.cn/img_convert/3732c6cc93f23e80b3dcb4f704b64fac.png)
![604451d4f7519b887873517d55301966.png](https://img-blog.csdnimg.cn/img_convert/604451d4f7519b887873517d55301966.png)
![c0534562f20aaea92c51fd962d203a4b.png](https://img-blog.csdnimg.cn/img_convert/c0534562f20aaea92c51fd962d203a4b.png)
问题分析与计算结果:
1.连续体模型
连续体模型仍然切割成2块板,两块板几何面之间的关键点和线是共用的,划分单元后,两块板之间的节点也是共用的。有限元参与计算的是单元和节点,与多少个几何体无关,只要模型的节点和单元连续,建立的仍然是连续体模型。
下图中,面编号以A开头,线编号以L开头,关键点编号是纯数字。上下几何面A1和A2之间共用关键点和线。
![cd13da82ed883c5162a4bf2819e9ea95.png](https://img-blog.csdnimg.cn/img_convert/cd13da82ed883c5162a4bf2819e9ea95.png)
下图是节点编号,在红线位置处为上下几何面之间共用节点。
![be7947631e138058dc9be94c8caa9162.png](https://img-blog.csdnimg.cn/img_convert/be7947631e138058dc9be94c8caa9162.png)
下图是边界条件,施加压力100MPa
![7b474c947843944dd75be9d422e1b6e1.png](https://img-blog.csdnimg.cn/img_convert/7b474c947843944dd75be9d422e1b6e1.png)
下图是整个板Y方向位移云图,往下压缩了0.0455 mm,假想的下板压缩了0.01140 mm。
![f07c833b5731b5d5cdf0829fddce558e.png](https://img-blog.csdnimg.cn/img_convert/f07c833b5731b5d5cdf0829fddce558e.png)
下图是整个板X和Y方向应力云图,分别为-30 MPa和-100 MPa。
![aadcae54bfd93822012ad8dea9f0db01.png](https://img-blog.csdnimg.cn/img_convert/aadcae54bfd93822012ad8dea9f0db01.png)
2.接触非线性模型
切割成2块板,两块板几何面之间的关键点和线是不共用的,划分单元后,两块板之间的节点也是不共用的。有限元参与计算的单元和节点是不连续的,需要创建接触关系,采用接触非线性模型才能计算。
![632d437580794c122e1b44fc78245c4a.png](https://img-blog.csdnimg.cn/img_convert/632d437580794c122e1b44fc78245c4a.png)
上下板的节点不匹配,上下板之间在红色位置处是初始无间隙也没有嵌入。下板单元n1取11份,上板单元n2取22份