Question

各位ANSYS大侠，下面是论坛里面一个二维桩土接触分析实例，为什么结果X方向云图不对称？

/PREP7
!**有限元单元类型**
ET,1,PLANE42
KEYOPT,1,3,2 !打开PLANE42平面应变选项
ET,2,TARGE169
ET,3,CONTA171
KEYOPT,3,12,0 !选择接触类型(CONTA171)为STANDARD
R,1,,,1,0.0001,, !设置实常数,法向接触刚度因子1，最大渗透范围0.001
RMORE,,,2.6E4 !设置实常数,最大接触摩擦26kPa
!**地层材料模型**
!***粉质粘土***
MP,DENS,1,1950
MP,EX,1,1E7
MP,NUXY,1,0.3
MP,MU,1,0.3 !设置摩擦系数为0.3
TB,DP,1
TBDATA,1,3.0E4,21
!***中粗沙***
MP,DENS,2,2100
MP,EX,2,4E7
MP,NUXY,2,0.25
MP,MU,2,0.9 !设置摩擦系数为0.9
TB,DP,2
TBDATA,1,0,38
!***卵石***
MP,DENS,3,2180
MP,EX,3,5.4E7
MP,NUXY,3,0.2
MP,MU,3,0.1 !设置摩擦系数为0.1
TB,DP,3
TBDATA,1,0,45
!***800钻孔桩C30钢筋砼***
MP,DENS,4,2500
MP,EX,4,3e10
MP,NUXY,4,0.2
!**地层几何模型**
!***粉质粘土层***
BLC4,-20,0,19.6,-5
BLC4,0.4,0,19.6,-5
!***中粗沙层***
BLC4,-20,-5,19.6,-5
BLC4,0.4,-5,19.6,-5
!***卵石层***
BLC4,-20,-10,19.6,-56
BLC4,0.4,-10,19.6,-56
BLC4,-0.4,-23,0.8,-43
ALLSEL,ALL !合并重复节点并重新排序
AGLUE,ALL
NUMCMP,ALL
!***桩模型，桩径0.8m***
BLC4,-0.4,0,0.8,-23
!**有限元网格划分**
!***粉质粘土层***
ASEL,S,,,1
ASEL,A,,,2
ESIZE,1
TYPE,1
MAT,1
AMESH,ALL
!***中粗沙层***
ASEL,S,,,3
ASEL,A,,,4
ESIZE,1
TYPE,1
MAT,2
AMESH,ALL
!***卵石层***
ASEL,S,,,5
ESIZE,1
TYPE,1
MAT,3
AMESH,ALL
LCCAT,22,19
LCCAT,20,23
ASEL,S,,,6
ASEL,A,,,7
ESIZE,1
TYPE,1
MAT,3
AMESH,ALL
!***桩模型***
ASEL,S,,,8
ESIZE,1
TYPE,1
MAT,4
AMESH,ALL
!***接触部分***
!****桩部分接触面****
LSEL,S,,,28
LSEL,A,,,26
NSLL,S,1 !选出需要加接触的节点
TYPE,2 !桩为刚性面TARGET169单元
REAL,1
ESURF,ALL
!****土体部分接触面****
LSEL,S,,,2
LSEL,A,,,8
LSEL,A,,,15
LSEL,A,,,16
LSEL,A,,,22
LSEL,A,,,23
NSLL,S,1 !选出需要加接触的节点
TYPE,3 !土为柔性面CONTAC171单元
REAL,1
ESURF,ALL
ALLSEL,ALL !全部选定，避免前面操作的漏选
FINISH
!*求解*
/SOLU
!**施加约束**
LSEL,S,LOC,X,-20 !选出坐标X=-20的所有线
DL,ALL,,UX !将X=-20的所有线进行X方向约束
LSEL,S,LOC,X,20 !选出坐标X=20的所有线
DL,ALL,,UX !将X=-20的所有线进行X方向约束
LSEL,S,LOC,Y,-66 !选出坐标Y=-66的所有线
DL,ALL,,UY !将Y=-66的所有线进行Y方向约束
ALLSEL,ALL !全部选定，避免前面操作的漏选
!**求解参数设置**
ANTYPE,STATIC !分析类型为静力
ACEL,,9.8 !施加重力加速度9.8
NLGEOM,ON !打开大变形选项
NROPT,FULL,,ON !打开牛顿-辛普生选项，加快收敛速度
SOLVE

Answer 1

天知，地知，你知，Ansys知。