"contact cmat default model linear"

1 引言

由于近日遭遇到了暂时无法解决的生活难题, 打破了原有的平静心情, 再加上下周就要开学, 因此彻彻底底地放纵了几天. Anyway, 现在开始逐渐恢复到以前的生活节奏, 休息休息, 先用线性模型热热身.

线性接触模型"contact cmat default model linear" 是PFC最简单也是最常用的模型, 它是PFC最早发展的一个线性模型. 在这个模型中,  不活动活塞和参考间隙为零，可以安装在ball-ball和ball-facet接触上，在FISH中称为linear。 基于线性模型的表面间隙如下图所示.

如果且仅当表面间隙小于或等于零时，线性模型的接触才是活动的。对于不活动的触点，不计算力-位移。当表面间隙为零时，名义表面与实际表面重合。

如果对应的接触类型不能使用指定的接触模型，则不会修改默认槽。线性模型提供了相互平行作用的线性和活塞组件。线性组件提供线性弹性（无张力）摩擦行为，而活塞组件提供粘性行为。两种分量都在一个很小的区域内作用，因此只传递一个力。本笔记简要总结了这个模型最常用的参数设置.

2 Method

method用来设置linear模型的变形参数, 只有两个参数: emod是颗粒的有效模量; kratio是颗粒接触法向与切向刚度之比. 如果不考虑变形, 在模拟中可以不设置method.

method deformability emod 1.0e9 kratio 0.0

method deformability emod [materialMod] kratio 1.0 

3 Property

property用来设置接触的属性, 即使是线性模型,这部分设置也相当复杂. 最常用的参数包括:

(1) kn: 法向刚度

(2) ks: 切向刚度

(3) fric: 摩擦系数

(4) dp_nratio: 法向临界阻尼比

(5) dp_mode: 阻尼模式

 0: 全法向和全剪切

 1: 非拉伸的法向和全剪切

 2: 全法向和滑落-切断剪切

 3: 非拉伸的法向和滑落-切断剪切

(6) lin_mode: 法向力的更新模式, 默认值是0

 0: 更新是绝对的

 1: 更新是增量式的

property fric 0.3 dp_nratio 0.5 dp_mode 1 lin_mode 1

property kn 1.0e6 ks 1.0e6 fric 0.25 dp_nratio 0.1

property kn 1e4 ks 5e3 fric @cfric dp_nratio 0.2

prop dp_nratio 0.8 dp_mode 1