接触分析实例---包含初始间隙(include notation)

zq0046

fini
/clear,nostart

/prep7
et,1,82
KEYOPT,1,3,3
r,1,0.5
mp,ex,1,1e9
mp,prxy,1,0.3

k,1,0,0
k,2,10,0
k,3,10,5
k,4,6.2,5
k,5,7.5,3.4
k,6,2.5,3.4
k,7,3.8,5
k,8,0,5
a,1,2,3,4,5,6,7,8
LFILLT,6,5,0.18, ,   
LFILLT,5,4,0.18, ,   
FLST,2,3,4  
FITEM,2,9   
FITEM,2,11  
FITEM,2,10  
AL,P51X
FLST,2,3,4  
FITEM,2,13  
FITEM,2,14  
FITEM,2,12  
AL,P51X
FLST,2,3,5,ORDE,2   
FITEM,2,1   
FITEM,2,-3  
AADD,P51X

rect,0,10,4.8,5
ASBA,       4,       1  

gap=0.02
k,24,6.2-gap,5
k,25,7.5-gap,3.4
k,26,2.5+gap,3.4
k,27,3.8+gap,5
a,24,25,26,27

LFILLT,4,3,0.2, ,   
LFILLT,3,2,0.2, ,   
FLST,2,3,4  
FITEM,2,7   
FITEM,2,10  
FITEM,2,8   
AL,P51X
FLST,2,3,4  
FITEM,2,13  
FITEM,2,14  
FITEM,2,11  
AL,P51X
FLST,3,2,5,ORDE,2   
FITEM,3,3   
FITEM,3,-4  
ASBA,       1,P51X  

rect,3.8+gap,6.2-gap,5,10
rect,3.8+gap,3.8+gap+8,10,12
FLST,2,3,5,ORDE,3   
FITEM,2,1   
FITEM,2,3   
FITEM,2,5   
AADD,P51X   
rect,3.8+gap+8,3.8+gap+8+2,10,12
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,4   
AGLUE,P51X  

CYL4,2.0,1.8,0.6  
CYL4,7.0,1.8,0.6
FLST,2,3,5,ORDE,3   
FITEM,2,2   
FITEM,2,4   
FITEM,2,-5  
AOVLAP,P51X

esize,0.2
amesh,all

FLST,5,135,2,ORDE,32
FITEM,5,485
FITEM,5,576
FITEM,5,-577
FITEM,5,621
FITEM,5,-625
FITEM,5,707
FITEM,5,-711
FITEM,5,716
FITEM,5,741
FITEM,5,-745
FITEM,5,750
FITEM,5,-751
FITEM,5,766
FITEM,5,797
FITEM,5,-798
FITEM,5,854
FITEM,5,888
FITEM,5,-938
FITEM,5,1101
FITEM,5,1103
FITEM,5,1420
FITEM,5,1628
FITEM,5,1653
FITEM,5,1696
FITEM,5,1699
FITEM,5,-1702   
FITEM,5,1726
FITEM,5,-1728   
FITEM,5,1852
FITEM,5,-1874   
FITEM,5,2044
FITEM,5,-2066   
CM,_Y,ELEM  
ESEL, , , ,P51X
CM,_Y1,ELEM
CMSEL,S,_Y  
CMDELE,_Y   
EREF,_Y1, , ,1,0,1,1
CMDELE,_Y1  


ET,2,TARGE169   
ET,3,CONTA172     
R,3, , ,
R,3,0,0,0.1, 10,0,0
R,4, , ,
R,4,0,0,0.1, 10,-0.02,0  


lsel,s,,,9
lsel,a,,,5
lsel,a,,,12
nsll,s,1
type,3
real,3
esurf,all
alls,

lsel,s,,,19
lsel,a,,,20
nsll,s,1
type,3
real,4
esurf,all
alls,

lsel,s,,,7
lsel,a,,,3
lsel,a,,,11
nsll,s,1
type,2
real,3
esurf,all
alls,

lsel,s,,,25
lsel,a,,,26
nsll,s,1
type,2
real,4
esurf,all
alls,

FLST,2,2,5,ORDE,2   
FITEM,2,4   
FITEM,2,-5  
DA,P51X,ALL,
FLST,2,1,4,ORDE,1   
FITEM,2,6   
SFL,P51X,PRES,500,  

/solu
antype,0
nlgeom,on
outres,all,all
nsubst,200,200,2
neqit,1000
solve

[ 本帖最后由 zq0046 于 2006-8-4 16:37 编辑 ]

