幫忙看看,如何定接觸關系?
@ANSYS專家
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
*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
!!!!!!!!?????
/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

寧博士CAE團隊 


















