含裂纹的模型,用ANSYS网格划分

/prep7

et,1,plane183 !!!!!八节点平面奇异性单元

mp,ex,1,236.036e4 !!!!!弹性模量
mp,nuxy,1,0.267 !!!!!泊松比

k,1,0,0
k,2,0.006E-3,0
k,3,0,0.006E-3
k,4,-0.006E-3,0

csys,1
l,1,2
l,2,3
l,3,4
l,4,1

a,1,2,3,4

KSCON,1,0.002E-3,1,8,1,

LESIZE,2, , ,4, , , , ,1
LESIZE,3, , ,4, , , , ,1

amesh,1

csys,0

k,5,0.02E-3,0
k,6,0.02E-3,0.02E-3
k,7,0,0.02E-3
k,8,-0.02E-3,0.02E-3
k,9,-0.02E-3,0

l,2,5
l,5,6
l,6,7
l,7,8
l,8,9
l,9,4
l,7,3

al,5,6,7,11,2
al,8,9,10,3,11

LESIZE,6, , ,2, , , , ,1
LESIZE,7, , ,2, , , , ,1
LESIZE,8, , ,2, , , , ,1
LESIZE,9, , ,2, , , , ,1

LESIZE,5, , ,4, , , , ,1
LESIZE,10, , ,4, , , , ,1
LESIZE,11, , ,4, , , , ,1

AMAP,2,3,2,5,7
AMAP,3,4,3,7,9

k,10,0.02E-3,0.04E-3
k,11,-0.02E-3,0.04E-3

a,6,10,11,8,7

LESIZE,12, , ,2, , , , ,1
LESIZE,13, , ,4, , , , ,1
LESIZE,14, , ,2, , , , ,1
amesh,4

k,12,0.38E-3,0
k,13,0.38E-3,0.4E-3
k,14,-0.02E-3,0.4E-3

a,5,6,10,13,12
a,10,13,14,11
LESIZE,16, , ,4, , , , ,1
LESIZE,18, , ,4, , , , ,1

LESIZE,17, , ,10,1/20 , , , ,1
LESIZE,15, , ,10,20 , , , ,1
LESIZE,19, , ,10,1/20 , , , ,1

AMAP,5,10,5,12,13
AMAP,6,11,10,13,14











!!!!!!!!面7
k,15,0.98E-3,0
k,16,0.98E-3,1E-3
k,17,-0.02E-3,1E-3
a,12,15,16,17,14,13

LESIZE,20,0.1E-3
LESIZE,21,0.1E-3
LESIZE,22,0.1E-3
LESIZE,23,0.1E-3
amesh,7

!!!!!!!!!!!!!!!以上为网格划分，下面为计算选项设置，包括施加约束、荷载、以及求解选项设置

/SOLU

!!!!!!左边约束x方向位移，等效对称约束
DL,9, ,UX,
DL,14, ,UX,
DL,19, ,UX,
DL,23, ,UX,
!!!!!!下边除了裂纹面外，约束y方向位移，等效对称约束
DL,1, ,UY,
DL,5, ,UY,
DL,17, ,UY,
DL,20, ,UY,

!!!!!!上边约束x，y两方向位移；等效固定约束
DL,22, ,UX,
DL,22, ,UY,

!!!!!!!右边施加均布荷载，大小只需修改1000E3的地方即可
SFL,21,PRES,1000E3,

ANTYPE,0
DELTIM,0.1,0,0
OUTRES,ERASE
OUTRES,ALL,ALL
TIME,1

solve
save
finish