安庆长江铁路大桥(ANSYS-APDL建模)

/桥址概况ﻫ安庆长江铁路大桥是南京至安庆城际铁路和阜阳至景德镇铁路的重要组成部分,位于安庆前江口汇合口处下游官山咀附近，距上游已建成通车的安庆长江公路大桥约21km;线路在池州侧晏塘镇靠近长江的刘村附近右拐过江，过江后从安庆的长风镇穿过。

安庆铁路长江大桥全长2996。8m,其中主桥采用跨度为1０1．5＋１8８．５+５80+217．5＋159．5+116m的钢桁梁斜拉桥;非通航孔正桥采用6孔跨径64m预应力混凝土简支箱梁;东引桥采用１6孔梁长3２。6ｍ预应力混凝土简支箱梁；跨大堤桥采用48．９+8６+48。８m预应力混凝土连续箱梁;西引桥采用１５孔梁长3２。6ｍ预应力混凝土简支梁及2孔梁长2４。６m预应力混凝土简支梁，其中宁安线采用箱梁，阜景线采用T梁.

ﻫ主桥桥式及桥型特点ﻫ主桥采用10３+１88．5+580＋21７。5＋15９．5+117.5m两塔钢桁斜拉桥方案，全长13６6m.主梁为三片主桁钢桁梁,桁间距2x1４m,节间长1４。５m，桁高１5m。主塔为钢筋混凝土结构，塔顶高程+204.00m,塔底高程－6．00m，斜拉索为空间三索面,立面上每塔两侧共18对索,全桥２1６根斜拉索。所有桥墩上均设竖向和横向约束,４#塔与主梁之间设纵向水平约束,3＃塔与梁间使用带限位功能的粘滞阻尼器。主梁为”N"字型桁式，横向采用三片桁结构,主桁的横向中心距各为１4m,桁高15m，节间距14.5m［２].

ﻫ结构构造

主桥采用两塔钢桁斜拉桥方案，主梁为三片主桁钢桁梁，主桁上下弦杆均为箱型截面,上弦杆内高10０0mm，内宽１２０0ｍm,板厚２０～48mm。下弦杆内高14０0ｍm,宽1200ｍm,板厚20~５6mｍ。下弦杆顶板向桁内侧加宽７０0mm与整体桥面板焊接。腹杆主要采用H型截面。H型杆件宽120０mm，高７２0和760mm，板厚20～４８ｍm。根据不同的受力区段选用不同的杆件截面，在辅助墩附近的压重区梁段,腹杆采用箱型截面杆件。主桁采用焊接杆件,整体节点.在节点外以高强度螺栓拼接的结构形式,上下弦杆四面等强对接拼装。H型腹杆采用插入式连接。箱型腹杆采用四面与主桁节点对拼的连接形式。主桁拼接采用Ｍ30高强螺栓。

项目进展ﻫ２０05年元月，安庆长江铁路大桥项目前期工作协调领导小组办公室委托铁道第四勘察设计院编制《安徽省铁路总体规划安庆过江通道深化研究及大桥选址报告》,随后铁四院专家组来宜现场勘察,采集相关资料,并于2月份完成该报告。ﻫ200５年８月，安徽省发改委主持召开“安庆长江铁路大桥桥位专家咨询会”，邀请中国工程院院士陈新等知名专家对大桥桥位进行咨询研究,并对选址报告进行评审。ﻫ2０05年10月，经省部协商，铁道部将沿江城际铁路及安庆长江铁路大桥项目补列入国家“十一五”规划，并向铁四院下达前期工作任务书。

200５年12月，铁道部主持召开了宁安城际铁路及安庆长江铁路大桥项目预可研报告审查会。

2006年元月,宁安城际铁路及安庆长江铁路大桥项目列入国家“十一五”规划.ﻫ2006年７月5日，安庆长江铁路大桥设计竞标工作会议在北京举行,经研究、审核、竞标，铁四院与中铁大桥院联合体中标。

2０06年９月,中铁大桥院和铁四院正式启动工程可行性研究报告的编制程序，同时编制了桥梁建设对长江航道的影响书，河势分析,桥位所在枯、中水流影响、流速及航迹图。ﻫ2006年１0月9日,安庆市与铁四院、中铁大桥院联合举行了宁安城际铁路及安庆长江铁路大桥工可勘察开钻仪式；1０月份武汉水文站和中铁大桥院近1０0人的设计、水文、勘测、物探等四个工作组进入安庆进行实地勘测。

20０９年3月2０日,安庆长江铁路大桥开工建设表１．２安庆长江大桥主要工程数量表项目材料规格单位数量备注混凝土主塔Ｃ50m351２00塔座Ｃ40ｍ3７９20预应力钢筋钢绞线T398ｆｐk＝1８60MPa粗钢筋T3８0fpｋ=９３０MPa普通钢筋Ｈrb３35Ｔ11２70钢料索导管2０号热轧无缝T90垫板Q235BT70检查设施Q２35ＢＴ187劲性骨架Q235ＢT12８0锚具群锚15—19套8４0JＬ—３２粗钢筋锚具套１1160波纹管100mmＭ1８1６24５mｍＭ57100金属扩张网ＧＭ50075型Ｍ2１550表2．1设计构件节点和单元数全桥节点数16９403全桥单元数２94６02主桁单元数1７66８0桥面板单元数1174９9斜拉索单元数２16桥塔单元数２07表2。2安庆长江铁路大桥模型主要参数列表材料弹性模量(Pa）泊松比密度（ｋg/m3）结构构件Q37０qＤ钢2.1e１10。３7850上弦杆、腹杆Ｑ420qD钢2。1ｅ110。37850下弦杆C50混凝土3.5e70．1625０0主塔C40混凝土３.５e70.16２500塔座斜拉索2．１e110．3７850斜拉索表2。３模型实常数结构部分实常数截面性质梁单元(beａｍ44）R1未定义参数（由截面文件直接读入）板壳单元（sｈell63)R２厚度16ｍm拉索单元（ｌink８)R3，R4，Ｒ5,R6（拉索面积)Ｒ,3，１３４.３1ｅ-4，Ｒ，４,１15。８4e-4,R，5，1０８.９１e-4，Ｒ，６，100。000e-4，/neｒr，0

！将参数读入数据库

fｉnish

/cｌearﻫ/PREＰ7

！单元类型编号

et，１，bｅam44

eｔ,2,sheｌl６3ﻫet，３,ｌink8ﻫ！材料类型编号Q4２0／Q３70/拉索－1，混凝土C40/C50—2,刚臂—3ﻫｍp，eｘ，1,2.1ｅ11ﻫmｐ,prxy，1,0．3

mp，denｓ，1,7850ﻫmp,aｌpx，１，１e-5

ｍp，ex，2,3。5e７ﻫmｐ，prｘｙ，2,0.１6ﻫmｐ,deｎs，2,250０ﻫmp,alpｘ,2，1e—５

ﻫｍp，ex，3,1e1５

mp,pｒxy,3，1

mp，dens,3，0

mp,alpx，3,0ﻫﻫ！读取自定义截面并编号

＊dｉm，namｅ，char，２4

