[转帖]钢管混凝土拱桥建模实例
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
!!!!!!!!!!*****************************以上已经校核过!
回复：（AaronSpark）[转帖]钢管混凝土拱桥建模实例
secnum,60n,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
谢谢楼主，请问楼主有没有钢管相贯节点的命令流实例，谢谢 命令流没有，gui得到是有一个，看帖子
http://forum.vibunion.com/forum/viewthread.php?tid=17479 谢谢楼主 谢谢！！非常
页:
[1]