jianghaibo

为啥没人顶呢？我来

jinqinghai

把本例加上注释, 可以发到接触版领加分.

lxrswf

顶,很好啊,但如果要在裂隙中加入水压力该怎么加啊?谢谢!!

njweiwei2

是个什么样的实际问题，能说说吗？

lxrswf

现实是我要模拟的地下洞室长期处于地下水中,而岩石有许多的裂隙存在,要模拟地下水压力对洞室稳定性影响,但我现在不知道水压力该用什么命令加上去,如何加?真急啊!

zq0046

The first question I found is the following net.
If you want to know about more information, pls refer to it.

http://www.simwe.com/forum/viewt ... ghlight=%2Bjstj2000

majing613

先顶顶，在好好看看

leiy

1# zq0046 谢谢了 呵呵

dyoulinghot

为什么总是要用命令流的形式啊，哎

1143942969

学习了！！！

0418142

学习参考！！！！

firemanxll

非常不错的接触分析

ningfengping

命令看着很麻烦，不知什么意思，谁能帮忙解读下，谢谢

shaneluv

先收下了，接触分析不会做。。。

xss

/prep7
et,1,82
keyopt,1,3,3
r,1,0.5
mp,ex,1,1e9
mp,prxy,1,0.3

k,1,0,0
k,2,10,0
k,3,10,5
k,4,6.2,5
k,5,7.5,3.4
k,6,2.5,3.4
k,7,3.8,5
k,8,0,5
a,1,2,3,4,5,6,7,8
lfillt,6,5,0.18,
lfillt,5,4,0.18,
lsel,s,line,,9
lsel,a,line,,11
lsel,a,line,,10
al,all
lsel,s,line,,13
lsel,a,line,,14
lsel,a,line,,12
al,all
asel,s,area,,1,3
aadd,all

rect,0,10,4.8,5
asba,4,1

gap=0.02
k,24,6.2-gap,5
k,25,7.5-gap,3.4
k,26,2.5+gap,3.4
k,27,3.8+gap,5
a,24,25,26,27

lfillt,4,3,0.2,
lfillt,3,2,0.2,
lsel,s,line,,7
lsel,a,line,,10
lsel,a,line,,8
al,all
lsel,s,line,,13
lsel,a,line,,14
lsel,a,line,,11
al,all
asel,s,area,,3,4
asba,1,all

rect,3.8+gap,6.2-gap,5,10
rect,3.8+gap,3.8+gap+8,10,12
asel,s,area,,1
asel,a,area,,3
asel,a,area,,5
aadd,all
rect,3.8+gap+8,3.8+gap+8+2,10,12
asel,s,area,,1
asel,a,area,,4
aglue,all

cyl4,2.0,1.8,0.6
cyl4,7.0,1.8,0.6
asel,s,area,,2
asel,a,area,,4,5
aovlap,all

esize,0.2
amesh,all

esel,s,elem,,485
esel,a,elem,,576,577
esel,a,elem,,621,625
esel,a,elem,,707,711
esel,a,elem,,716
esel,a,elem,,741,745
esel,a,elem,,750,751
esel,a,elem,,766
esel,a,elem,,797,798
esel,a,elem,,854
esel,a,elem,,888,938
esel,a,elem,,1101
esel,a,elem,,1103
esel,a,elem,,1420
esel,a,elem,,1628
esel,a,elem,,1653
esel,a,elem,,1696
esel,a,elem,,1699,1702
esel,a,elem,,1726,1728
esel,a,elem,,1852,1874
esel,a,elem,,2044,2066
eref,all,,,1,0,1,1
et,2,targe169
et,3,conta172
r,3,,
r,3,0,0,0.1,10,0,0
r,4,,
r,4,0,0,0.1,10,-0.02,0

lsel,s,,,9
lsel,a,,,5
lsel,a,,,12
nsll,s,1
type,3
real,3
esurf,all
alls

lsel,s,,,19
lsel,a,,,20
nsll,s,1
type,3
real,4
esurf,all
alls

lsel,s,,,7
lsel,a,,,3
lsel,a,,,11
nsll,s,1
type,2
real,3
esurf,all
alls

lsel,s,,,25
lsel,a,,,26
nsll,s,1
type,2
real,4
esurf,all
alls

asel,s,area,,4,5
da,all,all
lsel,s,line,,6
sfl,all,pres,500

/solu
antype,0
nlgeom,on
outres,all,all
nsubst,200,200,2
neqit,1000
solve