固体力学中怎么解决接触问题啊

sky007

求助各位看看受压下的接触问题啊，谢谢啊

问题描述：利用固体力学模块建立带有圆柱形空腔的圆柱体（内外圆柱中心相对应）二维轴对称模型（材料使用线弹性材料）。圆柱体的半径为25mm，圆柱腔的半径为15mm，材料的杨氏模量为10MPa，泊松比为0.49，密度为1300Kg/m3。圆柱体外边界指定法向方向上的位移为0，在外圆柱体最上边界施加压力，范围从0MPa～5MPa，下边界固定，要求计算圆柱体的稳态变形。在压力增大到一定程度时，圆柱形空腔的水平边界和垂直边界会在相交的地方发生接触。我们在模型中加入接触边界条件后，边界位移会出现突变的凸起，甚至是交叉（两边界彼此渗入到对方里面），接触面不是光滑过渡。请问如何求解包含接触的稳态变形或在求解接触问题过程中要注意什么问题。如有时间，希望各位大侠求解一下该模型并作说明，或给出建议。谢谢！