钢管桩结构分析，特别是基底接触与桩顶接触

zgq0218

小弟做了一个海上钢管桩的结构分析，其中有些不明白之处，请给位大侠给予指正。

钢管桩底部我直接定义在了基岩处。而实际上基岩上部还有4米的不比较松散的淤泥等，这块小弟不明白该怎么做，请给予帮助。谢谢

/prep7
et,1,beam4!截面对称
et,2,beam44!截面不对称
et,3,CONTAC52 !接触单元
mp,ex,1,2e11!钢结构
mp,dens,1,7850
mp,ex,2,3.5e10!混凝土结构
mp,dens,2,2500
mp,ex,3,2e15!钢臂
KEYOPT,2,2,1
KEYOPT,2,6,0
KEYOPT,2,9,0
KEYOPT,2,10,0  
KEYOPT,2,7,11  !i结点释放约束
KEYOPT,2,8,11  
!定义截面特性
!钢管φ820 δ10
r,1,0.0254469,2.08728e-3,2.08729e-3,0.820,0.820
!钢管φ273 δ6
r,2,0.005032831,4.48708e-5,4.48708e-5,0.273,0.273
![22a槽钢
!r,3,0.003183413,2.39365e-5,1.577257e-6,0.22,0.077
!r,3,31.83/100/100,2393.65/100/100/100/100,157.73/100/100/100/100,0.11,2.1/100
!RMODIF,3,21,0.11,(7.7-2.1)/100
!RMODIF,3,7,31.83/100/100,2393.65/100/100/100/100,157.73/100/100/100/100,0.1,2.01/100
!RMODIF,3,23,0.11,(7.7-2.1)/100
!A,Izz,Iyy,hz,hy
!2-[20a槽钢
!r,3,0.0057655,3.56074e-5,1.871944e-5,0.2,0.146
!2-[16a槽钢
!r,3,0.004390827,1.732438e-5,1.037914e-5,0.16,0.126

!钢管混凝土φ800
r,4,0.5026548,2.01061e-2,2.01061e-2,0.800,0.800
!2-[22a槽钢
r,5,0.006366827,4.7873e-5,2.31224e-5,0.22,0.154
!垫梁
r,6,0.026688,3.37289e-3,1.709629e-4,0.9,0.4

!钢管，顶坐标为0
h1=15!钢管总长
n,1
n,31,0,-h1
fill
*get,nmax,node,,num,max!获得最大结点号
*do,i,1,nmax,1
n,100+i,nx(i)+2.5,ny(i)
n,200+i,nx(i),ny(i),-2.5
n,300+i,nx(i)+2.5,ny(i),-2.5
*enddo

!钢管单元
*do,i,1,nmax-1,1
type,1
mat,1
real,1
e,i,i+1,101
e,i+100,i+1+100,1
e,i+200,i+1+200,301
e,i+300,i+1+300,201
*enddo

!管内混凝土单元
*do,i,1,nmax-1,1
type,1
mat,2
real,4
e,i,i+1,101
e,i+100,i+1+100,1  
e,i+200,i+1+200,301
e,i+300,i+1+300,201
*enddo

!钢管间顺桥向联系   
type,1
mat,1
real,5
!定义参考点
n,100000,0,1,0
n,100001,0,1,-2.5
NUMSTR,ELEM,1001,  !斜杆
n,1001,0.41,-1
n,1011,2.5-0.41,-4
fill
n,1021,2.5-0.41,-1
n,1031,0.41,-4
fill
*do,i,1,10,1
e,1000+i,1000+i+1,100000
e,1020+i,1020+i+1,100000
*enddo
edele,1010,1012,2
ndele,1026
e,1025,1006,100000
e,1006,1027,100000
type,1
mat,1
real,2
!平杆
n,1041,0.41+0.168,-1
n,1048,2.5-0.41-0.168,-1
fill
n,1051,0.41+0.168,-4
n,1058,2.5-0.41-0.168,-4
fill
*do,i,1,7,1
e,1040+i,1040+i+1,100000
e,1050+i,1050+i+1,100000
*enddo
e,1001,1041,100000
e,1021,1048,100000
e,1031,1051,100000
e,1058,1011,100000
!斜杆
type,1
mat,1
real,5
n,1101,0,-1,-0.41
n,1111,0,-4,-2.5+0.41
fill
n,1121,0,-1,-2.5+0.41
n,1131,0,-4,-0.41
fill
*do,i,1,10,1
e,1100+i,1100+i+1,100001
e,1120+i,1120+i+1,100001
*enddo
edele,1050,1052,2
ndele,1126
e,1125,1106,100001
e,1106,1127,100001
type,1
mat,1
real,2
!平杆
n,1141,0,-1,-0.41-0.168
n,1148,0,-1,-2.5+0.41+0.168
fill
n,1151,0,-4,-0.41-0.168
n,1158,0,-4,-2.5+0.41+0.168
fill
*do,i,1,7,1
e,1140+i,1140+i+1,100001
e,1150+i,1150+i+1,100001
*enddo
e,1101,1141,100001
e,1121,1148,100001
e,1131,1151,100001
e,1158,1111,100001
allsel
!单元数,1001~1040
egen,2,200,1001,1040,1,,,,,,,,-2.5  !向-z方向复制
!单元数,1041~1080
egen,2,200,1041,1080,1,,,,,,2.5,,  !向+x方向复制

!刚臂
type,1
mat,3
real,1
e,3,1001
e,3,1101
e,9,1031
e,9,1131
e,103,1021
e,103,1301
e,109,1011
e,109,1331
e,203,1121
e,203,1201
e,209,1111
e,209,1231
e,303,1221
e,303,1321
e,309,1211
e,309,1311
allsel
!将该部分定义成组
nsel,s,loc,x,0,2.5
nsel,a,loc,z,-2.5,0
esln,s,1,all
cm,kuai1,elem
allsel

egen,2,2000,kuai1,,,,,,,,,,-7.2 !向-z方向复制
!egen,2,4000,kuai1,,,,,,,,11.25,,  !向+x方向复制
!egen,2,6000,kuai1,,,,,,,,11.25,,-7.2  !向对角方向复制
allsel

!钢管组横桥向联系   
NUMSTR,ELEM,4001,  !斜杆
type,1
mat,1
real,5
n,4001,0,-1,-2.5-0.41
n,4011,0,-4,-2.5-4.7+0.41
fill
n,4021,0,-1,-2.5-4.7+0.41
n,4031,0,-4,-2.5-0.41
fill
*do,i,1,10,1
e,4000+i,4000+i+1,100000
e,4020+i,4020+i+1,100000
*enddo
edele,4010,4012,2
ndele,4026
e,4025,4006
e,4006,4027
type,1
mat,1
real,2
!平杆
n,4041,0,-1,-2.5-0.41-0.388
n,4048,0,-1,-2.5-4.7+0.41+0.388
fill
n,4051,0,-4,-2.5-0.41-0.388
n,4058,0,-4,-2.5-4.7+0.41+0.388
fill
*do,i,1,7,1
e,4040+i,4040+i+1,100000
e,4050+i,4050+i+1,100000
*enddo
e,4001,4041
e,4021,4048
e,4031,4051
e,4058,4011
allsel
!单元数4001~4040
egen,2,100,4001,4040,1,,,,,,2.5,,  !向+x方向复制

!刚臂
type,1
mat,3
real,1
e,203,4001
e,209,4031
e,2003,4021
e,2009,4011
e,303,4101
e,309,4131
e,2103,4121
e,2109,4111
allsel
eplot
!将该部分定义成组
nsel,s,loc,x,0,2.5
nsel,a,loc,z,-13.75,0
esln,s,1,all
cm,kuai2,elem
allsel
egen,2,5000,kuai2,,,,,,,,11.25,,  !向+x方向复制
allsel

!钢管组顺桥向联系
!平杆
NUMSTR,ELEM,10001,  
type,1
mat,1
real,2
n,10001,2.5+0.41,-1
n,10021,2.5+8.75-0.41,-1
fill
n,10031,2.5+0.41,-4
n,10051,2.5+8.75-0.41,-4
fill
*do,i,1,20,1
e,10000+i,10000+i+1
e,10030+i,10030+i+1
*enddo
type,1
mat,1
real,5
!斜杆
n,10061,2.5+0.41+0.3965,-1-0.3
n,10069,2.5+4.375-0.3965,-4+0.3
fill
n,10071,2.5+4.375-0.3965,-1-0.3
n,10079,2.5+0.41+0.3965,-4+0.3
fill
*do,i,1,8,1
e,10060+i,10060+i+1
e,10070+i,10070+i+1
*enddo
edele,10048,10050,2
ndele,10075
e,10074,10065
e,10065,10076
e,10001,10061
e,10069,10041
e,10031,10079
e,10071,10011
n,10081,2.5+4.375+0.3965,-1-0.3
n,10089,2.5+8.75-0.41-0.3965,-4+0.3
fill
n,10091,2.5+8.75-0.41-0.3965,-1-0.3
n,10099,2.5+4.375+0.3965,-4+0.3
fill
*do,i,1,8,1
e,10080+i,10080+i+1
e,10090+i,10090+i+1
*enddo
edele,10070,10072,2
ndele,10095
e,10094,10085
e,10085,10096
e,10021,10091
e,10099,10041
e,10051,10089
e,10081,10011
n,10053,2.5+4.375,-1-0.5
n,10057,2.5+4.375,-4+0.5
fill
*do,i,1,4,1
e,10052+i,10052+i+1
*enddo
e,10011,10053
e,10041,10057
!单元数10001~10090
egen,2,100,10001,10090,1,,,,,,,,-2.5  !向-z方向复制
n,10201,2.5+4.375,-1,-0.41
n,10211,2.5+4.375,-1,-2.5+0.41
fill
n,10221,2.5+4.375,-4,-0.41
n,10231,2.5+4.375,-4,-2.5+0.41
fill
*do,i,1,10,1
e,10200+i,10200+i+1
e,10220+i,10220+i+1
*enddo
e,10201,10011
e,10211,10111
e,10221,10041
e,10231,10141
allsel
eplot
!单元数10001~10200
egen,2,6000,10001,10200,1,,,,,,,,-7.2  !向-z方向复制

!刚臂
type,1
mat,3
real,1
e,103,10001
e,109,10031
e,5003,10021
e,5009,10051
e,303,10101
e,309,10131
e,5203,10121
e,5209,10151
e,2103,16001
e,2109,16031
e,7003,16021
e,7009,16051
e,2303,16101
e,2309,16131
e,7203,16121
e,7209,16151
allsel
eplot

!钢管顶垫梁
type,1
mat,1
real,6
n,501,-0.341,0.01,0.341
n,502,,0.01,
n,506,1.25,0.01,-1.25
fill
n,510,2.5,0.01,-2.5
fill
n,511,2.5+0.341,0.01,-2.5-0.341
n,521,-0.341,0.01,-2.5-0.341
n,522,,0.01,-2.5
n,526,1.25,0.01,-1.25
fill
n,530,2.5,0.01,
fill
n,531,2.5+0.341,0.01,0.341
*do,i,1,10,1
e,500+i,500+i+1,100000
e,520+i,520+i+1,100001
*enddo
edele,4996,4998,2
ndele,526
e,525,506,100001
e,506,527,100001
allsel
n,601,-0.341,0.01,-7.2+0.341
n,602,,0.01,-7.2
n,606,1.25,0.01,-7.2-1.25
fill
n,610,2.5,0.01,-7.2-2.5
fill
n,611,2.5+0.341,0.01,-7.2-2.5-0.341
n,621,-0.341,0.01,-7.2-2.5-0.341
n,622,,0.01,-7.2-2.5
n,626,1.25,0.01,-7.2-1.25
fill
n,630,2.5,0.01,-7.2
fill
n,631,2.5+0.341,0.01,-7.2+0.341
*do,i,1,10,1
e,600+i,600+i+1,102000
e,620+i,620+i+1,102001
*enddo
edele,5018,5020,2
ndele,626
e,625,606,102001
e,606,627,102001
allsel

n,701,11.25-0.341,0.01,0.341
n,702,11.25,0.01,
n,706,11.25+1.25,0.01,-1.25
fill
n,710,11.25+2.5,0.01,-2.5
fill
n,711,11.25+2.5+0.341,0.01,-2.5-0.341
n,721,11.25-0.341,0.01,-2.5-0.341
n,722,11.25,0.01,-2.5
n,726,11.25+1.25,0.01,-1.25
fill
n,730,11.25+2.5,0.01,
fill
n,731,11.25+2.5+0.341,0.01,0.341
*do,i,1,10,1
e,700+i,700+i+1,105000
e,720+i,720+i+1,105001
*enddo
edele,5040,5042,2
ndele,726
e,725,706,105001
e,706,727,105001
allsel
n,801,11.25-0.341,0.01,-7.2+0.341
n,802,11.25,0.01,-7.2
n,806,11.25+1.25,0.01,-7.2-1.25
fill
n,810,11.25+2.5,0.01,-7.2-2.5
fill
n,811,11.25+2.5+0.341,0.01,-7.2-2.5-0.341
n,821,11.25-0.341,0.01,-7.2-2.5-0.341
n,822,11.25,0.01,-7.2-2.5
n,826,11.25+1.25,0.01,-7.2-1.25
fill
n,830,11.25+2.5,0.01,-7.2
fill
n,831,11.25+2.5+0.341,0.01,-7.2+0.341
*do,i,1,10,1
e,800+i,800+i+1,107000
e,820+i,820+i+1,107001
*enddo
edele,5062,5064,2
ndele,826
e,825,806,107001
e,806,827,107001
allsel


！此处需要钢管顶垫梁与钢管顶之间做接触分析，小弟在此做的不对，需指正

type,3
mat,3
real,1
e,1,502
e,101,530
e,201,522
e,301,510
e,2001,602
e,2101,630
e,2201,622
e,2301,610
e,5001,702
e,5101,730
e,5201,722
e,5301,710
e,7001,802
e,7101,830
e,7201,822
e,7301,810
allsel
eplot


/solu           
ACEL,0,9.8,0,
nSEL,S,,,31
nSEL,a,,,131
nSEL,a,,,231
nSEL,a,,,331
nSEL,a,,,2031
nSEL,a,,,2131
nSEL,a,,,2231
nSEL,a,,,2331
nSEL,a,,,5031
nSEL,a,,,5131
nSEL,a,,,5231
nSEL,a,,,5331
nSEL,a,,,7031
nSEL,a,,,7131
nSEL,a,,,7231
nSEL,a,,,7331
d,all,ux
d,all,uy
d,all,uz
!d,all,all
allsel

!f,1,fz,-10000*0.05*1e4/16
!f,101,fz,-10000*0.05*1e4/16
!f,201,fz,-10000*0.05*1e4/16
!f,301,fz,-10000*0.05*1e4/16
!f,2001,fz,-10000*0.05*1e4/16
!f,2101,fz,-10000*0.05*1e4/16
!f,2201,fz,-10000*0.05*1e4/16
!f,2301,fz,-10000*0.05*1e4/16
!f,5001,fz,-10000*0.05*1e4/16
!f,5101,fz,-10000*0.05*1e4/16
!f,5201,fz,-10000*0.05*1e4/16
!f,5301,fz,-10000*0.05*1e4/16
!f,7001,fz,-10000*0.05*1e4/16
!f,7101,fz,-10000*0.05*1e4/16
!f,7201,fz,-10000*0.05*1e4/16
!f,7301,fz,-10000*0.05*1e4/16
f,1,fx,-10000*0.05*1e4/16
f,101,fx,-10000*0.05*1e4/16
f,201,fx,-10000*0.05*1e4/16
f,301,fx,-10000*0.05*1e4/16
f,2001,fx,-10000*0.05*1e4/16
f,2101,fx,-10000*0.05*1e4/16
f,2201,fx,-10000*0.05*1e4/16
f,2301,fx,-10000*0.05*1e4/16
f,5001,fx,-10000*0.05*1e4/16
f,5101,fx,-10000*0.05*1e4/16
f,5201,fx,-10000*0.05*1e4/16
f,5301,fx,-10000*0.05*1e4/16
f,7001,fx,-10000*0.05*1e4/16
f,7101,fx,-10000*0.05*1e4/16
f,7201,fx,-10000*0.05*1e4/16
f,7301,fx,-10000*0.05*1e4/16              

f,1,fy,-10000*1e4/16
f,101,fy,-10000*1e4/16
f,201,fy,-10000*1e4/16
f,301,fy,-10000*1e4/16
f,2001,fy,-10000*1e4/16
f,2101,fy,-10000*1e4/16
f,2201,fy,-10000*1e4/16
f,2301,fy,-10000*1e4/16
f,5001,fy,-10000*1e4/16
f,5101,fy,-10000*1e4/16
f,5201,fy,-10000*1e4/16
f,5301,fy,-10000*1e4/16
f,7001,fy,-10000*1e4/16
f,7101,fy,-10000*1e4/16
f,7201,fy,-10000*1e4/16
f,7301,fy,-10000*1e4/16  


allsel
eplot
/SOLU
ANTYPE,STATIC  !定义分析类型为静力分析
PSTRES,ON    !打开预应力开关
ALLSEL,ALL     ！选择整个模型
SOLVE     ！求静力解

/ESHAPE,1.0

/eof
FINISH     ！退出求解器
/SOLU     ！重新进入求解器
ANTYPE,BUCKLE   ！定义分析类型为特征值屈曲分析
BUCOPT,SUBSPACE,5 ！选择特征值的提取方式为子空间迭代法，并且提取前5阶特征值
MXPAND,1    ！扩展第一阶特征值，即求对应的特征向量
MXPAND,2    ！扩展第二阶特征值，即求对应的特征向量
SOLVE     ！求解

/eof
/post1
beam4
beam44
SMAX NMISC -- 1 3
SMIN NMISC -- 2 4
PLNSOL,U,x,0,1  
   
ETABLE, ,NMISC, 1   
ETABLE, ,NMISC, 2   
ETABLE, ,NMISC, 3   
ETABLE, ,NMISC,4   
PLLS,NMIS1,NMIS1,1,0