name(１）＝’Sｅc_1ａ＇ﻫｎame(2)=＇Sｅc_2a＇ﻫname(3）＝'Seｃ＿3a'ﻫnａme(4）=＇Ｓec_4a'ﻫnamｅ(5)＝＇Seｃ＿5ａ’

name(6）＝＇Sec_6a＇ﻫnａｍｅ（７）=’Sec_7a’

nａme(8)＝’Ｓeｃ_1’ﻫｎame(9）=＇Sec＿２’ﻫｎame（１0)='Ｓｅc_３’ﻫname（11）=’Sec_４’

nａme（12）='Ｓec_5'

ｎａme（1３）=’Ｓec＿６'ﻫnａmｅ(１4）=’Sec_7’

name(15)='Sec＿8'

ｎaｍe(1６）=＇Sｅc_9’

name(17）＝'Seｃ_10'

name（1８)=＇Sec_11’

naｍe（１9)=’Sec＿1２’ﻫｎａme（20）=＇Sec_13＇

naｍe(21）=’Sec_14'ﻫnａme（22)=’Seｃ_1５’

name(23）=’Sec_1６'

name（24)＝'Sｅc_ｘ2'ﻫ

*ｄｏ，ｉ,1,2４,1

sectyｐe,i，ｂｅam，mesｈ，

seｃreａｄ，name（ｉ），，,meshﻫ＊ｅnddo

ﻫ！将tａ_ｓec参数读入数据库ﻫ!读取自定义截面并编号

*dｉm，nａme2，chａr，15ﻫｎamｅ2(１）='ta＿1'

naｍe2(2）=＇ｔa_2’

name2（3）='tａ_3＇

name2（4）='tａ_4'ﻫnａme2（5）='ta_5＇

ｎａme2(６)＝＇ta_６'ﻫｎａmｅ２（7）＝'ta_7'ﻫｎａｍe2（8)=’ｔａ＿8＇

name2(９）=＇ta_９＇ﻫname2(１０）＝'ｔa＿１0’ﻫｎａｍe2(11)=’ta＿11＇ﻫnamｅ2(12）＝'tａ_１2'

ｎamｅ２（1３)='ta_１3'ﻫnａｍe2（１4）=’ta＿１４'

ﻫ＊ｄo,i,1，14，1ﻫ

secｔyｐe，４０+i,beａm,mesh,

seｃread，naｍｅ2(i），,,ｍｅｓhﻫ＊enｄｄo

SECTYPE,

26，BEAM,HREC,，0

ＳECＯＦFSＥＴ,ＵＳER,0.1,0．２６6

ＳＥCＤＡTA，0.2，0。2６6,0.006,0。0０６，０．００6，0．0０６,0,0，0,0

ﻫSECTYPE,

２7，BEAＭ,Ｉ，，0ﻫSECＯＦＦＳEＴ,ＵＳEＲ,0，0.５84

ﻫSECDＡＴA,0.２40，0．２40,０．６，0。016，0。016,０.０12，0,0，0，0ﻫSＥCTYPE,

28,BＥAM，T，T_nin,0ﻫSECＯFFSET,ＵSER,0，1．4

SEＣDATA,０。6，１．4，0．032,0.0１8,０，0，0，０，0，0ﻫSECＴYPE,

29，BEAM,T，T_max，0ﻫSＥCOFFSET，USER,0，2.1ﻫSECDAＴＡ，0.6，２．１，0.032，0．０16，0,0,0,0，0，0

SＥＣTYPE，

31,BＥＡＭ,ＨＲEC,，0

SECＯFFSET，CＥNＴﻫＳECDATA，0.5４８，0.500，0。０１6，0.016,０.01６，0.016,0,0，0,0

SＥＣTYPE，

32，BＥAＭ，I,,０ﻫSＥＣOFＦSET，CENT

SEＣＤATＡ,0．560,0。５6０,0．500，0。028，0.028，0.0１６，0，０，0,０ﻫ!刚臂的截面ﻫSECTＹPE，

１00,ＢEAM，CＳＯLIＤ,,0

ＳECＯFFSEＴ，ＣENTﻫSＥCDATA，0．06，０,0，0,0，０，0,０，0，０

ﻫr,1，ﻫr，２，０.016,ﻫ！拉索截面控制

！７-3４9ﻫr,3,134．31e－4，ﻫﻫ！7－3０1ﻫr，4,115．８4e—4，ﻫﻫ！7-２83

ｒ，５，１08.91e－4,ﻫ

!7-24１

r，6,1０0．0０0ｅ-4，

ﻫｆｉｎish

!开始建立模型

/pｒep７ﻫ！＊*＊*＊**＊＊＊＊＊*＊＊*＊＊*＊*＊＊**＊＊**＊*＊＊＊＊**＊***＊＊＊*

！ﻫ！

！xｉaXiaｎGａn

！ﻫ!

！**＊＊＊＊**＊＊＊*＊＊＊＊**＊*＊＊＊＊＊＊＊＊＊*＊**＊＊＊＊＊**＊＊*＊＊

＊SET，j，-14

＊do，i,1，(94＊25+1）＊2，94*2５+1ﻫn，i，0，ｊ，0

n，94*25+i，136３．０00,ｊ，０

FILL，i，94＊25＋i，94＊25-1，，1,，，ﻫ!中桁ﻫ＊ｉf,ｉ，ｅq，９４＊25＋2,tｈｅnﻫ！7截面ﻫtype，１

ﻫｍaｔ，１

secnum，１4

ﻫ*dｏ,a，1，３*2５，1

E，i—1+a，i+a

＊enddo

ﻫ＊ｄo,ｂ,１，１1*25,1ﻫE，３5＊２5－1+i+b，35＊２５+i+ｂ

*enｄdｏﻫ*dｏ，ｃ，1，8＊25,1

Ｅ,７７*25－1+ｉ+ｃ,７7＊25＋i+c

＊enddｏﻫ!６截面

type,１

ﻫmａt,1

secｎum，13

ﻫ＊do,a，1,３＊25，1ﻫE,i+3*２５－1+a，i+３*25+aﻫ＊endｄｏ

ﻫ*do,b,1,3＊２5,1

E，32＊2５-1+i+ｂ,3２＊2５＋i＋b

ﻫ＊enｄdo

ﻫ＊do，c，１，3*2５,１ﻫE，4６＊２5-１+i+c，46*25+ｉ+c

ﻫ＊endｄo

*ｄｏ,ｄ，１,１＊25，1ﻫE，i+７6＊25－1+d,i+76＊2５+d

*ｅｎddoﻫ＊do，e，1，7*25,1ﻫE，8７*2５－１+i+e，８7＊2５+ｉ+e

＊enddo

!５截面

type，1

mat,1

ﻫsｅｃnuｍ,１2

*dｏ,ａ,１，２5,1

E，8*25-1+i+ａ,8＊２5+i＋aﻫ*eｎddｏ

ﻫ＊do，b，1,2＊25,1ﻫE，30＊25-1+i+b，30＊25＋i+b

*eｎｄdo

ﻫ＊do,ｃ,1，2＊25，１ﻫE，49＊２５－１+i+c,49＊25+ｉ+cﻫ＊enddo

*do,d，1,２＊25，1ﻫE，72＊2５-1+i+ｄ，７2*2５＋ｉ+ｄﻫ*endｄoﻫ！4截面ﻫｔｙｐe,１

ﻫmａｔ,1

sｅcnum，１１

ﻫ*do，a,１,3＊25,１ﻫE,9*25—1+i+ａ，9＊２５+i+ａﻫ＊enddo

ﻫ＊do,b，1，４*２５,1

E,2６＊2５－１+ｉ＋b,２６＊25+i＋b

ﻫ＊enｄdo

＊dｏ，c,1,3＊25,1

E，51＊25—1+i+c，5１＊2５＋i+c

*ｅnddｏﻫ＊dｏ,d,1，4＊25，1

E，68*25—1+i+d,68*２5+ｉ+dﻫ＊enｄdo

!３截面

tyｐe，1

mat，1

secnum,１0

ﻫ＊do，ａ，１,1＊25，1ﻫE,7*25—１＋i+ａ,7＊25+i＋aﻫ＊ｅnｄdｏ

ﻫ*dｏ，ｂ,1，５*25,1ﻫE，12＊2５—1+i+b，12＊25＋i＋b

＊eｎddｏ

ﻫ＊ｄo，ｃ，1，4*25,1

Ｅ，22＊2５－１+i+c,2２*25＋i+c

＊enddoﻫ*do,d，1,4＊25,１

Ｅ,54*25-１+i＋ｄ，54*25＋i+dﻫ＊endｄo

＊ｄo，ｅ,1,５*25，1ﻫＥ，63*25—1＋i+ｅ,63＊２5＋i+eﻫ＊enddo

！2截面ﻫtyｐｅ，1

ﻫmat,1

ﻫseｃnum，9

*do，a,1，1＊２5,1ﻫE，６＊２5－1+i＋a，6＊2５+i+aﻫ*ｅnddo

ﻫ＊dｏ，b,１，2*２5，1ﻫＥ，17*25—1+ｉ+b，17*25+i+b

ﻫ*enddo

ﻫ＊dｏ，c,1,1＊25,1ﻫＥ，21*２5－1+ｉ+ｃ，2１＊２5＋i+cﻫ＊enddoﻫ*do，d,1，１*25，1ﻫE，58＊25－１+i+d，58*25+i+d

＊eｎddoﻫ＊dｏ,e，1，2*25,1ﻫE,６１*25—1+i＋e，61*25+ｉ+ｅﻫ*eｎddｏ

＊do,f,1，2＊25,1

E,7４*２5—１+ｉ+f,7４*２5+i+ｆﻫ＊eｎddo

*do,ｇ,１，1＊2５，1

E，8５*25-1+i+g，85＊２5+i+g

＊eｎdｄo

!1截面

type,１

mａt，1

ﻫsｅｃnum,8

＊do,a,１，２＊25,1ﻫE,19*25-１+i+ａ，19＊２5+i+ａﻫ＊eｎｄdo

ﻫ*ｄo,b，１，2＊2５,1

E，５９＊２5—１+i+b,59*2５+ｉ+b

*enddo

＊ｄo，c，1，1＊25，1

Ｅ，86＊2５－1＋ｉ+ｃ，86*25＋i+c

＊enddoﻫ

！左边桁

*elｓｅif，ｉ,eq，(9４*２5+1）＊２+１，ｔhｅnﻫ！7截面ﻫｔypｅ，1

ﻫmat，１

ﻫseｃnum，7

ﻫ＊dｏ,a，1，3*25,1

E，i+a，ｉ-1+a

*enｄdo

＊do,ｂ,１，11*25,１

E,35＊２５+i+b，３5＊25—１+i+b

*eｎｄｄo

＊ｄo，c，１，8＊25，1

E,77*２5+i＋c，77*2５－１+i+ｃ

＊eｎddoﻫ！6截面ﻫtyｐe,1

mａｔ,1

ｓecnum,6

＊do,a，1,3＊25，1ﻫＥ,i+３*25+ａ，i＋3*2５-１+a

*endｄｏ

＊dｏ,b，1,3＊25，1ﻫE，32*25+i+ｂ,32＊25—1+i+b

＊enddo

ﻫ*do，c,1，3＊25,1

E，46*25+ｉ＋c，4６*２5—１＋i+c

＊eｎdｄoﻫ*do，d,1,１*２5，1ﻫE，i+76*2５+d,i+７６＊25-1＋d

*ｅｎddｏﻫ*do，e，1，７*25,1

Ｅ，87＊25+i＋e，８7＊2５-1+ｉ＋ｅ

*endｄo

ﻫ!5截面ﻫtypｅ,１

mat，1

ﻫｓｅcnｕm，5

＊dｏ，a,１,25,1

Ｅ,８＊2５+ｉ+ａ,8＊2５—1＋ｉ＋aﻫ＊ｅｎｄｄo

ﻫ*dｏ,b，1,２＊25,1

E，3０＊2５+i+b，３0*２5—1＋i+ｂ

ﻫ＊enddo

＊ｄo,ｃ，1，2*25，１

E，4９*２5+i+c，49＊25－1＋i+c

＊enddo

*ｄo，ｄ,1,2＊25,1

E,72*２５+i+ｄ,72＊25-1＋i+d

*eｎddo

！4截面

type,１

ﻫmａｔ，1

ｓecnuｍ，４

ﻫ＊do，a,１,3*２5，1ﻫＥ，9*25＋i+ａ，9＊２5-1＋i+a

＊ｅnddo

ﻫ＊do，b，1，4＊25，１

E,２6＊25＋i+ｂ，26＊２5—１＋i+b

*eｎddo

ﻫ＊do,c，１,3＊２5，1

Ｅ，５1＊25+ｉ+c，51＊２5-1+ｉ+c

＊endｄoﻫ＊do，d，1，４*２5，1ﻫＥ，68*25+i+d,6８＊25—1＋i+dﻫ*enddoﻫ！3截面

tyｐe,1

mat,１

ﻫsｅｃｎum,３

ﻫ*do，a,1，1＊25,1ﻫE，7＊２5+i+a,7＊２5-1＋i+a

＊ｅnｄdo

＊ｄo，ｂ，1，5＊25,1

Ｅ，1２*25+i+b，12*２5－1+i+b

ﻫ＊enddo

*ｄｏ,c，1，4＊２5，1

Ｅ，２2＊25+i+c，２2＊2５—１+ｉ+c

＊eｎｄｄｏﻫ＊do,d,１,4*2５，１ﻫE，54*２5+i＋d,5４＊25—1+i+d

＊eｎddoﻫ＊do,e，1，5＊25,1

Ｅ，６３＊２5＋i+ｅ，6３＊25－１+i+e

*ｅnddo

！2截面ﻫtype,1

mat，1

ﻫsecnum,２

＊do,a,1，1＊25，１

E,6*25+i+a，6*２5—1＋i+aﻫ＊enddo

*ｄｏ,ｂ,1,2＊25,1ﻫE,１7＊25+i+ｂ，17*2５—１+ｉ+ｂ

ﻫ*enddo

ﻫ＊ｄｏ，c,１，1*25,１ﻫE，21*25+i＋c，２1＊25－1+ｉ+cﻫ*enddoﻫ＊do,d，1,1*２5，1

E，５8＊２5+ｉ＋ｄ，5８*２5－1+i+d

＊eｎddo

＊do,e，1,2＊25,1ﻫE，6１*25+i+ｅ，61*25－1＋i+ｅﻫ*ｅnddo

＊do，f,1，2*25，1ﻫE，74*２5+i+f，74＊2５-１+ｉ+fﻫ*ｅndｄoﻫ＊do,g,1，1*2５,1

E,８5*25+i＋ｇ，８５＊2５－1+ｉ+gﻫ*enddｏ

！１截面

typｅ，１

ﻫｍａt，1

secnum，1

*do,ａ，１，2＊25，１

E，19＊25+i+a,19*25－1+i+a

＊enｄdｏ

ﻫ*ｄo,ｂ,1，2*25，１

E，59＊2５＋i＋b，５9*25-1+ｉ+b

*ｅndｄo

*ｄo，c，1，1＊２5,1

Ｅ，86＊25+i＋c,8６*25—１+i+c

＊enddoﻫ

!右边桁ﻫ*elsｅ

ﻫ！7截面

tｙpe，１

mat，1

secnｕm,7

ﻫ*ｄo，a,１，3*2５,１

E，i—１+a，i+aﻫ*ｅｎｄdo

＊do，b，1，1１*25，１

E，35*25-1+i+b,35＊２5+i+b

ﻫ＊enddo

＊ｄo,ｃ,1，８*25,1

E，７7*2５—1+i＋c,77*25＋i+ｃ

＊enddｏﻫ！6截面ﻫtyｐe，1

mat,1

ﻫsecnuｍ,6

*dｏ,a，1，３*2５，1

E，i+3*2５－1+a,ｉ+3＊2５+ａﻫ＊eｎｄdo

＊ｄo，b，１，3*25,1ﻫE，32＊2５—1+i＋ｂ，32*25+i+ｂ

ﻫ*enddo

＊do，c，１,3＊25，1

Ｅ，46*２5—1+i+ｃ，４6＊25+i+c

ﻫ＊ｅnddoﻫ*do，ｄ,１,１＊２５,1ﻫE,i+7６＊２5—1+ｄ,i+76*25+ｄ

＊enddoﻫ*do,ｅ，1，７＊25，1

Ｅ，８7＊25-1+ｉ+e，87*2５＋i+ｅ

ﻫ*endｄo

ﻫ！5截面

tｙpe，1

mat,1

ｓecnuｍ，5

ﻫ＊do,a，１,25,1ﻫＥ,８*25-１+i+a，8＊25＋i+ａ

*enddo

*ｄo，b，1,２＊２5,１

Ｅ，３0＊25-１+i+b，30＊２５＋ｉ＋b

＊enｄｄo

*do，c，1,2＊2５，1

E，49＊２5－１+i+ｃ，49＊2５＋ｉ＋cﻫ*enddo

*ｄo，ｄ，1，2＊２5,1ﻫE，72＊25-１+i＋d，7２＊2５+i+dﻫ＊enddoﻫ！4截面ﻫｔype，1

ｍat,1

secnum，4

ﻫ＊do,a,1，3*25,1

Ｅ,９＊25-１＋i＋ａ，9＊2５+ｉ+aﻫ*eｎddo

＊dｏ，b，1,4＊2５，1ﻫE，２６＊2５-１+i＋b,２６*25+i＋ｂ

ﻫ＊enddo

＊ｄo，c，1，３＊25,1

Ｅ,５１＊２5—1+i＋c，5１*25+i＋ｃﻫ＊enddo

*dｏ，d，1,4＊25,１ﻫE，６8*25-1+i+d，6８＊２5+i+dﻫ*ｅnddo

!3截面ﻫtypｅ，1

mat,１

secnum，3

ﻫ＊do,a，1，１*25，1ﻫE，７＊2５－1+i+ａ,7＊25+ｉ+aﻫ＊ｅnｄｄo

ﻫ＊dｏ，b,1，5＊25,１ﻫE，１2*２5－１+ｉ＋b，12*25＋i＋b

＊ｅｎｄdo

＊do，ｃ,１，４＊25,１

E，2２*２5-１＋i+c，22＊25+i+cﻫ*enddｏﻫ＊do,d,1，４＊25,１ﻫE，５4*25—1+i＋ｄ,54*2５+i+d

＊ｅnddo

＊dｏ，e，1，5*25，1

E，63*25-１+i+e，６3*25＋ｉ+eﻫ＊enddoﻫ！２截面ﻫtype，1

ﻫmａｔ，1

ﻫseｃnum，２

ﻫ＊do，a，１,1*2５，1

E，6*2５－１+ｉ＋a,６*25+i+aﻫ*endｄｏ

＊do，b，1,2＊25，1

E，1７＊25－1+i+ｂ,1７＊25+ｉ＋b

ﻫ＊enddo

*do，c，1,１＊25，1

E，21*２5－1＋ｉ+ｃ，2１＊25+i+c

*ｅnddｏﻫ＊do，d,1，1*2５,1

Ｅ，58＊25-1+i+d,５８*２5+i+d

＊eｎddoﻫ＊do，e，１，2*25，1

E,６１＊2５—1＋ｉ+e，61＊25+i+e

＊endｄo

*ｄo，f，１,2*2５，1

E，74＊2５-1+ｉ+ｆ，７4＊2５+i+f

*eｎddoﻫ*dｏ,g,1,１＊2５，1ﻫE,85*25—1＋i+g，８5＊2５＋i+ｇ

*endｄo

!1截面ﻫtｙｐe,1

maｔ,１

ｓecnum，1

*ｄo,a,1，2*25,1

E,1９*25－1＋ｉ+a，19＊２５+i+a

＊enddo

*do，b，１,２＊2５,1

Ｅ,59＊２5－1+i+b,59＊25＋ｉ+b

ﻫ*enddo

ﻫ*ｄo，c,1，１*25，1ﻫE，86*25－1+i+ｃ，８6＊25+i＋c

*ｅnddo

*endifﻫ

ﻫ＊SET，ｊ，j+１4ﻫ＊endｄo!＊＊＊＊＊＊*＊＊＊＊＊＊＊＊＊＊***＊*＊＊＊＊＊**＊*＊*＊＊*＊＊***＊＊*＊ﻫ！

!

！T_beamAＮDshell６3'sｎodes

!

！ﻫ！＊＊*＊**＊**＊＊**＊＊*＊＊*＊＊＊＊＊＊＊**＊＊***＊***＊＊＊**＊＊*

/prｅｐ7

nｇen，2，705３，1，2３5１,１,0,0,0.7，ﻫngen,4,23５１,7０5４，7０54+2350，1，０，０。57，0,ﻫngen,4,２351,7054+3*2351，７0５４＋３*２351＋2３５０,1,0,０。60,0,ﻫｎgen,3，235１，７０54＋6*２３51,7054+6*２３5１＋2３5０,1，０，0．45，0,

ngen,2,2351,７05４+8＊2３５1，7０54+８＊2351+2３５０，1,0,0.6０,0，

ngen,3,235１，7０5４+9*2351，7054+9＊23５1+2３５０,1，０,0。45,0,ﻫngen,2，2３５１,7054+１1*2３５1，7０54＋１1＊23５1+235０，1,0，0.7０,0，

ｌocal,11,０,0,—７，0ﻫﻫnsym,y,2351*12，7054+2351，３29１４+2*23５1，1ﻫﻫcsys，0

nsｙｍ,y,5877５,7054,65828,1ﻫ

nｇen，2,1222５2,2352，２３52+2350，1，0，0,0。7,ﻫ

！生成边上４条Ｕ肋

type,1

ｍａｔ，1

secnum，2６ﻫ*ｄo，h,6，9，１ﻫ＊do，ｋ,1,235０，1ﻫe,２３５1＊ｈ＋k,23５1＊h+１+kﻫ*enｄdo

＊ｅｎｄdoﻫ

!生成T梁1

typｅ，1

mａt,1

secｎum，27ﻫ＊ｄo，k，1，２３50，1

e,23５1＊10＋k,2３51＊10+１+ｋ

＊enｄdｏﻫ

！生成U肋２条

typｅ，1

ﻫｍat，１

ﻫsecnuｍ，2６ﻫ＊ｄo,h，11，12，１

*ｄo,k，1,2３50，１

ｅ，２351＊h＋ｋ，2351＊h+1+k

*endｄo

*ｅnddｏ

！生成Ｔ梁2ﻫtypｅ，1

ﻫmat,1

ﻫｓecnum，27ﻫ*ｄo，k,1,2350，1

e，2３51*13+k,23５1*13+1+k

＊enddｏ

ﻫ！生成U肋另外2条ﻫtｙpe,１

ﻫｍaｔ，1

ﻫsｅcnｕm，26ﻫ＊ｄｏ,ｈ,14，１5,1ﻫ*do，k，1,２35０，1ﻫｅ，2３51*h+ｋ，２351＊ｈ＋1+kﻫ*enddｏ

＊endｄｏ

ﻫ！镜像四分之一纵向的单元

csｙs,11

!（最好用*ｇet最大节点号单元号)

ｅsym，，2３5１*12,７051,3055０,1ﻫﻫ！生成横梁

typｅ,1

maｔ，1ﻫsecｎuｍ，29

＊do,h，7054,9４04，25ﻫ*do,ｋ,h，h＋1１＊2351,２351

e,k，k＋2351

*enｄｄo

＊enddｏ

ﻫ＊dｏ，h，37617,399６７，２５

*do,k，ｈ，h＋10＊2351，235１ﻫe,ｋ,k＋23５1ﻫ＊ｅnddoﻫ＊enddｏﻫﻫ＊do，h，3５266，37616，25

ｅ，h，h+2３５1*１２ﻫ＊enddo

ﻫ！生成横肋

type,1

maｔ，1

ｓｅcｎum，2８

＊do，ah,7059,93８４,25

*do,ａc，0，15,５ﻫ＊dｏ,ak,aｈ+ac，ah+ac＋１1＊235１,2351

e，ak，ａk+235１ﻫ＊enｄｄoﻫ＊enｄdo

*ｅnｄｄoﻫﻫ*do,ah,37６22，3９94７，２５

＊do,ａｃ,０，1５，５

*do，ａk，ah+ac,ah＋ac+1０＊2３5１,２351ﻫe,aｋ，aｋ+2351

＊enddoﻫ*enddoﻫ＊ｅnddo

＊do，ａｈ,3５27１,37596，25ﻫ＊do,ac，０,15,5ﻫe，ａh＋ａc，ah+ac＋２3５１＊1２ﻫ＊enddo

＊enddoﻫﻫ!镜像一半单元ﻫｃsｙs,0

ﻫeｓym,,58775,7051，653５4,1ﻫ

！补充横梁

ｔypｅ，１

mat，１

ﻫｓeｃｎum,２９

＊ｄo，ｈ，376１7，39９67，25ﻫe,ｈ，h＋869８７ﻫ＊ｅｎddo

*do，h,96３9２,9８７４2，2５ﻫe，h，h＋2821２

*ｅｎddoﻫﻫ!补充横肋ﻫtypｅ，1ﻫmａt,1

sｅcｎｕm，2８ﻫ＊do，ah,37６22,3９947,2５ﻫ*do，ac,0,1５，５

e，ah+ac，ah+ａc＋86９8７ﻫ*enｄdｏﻫ＊enddoﻫ*do，ａh，96３９７,98722,25

*ｄo，ac，０，１5,5

e,aｈ+ac,ａｈ＋ac+28212ﻫ*enｄｄｏﻫ＊enｄdo

ﻫ!刚臂

ｔyｐｅ,1

mat,3

ﻫsecnum，100

＊do，h，1，２351,1

e,h,h+７053

＊enｄｄoﻫ＊do,h,4703,7053，１

e，h,h＋61126

*ｅndｄoﻫ*do,h,2352,4702，1ﻫe，h，ｈ＋１246０4－2３52

＊enddｏﻫﻫ！**＊＊＊**＊*＊**＊＊＊＊*＊＊＊＊＊＊＊＊＊＊＊**＊**＊＊＊*＊*＊＊***＊ﻫ！

！

!shell

!ﻫ!ﻫ！*＊＊＊***＊**＊＊＊＊*＊＊＊****＊＊*＊＊*＊**＊＊＊＊*＊＊*＊*＊＊*＊

typｅ，２

ﻫmａｔ，1

real,2ﻫ*do,i，0，11，1

ﻫ＊do，h，7054，7０５４+234９,1ﻫe,h+i*2351,ｈ+１+ｉ*2351,h＋2352+ｉ*2３５1,h＋２351+i＊23５1ﻫ*enddoﻫ＊ｅnｄdo

ﻫ＊do，ｈ，35266，35２6６+2349，1ﻫｅ,ｈ,h＋1,h＋6347９－35266,ｈ+63478—３5２66

*ｅnddo

ﻫ*do,i,0，－10，—1

ﻫ*ｄo，ｈ，634７８,６3４78+２３4９,1

ｅ，h+i＊2３51，h＋1+i＊２351,h＋ｉ*2351－２350,ｈ+i*2351－2351ﻫ＊enddo

*enddo

＊dｏ，ｈ,3７617，37６17+23４9，１

ｅ,h，h+１，h+1+124604—376１7，h+1２4６04－3７6１7ﻫ*enddo

ﻫcｓys,０

esｙｍ，，5８７75，131６５4，190403,１ﻫ

＊ｄo，h，1２４60４，１２4６04+2３４9,1ﻫe,ｈ,h＋1,h+1—（12４604-9６392），ｈ－(12４６04－９639２）ﻫ*enddoﻫﻫ

!＊＊＊**＊***＊＊＊*＊*＊＊＊*＊＊***＊＊＊＊*＊＊*＊＊*****＊**＊＊＊ﻫ！

!

！SｈangXｉanGanﻫ！ﻫ!

!＊＊****＊＊*＊****＊*＊＊*＊＊＊*＊＊**＊***＊＊*＊＊**＊***＊*＊

/prｅp７

ｒeal，1ﻫ＊ＳEＴ，Nmax1，３0０0００

＊ＳET，j,-１4

＊do,ｉ,Nmａｘ1+1，Nmａx1+（92＊25+1)＊２，92＊25+1

ｎ，i，14。5０，j，15.０0

n，92＊25＋i,134８.５0,j，15.00ﻫFILＬ，ｉ，92＊25＋i,，,，,，

！中桁ﻫ*iｆ，i,eq,Nmax1+９2*25+2，tｈenﻫ！13截面ﻫtyｐe，1

ﻫｍａt，1

ﻫｓecnuｍ,２０

*do，a，1，2*25,１ﻫE，i—1+ａ,i+a

*eｎddo

＊do，b,1，６*２5，1

E,3＊２5－１＋i＋b，3*25+i＋ｂ

*eｎdｄｏ

*ｄo，c，1，３*２5,1ﻫＥ,７0＊25－１+ｉ+ｃ，70*25+i+c

＊ｅndｄｏﻫ*do,d,1,１0＊25,1

Ｅ，74＊25-１＋ｉ＋d，７4*25+i+ｄ

*enddｏﻫ＊do，e，1,３*２５，１ﻫE,８5*25-１+ｉ+e，８5*25+ｉ＋e

＊eｎddo

＊ｄo，f，1,1＊25,1ﻫE,91*25－1+i＋ｆ，9１＊25+i+f

ﻫ*eｎddoﻫ!１2截面ﻫｔyｐe,１

mat,1

ﻫｓｅcnum，19

＊dｏ，a，１，１＊25，1

E，i+２＊25-1+a,ｉ+2*2５+ａ

*enｄdo

＊do,b，1,1*2５，１ﻫE，9＊25—1+ｉ＋ｂ，９＊２５＋ｉ+b

＊endｄo

ﻫ＊do,ｃ,1，8＊２5,1ﻫE，35*25—1+i+ｃ，35*２5+ｉ＋ｃ

*enｄdｏﻫ＊do，ｄ，1，１*25,１ﻫＥ，i＋６9＊25－1＋d,ｉ＋6９*25+dﻫ＊ｅｎdｄoﻫ*ｄo，e，１,1＊25,1ﻫＥ，73*25－1＋i+e，73*25+i+e

ﻫ*eｎddoﻫ＊do,f,1，3*25,1ﻫＥ,88＊25—１+ｉ+f,88＊２５+i+f

*enｄdo

！１1截面

ｔyｐe，1

mat，1

seｃnuｍ,18

ﻫ*do，a，1,25,1ﻫE，10*25-１+i+a，１0*2５+i＋ａﻫ＊enddo

*do,b,１，1＊２5，1

E,1８*25—１＋ｉ+b，18*25+i+b

＊eｎｄｄo

ﻫ＊dｏ，c,1,1*25，１ﻫE,26＊25-1＋i+ｃ，２6*25+ｉ＋cﻫ*enddo

＊do,ｄ，1，４*2５,1ﻫE，3１＊２5—1+i+d，31*25+i+d

＊eｎddｏ

＊ｄo,ｅ,1，1*25,1ﻫE,43＊２5-１＋i+ｅ，43＊25+i+ｅﻫ*enddｏ

*ｄo,f,1，５*25,1ﻫE，45＊25—1+i＋f,４５*２5＋i+ｆ

＊eｎdｄｏ

＊do，g,１，1＊25，1ﻫＥ,59＊25－1+i+g,59＊2５+i＋gﻫ*enddoﻫ＊dｏ，ｈ,1，１＊２5,１ﻫE,68＊25-1＋i＋h，68＊2５+ｉ+ｈ

*ｅndｄoﻫ！１0截面

type，1

mat，1

sｅｃnｕｍ，1７

ﻫ*ｄo,a，1,１＊25,1

E，１1＊25—１+i+ａ,11*25＋i+a

*ｅnｄdｏ

＊do，b,1，１＊25，1

Ｅ，17*25—1+i+ｂ，1７＊２5+i+ｂ

ﻫ*endｄo

＊do，ｃ,1,1＊２5，１ﻫＥ,20*25—1+i＋c，２0＊25＋ｉ＋c

＊eｎｄdoﻫ＊ｄo,ｄ,1,1*25,1ﻫE,25＊25-1+ｉ+ｄ,２5＊２5+i+dﻫ*ｅnｄdｏ

＊dｏ，e,１，4*25,1ﻫE,27＊2５-1+i＋e,2７＊25+i+e

*eｎｄｄoﻫ＊dｏ，f，１，3＊25,1

E，５0*25—1+i+f，50＊25+i+ｆ

＊ｅｎｄdｏ

＊ｄｏ，g,１,1＊2５,1ﻫE,5７＊２5－1＋i+g，57＊２5+ｉ+g

＊ｅndｄｏﻫ＊do，h，1，1＊25，1ﻫE,60＊25—1+i+h，6０＊２5+i+ｈﻫ＊endｄo

*dｏ，hh,1，１*2５,1

E，６7＊25-１+ｉ+hｈ,６7＊25+i+hh

＊ｅnｄdoﻫ!9截面ﻫtｙｐe，1

mat,１

secnum，16

*do,a,１,５＊2５,1ﻫE，12＊2５—１＋i＋a，12＊25+i+aﻫ*enｄｄｏ

ﻫ＊dｏ，b，１,1＊２5,１

E,19*２5—1+i＋b,19＊25+i＋b

＊enddｏ

*ｄo，ｃ，1，4＊2５，１ﻫＥ，21*25-１+ｉ+c，21＊２5＋i+cﻫ＊enddoﻫ＊do，d，1，4＊２5，1

E，53＊２5—1+i＋d，53＊２5+i+d

＊enddo

＊dｏ,e,1,1＊25，１ﻫE，５8＊2５—1+ｉ+e,５8*2５＋i+eﻫ*enddoﻫ＊ｄo,g,１，6*25,1

E，61*25-1+i+g，61＊25+ｉ+gﻫ＊enddoﻫ＊dｏ,h，1，1＊2５，１ﻫE,８４＊25-1+i+h,84*25＋ｉ+hﻫ*enｄdo

!8截面ﻫｔype,１

ﻫｍat，1

ﻫsｅcnum，１5

ﻫ＊dｏ,a，1，１*25,1ﻫE,４4＊25-1+ｉ＋a,４4＊25+ｉ+ａﻫ＊enｄｄoﻫ

!左边桁ﻫ＊elsｅif,i，ｅｑ，Nmax1+（９2*25+1）*2+1，theｎ

！1３截面ﻫｔyｐe,1

ｍat，1

ﻫsecnum,20

ﻫ*do，a,１，２*25,１

E，i—1+ａ,i＋aﻫ*enｄdo

ﻫ＊do,b，１,6*25，1ﻫE，３*２5-１+i+b，3＊２5+i＋ｂ

ﻫ*endｄoﻫ＊do,ｃ，1，3*25,1

E，70＊25—１+i＋c，70*2５+i+ｃ

ﻫ*ｅnｄdo

＊do，d,1,10*25，1

E,7４＊25－１+i+d，7４＊2５+i＋d

*enddｏ

*ｄｏ，e,1,3*２5，1

E,８5＊25－1+i+e,８５＊2５+i＋e

＊enｄdoﻫ*do,f，1，1＊2５，1

E,91*25-１+i+f,９1＊25+ｉ+f

ﻫ＊ｅｎddo

！12截面ﻫｔype，１

ﻫmaｔ,1

sｅcnum，１９

ﻫ＊dｏ,a，1,1＊２5,１ﻫE，i+2*25－1＋a，ｉ+２*25＋ａﻫ*endｄo

ﻫ*dｏ,ｂ,1，1＊２5,1

E,9*2５—1+i+ｂ，９*25+i+b

＊enｄdo

ﻫ＊do，c,１，8*２5,1ﻫＥ,３5＊25—１＋i+c,３5＊25+i+c

ﻫ＊eｎddo

＊ｄo，d，1,1*２5,1

Ｅ,ｉ+69＊2５－1+d，i+６9*25+d

*eｎddoﻫ＊do，e，１，１*２5，1

Ｅ，73*2５—1+ｉ+ｅ，73*25+i＋e

＊enddoﻫ*ｄo，f，1，3＊25，1

E，８8＊25－1+ｉ+f，88＊25+ｉ+ｆ

ﻫ*enｄdo

ﻫ！11截面ﻫtype,1

ﻫmaｔ，1

secnuｍ,18

＊ｄｏ,a，1,25,1ﻫE，10＊2５—1+i＋a,10*25＋i+aﻫ*endｄo

*dｏ,b，１,１＊2５，1ﻫE，18*25－１＋ｉ+b，18＊2５＋ｉ+ｂ

ﻫ＊enｄdo

*ｄo,c，1，１＊2５,1ﻫE,2６*２5—1＋i＋ｃ,２６*２5+i+c

＊ｅnddｏ

＊dｏ,d,1，4＊2５,1ﻫＥ，３1*25—１+i+d，31＊25＋i＋ｄﻫ*eｎddｏ

*do,e,１，1＊25，1

E，43＊25-1+i+e,43＊25+i+ｅ

*enｄｄoﻫ＊do,f,1，５*25,１ﻫE，45＊25—１+ｉ+f,4５＊2５+i＋ｆ

＊enddoﻫ＊ｄo,g,１，１*2５,1ﻫＥ,59*2５-1+ｉ+g，59＊25+i+ｇ

＊ｅｎddo

*dｏ,h，1,1＊25，1ﻫE，68*25－1＋ｉ+h，68＊25＋i＋h

＊enddoﻫ!10截面ﻫtype，１

mat，1

ﻫｓecｎum,17

＊do，ａ,1,１＊25，1ﻫE，11*2５-1+i+a，11*25＋ｉ+aﻫ*endｄo

ﻫ*do，b,1,1*2５，1

E，17＊2５—1+ｉ+b,17*２5+ｉ+ｂ

＊enddｏ

＊ｄｏ，c，１，1＊25,1

E,２０*25-1＋i+c,20＊２５+i+ｃﻫ＊ｅnddｏ

＊dｏ，ｄ,１,1＊２5，1

E,25＊25－1+ｉ＋d,2５*２5+ｉ+dﻫ＊enddｏﻫ＊do，e，1,4*２5,1ﻫＥ，27*25－１+i＋e，２7*25+i+eﻫ＊ｅndｄo

*ｄo，f,1，３＊25，1ﻫE,50*25—１＋i+f，50＊25+i+f

*enｄdoﻫ＊do，g,1，1＊25，1

E，57*25-１+i+g，57＊2５+ｉ+gﻫ＊ｅnddo

*ｄｏ，ｈ,１,1*2５，１

E，６0＊25－1＋i+h,6０*2５＋i＋hﻫ＊eｎdｄｏ

＊do，hh，１，1＊25，1

Ｅ，６7*2５—1＋i+hh，6７＊2５+ｉ+hh

*ｅnddｏ

!9截面ﻫtｙｐe，1

ﻫmat，1

ﻫsｅcｎｕm,1６

＊do，a，1,５*２5,1

E,12＊２5-1+i＋a,1２＊２5+i+ａ

*enｄｄo

ﻫ＊do,b，1,1＊25，1

E，１9*2５—１+i＋b,１9＊25＋ｉ＋b

＊enddo

ﻫ＊do,c，1，４*２５，１ﻫＥ,２１*25－1+i+c，21＊25+i＋ｃ

＊enddoﻫ＊do，d，１，４*2５，1ﻫE，５3＊２5—1+i+d，53＊25+i+ｄﻫ＊ｅnｄdo

*do,e,１，1＊2５，１ﻫE,5８*25－１＋i+e，58＊2５+i＋eﻫ＊endｄｏﻫ＊do,g,１，６*25，1ﻫE,6１＊２5—1+ｉ+g，61*2５+i+g

*enddoﻫ＊dｏ，h,1,1＊２5,1

E，84＊２５—１+i+h,84*2５+i+ｈ

＊enｄdoﻫ!8截面

ｔｙpe，１

mat，1

ﻫsｅｃnum，1５

ﻫ＊do，a,１,1*25,1

Ｅ,44＊2５—1+i+a,４4＊２５+i+ａﻫ*ｅndｄo

ﻫ！右边桁

*ｅｌse

！13截面

type，1

ｍaｔ，１

ﻫｓｅｃｎｕm，2０

ﻫ*ｄo，a,1,2*25，1ﻫE，i-1+a，i+a

＊enｄｄo

ﻫ*ｄｏ,b，1,6＊25，1ﻫE,3*25－１+i+b，3*２5+i+b

＊enddｏﻫ＊ｄo，c,1，3＊25,1ﻫE，７0＊25－1+i+c,70＊25+i＋c

＊enｄdoﻫ*do，d，１,10＊25,1

Ｅ，74*25-１＋i＋ｄ,74＊2５+i+d

＊enddｏ

*do，e，1,3*25，１ﻫＥ，８5*25—１＋i＋e,８5*25＋i+e

ﻫ*enddoﻫ*do,f,1,1*２5，1ﻫE，９1＊25—1+i+f，91＊２5+i+ｆ

ﻫ＊enddｏﻫ！12截面

ｔypｅ,1

ﻫmａt，１

sｅcｎum，1９

＊ｄo，a，1，1*25，1

Ｅ，i+2*２5—１+ａ，ｉ+2＊2５+aﻫ＊eｎdｄo

＊do,ｂ，1,１*25，１ﻫE，9＊２5-1+ｉ+b，９＊25＋ｉ+ｂ

*enddo

ﻫ＊dｏ,c，1，8*25，１ﻫE,35＊2５—1+i+ｃ，35＊25+i+c

ﻫ＊endｄoﻫ＊dｏ,d，1,1＊25，１

E,ｉ+6９＊25－1+d，i+69＊25+ｄﻫ＊ｅnddoﻫ*dｏ,e,1，1*２5,1ﻫE，73*25—１+i+ｅ，73＊25＋i+e

ﻫ*eｎddo

*dｏ,f，1，3*２5，1ﻫE，88＊２5—１+i+f，8８*２５+i+ｆ

ﻫ＊enｄｄo

ﻫ！11截面

tyｐｅ,1

ﻫmat，1

ﻫsecnum，18

＊do，a,1，25，1

E,10＊２5－１+i＋a，10＊2５+i+aﻫ*enddo

*do,b,1，1＊２5，１

E，18*25－１＋i+ｂ，1８＊2５+i+b

ﻫ*enddo

*do，ｃ,1，1*2５,1ﻫE，2６＊25－1＋i+c，26*25+i+cﻫ*eｎddoﻫ＊do,d，１,4＊２５,１

Ｅ,3１*25-１＋i+d,３１*25+i+dﻫ*enddoﻫ＊ｄo，ｅ,1,1＊25，1ﻫE,43＊２5—１+i+e,43*25＋i+eﻫ*ｅｎddｏﻫ*ｄｏ，f，1，5*25,1ﻫE，４5＊25—1+ｉ＋ｆ,45＊25+ｉ＋f

＊ｅｎdｄｏ

*ｄo，g,1，1＊25，1

E，59＊25—1+i＋g，59＊25+i+g

＊eｎdｄo

＊ｄo，h，１,1*２5，1ﻫE,6８＊25－1+i＋h,68＊２５+ｉ+hﻫ＊enｄdo

!10截面

ｔyｐｅ,1

mat，1

ﻫsecnum，１7

ﻫ＊dｏ，a,1，１＊25，1

E，11＊25—1+ｉ＋a，1１*２5+ｉ+aﻫ*endｄo

＊dｏ，b，1,1*2５，1ﻫE,１7＊2５-1＋i＋b，17*25+i+b

＊ｅnddｏ

＊do，c，1，１＊２5，１

E,20＊25-1+i＋ｃ，２0＊25+i+ｃﻫ*eｎddoﻫ＊do，d，1，１＊25，1ﻫE，25＊25-１+i+ｄ,25＊25＋i＋dﻫ*enddoﻫ*ｄo，e,1，4＊25,1

Ｅ，27＊25—１+ｉ+e,27＊25+ｉ+eﻫ＊enddｏ

*do，f,1,3*25,1ﻫE，5０*２5－1＋i+f，５0＊25＋ｉ+ｆﻫ*ｅnddo

*ｄo，ｇ,1，1*２５，１

E，57*2５—1+ｉ+g，57*25＋i+gﻫ＊ｅndｄoﻫ＊ｄo,h，1,１*25,1

E，60*25—1+i+h,60＊2５＋i+ｈﻫ＊enddo

＊do,hｈ,1,1＊25,1

Ｅ，６７*25—1+i+hh，67*25+ｉ+hh

＊enｄdo

！9截面ﻫtｙｐe，1

mat，1

ｓecnum，16

＊dｏ,a,1,５*25,1ﻫE，12＊25－１+i+a，12*25+i+a

＊endｄo

ﻫ＊dｏ，b，1,1＊2５，１

Ｅ，１9*２5－1+ｉ+b,19＊25＋i＋ｂ

ﻫ*enｄdo

*dｏ，c,1,4*2５,1ﻫＥ,21＊２5—1+i+c,21＊２５+i＋c

＊enddｏﻫ＊do,d，1，4＊25,１

E，53＊2５-１+ｉ+d,53*25+i+dﻫ＊ｅnddｏﻫ＊dｏ，e，１,１＊25，１

E,58＊２5—1+i+e,58＊25＋i＋e

*endｄo

＊ｄo，g,１,6＊25，１

Ｅ，61＊２５-1+ｉ+g,６1*2５＋i＋gﻫ*enddoﻫ*do，h，1，1*2５，1ﻫＥ，84＊25—1＋i＋ｈ，84＊25+i+ｈ

*enddｏ

！8截面

tｙｐｅ,1

ｍａt,1

ﻫsecnuｍ，1５

ﻫ*do，a，1,1*25,1

Ｅ,44＊25－1+ｉ＋a，44*25+ｉ+ａ

＊enddｏ

ﻫ*ｅｎdiｆ

ﻫ＊ＳEＴ，ｊ，j+14ﻫ＊enddo！*＊***＊＊*＊＊*＊**＊＊*＊*＊*＊*＊＊**＊＊*＊*＊***＊＊***＊*＊＊

！ﻫ!ﻫ!ＦuGanﻫ！ﻫ!

！＊*＊＊＊**＊*＊＊*＊*＊＊*＊＊***＊*＊***＊＊＊＊*＊＊＊***＊*＊＊*＊ﻫ/prep7

＊SET，j，－14

ﻫ＊dｏ，i,1，（９4＊25+１）＊２,9４＊2５+１

k,i,０,j,0

ｋ,9４＊25＋i,136３.000，j，0

kFILＬ，i，9４*2５+i，9３,，，，，

ｊ=j＋1４

*enｄｄoﻫ＊SＥT,Ｎmax1,10０00

＊SET，ｊ，-14

ﻫ*ｄo,i,Nmax1+1，Nmaｘ1+（92*25＋1）*２，92＊25+1

k，i，14．50,j,１5.00

ﻫk，9２＊2５＋ｉ，1348。５0,j，１5。０0ﻫkFILL，ｉ,92*２５+i，91，，,,，ﻫj=j＋1４ﻫ＊enddｏﻫ！连接右斜线1ﻫj=Ｎmax1+1

*do，n,１,（９4*2５+1）＊２，９4*25+１

kｋ=0ﻫ*dｏ，ｉ,n，n－1＋20＊２５,２5ﻫｌsｔr，i,ｊ

kk＝ｋｋ+1ﻫ*if，kｋ，eｑ,１,tｈen

ｌatｔ,1,1，1,，,,2１

esize,，２5

ｌmesh，ａll

＊elseif,ｋk,ｇe，4，aｎd，kk,le，６，then

latｔ，1,1,1,,,，2１ﻫesｉzｅ，，25

lmesｈ，aｌl

＊eｌseif，kｋ，gｅ,8,and，kk，le，９,thenﻫlatt,1,1,1，，,，２１ﻫesｉｚｅ，,２5

ｌmesｈ,ａllﻫ＊elsｅiｆ,kｋ，eｑ,7,thｅnﻫｌatt,1,1，１，,，，２3ﻫeｓiｚｅ，,2５

lmeｓｈ,all

＊ｅlseﻫlaｔｔ,1,1,1,,，,２3

esiｚe，，25ﻫｌｍｅsh，all

*endiｆﻫj＝j＋25ﻫ＊enddoﻫj=ｊ＋2301-（20-0）＊２５ﻫ*enddｏﻫ!连接右斜线2

j=Ｎmax１＋1+４0＊2５

*do,n，1，(９4*25＋1)＊2，9４*２５+1

kk=0ﻫ*do，i,n＋40＊25，ｎ-1＋６0*25,２5

lstr,ｉ,ｊ

ｋｋ=kk+1

*if,kk，eｑ,２0，ｔhen

lａtt，1,１,1，，，，21

esiｚe，，25ﻫlmｅsh,aｌlﻫ＊ｅlseif，kk,ge,3,anｄ,kk，le，9,tｈenﻫlaｔt，1,１，1,，,,21

esiｚe，,25ﻫlmeｓｈ,alｌ

*elsｅﻫlatt，1,１，１,,,，2２

esiｚe,，２5

lｍesｈ，all

＊endiｆ

j=ｊ+2５

*ｅnddoﻫj＝ｊ+２30１－(60-40)＊25ﻫ*eｎddoﻫ!连接左斜线１ﻫj=Nmax１+１９＊25＋1ﻫ＊ｄo,n，1,（94＊25＋1）＊2,9４＊25+1

kｋ=0ﻫ＊dｏ,i，ｎ+2１＊2５，ｎ－１+41*2５，25ﻫlstr，i，jﻫkk＝kk＋１

＊ｉｆ，kk,ｅq,1,thｅn

laｔｔ,1,1,１，,，，２１

eｓｉze，，2５ﻫｌmesｈ，allﻫ＊elseif，kｋ,ｇe,13,and，kk，le，18,then

latt,1,1,1,，，，２1

esｉze,,2５

lmｅsh，allﻫ*eｌｓeﻫlａtｔ,1,1，1,，，,22ﻫｅｓiｚe,,２5

lmesh,all

*ｅndｉfﻫj=j＋2５ﻫ*eｎddoﻫｊ=j+2301-(４１-２1)*25ﻫ＊enddoﻫ！连接左斜线２

ｊ=Nｍax1＋59*2５＋1ﻫ＊ｄo,ｎ,１，（９４＊２5+1)*2,９４*2５+1ﻫkk＝０ﻫ＊do，i,n+６1＊25，n-1+９5＊25,２5

lstr，i，jﻫkｋ=kk＋1ﻫ*if,kk，gｅ,1,aｎd,kk，le，12，then

latt,1,1，1,，，，22ﻫesizｅ，，25ﻫlmesh,ａll

＊ｅlsｅiｆ，ｋk,ge，２1，ａｎd,kk，le，２２,thenﻫｌatt,1,1，１，,，，２2

esizｅ，,25ﻫlmesh，allﻫ*ｅlsｅiｆ,kｋ，ge，15，aｎd，kk，le，１6,tｈenﻫlaｔｔ,1，1，１,，，，23ﻫesｉze，,25

lｍesh,ａll

＊elsｅif,kk,ge，24,and,kk，ｌe,2８，then

ｌａｔｔ，１,1,1,,,，２3

ｅsize,，２５

ｌmesh,allﻫ*elｓeif，ｋk，ge,３3,ａnｄ，ｋk,lｅ，３4,tｈenﻫlatt,1，1，1，,，,23

eｓizｅ，，２5ﻫlｍesｈ，ａll

*ｅlｓe

latt,１，１，1，,,，２１

eｓｉｚe，,2５

ｌmesh,alｌﻫ*endifﻫｊ=ｊ+2５

＊enddo

j＝ｊ+23０1－(９5-６１)*25

*enddo

！连接ｓｈｕ线ﻫj＝Nmａx1＋1ﻫ*do，n,１，（94＊2５＋1)*2,94＊25+1ﻫkk=0

*dｏ，ｉ,n+1*25，ｎ-1＋94*2５，２5ﻫlstｒ，i，ｊ

kｋ＝ｋk+１

＊if,kk，gｅ,1,and,ｋk，le,3，ｔhenﻫlａｔt,1,１，1，,,,22ﻫｅsｉze，，2５

ｌmｅsh,allﻫ*elseiｆ,ｋk，ge，22，and,kｋ，ｌe,58,thenﻫlａtt，１,1,1，，，，22

ｅｓize，，2５ﻫlｍｅsh，alｌﻫ*elｓeiｆ,kk，ge，62，and，kｋ,ｌｅ,74,tｈen

ｌａtｔ，1，１，1，，，，２2ﻫeｓize，,２5ﻫlmesh，ａll

＊elsｅiｆ,ｋk，ｇe，78，anｄ，ｋk，le，８3,thｅｎﻫｌａtｔ，1，1，1，,,,22

esizｅ,，25ﻫlmｅsｈ,ａlｌﻫ＊elｓeif，kk，gｅ，89,ａnｄ，kｋ，ｌe,9３，theｎﻫlatｔ，1,1，１，,，，22

esize,，25

ｌmesh,allﻫ＊eｌsｅif,kk,eq，7，tｈｅｎﻫｌａtt,1，，1,,,,23ﻫesiｚe，,２５ﻫlｍeｓh，allﻫ*ｅｌsｅif，kｋ,eq，60,ｔhｅn

latt，１,1，１,，,，23ﻫeｓize，，2５

ｌmｅsh，allﻫ*eｌseif，ｋk,ｅq，75，ｔhen

latt,1，1,1,,,,2３ﻫｅsize，，2５

lmesh，ａｌlﻫ＊ｅlsｅif，ｋk，eq,8６,thenﻫlaｔt，1，1,１，，,，23

eｓiｚe，，25ﻫｌmｅｓh，ａlｌﻫ*elseif,kk，eｑ，８７,tｈenﻫlａｔt，１，1,1，，，,23

ｅsize，,25

lmesh,aｌl

*ｅlse

latt，1，1,1，，,，21ﻫｅsizｅ,，25ﻫlmesh,allﻫ*enｄifﻫj＝ｊ+2５

*endｄｏﻫｊ＝j+２301－９3*25

＊enddoﻫｍｏdmsｈ,dｅtack

ｌdele，all

kdeｌe,all

！*＊＊**＊*＊**＊**＊*＊＊*＊＊＊*＊****＊*＊*＊＊＊＊＊**＊＊**＊*＊

！

!

!2X+上平联ﻫ!

！ﻫ!***＊＊＊*＊**＊＊****＊＊*＊*＊＊*＊*＊＊＊*****＊＊＊＊***＊＊＊*ﻫ／prep７

k,1，１4.5，-14,１5.0０

k，2，14.5,—7，１5.００ﻫk，3，14。５，—１４，10.2

ｋ,4,14。5,－７，1０.2

lstｒ,１,4

lstｒ,2,３

lｓtr，３，4ﻫlｇen，２,all，，，，７,

ｎummrｇ，kp

numcmp，kpﻫｌstr,2，4ﻫlgen，９３，all，,，14。5，,

lｇeｎ，2,ａll,，,,1４，，

latt,1,，１,,,,24ﻫeｓize，，１0ﻫｌmesｈ,ａｌlﻫｍodｍsｈ，dｅtaｃk

ldele，aｌl

ｋｄele,ａll

＊SEＴ,ｊ,—１4

ﻫ*do,i，1,(92＊25＋1）＊2,92＊２5+1ﻫｋ,i，１4.50，ｊ，1５。00

ﻫk，９2*25+i，１3４8．50,j，1５.00

kFILL，ｉ,92*2５+i,９1，,，，，ﻫj=j+14

＊endｄo

*do,i，1，9２＊25+1,2５

lsｔr,i,ｉ+９２*２5+1

lstr，ｉ+9２＊2５+1，i+（９2*25+1）＊２

*enddo

lａｔｔ，１，,1,，,，31ﻫesiｚe,，２0

lｍesh，ａｌｌﻫmodmsｈ,ｄetacｋ

ｌdeｌe，all

＊do，i，１,92＊25，25

ｌｓｔr,i，i+92＊25+1+1＊25

lｓtr,i＋92*25+1，i＋（92*25+1)*2+1*25ﻫ*ｅnｄdｏ

＊dｏ,i,1+1*2５,９2＊25+1，25ﻫlsｔｒ，i，ｉ+9２＊25+1－25

lsｔr，i+９2＊2５+1，ｉ+(92*２５+１）*２－25ﻫ*enddo

nuｍmrg,kpﻫnuｍｃmp，kp

laｔt，1,,１，，，，32ﻫｅsize，,20

ｌmesh,allﻫmoｄｍsh，ｄeｔack

ｌdelｅ，all

kｄｅｌｅ,allﻫ!＊＊*＊＊＊＊＊*＊**＊＊***＊＊**＊＊**＊＊＊*＊＊***＊＊**＊*＊＊＊＊＊*＊*！

！ﻫ！1#ta

！ﻫ!ﻫ!＊*＊＊＊＊**＊＊＊＊＊*＊*＊＊*＊＊＊*＊*＊**＊＊*＊＊＊***＊＊＊**＊***＊＊ﻫ/ｐrep7

nn＝４０0０0０

ｎ，1+nｎ，2９０，0，158．９

ｎ，１3＋nn,２９０，0,132．５ﻫn，２0+ｎｎ,2９0,0,115

n，22＋nｎ，２９0，0,1１0。５ﻫn,６3+nn，2９0，0,１1０.５

n，５２＋nn，290,-19，－6

n,６2＋ｎn,2９0，-11.5，－34．5

n，9３＋nn，２９0,１９，-6ﻫn，103＋nn,29０，11.5,—3４.5ﻫn，1０4+nn，290,—19，-６ﻫn,107+nn,290，19，－6

fiｌl，1+ｎｎ,13+nn，11，ﻫｆill，13+nn，20+nn,６，ﻫｆiｌl，2０+ｎn,2２+ｎn,1，ﻫfill,22+nn,52+nn，29，ﻫfiｌl，52+nn，62+ｎn,9，ﻫｆｉll，6３＋nn,9３+ｎn,29,ﻫfiｌl，9３＋ｎn，10３+nｎ，9，ﻫfｉll,104＋ｎｎ,1０７+nn，２,

ｒeal，1

tｙpe，1

ﻫmａt,２

ﻫsecｎum,53

＊ｄo,a,1+ｎｎ，12+nn，1ﻫE,a,a+1ﻫ＊endｄo

tｙpe,1

ﻫmat，2

ﻫｓecnum，5２

＊do,ｂ，1３+nn，１９＋nn,１

E,b，b＋1

*ｅnｄdoﻫtｙpｅ，１

mａｔ,２

ﻫsecｎｕm，５1ﻫ*dｏ，ｂ,20+ｎn,21+nn，1ﻫE,b,b＋１

＊ｅｎddoﻫtype，１

ﻫｍaｔ,2

ﻫseｃnum，46

＊do，b,22+nn,51+nn，1ﻫE，b,ｂ+1

ﻫ*enddoﻫtyｐｅ,1

ﻫmat，2

ﻫsｅcnｕｍ,４2

＊do，b,５2+nn，6１＋nn,1

E,ｂ，b+１

＊enｄdo

type，１

ﻫmａｔ，2

ﻫsecnｕm,46ﻫ＊do,b,63+ｎn,92＋nn，1ﻫＥ，b，b+１

ﻫ＊enddｏﻫｔype，１

mａt,2

ｓecｎum，42

＊dｏ,b，9３+nn，102+nｎ，1

Ｅ,ｂ,b+1

＊enddoﻫtyｐe,1

ﻫmaｔ，２

ﻫsecｎum，5４ﻫ＊do,b，104＋nn，106＋nn，1ﻫE,b,ｂ+1

ﻫ*enddo

！***＊**＊＊*＊*＊*＊＊＊*＊****＊*****＊***＊*＊＊＊*＊＊*＊*＊*＊＊＊ﻫ！

！2＃tａ

!

！**＊＊＊＊***＊＊**＊*＊＊＊*＊*＊*＊＊*＊＊＊＊*＊*＊*＊＊＊*＊*****＊*＊＊ﻫｎn＝4５0０00ﻫn,1+ｎn,870，0，１58。９

n，１3＋ｎｎ，８7０,０,132．5ﻫn,2０＋nｎ，870,0，１15

n,２２+nn，８70,0，11０。5

n，6３+ｎn，87０,0,11０．5

n，52+nn，8７0,—19，—６

ｎ，62+nn,87０，－１1。5，-3４.5ﻫn,93+nｎ,8７0,19，-６ﻫn，103+nn，8７0,１1。５，—34.5

n,1０4＋nn，870，-19,-６

ｎ,1０７+nn，870,１9,—6

ｆiｌl,１+ｎn，1３＋nn，11，ﻫfill，1３+nｎ，20+nｎ,6,

ｆill,2０+ｎn,22＋ｎn,１,

fｉll,2２+nn,52＋nn,29,ﻫｆill,52＋nn,62＋nｎ，9，ﻫfill，63＋nn，９3+nn，２９，

fｉll，9３+ｎn，１03+nn，9，ﻫfｉll,1０4+ｎn，107＋nｎ，２,

real，1ﻫtyｐe,1

ﻫmat,２

secｎum，53

ﻫ＊ｄｏ，a,1+nn，12+nn，1ﻫE,a，a+1

＊endｄoﻫtypｅ，1

ﻫmａt，2

secｎum，５2

＊ｄo,b，13+ｎn，１9+nn，1ﻫE,b，ｂ+１

＊enｄｄo

type，１

ｍat，２

secnum，5１

＊do，b，20+nn,２1+nn,1

E，b，b+１

*eｎddo

tｙpe，1

ﻫmaｔ,2

ｓecｎum,４6ﻫ＊do，b,２２+nn，５1＋nn，１

Ｅ,b，b+1

＊enｄdoﻫtyｐe，1

ﻫmat，2

ﻫsｅcｎum，42ﻫ＊do，b，５２+ｎｎ，６1+nｎ，1

E,b，b+１

ﻫ*enddoﻫtypｅ,1

mat,2

secnum，46

＊dｏ，b,63+nn,92＋nn，1

E,ｂ，b＋1

ﻫ＊endｄｏ

ｔyｐe,1

maｔ，2

secnｕｍ,42

＊ｄｏ,b，９3+nn,1０2＋ｎn,1

E，b，b+１

ﻫ*ｅnｄｄｏﻫtype,１

mat，2

secnuｍ，５４ﻫ*do,b,1０4＋nn，１06+nn，1ﻫE，b，ｂ+1

*enddo

!*＊＊****＊*＊*＊＊＊＊＊＊＊＊＊＊*＊＊**＊＊＊＊**＊*＊*＊**＊*