如何模拟钢管混凝土
如何模拟钢管和混凝土之间的接触，用什么单元比较好呢！网格划分的时候应该注意些什么？ 这个问题比较复杂，很多参数选择比较困难，除非能够做实验来获取，比如摩擦系数等
下面转两个例子你可以参考一下
[ 本帖最后由 AaronSpark 于 2006-11-11 07:37 编辑 ] 《Algor、Ansys在桥梁工程中的应用方法与实例》中最后一个例子
finish
/clear
/prep7
/title,the analyse of steel arch bridge
et,1,82
cyl4,0.6,0.6,0.6,,0.586
cyl4,0.6,2.4,0.6,,0.586
rectng,0.35,0.364,0.8,2.2
rectng,0.85,0.836,0.8,2.2
asel,all
aptn,all
adele,5,8,1
allsel
aadd,all
smrtsize,5
amesh,all
secwrite,gg,sect,,1
sectype,1,beam,mesh,
secoffset,cent,,,
secread,'gg','sect',,mesh
asel,all
aclear,all
adele,all,,,1
/replot
cyl4,0.6,0.6,0.586
cyl4,0.6,2.4,0.586
rectng,0.364,0.836,0.8,2.2
allsel
aadd,all
smrtsize,5
amesh,all
secwrite,hnt,sect,,1
sectype,2,beam,mesh,
secoffset,cent,,,
secread,'hnt','sect',,mesh
asel,all
aclear,all
adele,all,,,1
/replot
k,1,-0.5,0
k,2,-0.75,0.25
k,3,-0.75,0.5
k,4,-0.55,0.5
k,5,-0.55,0.85
k,6,0.55,0.85
k,7,0.55,0.5
k,8,0.75,0.5
k,9,0.75,0.25
k,10,0.5,0
a,1,2,3,4,5,6,7,8,9,10
smrtsize,5
amesh,all
secwrite,hl1,sect,,1
sectype,4,beam,mesh,
secoffset,cent,,,
secread,'hl1','sect',,mesh
asel,all
aclear,all
adele,all,,,1
/replot
k,1,-0.5,0
k,2,-0.5,1
k,3,-0.75,1.25
k,4,-0.75,1.5
k,5,0.75,1.5
k,6,0.75,1.25
k,7,0.5,1
k,8,0.5
a,1,2,3,4,5,6,7,8
smrtsize,5
amesh,all
secwrite,hl2,sect,,1
sectype,5,beam,mesh,
secoffset,cent,,,
secread,'hl2','sect',,mesh
asel,all
aclear,all
adele,all,,,1
/replot
k,1,-0.5,0
k,2,-0.5,1.04
k,3,-0.75,1.04+0.25
k,4,-0.75,1.04+0.5
k,5,-0.55,1.04+0.5
k,6,-0.55,1.04+0.5+0.35
k,7,0.55,1.04+0.5+0.35
k,8,0.55,1.04+0.5
k,9,0.75,1.04+0.5
k,10,0.75,1.04+0.25
k,11,0.5,1.04
k,12,0.5,0
k,13,-0.2,0.58
k,14,-0.2,1.04+0.5-0.3
k,15,0.2,1.04+0.5-0.3
k,16,0.2,0.58
a,1,2,3,4,5,6,7,8,9,10,11,12
a,13,14,15,16
asba,1,2
smrtsize,5
amesh,all
secwrite,hl3,sect,,1
sectype,6,beam,mesh,
secoffset,cent,,,
secread,'hl3','sect',,mesh
asel,all
aclear,all
adele,all,,,1
/replot
ldele,all
kdele,all
k,1,-0.5,0
k,2,-0.5,1.23
k,3,-0.75,1.23+0.25
k,4,-0.75,1.23+0.5
k,5,-0.55,1.23+0.5
k,6,-0.55,1.23+0.5+0.35
k,7,0.55,1.23+0.5+0.35
k,8,0.55,1.23+0.5
k,9,0.75,1.23+0.5
k,10,0.75,1.23+0.25
k,11,0.5,1.23
k,12,0.5,0
k,13,-0.2,0.38
k,14,-0.3,0.38+0.2
k,15,-0.3,1.23+0.5-0.3
k,16,-0.2,1.23+0.5-0.1
k,17,0.2,1.23+0.5-0.1
k,18,0.3,1.23+0.5-0.1-0.2
k,19,0.3,0.38+0.2
k,20,0.2,0.38
a,1,2,3,4,5,6,7,8,9,10,11,12
a,13,14,15,16,17,18,19,20
asba,1,2
smrtsize,5
amesh,all
secwrite,hl4,sect,,1
sectype,7,beam,mesh,
secoffset,cent,,,
secread,'hl4','sect',,mesh
asel,all
aclear,all
adele,all,,,1
/replot
k,1,-0.5,0
k,2,-0.5,1.23
k,3,-0.75,1.23+0.25
k,4,-0.75,1.23+0.5
k,5,-0.55,1.23+0.5
k,6,-0.55,1.23+0.5+0.35
k,7,0.55,1.23+0.5+0.35
k,8,0.55,1.23+0.5
k,9,0.75,1.23+0.5
k,10,0.75,1.23+0.25
k,11,0.5,1.23
k,12,0.5,0
a,1,2,3,4,5,6,7,8,9,10,11,12
smrtsize,5
amesh,all
secwrite,hl5,sect,,1
sectype,8,beam,mesh,
secoffset,cent,,,
secread,'hl5','sect',,mesh
asel,all
aclear,all
adele,all,,,1
/replot
etdele,1
/replot
et,1,beam44
mp,ex,1,2.1e11
mp,dens,1,7800
mp,prxy,1,0.3
n,90000,0,0,30
et,2,beam44
mp,ex,2,3.5e10
mp,dens,2,2600
mp,prxy,2,0.1667
et,3,beam44
mp,ex,3,3.0e10
mp,dens,3,2600
mp,prxy,3,0.1667
et,4,beam44
mp,ex,4,3.0e10
mp,dens,4,2600
mp,prxy,4,0.1667
et,5,beam44
mp,ex,5,2.1e11
mp,dens,5,7800
mp,prxy,5,0.3
et,6,link10
mp,ex,6,1.9e11
mp,dens,6,7800
mp,prxy,6,0.3
keyopt,6,3,0
et,7,beam44
mp,ex,7,3.0e10
mp,dens,7,2600
mp,prxy,7,0.1667
et,8,beam44
mp,ex,8,3.3e10
mp,dens,8,2600
mp,prxy,8,0.1667
et,9,beam44
mp,ex,9,3.0e10
mp,dens,9,2600
mp,prxy,9,0.1667
et,10,link10
mp,ex,10,2.0e11
mp,dens,10,7800
mp,prxy,10,0.3
keyopt,10,3,0
et,11,beam44
mp,ex,11,2.85e10
mp,dens,11,2500
mp,prxy,11,0.1667
n,90000,0,0,30
n,1,0,2.025
n,2,6.25,2.025
n,3,6.25+5.1,2.025
n,4,6.25+5.1+0.2+1.9/2,2.025
n,5,6.25+5.1+2.3,2.025
n,6,6.25+5.1+2.3+2.23,2.025
nsym,x,10,2,6,1
type,4
r,1
real,1
mat,4
secnum,7
e,1,2
secnum,6
e,2,3
secnum,5
e,3,4
e,4,5
secnum,4
e,5,6
secnum,7
e,1,12
secnum,6
e,12,13
secnum,5
e,13,14
e,14,15
secnum,4
e,15,16
egen,2,100,1,10,1,,,,,,0,4.7+1.25/2+1.5/2
egen,17,100,11,20,1,,,,,,0,5.3+1.5
egen,2,100,171,180,1,,,,,,0,4.7+1.5/2+1.25/2
cm,hl,elem
esel,none
type,3
r,2
real,2
mat,3
sectype,10,beam,rect,zl
secdata,0.4,1.2
secnum,10
e,1,101
e,2,102
e,3,103
e,5,105
e,12,112
e,13,113
e,15,115
*do,i,0,15,1
e,101+i*100,201+i*100
e,102+i*100,202+i*100
e,104+i*100,204+i*100
e,112+i*100,212+i*100
e,114+i*100,214+i*100
*enddo
e,1701,1801
e,1702,1802
e,1703,1803
e,1705,1805
e,1712,1812
e,1713,1813
e,1715,1815
cm,zl,elem
esel,none
*dim,y,array,16
*dim,cs,array,16
*dim,z,array,16
y(1)=8.1
y(2)=14.9
y(3)=21.7
y(4)=28.5
y(5)=35.3
y(6)=42.1
y(7)=48.9
y(8)=55.7
y(9)=62.5
y(11)=21.7
y(12)=28.5
y(13)=42.1
y(14)=48.9
y(15)=59.1
y(16)=2.025
f1=25
m1=1.1
k=log(m1+sqrt(m1*m1-1))
*do,i,1,16,1
cs(i)=1-y(i)*2/125
z(i)=25-f1/(m1-1)*(cosh(k*cs(i))-1)
*enddo
*dim,ygd,array,23
ygd(1) =0.0040004 +0.0129
ygd(2) =0.010685 +0.0222
ygd(3) =0.019483 +0.0303
ygd(4) =0.029219 +0.0369
ygd(5) =0.039022 +0.0423
ygd(6) =0.048243 +0.0465
ygd(7) =0.05639 +0.0495
ygd(8) =0.06307 +0.0515
ygd(9) =0.068014 +0.0521
ygd(10) =0.071003 +0.0515
ygd(11) =0.071739 +0.0495
ygd(12) =0.069962 +0.0465
ygd(13) =0.065514 +0.423
ygd(14) =0.058289 +0.0369
ygd(15) =0.048248 +0.0304
ygd(16) =0.035507 +0.0224
ygd(17) =0.20415 +0.0123
ygd(18) =0.06577 +0.040234
ygd(19) =0.069769 +0.041215
ngen,2,10000,104 ,114 ,10,0,0,z(1)+ygd(1 )
ngen,2,10000,204 ,214 ,10,0,0,z(2)+ygd(2 )
ngen,2,10000,304 ,314 ,10,0,0,z(3)+ygd(3 )
ngen,2,10000,404 ,414 ,10,0,0,z(4)+ygd(4 )
ngen,2,10000,504 ,514 ,10,0,0,z(5)+ygd(5 )
ngen,2,10000,604 ,614 ,10,0,0,z(6)+ygd(6 )
ngen,2,10000,704 ,714 ,10,0,0,z(7)+ygd(7 )
ngen,2,10000,804 ,814 ,10,0,0,z(8)+ygd(8 )
ngen,2,10000,904 ,914 ,10,0,0,z(9)+ygd(9 )
ngen,2,10000,1004 ,1014 ,10,0,0,z(8)+ygd(10)
ngen,2,10000,1104 ,1114 ,10,0,0,z(7)+ygd(11)
ngen,2,10000,1204 ,1214 ,10,0,0,z(6)+ygd(12)
ngen,2,10000,1304 ,1314 ,10,0,0,z(5)+ygd(13)
ngen,2,10000,1404 ,1414 ,10,0,0,z(4)+ygd(14)
ngen,2,10000,1504 ,1514 ,10,0,0,z(3)+ygd(15)
ngen,2,10000,1604 ,1614 ,10,0,0,z(2)+ygd(16)
ngen,2,10000,1704 ,1714 ,10,0,0,z(1)+ygd(17)
ngen,2,3000,104 ,114 ,10,0,0,0
ngen,2,3000,204 ,214 ,10,0,0,0
ngen,2,3000,304 ,314 ,10,0,0,0
ngen,2,3000,404 ,414 ,10,0,0,0
ngen,2,3000,504 ,514 ,10,0,0,0
ngen,2,3000,604 ,614 ,10,0,0,0
ngen,2,3000,704 ,714 ,10,0,0,0
ngen,2,3000,804 ,814 ,10,0,0,0
ngen,2,3000,904 ,914 ,10,0,0,0
ngen,2,3000,1004,1014 ,10,0,0,0
ngen,2,3000,1104,1114 ,10,0,0,0
ngen,2,3000,1204,1214 ,10,0,0,0
ngen,2,3000,1304,1314 ,10,0,0,0
ngen,2,3000,1404,1414 ,10,0,0,0
ngen,2,3000,1504,1514 ,10,0,0,0
ngen,2,3000,1604,1614 ,10,0,0,0
ngen,2,3000,1704,1714 ,10,0,0,0
ngen,2,6000,104 ,114 ,10,0,0,z(1)+ygd(1 )
ngen,2,6000,204 ,214 ,10,0,0,z(2)+ygd(2 )
ngen,2,6000,304 ,314 ,10,0,0,z(3)+ygd(3 )
ngen,2,6000,404 ,414 ,10,0,0,z(4)+ygd(4 )
ngen,2,6000,504 ,514 ,10,0,0,z(5)+ygd(5 )
ngen,2,6000,604 ,614 ,10,0,0,z(6)+ygd(6 )
ngen,2,6000,704 ,714 ,10,0,0,z(7)+ygd(7 )
ngen,2,6000,804 ,814 ,10,0,0,z(8)+ygd(8 )
ngen,2,6000,904 ,914 ,10,0,0,z(9)+ygd(9 )
ngen,2,6000,1004,1014 ,10,0,0,z(8)+ygd(10)
ngen,2,6000,1104,1114 ,10,0,0,z(7)+ygd(11)
ngen,2,6000,1204,1214 ,10,0,0,z(6)+ygd(12)
ngen,2,6000,1304,1314 ,10,0,0,z(5)+ygd(13)
ngen,2,6000,1404,1414 ,10,0,0,z(4)+ygd(14)
ngen,2,6000,1504,1514 ,10,0,0,z(3)+ygd(15)
ngen,2,6000,1604,1614 ,10,0,0,z(2)+ygd(16)
ngen,2,6000,1704,1714 ,10,0,0,z(1)+ygd(17)
dgarea=109*3.1415926*0.25*0.007*0.007
dgyl=0.3*1670*1e6
dgl=dgyl*dgarea
edg=1.95e11
*dim,L0,array,17
*dim,L,array,17
*dim,deltL,array,17
*dim,istrn,array,17
*dim,fdg,array,17
*do,i,1,9,1
L0(i)=z(i)
L0(18-i)=z(i)
*enddo
fdg(1 )=dgl
fdg(2 )=dgl
fdg(3 )=dgl
fdg(4 )=dgl
fdg(5 )=dgl
fdg(6 )=dgl
fdg(7 )=dgl
fdg(8 )=dgl
fdg(9 )=dgl
fdg(10)=dgl
fdg(11)=dgl
fdg(12)=dgl
fdg(13)=dgl
fdg(14)=dgl
fdg(15)=dgl
fdg(16)=dgl
fdg(17)=dgl
type,6
mat,6
*do,i,1,17,1
deltL(i)=fdg(i)*L0(i)/(edg*dgarea)
L(i)=L0(i)+deltL(i)
istrn(i)=deltL(i)/L(i)
r,i+1000,dgarea,istrn(i)
*enddo
*do,i,1,17,1
real,1+1000
e,i*100+4+3000,i*100+4+6000
e,i*100+10+4+3000,i*100+10+4+6000
*enddo
cm,dg,elem
esel,none
allsel
real,1018
n,20001, 0 ,21.7 ,z(11)
n,20002, 0 ,28.5 ,z(12)
n,20003, 0 ,42.1 ,z(13)
n,20004, 0 ,48.9 ,z(14)
n,20005, 0 ,59.1 ,z(15)
n,20006, 0 ,65.9 ,z(15)
n,20007, 0 ,76.1 ,z(14)
n,20008, 0 ,82.9 ,z(13)
n,20009, 0 ,96.5 ,z(12)
n,20010, 0 ,103.3 ,z(11)
n,20011, 12.5 ,21.7 ,z(11)
n,20012, 12.5 ,28.5 ,z(12)
n,20013, 12.5 ,42.1 ,z(13)
n,20014, 12.5 ,48.9 ,z(14)
n,20015, 12.5 ,59.1 ,z(15)+ygd(18)
n,20016, 12.5 ,65.9 ,z(15)+ygd(19)
n,20017, 12.5 ,76.1 ,z(14)
n,20018, 12.5 ,82.9 ,z(13)
n,20019, 12.5 ,96.5 ,z(12)
n,20020, 12.5 ,103.3 ,z(11)
n,20021, -12.5 ,21.7 ,z(11)
n,20022, -12.5 ,28.5 ,z(12)
n,20023, -12.5 ,42.1 ,z(13)
n,20024, -12.5 ,48.9 ,z(14)
n,20025, -12.5 ,59.1 ,z(15)+ygd(18)
n,20026, -12.5 ,65.9 ,z(15)+ygd(19)
n,20027, -12.5 ,76.1 ,z(14)
n,20028, -12.5 ,82.9 ,z(13)
n,20029, -12.5 ,96.5 ,z(12)
n,20030, -12.5 ,103.3 ,z(11)
n,20000, 0 ,62.5 ,z(15)
type,5
mat,5
sectype,500,beam,ctube,windc
secdata,0.486,0.50
secnum,500
e,20002,10404,90000
e,20002,10414,90000
e,20002,10304,90000
e,20002,10314,90000
e,20004,10704,90000
e,20004,10714,90000
e,20004,10604,90000
e,20004,10614,90000
e,20000,20015,90000
e,20000,20016,90000
e,20000,20025,90000
e,20000,20026,90000
e,20007,11104,90000
e,20007,11114,90000
e,20007,11204,90000
e,20007,11214,90000
e,20009,11404,90000
e,20009,11414,90000
e,20009,11504,90000
e,20002,11514,90000
cm,fc,elem
esel,none
n,30000,12.5,0,0
n,30001,12.5,125,0
n,30002,12.5,2.025,z(16)
n,30003,12.5,125-2.025,z(16)
n,40000,12.5,0,0
n,40001,12.5,125,0
n,40002,12.5,2.025,z(16)
n,40003,12.5,125-2.025,z(16)
*dim,zgzx,array,23
*dim,ygzx,array,23
zgzx(1 )=30000
zgzx(2 )=30002
zgzx(3 )=10104
zgzx(4 )=10204
zgzx(5 )=10304
zgzx(6 )=10404
zgzx(7 )=10504
zgzx(8 )=10604
zgzx(9 )=10704
zgzx(10 )=10804
zgzx(11 )=20015
zgzx(12 )=10904
zgzx(13 )=20016
zgzx(14 )=11004
zgzx(15 )=11104
zgzx(16 )=11204
zgzx(17 )=11304
zgzx(18 )=11404
zgzx(19 )=11504
zgzx(20 )=11604
zgzx(21 )=11704
zgzx(22 )=30003
zgzx(23 )=30001
zgzx(1 )=40000
zgzx(2 )=40002
zgzx(3 )=10114
zgzx(4 )=10214
zgzx(5 )=10314
zgzx(6 )=10414
zgzx(7 )=10514
zgzx(8 )=10614
zgzx(9 )=10714
zgzx(10 )=10814
zgzx(11 )=20025
zgzx(12 )=10914
zgzx(13 )=20026
zgzx(14 )=11014
zgzx(15 )=11114
zgzx(16 )=11214
zgzx(17 )=11314
zgzx(18 )=11414
zgzx(19 )=11514
zgzx(20 )=11614
zgzx(21 )=11714
zgzx(22 )=40003
zgzx(23 )=40001
n,90001,12.5,0,30
n,90002,-12.5,0,30
type,1
mat,1
r,2000
secnum,1
*do,i,1,22
e,zgzx(i),zgzx(i+1),90001
e,ygzx(i),ygzx(i+1),90002
*enddo
cm,gzgg,elem
esel,none
type,2
mat,2
r,2001
secnum,2
*do,i,1,22
e,zgzx(i),zgzx(i+1),90001
e,ygzx(i),ygzx(i+1),90002
*enddo
cm,gzhnt,elem
esel,none
type,9
r,6
real,6
mat,9
sectype,70,beam,rect,dizuo
secdata,1.4,4
secnum,70
e,30001,30003,90001
e,40001,40003,90002
e,30000,30002,90001
e,30000,40002,90002
cm,gzdz,elem
esel,none
ngen,2,2000,1,16,1,0,-2.025,0
ngen,2,100,1801,1816,1,0,2.025,0
type,7
mat,7
r,7
real,7
sectype,60,beam,rect,gail
secdata,4,1.6
!!!!!!!!!!*****************************以上已经校核过!
secnum,60
n,90005,0,0,20
n,90006,0,125,20
e,1901, 1902, 90006
e,1902, 1903, 90006
e,1903, 1904, 90006
e,1904, 1905, 90006
e,1905, 1906, 90006
e,1906, 1907, 90006
e,1907, 1908, 90006
e,1908, 1909, 90006
e,1909, 1910, 90006
e,1910, 1911, 90006
e,1911, 1912, 90006
e,1912, 1913, 90006
e,1913, 1914, 90006
e,1914, 1915, 90006
e,1915, 1916, 90006
e,2001, 2002, 90005
e,2002, 2003, 90005
e,2003, 2004, 90005
e,2004, 2005, 90005
e,2005, 2006, 90005
e,2006, 2007, 90005
e,2007, 2008, 90005
e,2008, 2009, 90005
e,2009, 2010, 90005
e,2010, 2011, 90005
e,2011, 2012, 90005
e,2012, 2013, 90005
e,2013, 2014, 90005
e,2014, 2015, 90005
e,2015, 2016, 90005
cm,gailiang,elem
esel,none
n,90003,14.5,0,0
n,90004,-14,5,0,0
type,8
r,8
real,8
sectype,50,beam,hrec,dun
secdata,3,4,0.8,0.8,0.8,0.8
ngen,2,10000,2004,2014,10,0,0,-10
ngen,2,10000,1904,1914,10,0,0,-10
secnum,50
e,2004,12004,90003
e,2014,12014,90004
e,1904,11904,90003
e,1914,11914,90004
cm,dunzhu,elem
esel,none
type,11
mat,11
r,1
sectype,40,beam,rect,chengtai
secdata,12.7,9
r,11
real,11
secnum,40
ngen,2,20000,2004,2014,10,0,0,-11.5
ngen,2,20000,1904,1914,10,0,0,-11.5
ngen,2,30000,2004,2014,10,0,0,-13
ngen,2,30000,1904,1914,10,0,0,-13
e,12004,22004,90003
e,11904,21904,90003
e,12014,22014,90004
e,11904,21904,90004
e,22004,32004,90003
e,21904,31904,90003
e,22014,32014,90004
e,21904,31904,90004
cm,chtai,elem
esel,none
type,9
mat,9
r,9
real,9
n,90007,12.5,0,10
n,90008,-12.5,0,10
cm,hxl,elem
esel,none
n,50000,6.15+1.85,2.7,-13
n,50001,6.15+1.85,-2.7,-13
n,50002,6.15+1.85+4.5,2.7,-13
n,50003,6.15+1.85+4.5,-2.7,-13
n,50004,6.15+1.27-1.85,2.7,-13
n,50005,6.15+1.27-1.85,-2.7,-13
nsel,s,,,50000,50005,1
nsym,x,10,50000,50005,1
allsel
ngen,2,100,50000,50015,1,0,125,-2.5
ngen,10,1000,50000,50015,1,0,0,-2.5
esel,none
type,11
mat,11
r,11
real,11
sectype,20,beam,csolid,pier
secdata,0.9
secnum,20
*do,i,0,8,1
*do,k,0,5,1
e,50000+k+i*1000,50000+k+(i+1)*1000,90003
e,50000+k+10+i*1000,50000+k+10+(i+1)*1000,90003
e,50000+k+100+i*1000,50000+k+100+(i+1)*1000,90003
e,50000+k+110+i*1000,50000+k+110+(i+1)*1000,90003
*enddo
*enddo
cm,zhuang,elem
esel,none
d=1.8
kf=0.9
k0=1+1/d
b1=0.6
L1=3.6
h1=3*(d+1)
k=b1+(1-b1)/0.6*L1/h1
B=Kf*k0*k*d
*dim,m,array,5
m(1)=8
m(2)=15
m(3)=20
m(4)=25
m(5)=4000
*dim,h,array,5
h(1)=5
h(2)=5
h(3)=5
h(4)=5
h(5)=10
*dim,c,array,4
c(1)=m(1)*h(1)
c(2)=m(2)*(h(1)+h(2))
c(3)=m(3)*(h(1)+h(2)+h(3))
c(4)=m(4)*(h(1)+h(2)+h(3)+h(4))
*dim,e,array,4
e(1)=0.5*c(1)*h(1)*B*1e6
*do,i,2,4,1
e(i)=0.5*(c(i-1)+c(i))*h(1)*B*1e6
*enddo
et,20,combin14
r,20,e(1)
r,30,e(2)
r,40,e(3)
r,50,e(4)
r,60,m(5)*3.1415926*B*B/4*1e6
ngen,2,10000,51000,51115,1,1,0,0
ngen,2,20000,51000,51115,1,0,1,0
ngen,2,30000,51000,51115,1,0,0,1
type,20
real,20
*do,i,0,5,1
e,51000+i,51000+i+10000
e,51000+i+10,51000+i+10+10000
e,51100+i,51100+i+10000
e,51100+i+10,51100+i+10+10000
e,51000+i,51000+i+20000
e,51000+i+10,51000+i+10+20000
e,51100+i,51100+i+20000
e,51100+i+10,51100+i+10+20000
*enddo
cm,spring1,elem
esel,none
ngen,2,10000,53000,53115,1,1,0,0
ngen,2,20000,53000,53115,1,0,1,0
ngen,2,30000,53000,53115,1,0,0,1
type,20
real,30
*do,i,0,5,1
e,53000+i,53000+i+10000
e,53000+i+10,53000+i+10+10000
e,53100+i,53100+i+10000
e,53100+i+10,53100+i+10+10000
e,53000+i,53000+i+20000
e,53000+i+10,53000+i+10+20000
e,53100+i,53100+i+20000
e,53100+i+10,53100+i+10+20000
*enddo
cm,spring2,elem
esel,none
ngen,2,10000,55000,55115,1,1,0,0
ngen,2,20000,55000,55115,1,0,1,0
type,20
real,40
*do,i,0,5,1
e,55000+i,55000+i+10000
e,55000+i+10,55000+i+10+10000
e,55100+i,55100+i+10000
e,55100+i+10,55100+i+10+10000
e,55000+i,55000+i+20000
e,55000+i+10,55000+i+10+20000
e,55100+i,55100+i+20000
e,55100+i+10,55100+i+10+20000
*enddo
cm,spring3,elem
esel,none
ngen,2,10000,57000,57115,1,1,0,0
ngen,2,20000,57000,57115,1,0,1,0
未完...... 续上....
type,20
real,50
*do,i,0,5,1
e,57000+i,57000+i+10000
e,57000+i+10,57000+i+10+10000
e,57100+i,57100+i+10000
e,57100+i+10,57100+i+10+10000
e,57000+i,57000+i+20000
e,57000+i+10,57000+i+10+20000
e,57100+i,57100+i+20000
e,57100+i+10,57100+i+10+20000
*enddo
cm,spring4,elem
esel,none
ngen,2,10000,58000,58115,1,1,0,0
ngen,2,20000,58000,58115,1,0,1,0
type,20
real,60
*do,i,0,5,1
e,58000+i,58000+i+10000
e,58000+i+10,58000+i+10+10000
e,58100+i,58100+i+10000
e,58100+i+10,58100+i+10+10000
e,58000+i,58000+i+20000
e,58000+i+10,58000+i+10+20000
e,58100+i,58100+i+20000
e,58100+i+10,58100+i+10+20000
*enddo
cm,spring5,elem
esel,none
ngen,2,10000,59000,59115,1,1,0,0
ngen,2,10000,59000,59115,1,1,0,0
ngen,2,20000,59000,59115,1,0,1,0
ngen,2,20000,59000,59115,1,0,1,0
ngen,2,30000,59000,59115,1,0,0,-1
ngen,2,30000,59000,59115,1,0,0,-1
type,20
real,60
*do,i,0,5,1
e,59000+i,59000+i+10000
e,59000+i+10,59000+i+10+10000
e,59100+i,59100+i+10000
e,59100+i+10,59100+i+10+10000
e,59000+i,59000+i+20000
e,59000+i+10,59000+i+10+20000
e,59100+i,59100+i+20000
e,59100+i+10,59100+i+10+20000
e,59000+i,59000+i+30000
e,59000+i+10,59000+i+10+30000
e,59100+i,59100+i+30000
e,59100+i+10,59100+i+10+30000
*enddo
cm,spring6,elem
allsel
allsel
nsel,s,loc,z,-50,-0.01
nsel,r,loc,y,-10,10
cm,zdnode,node
esln,s,,all
cm,zdelem,elem
allsel
nsel,u,,,zdnode
nsym,y,100000,all
allsel
esel,u,,,zdelem
esym,,100000,all
allsel
xgarea=8*27*7*3.1415926*0.25*0.005*0.005
exg=1.95e11
*dim,fxg,array,2
*dim,Llast,array,2
*dim,Listrn,array,2
*dim,detL,array,2
xgyl=0.3*1860*1e6
xgl=xgyl*xgarea
fxg(1)=xgl
fxg(2)=xgl
Lzero=250
*do,i,1,2,1
detL(i)=Lzero*Fxg(i)/(exg*xgarea)
Llast(i)=Lzero+detL(i)
Listrn(i)=detl(i)/Llast(i)
r,i+100,xgarea,Listrn(i)
*enddo
type,10
mat,10
real,101
e,30001,130001
real,102
e,40001,140001
cm,xigan,elem
esel,none
allesl
cp, 1, ux ,2004 ,30000
cp, 1, uy ,2004 ,30000
cp, next, uz ,2004 ,30000
cp, next, roty ,2004 ,30000
cp, next, rotz ,2004 ,30000
cp, next, ux ,1904 ,30001
cp, next, uy ,1904 ,30001
cp, next, uz ,1904 ,30001
cp, next, roty ,1904 ,30001
cp, next, rotz ,1904 ,30001
cp, next, ux ,2014 ,40000
cp, next, uy ,2014 ,40000
cp, next, uz ,2014 ,40000
cp, next, roty ,2014 ,40000
cp, next, rotz ,2014 ,40000
cp, next, ux ,1914 ,40001
cp, next, uy ,1914 ,40001
cp, next, uz ,1914 ,40001
cp, next, roty ,1914 ,40001
cp, next, rotz ,1914 ,40001
cp, next, ux ,101914 ,140001
cp, next, uy ,101914 ,140001
cp, next, uz ,101914 ,140001
cp, next, roty ,101914 ,140001
cp, next, rotz ,101914 ,140001
cp, next, ux ,101904 ,130001
cp, next, uy ,101904 ,130001
cp, next, uz ,101904 ,130001
cp, next, roty ,101904 ,130001
cp, next, rotz ,101904 ,130001
allsel
nsel,s,,,30000,130000,100000
nsel,a,,,40000,140000,100000
nummrg,node
allsel
*do,i,1,3,1
cp,next,all,1800+i,1900+i
*enddo
allsel
*do,i,5,6,1
cp,next,all,1800+i,1900+i
*enddo
allsel
*do,i,1,2,1
cp,next,all,1800+i+11,1900+i+11
cp,next,all,i+11,2000+i+11,i+100000+11
*enddo
allsel
*do,i,4,5,1
cp,next,all,1800+i+11,1900+i+11
cp,next,all,i+11,2000+i+11,i+100000+11
*enddo
allsel
*do,i,1,3,1
cp,next,all,i,2000+i,100000+i
*enddo
allsel
*do,i,5,6,1
cp,next,all,i,2000+i,100000+i
*enddo
allsel
*do,i,1,3,1
cp,next,all,101800+i,101900+i
*enddo
allsel
*do,i,5,6,1
cp,next,all,101800+i,101900+i
*enddo
allsel
*do,i,1,2,1
cp,next,all,101800+i+11,101900+i+11
*enddo
allsel
*do,i,4,5,1
cp,next,all,101800+i+11,101900+i+11
*enddo
allsel
nsel,s,loc,z,-13
nsel,r,loc,y,-3,3
nsel,r,loc,x,8,20
cp,next,all,all
allsel
nsel,s,loc,z,-13
nsel,r,loc,y,-3,3
nsel,r,loc,x,-8,-20
cp,next,all,all
allsel
nsel,s,loc,z,-13
nsel,r,loc,y,122,128
nsel,r,loc,x,-8,-20
cp,next,all,all
allsel
nsel,s,loc,z,-13
nsel,r,loc,y,122,128
nsel,r,loc,x,8,20
cp,next,all,all
allsel
nsel,s,loc,z,-13
nsel,r,loc,y,-122,-128
nsel,r,loc,x,8,20
cp,next,all,all
allsel
nsel,s,loc,z,-13
nsel,r,loc,y,-122,-128
nsel,r,loc,x,-8,-20
cp,next,all,all
allsel
*do,i,1,17,1
cp, next, ux, i*100+4 ,i*100+4+3000
cp, next, ux, i*100+4+6000 ,i*100+4+10000
cp, next, uy, i*100+4 ,i*100+4+3000
cp, next, uy, i*100+4+6000 ,i*100+4+10000
cp, next, uz, i*100+4 , i*100+4+3000
cp, next, uz, i*100+4+6000 ,i*100+4+10000
cp, next, ux, i*100+4+10 ,i*100+4+3000+10
cp, next, ux, i*100+4+6000+10 ,i*100+4+10000+10
cp, next, uy, i*100+4+10 ,i*100+4+3000+10
cp, next, uy, i*100+4+6000+10 ,i*100+4+10000+10
cp, next, uz, i*100+4+10 ,i*100+4+3000+10
cp, next, uz, i*100+4+6000+10 ,i*100+4+10000+10
cp, next, ux, i*100+4+100000 ,i*100+4+3000+100000
cp, next, ux, i*100+4+6000+100000 ,i*100+4+10000+100000
cp, next, uy, i*100+4+100000 , i*100+4+3000+100000
cp, next, uy, i*100+4+6000+100000 ,i*100+4+10000+100000
cp, next, uz, i*100+4+100000 ,i*100+4+3000+100000
cp, next, uz, i*100+4+6000+100000 ,i*100+4+10000+100000
cp, next, ux, i*100+4+100000+10 ,i*100+4+3000+100000+10
cp, next, ux, i*100+4+6000+100000 ,i*100+4+10000+100000+10
cp, next, uy, i*100+4+100000+10 ,i*100+4+3000+100000+10
cp, next, uy, i*100+4+6000+100000 ,i*100+4+10000+100000+10
cp, next, uz, i*100+4+100000+10 ,i*100+4+3000+100000+10
cp, next, uz, i*100+4+6000+100000 ,i*100+4+10000+100000+10
*enddo
allsel
allsel
nsel,u,,,erase
nsel,s,loc,x,13.5
nsel,a,loc,x,9
nsel,a,loc,x,18
nsel,a,loc,x,-11.5
nsel,a,loc,x,-7
nsel,a,loc,x,-16
nsel,a,loc,y,3.7
nsel,a,loc,y,-1.7
nsel,a,loc,y,128.7
nsel,a,loc,y,123.3
nsel,a,loc,y,-128.7
nsel,a,loc,y,-123.3
nsel,a,loc,z,-40,-13
d,all,all
allsel
nsel,s,loc,z,-36,5
d,all,all
/color,pbak,off
/replot
allsel
nsel,s,loc,y,0
nsel,r,loc,z,0
nummrg,node,all
esel,none
allsel
/solu
allsel
mpdele,dens,all
mp,dens,3,2600
mp,dens,4,2600
acel,,,9.8
pstres,on
solve
save
/solu
allsel
mpdele,dens,3,11,1
mp,dens,1,7800
mp,dens,2,2600
*do,i,1,17,1
r,1000+i,dgarea,0
*enddo
r,101,xgarea,0
r,102,xgarea,0
acel,,,9.8
solve
/solu
allsel
*do,i,1,17,1
r,1000+1,dgarea,0
*enddo
r,101,xgarea,0
r,102,xgarea,0
r,101,xgarea,0
r,102,xgarea,0
qmb=(250/1.5)*32/6.8*5*1000*9.8
acel,,,9.8
mpdele,dens,3,4,1
mp,dens,1,7800
mp,dens,2,2600
mp,dens,3,4167.7
mp,dens,4,4167.7
mp,dens,5,7800
mp,dens,6,7800
mp,dens,7,2600
mp,dens,8,2600
mp,dens,9,2600
mp,dens,10,7800
mp,dens,11,2600
solve
/solu
allsel
acel,,,9.8
r,101,xgarea,Listrn(1)*1.3
r,102,xgarea,Listrn(2)*1.3
pstres,on
solve
/solu
allsel
r, 1001, dgarea,1.6787e-03+7.64e-05
r, 1002, dgarea,1.7538e-03+8.48e-07
r, 1003, dgarea,1.7547e-03+1.03e-06
r, 1004, dgarea,1.7400e-03+2.71e-06
r, 1005, dgarea,1.7292e-03+6.17e-06
r, 1006, dgarea,1.7303e-03+6.37e-06
r, 1007, dgarea,1.7267e-03+6.47e-06
r, 1008, dgarea,1.7279e-03+7.69e-06
r, 1009, dgarea,1.7365e-03+7.91e-06
r, 1010, dgarea,1.7277e-03+7.81e-06
r, 1011, dgarea,1.7263e-03+6.56e-06
r, 1012, dgarea,1.7299e-03+6.24e-06
r, 1013, dgarea,1.7290e-03+5.53e-06
r, 1014, dgarea,1.7412e-03+2.03e-06
r, 1015, dgarea,1.7504e-03+3.03e-06
r, 1016, dgarea,1.7587e-03+1.28e-05
r, 1017, dgarea,1.6143e-03+4.90e-05
pstres,on
solve
finish
/solu
antype,buckle
bucopt,subsp,4
outpr,nsol,all
save
solve
finish
/solu
expass,on
mxpand,4,,,
outpr,all
outres
solve
finish
save
finish
代码提供者：maple520 某铁路钢管混凝土拱桥
/prep7
/title,analysis of a cftb
!主梁节点
/prep7
z=-3.0124 !(=-6.0248/2)
y=-22.9207! (=-45.3653/2-0.2380)
b=0
n,b+1,-84,y,z$n,b+2,-69,y,z
n,b+3,-57,y,z $n,b+4,-45,y,z
*do,i,b+5,b+14,1
n,i,-45+9*(i-4),y,z
*enddo
n,b+15,56,y,z
ngen,2,15,all,,,,,5.72
!主梁左侧节点由对称产生，这里不写
!吊杆底部节点放在后面
！拱肋关键点
c=30!c的值可变
n,c+1,0,1.3454,-0.2130 $n,c+26,0,-1.3454,0.2130
n,c+2,-3,1.2736,-0.2016 $n,c+27,-3,-1.4182,0.2245
n,c+3,-6,1.0591,-0.1677 $n,c+28,-6,-1.6374,0.2592
n,c+4,-9,0.7014,-0.1110 $n,c+29,-9,-2.0027,0.3170
n,c+5,-12,0.2006,-0.0316$n,c+30,-12,-2.5146,0.3981
n,c+6,-15,-0.4438,0.0702 $n,c+31,-15,-3.1728,0.5023
n,c+7,-18,-1.2311,0.1949$n,c+32,-18,-3.9785,0.6298
n,c+8,-21,-2.1621,0.3423$n,c+33,-21,-4.9315,0.7807
n,c+9,-24,-3.2372,0.5125$n,c+34,-24,-6.0321,0.9549
n,c+10,-27,-4.4564,0.7055$n,c+35,-27,-7.2811,1.1526
n,c+11,-30,-5.8205,0.9214$n,c+36,-30,-8.6785,1.3738
n,c+12,-33,-7.3287,1.1602$n,c+37,-33,-10.2262,1.6188
n,c+13,-36,-8.9825,1.4219$n,c+38,-36,-11.9238,1.8876
n,c+14,-39,-10.7820,1.7068$n,c+39,-39,-13.7722,2.1802
n,c+15,-42,-12.7303,2.0152$n,c+40,-42,-15.7699,2.4964
n,c+16,-45,-14.8188,2.3459$n,c+41,-45,-17.9269,2.8379
n,c+17,-48,-17.0570,2.7002$n,c+42,-48,-20.2351,3.2033
n,c+18,-51,-19.4426,3.0778$n,c+43,-51,-22.6984,3.5932
n,c+19,-54,-21.9766,3.4790$n,c+44,-54,-25.3178,4.0079
n,c+20,-57,-24.6590,3.9036$n,c+45,-57,-28.0949,4.4475
n,c+21,-60,-27.4903,4.3518$n,c+46,-60,-31.0313,4.9123
n,c+22,-63,-30.4717,4.8238$n,c+47,-63,-34.1279,5.4025
n,c+23,-66,-33.6029,5.3194$n,c+48,-66,-37.3871,5.9185
n,c+24,-69,-36.8825,5.8386$n,c+49,-69,-40.8127,6.4608
n,c+25,-71,-39.1646,6.1998
nsel,s,,,31,79,1
ngen,2,49,all,,,,0.2345/2,1.4815/2
!（49＝79+1-31）
nsel,s,,,80,128,1
ngen,2,49,all,,,,0.2345/2,1.4815/2
nsel,s,,,80,128,1
ngen,2,49,all,,,,0.2345/2,1.4815/2
nsel,all
nsel,s,,,32,49,1
nsym,x,146 ,all,,,
nsel,s,,,57,74,1
nsym,x,139 ,all,,,
nsel,s,,,81,98,1
nsym,x,133 ,all,,,
nsel,s,,,106,123,1
nsym,x,126 ,all,,,
nsel,s,,,130,147,1
nsym,x,120 ,all,,,
nsel,s,,,155,172,1
nsym,x,113,all,,,
nsel,all
nplot
!补右拱脚节点
n,286,55,-26.221,4.1509
n,287,55,-26.1038,4.8917
n,288,55,-25.9865,5.6324
n,289,57,-24.6590,3.9036
n,290,57,-24.4245,4.64435
n,291,57,-24.4245,5.3851
！横撑中心节点
nsel,s,,,5
ngen,2,287,all,,,0,10.9969,0
nsel,s,,,5
ngen,2,288,all,,,0,13.9382,0
nsel,s,,,7
ngen,2,287,all,,,0,18.9422,0
nsel,s,,,7
ngen,2,288,all,,,0,21.6896,0
nsel,s,,,9
ngen,2,287,all,,,0,21.5753,0
nsel,s,,,9
ngen,2,288,all,,,0,24.2661,0
nsel,s,,,11
ngen,2,287,all,,,0,18.9422,0
nsel,s,,,11
ngen,2,288,all,,,0,21.6896,0
nsel,s,,,13
ngen,2,287,all,,,0,10.9969,0
nsel,s,,,13
ngen,2,288,all,,,0,13.9382,0
nsel,all
nplot
*do,i,1,11,1
n,301+i,-45+9*(i-1),-22.6827,2.6876
*enddo
*do,i,1,6,1
nsel,s,,,105+3*(i-1)
ngen,2,208-2*(i-1),all,,,0,-0.4981,0.01358
*enddo
*do,i,1,5,1
nsel,s,,,234+3*(i-1)
ngen,2,85-2*(i-1),all,,,0,-0.4981,0.01358
*enddo
nsel,all
nplot
!立柱
!324，2要耦合
!325，3要耦合
!326，17要耦合
!327，18要耦合
nsel,s,,,2,3,1
ngen,2,322,all,,,0,-(1.6-0.238),0
nsel,s,,,17,18,1
ngen,2,309,all,,,0,-(1.6-0.238),0
nsel,s,,,326
ngen,2,2,all,,,0,0,3.131
nsel,s,,,327
ngen,2,2,all,,,0,0,1.196
nsel,all
nplot
!定义单元类型
et,1,beam4
et,2,link10
keyopt,2,2,2
keyopt,2,3,0
et,3,shell63
!定义材料
!拱肋弦杆
MP,EX,1,0.527e11
MP,nuxy,1,0.25
MP,dens,1,3188
!拱肋腹杆，横撑强、弱构件
MP,EX,2,2.06e11
MP,nuxy,2,0.3
MP,dens,2,7850
!拱肋上下平联
MP,EX,3,0.355e11
MP,nuxy,3,0.2
MP,dens,3,2600
!吊杆
MP,EX,4,1.95e11
MP,nuxy,4,0.3
MP,dens,4,7850
!立柱、立柱横梁
MP,EX,5,0.3e11
MP,nuxy,5,0.2
MP,dens,5,2500
!鱼骨梁主梁
MP,EX,6,0.355e11
MP,nuxy,6,0.2
MP,dens,6,4119
!鱼骨梁刺梁
MP,EX,7,3.55e11 !弹性模量输入主梁的十倍
MP,nuxy,7,0.2
MP,dens,7,0.1
!定义实常数（面积除了拱肋弦杆外，其他输入实际面积以和密度对应）
!拱肋弦杆
r,1,0.37,0.03956,0.03956,0.7,0.7
rmore,,0.02357
!拱肋腹杆、横撑强构件
r,2,0.01225,0.000233,0.000233,0.4,0.4
!拱肋横撑弱构件
r,3,0.004825,0.0000233,0.0000233,0.2,0.2
!立柱
r,4,1,0.083333,0.083333,1,1
rmore,,0.1426
!立柱横梁
r,5,1.44,0.1728,0.1728,1.2,1.2
rmore,,0.2957
!鱼骨梁主梁
r,6,6.718,2.061,103.4,2.5,1.6!注意输入
rmore,0.02357
!鱼骨梁刺梁
r,7,0.01,1030.4,20.016,0.1,0.1 !抗弯惯性矩为主梁的十倍
!吊杆
r,8,0.00535,1e-8
!上下平联板
r,9,0.5，0.5，0.5，0.5，
!主梁
！刺梁
type,1
mat,7
real,7
*do,i,1,15,1
e,i,i+15
*enddo
！主梁
type,1
mat,6
real,6
*do,i,1,14,1
e,i,i+1
!e,e+15,i+16!三主梁的边梁
*enddo
!拱肋
！弦杆
type,1
mat,1
real,1
*do,i,1,24,1
e,i+30,i+31
e,i+128,i+129
*enddo
*do,i,1,23,1
e,i+55,i+56
e,i+153,i+154
*enddo
*do,i,1,17,1
e,i+177,i+178
e,i+195,i+196
e,i+249,i+250
e,i+267,i+268
*enddo
e,31,178$ e,129,250
e,56,196$ e,154,268
e,213,286$ e,285,288
e,195,289$ e,267,291
！竖腹杆
type,1
mat,2
real,2
*do,i,1,24,1
e,i+30,i+55
e,i+128,i+153
*enddo
*do,i,1,18,1
e,i+177,i+195
e,i+249,i+267
*enddo
！斜腹杆
type,1
mat,2
real,2
*do,i,1,23,1
e,i+31,i+55
e,i+129,i+153
*enddo
*do,i,1,17,1
e,i+178,i+195
e,i+250,i+267
*enddo
e,178,56$e,250,154
e,213,289$ e,285,291
!横撑
！横撑强构件
type,1
mat,2
real,2
e,292,293$ e,292,43
e,294,295$ e,294,37
e,296,297$e,296,31
e,298,299$e,298,183
e,300,301$ e,300,189
！横撑弱构件
type,1
mat,2
real,3
e,293,43 $ e,292,68
e,295,37 $e,294,62
e,297,31$e,296,56
e,299,183$ e,298,201
e,301,189$ e,300,207
!吊杆
type,2
mat,4
real,8
e,302,318$ e,303,317$ e,304,316
e,305,315$ e,306,314$ e,307,313
e,308,319$ e,309,320$ e,310,321
e,311,322$ e,312,323
e,95,120$ e,92,117$ e,89,114
e,86,111$ e,83,108$ e,80,105
e,216,234$ e,219,237$ e,222,240
e,225,243$ e,228,246
!立柱
！立柱横梁
type,1
mat,5
real,5
e,324,326$ e,326,328
e,325,327$ e,327,329
！竖立柱
type,1
mat,5
real,4
e,328,54 $ e,329,50
!上下平联板
type,3
mat,3
real,9
！左平联
*do,i,1,24,1
e,i+30,i+31,i+80,i+79
e,i+79,i+80,i+129,i+128
*enddo
*do,i,1,23,1
e,i+55,i+56,i+105,i+104
e,i+104,i+105,i+154,i+153
*enddo
e,31,178,214,80$ e,80,214,250,129
e,56,196,232,105$ e,105,232,268,154
！右平联
*do,i,1,17,1
e,i+177,i+178,i+214,i+213
e,i+213,i+214,i+250,i+249
e,i+195,i+196,i+232,i+231
e,i+231,i+232,i+268,i+267
*enddo
e,195,289,290,231$ e,231,290,291,267
e,213,286,287,249$ e,249,287,288,285
/eshape,1
eplot
cs,11,0,297,296,298,1,1
csys,11
dsys,11
nsym,z,329,all!前面一共three two nine 个节点
ensym,585,,329,all !前面一共five eight five个单元
eplot
csys,0
dsys,0
/replo
nummrg,all, , , ,low
numcmp,all
*do,i,1,6,1
cp,next,ux,312+i,102+3*i
cp,next,uy,312+i,102+3*i
cp,next,uz,312+i,102+3*i
*enddo
*do,i,1,5,1
cp,next,ux,318+i,231+3*i
cp,next,uy,318+i,231+3*i
cp,next,uz,318+i,231+3*i
*enddo
*do,i,1,11,1
cp,next,ux,301+i,18+i
cp,next,uy,301+i,18+i
cp,next,uz,301+i,18+i
*enddo
*do,i,1,6,1
cp,next,ux,616+i,416+3*i
cp,next,uy,616+i,416+3*i
cp,next,uz,616+i,416+3*i
*enddo
*do,i,1,5,1
cp,next,ux,622+i,545+3*i
cp,next,uy,622+i,545+3*i
cp,next,uz,622+i,545+3*i
*enddo
*do,i,1,11,1
cp,next,ux,605+i,332+i
cp,next,uy,605+i,332+i
cp,next,uz,605+i,332+i
*enddo
*do,i,1,2,1
cp,next,ux,16+i,325+i
cp,next,uy,16+i,325+i
cp,next,uz,16+i,325+i
cp,next,rotx,16+i,325+i
cp,next,roty,16+i,325+i
cp,next,rotz,16+i,325+i
*enddo
*do,i,1,2,1
cp,next,ux,1+i,323+i
cp,next,uy,1+i,323+i
cp,next,uz,1+i,323+i
cp,next,rotx,1+i,323+i
cp,next,roty,1+i,323+i
cp,next,rotz,1+i,323+i
*enddo
*do,i,1,2,1
cp,next,ux,330+i,627+i
cp,next,uy,330+i,627+i
cp,next,uz,330+i,627+i
cp,next,rotx,330+i,627+i
cp,next,roty,330+i,627+i
cp,next,rotz,330+i,627+i
*enddo
/solu
antype,2
modopt,lanb,10
mxpand,10, , ,0
lumpm,0
pstres,0
modopt,lanb,10,0,100, ,off
d,286, , , , 288,2 ,all
d,289, , , , 291,2 ,all
d,600, , , , 602,2 ,all
d,603, , , , 605,2 ,all
d,467, , , , , ,all
d,369, , , , , ,all
d,491, , , , , ,all
d,393, , , , , ,all
d,55, , , , , ,all
d,153, , , , , ,all
d,79, , , , , ,all
d,177, , , , , ,all
d,344,,,,,,ux,uy,uz
d,30,,,,,,ux,uy,uz
d,15,,,,,,ux,uy,uz
d,330,,,,,,uy,uz
d,1,,,,,,uy,uz
d,16,,,,,,uy,uz
solve
finish
代码提供者：西王赏功 回复 5 # AaronSpark 的帖子
这两个都是钢管混凝土拱桥的例子吧。。如果想要分别用壳、实体单元去模拟钢管混凝土柱。。。应该怎么模拟?特别是他们之间的接触应该怎么去设置？望高手不吝解答。。。耐心等待中。。感谢！ 看不懂，好复杂啊
页:
[1]