ansys/ls-dyna中接触刚度，惩罚因子如何定义，又如何选择？

以壳体为例           接触刚度k=alfa*area/对角线长度   alfa 为罚函数因子， 如果网格尺寸扩大一倍  如果要维持接触刚度不变 alfa要减小为原来的1/2;  实际调试的过程中，不需要计算这个接触刚度， 总的原则是： 尽可能的增大罚函数因子，只要不出现面向的穿透， 但是不能无限制增大否则出现接触不稳丁