procast如何实现接触单元法

lnavy

本人做砂型铸造热裂纹模拟，涉及到应力场的计算，大家都知道，计算应力场时处理铸件与砂型的边界条件是很重要的，在考虑铸件和铸型（芯）的相互作用时，通常有以下几种处理方法：
! a. p; y& |0 n0 [（1）假定型（芯）具有完全退让性，对铸件的收缩没有阻碍，在应力模拟时将型（芯）部分省略。这种处理方法多用在砂型铸造应力模拟中。
% j. n# |8 Y) F' v7 o' [（2）假定型（芯）具有完全刚性，在铸件与型（芯）接触的部位的法线方向施加约束，该方法多用于金属性铸造过程的应力模拟。
- h3 B! {5 T' T' ~2 b6 F& l（3）将型（芯）作用简化处理，简化成弹性体或采用等效边界条件来计算。
（4）假定铸件与型（芯）紧密相连，在模拟中永不分开，让型（芯）参与计算。
- q/ D% d7 l7 F: @. w1 \（5）采用接触力学中的接触单元法对铸件/型（芯）边界进行处理。
    （6）建立型（芯）的机械阻力模型，通过对试件的铸型阻力实验，建立阻力计算模型。
     现在本人打算应用接触单元法，应该如何在procast里实现呢？请高手指点。

szcast01

找一下柳百成的书看看