procast如何实现接触单元法

lnavy

本人做砂型铸造热裂纹模拟，涉及到应力场的计算，大家都知道，计算应力场时处理铸件与砂型的边界条件是很重要的，在考虑铸件和铸型（芯）的相互作用时，通常有以下几种处理方法：
0 x& d; U2 [- ^4 a- a$ M2 M# M（1）假定型（芯）具有完全退让性，对铸件的收缩没有阻碍，在应力模拟时将型（芯）部分省略。这种处理方法多用在砂型铸造应力模拟中。
（2）假定型（芯）具有完全刚性，在铸件与型（芯）接触的部位的法线方向施加约束，该方法多用于金属性铸造过程的应力模拟。
（3）将型（芯）作用简化处理，简化成弹性体或采用等效边界条件来计算。
* \. A2 P2 x' z1 m/ x; X. L/ r9 d（4）假定铸件与型（芯）紧密相连，在模拟中永不分开，让型（芯）参与计算。
5 x) D1 o5 C2 s/ M/ G( @: ?6 ~（5）采用接触力学中的接触单元法对铸件/型（芯）边界进行处理。
    （6）建立型（芯）的机械阻力模型，通过对试件的铸型阻力实验，建立阻力计算模型。
     现在本人打算应用接触单元法，应该如何在procast里实现呢？请高手指点。

szcast01

找一下柳百成的书看看