procast如何实现接触单元法

lnavy

本人做砂型铸造热裂纹模拟，涉及到应力场的计算，大家都知道，计算应力场时处理铸件与砂型的边界条件是很重要的，在考虑铸件和铸型（芯）的相互作用时，通常有以下几种处理方法：
) z; n! Y& \  N* I7 P/ c（1）假定型（芯）具有完全退让性，对铸件的收缩没有阻碍，在应力模拟时将型（芯）部分省略。这种处理方法多用在砂型铸造应力模拟中。
& Z/ N& f2 l; c5 b1 k) X（2）假定型（芯）具有完全刚性，在铸件与型（芯）接触的部位的法线方向施加约束，该方法多用于金属性铸造过程的应力模拟。
9 O- @2 b; ~' W# A3 H( O（3）将型（芯）作用简化处理，简化成弹性体或采用等效边界条件来计算。
. A2 c6 p. M) R6 P（4）假定铸件与型（芯）紧密相连，在模拟中永不分开，让型（芯）参与计算。
（5）采用接触力学中的接触单元法对铸件/型（芯）边界进行处理。
9 y% x6 n; v& a3 k! u    （6）建立型（芯）的机械阻力模型，通过对试件的铸型阻力实验，建立阻力计算模型。
     现在本人打算应用接触单元法，应该如何在procast里实现呢？请高手指点。

szcast01

找一下柳百成的书看看