曲面与平面接触分析探讨

borine

如下图实例。

载荷及约束：筒体内圈施加面压力，大小沿Y坐标单调递增。平板底面全约束。
分析经过：
1、直接添加两接触面间的接触对，加0.18的摩擦系数，计算无法进行，提示某节点Y向位移超标；
        似乎是因为划分网格后，接触对的曲面最下方节点离平板尚有0.034左右的距离，查资料，发现需要设置接触对相关的实常数。
2、设置接触环，ICONT=-0.1，等效应力结果约300MPa，感觉不错；
3、由于初始状态，曲面离平面距离仅有0.034，因此感觉ICONT值设置有些偏大，可能会导致更多节点被限制在接触环上，参与接触。
      重新设置常数，ICONT=-0.035，保证了初始接触节点的正确选取，但应力结果增至600MPa；
4、由于初始网格密度为5，感觉有些偏大，重新划分，改至密度2，并重新设置ICONT值，结果有800MPa。

几次尝试之下，头脑一片空白。
圆筒外径200有余，网格密度应该比较合理，计算收敛性也不错，为何会有如此大的差别？莫非我接触对设置还存在问题？接触刚度是否需要设置？。。。。。。
另：有论文提及接触单元的接触应力（是否就是Contact total stress ？），且比结构的等效应力大很多，似乎不太合理吧。

最近刚开始尝试接触分析，仅此例还有不少疑问，自己先翻翻书，再来请教各位。

仿勒个真

哥们这方面的问题解决了没  我刚好做有曲率地表情况下基础与地基的相互作用 希望共同探讨。我的QQ：250954424