@Jamesonl 您好！感谢您回答我的问题。就这个问题，目前的关键是不能很好的建立接触关系。Abqus我还不怎么熟悉，只会点皮毛，ansys使用的时间比较长一些。 关于O型圈的大变形，挤入间隙，我计算这个模型的目的就是看铜三角环能否将间隙闭合，还有就是如何调整铜三角环的参数，例如尺寸，弹性模量等，使得在加高压时（O型圈左侧和缸筒内壁加高压，这个结构可以承载100MPa左右，我目前已经应用成功，现在想用有限元验算下），铜三角环的径向变形能对缸筒的径向变形进行补偿。 您看这个模型的建立是不是需要改变一下呢？沟槽我建立的线，然后直接划分为刚性TARGE169单元，如果先建立节点，再划分是不是可以重复划分接触对呢？我做大型接触机器的分析时，有时将接触面划分两次接触对，即：对称接触对，两个面即是接触面，又是目标面，这是对节点进行划分的。但这个平面模型，在已划分过接触单元的面，如中间的O型圈，重复划分接触单元时，会弹出一个提示，说是该面已经划分，在划分需要删除已存在的单元，忽略，检查接触对时，也发现只有TARGE面，没有contact面。 高手，能否帮助解答，研究一下呢？ 附上部分命令流，以供验证，谢谢！ 其中接触部分未完成，计算不了，我用“！！！！？？？”隔开，需要调整。 FINISH /OUTPUT,OUTPUT LOG,TXT /CLEAR,START /TITLE, ULTRAHIGH PRESSURES SEAL ANLYSIS /FILNAME,ULTRAHIGH PRESSURES SEAL ANLYSIS,1 /PNUM,LINE,1 !定义参数 /PREP7 *AFUN,DEG R0=5.7/2 !!O型圈半径 R1=100.5 !!O型圈槽底径 R2=105 !!缸筒内径 R3=175 !!缸筒外径 T=1 !!O型圈压力侧活塞杆与缸筒间隙 T1=8.5 !!槽宽 B=10 P=80 !!液体压力，MPa ES=1 !!单元尺寸 !!定义单元 ET,1,PLANE82 KEYOPT,1,3,2 !!Plane strain !!橡胶材料单元 ET,2,PLANE183 KEYOPT,2,1,0 KEYOPT,2,3,2 !!Plane strain KEYOPT,2,6,1 !!u-P公式 !!定义橡胶材料参数 EX1=7.8 NU1=0.5 C10=1.84 C01=0.47 DD=2*(1-2*NU1)/(C10+C01) !!铜环材料参数 MP,EX,1,1.08E5 MP,NUXY,1,0.34 !!橡胶材料参数 MP,EX,2,EX1 MP,NUXY,2,NU1 TB,HYPER,2,1,2,MOONEY TBDATA,2,C10,C01,DD !!缸筒材料参数，45# MP,EX,3,2.06E5 MP,NUXY,3,0.3 !!接触摩擦系数 MP,MU,4,0.01 !!定义点的坐标 *DIM,X,ARRAY,14,1 *DIM,Y,ARRAY,14,1 X(1,1)=0 X(2,1)=B X(3,1)=B X(4,1)=B+T1 X(5,1)=B+T1+(R2-R1)*SQRT(2) X(6,1)=2*B+T1+(R2-R1)*SQRT(2) X(7,1)=B+T1 X(8,1)=B+T1 X(9,1)=B+T1+(R2-R1)*SQRT(2) X(10,1)=B+T1-R0 X(11,1)=0 X(12,1)=2*B+T1+(R2-R1)*SQRT(2) X(13,1)=2*B+T1+(R2-R1)*SQRT(2) X(14,1)=0 Y(1,1)=R2-T Y(2,1)=R2-T Y(3,1)=R1 Y(4,1)=R1 Y(5,1)=R2 Y(6,1)=R2 Y(7,1)=R1 Y(8,1)=R2 Y(9,1)=R2 Y(10,1)=R1+R0 Y(11,1)=R2+1.5*T Y(12,1)=R2+1.5*T Y(13,1)=R3 Y(14,1)=R3 !!创建点 *DO,I,1,14,1 K,I,X(I,1),Y(I,1), *ENDDO *SET,I !!创建沟槽 *DO,I,1,5,1 LSTR,I,I+1 *ENDDO CM,CAO,LINE !!创建沟槽组 *SET,I !!创建铜环截面 A,7,8,9 !!创建O型圈 CYL4,X(10,1),Y(10,1),R0 !!创建缸筒截面 A,11,12,13,14 !!划分铜环 ESIZE,ES TYPE,1 MAT,1 AMESH,1 !!划分O型圈橡胶 ESIZE,ES TYPE,2 MAT,2 AMESH,2 !!划分缸筒 ESIZE,ES TYPE,1 MAT,3 AMESH,3 !!细化接触单元 LSEL,S,,,6,13,1 LPLOT LSLL,S,1 NSLL,S,1 NREFINE,ALL,,,1 ALLSEL,ALL LPLOT SAVE !!!!!!!!!!!???????????????? K=3 ET,K,TARGE169 ET,K+1,TARGE169 ET,K+2,CONTA172 R,K,,,0.2 !!接触刚度 KEYOPT,K+2,2,1 !!接触算法：罚函数法 KEYOPT,K+2,5,4 KEYOPT,K+2,9,1 KEYOPT,K+2,10,1 KEYOPT,K+2,12,0 MAT,4 REAL,1 TYPE,K TSHAPE,LINE LSEL,S,,,1,5 LMESH,ALL ESURF,,REVERSE LSEL,S,,,6 LSEL,A,,,13 NSLL,S,1 TYPE,K+1 ESURF,ALL !!!!!!!!!!!???????????????? /SOLU ANTYPE,0 SOLCONTROL,ON NLGEOM,ON OUTRES,ALL,ALL LNSRCH,1 !!线性搜索打开 !!耦合轴边线UY自由度,及模拟安装压缩量 LSEL,S,,,1,5 NSLL,S,1 CP,1,UY,ALL D,ALL,UX,0 D,ALL,UY,2.5*T-0.01 LSEL,S,LINE,,16 NSLL,S,1 D,ALL,UX,0 ALLSEL,ALL EPLOT KBC,0 TIME,1 CNVTOV,U,,0.05,2,0.01 CNVTOV,F,,0.02,2,0.01 NSUBST,10 ALLSEL,ALL SOLVE FINISH /POST1 /EFACET,1 PLNSOL,U,Y PLNSOL,S,EQV