下载文档原格式
ANSYS刚桁架桥模态分析命令流/TITLE,Truss Bridge Static Analysis !指定标题/COM,Structural !选择分析类型为结构分析/PREP7 !进入前处理器ET,1,BEAM188 !定义1号单元类型KEYOPT,1,3,2ET,2,SHELL181 !定义2号单元类型KEYOPT,1,3,2SECTYPE, 1, BEAM, I, duan, 0 !定义1号工字形截面SECOFFSET, CENT !截面质心不偏移SECDATA,0.4,0.4,0.4,0.016,0.016,0.016,0,0,0,0 !1号截面参数SECTYPE, 2, BEAM, I, XIANHENG, 0 !定义2号工字形截面SECOFFSET, CENT ! 截面质心不偏移SECDATA,0.4,0.4,0.4,0.012,0.012,0.012,0,0,0,0 ! 2号截面参数SECTYPE, 3, BEAM, I, FU, 0 ! 定义3号工字形截面SECOFFSET, CENT ! 截面质心不偏移SECDATA,0.3,0.3,0.4,0.012,0.012,0.012,0,0,0,0 ! 3号截面参数R,1,0.0187,0.17e-3,0.54e-3,0.4,0.4,0, !设置1号实常数R,2,0.0141,0.128e-3,0.415e-3,0.4,0.4,0, !设置2号实常数R,3,0.0117,0.541e-4,0.324e-3,0.3,0.4,0, !设置3号实常数R,4,0.3, , , , , , !设置4号实常数MP,EX,1,2.1E11 !定义1号材料弹性模量MP,PRXY,1,0.3 !定义1号材料泊松比MP,DENS,1,7850 !定义1号材料密度MP,EX,2,3.5E10 !定义2号材料弹性模量MP,PRXY,2,0.1667 !定义2号材料泊松比MP,DENS,2,2500 !定义2号材料密度N, ,0,0,-5, , , , !建立节点NGEN,4,4,ALL, , ,12, , ,1, !复制节点NGEN,2,1,ALL, , , , ,10,1, !复制节点NGEN,2,1,2,10,4, ,16, ,1, !复制节点NGEN,2,1,3,11,4, , ,-10,1, !复制节点TYPE,1 !选择1号单元类型MAT,1 !选择1号材料REAL,1 ! 选择1号实常数ESYS,0 !单元坐标系SECNUM,1 !选择1号截面TSHAP,LINE !选择线形单元E,11,14 !建立单元E,12,13 !建立单元TYPE,1 ! 选择1号单元类型MAT,1 ! 选择1号材料REAL,2 ! 选择2号实常数ESYS,0 ! 单元坐标系SECNUM,2 ! 选择2号截面TSHAP,LINE ! 选择线形单元E,2,6 !建立单元E,6,10E,10,14E,1,5E,5,9E,9,13E,3,7E,7,11E,4,8E,8,12E,1,2E,3,4E,5,6E,7,8E,9,10E,11,12E,13,14TYPE,1 ! 选择1号单元类型MAT,1 ! 选择1号材料REAL,3 ! 选择3号实常数ESYS,0 ! 单元坐标系SECNUM,3 ! 选择3号截面TSHAP,LINE ! 选择线形单元E,3,6 !建立单元E,6,11E,4,5E,5,12E,2,3E,1,4E,6,7E,5,8E,10,11E,9,12TYPE,2 ! 选择2号单元类型MAT,2 ! 选择2号材料REAL,4 ! 选择4号实常数ESYS,0 ! 单元坐标系TSHAP,QUAD ! 选择四边形单元E,1,2,6,5 !建立单元E,5,6,10,9E,9,10,14,13NSYM,X,14,ALL !所有节点以y0z平面对称ESYM, ,14,ALL ! 所有单元以y0z平面对称NUMMRG,ALL, , , ,LOW !合并重复节点单元，编号取较小者NUMCMP,ALL !压缩节点单元等编号FINISH !结束前处理器/SOL !进入求解器NSEL,S, , ,23,24 !选择节点D,ALL, ,0, , , ,UX,UY,UZ, , , , !约束3个自由度NSEL,S, , ,13,14 !选择节点D,ALL, ,0, , , , ,UY,UZ, , , , !约束两个自由度ALLSEL,ALL !选择所有ANTYPE,2 !选择分析类型，模态分析MODOPT,SUBSP,6,0,100, ,OFF !选择字空间法，提取6阶模态SUBOPT,8,4,10,0,0,ALL !字空间法设置SOLVE !开始求解FINISH !结束求解器/POST1 !进入普通后处理器SET,LIST !列表格阶模态SET,FIRST !读取第一阶模态PLNSOL, U,SUM, 2,1.0 !显示振型云图SET,NEXT !读取下一阶模态PLNSOL, U,SUM, 2,1.0 !显示振型云图SET,NEXT !读取下一阶模态PLNSOL, U,SUM, 2,1.0 !显示振型云图SET,NEXT !读取下一阶模态PLNSOL, U,SUM, 2,1.0 !显示振型云图SET,NEXT !读取下一阶模态PLNSOL, U,SUM, 2,1.0 !显示振型云图SET,NEXT !读取下一阶模态PLNSOL, U,SUM, 2,1.0 !显示振型云图FINISH !结束普通后处理器! /EXIT,ALL !退出程序。
钢桁架桥静力受力分析对一架钢桁架桥进行具体静力受力分析，分别采用GUI方式和命令流方式。
A 问题描述图6-15 钢桁架桥简图已知下承式简支钢桁架桥桥长72米，每个节段12米，桥宽10米，高16米。
设桥面板为0.3米厚的混凝土板，当车辆行驶于桥梁上面时，轴重简化为一组集中力作用于梁上，来计算梁的受力情况。
桁架杆件规格有三种，见下表：杆件截面号形状规格端斜杆 1 工字形400×400×16×16上下弦 2 工字形400×400×12×12 横向连接梁 2 工字形400×400×12×12其他腹杆 3 工字形400×300×12×12所用材料属性如下表：表6-3 材料属性参数钢材混凝土弹性模量EX 2.1×1011 3.5×1010泊松比PRXY 0.3 0.1667密度DENS 7850 2500B GUI操作方法1.创建物理环境1）过滤图形界面：GUI：Main Menu> Preferences，弹出“Preferences for GUI Filtering”对话框，选中“Structural”来对后面的分析进行菜单及相应的图形界面过滤。
2）定义工作标题：GUI：Utility Menu> File> Change Title，在弹出的对话框中输入“Truss Bridge StaticAnalysis”，单击“OK”。
如图6-16（a）。
指定工作名：GUI：Utility Menu> File> Change Jobname，弹出一个对话框，在“Enter new Name”后面输入“Structural”，“New log and error files”选择yes，单击“OK”。
如图6-16（b）。
图6-16（a）定义工作标题图6-16（b）指定工作名3）定义单元类型和选项：GUI：Main Menu> Preprocessor> Element Type> Add/Edit/Delete，弹出“Element Types”单元类型对话框，单击“Add”按钮，弹出“Library of Element Types”单元类型库对话框。
（完整版）ANSYS最常用命令流+中文注释（超级大全）ANSYS最常用命令流+中文注释VSBV, NV1, NV2, SEPO, KEEP1, KEEP2 —Subtracts volumes from volumes，用于2个solid相减操作，最终目的是要nv1-nv2=?通过后面的参数设置，可以得到很多种情况：sepo项是2个体的边界情况，当缺省的时候，是表示2个体相减后，其边界是公用的，当为sepo的时候，表示相减后，2个体有各自的独立边界。
keep1与keep2是询问相减后，保留哪个体？当第一个为keep时，保留nv1,都缺省的时候，操作结果最终只有一个体，比如：vsbv,1,2,sepo,,keep,表示执行1-2的操作，结果是保留体2，体1被删除，还有一个1-2的结果体，现在一共是2个体（即1-2与2），且都各自有自己的边界。
如vsbv,1,2,,keep,,则为1-2后，剩下体1和体1-2，且2个体在边界处公用。
同理，将v换成a 及l是对面和线进行减操作！mp,lab, mat, co, c1,…….c4 定义材料号及特性lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens）ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）co: 材料特性值，或材料之特性，温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数定义DP材料：首先要定义EX和泊松比：MP，EX，MA T，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MA T进入单元表并编辑添加单元表：TBDATA，1，CTBDATA，2，ψTBDATA，3，……如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,degVSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Type，是选择的方式，有选择（s）,补选（a），不选(u)，全选(all)、反选（inv）等,其余方式不常用Item, Comp 是选取的原则以及下面的子项如volu 就是根据实体编号选择，loc 就是根据坐标选取，它的comp就可以是实体的某方向坐标！其余还有材料类型、实常数等MIN, VMAX, VINC，这个就不必说了吧！,例：vsel,s,volu,,14vsel,a,volu,,17,23,2上面的命令选中了实体编号为14，17，19，21，23的五个实体VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体nv1:初始体号nv2:最终的体号ninc:体号之间的间隔kswp=0:只删除体kswp=1:删除体及组成关键点,线面如果nv1=all,则nv2,ninc不起作用其后面常常跟着一条显示命令VPLO,或aplo,nplo,这个湿没有参数的命令，输入后直接回车，就可以显示刚刚选择了的体、面或节点，很实用的哦！Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0” 使用正负号“1”仅用绝对值下面是单元生死第一个载荷步中命令输入示例：!第一个载荷步TIME,... !设定时间值（静力分析选项）NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿－拉夫森选项ESTIF,... !设定非缺省缩减因子（可选）ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点（不与活动单元相连的结点）D,ALL,ALL,0 !约束所有不活动的结点自由度（可选）NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SA VESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D, F,SF和BF命令得到更详细的解释。
!!连续梁桥/prep7et,1,4 !!!!定义梁单元et,2,21 !!!!定义mass21单元!!定义粱材料!!泊松比!!密度mp,ex,2,3.45e10 !!直线段梁材料和1M段梁材料mp,nuxy,2,0.2mp,dens,2,3302.153125mp,ex,3,3.45e10mp,nuxy,3,0.2mp,dens,3,3301.658695mp,ex,4,3.45e10mp,nuxy,4,0.2mp,dens,4,3299.906778mp,ex,5,3.45e10mp,nuxy,5,0.2mp,dens,5,3298.327219mp,ex,6,3.45e10mp,nuxy,6,0.2mp,dens,6,3292.351605mp,ex,7,3.45e10mp,nuxy,7,0.2mp,dens,7,3284.137255mp,ex,8,3.45e10mp,nuxy,8,0.2mp,dens,8,3271.802136mp,ex,9,3.45e10mp,nuxy,9,0.2mp,dens,9,3260.41903mp,ex,10,3.45e10mp,nuxy,10,0.2mp,dens,10,3248.193657mp,ex,11,3.45e10mp,nuxy,11,0.2mp,dens,11,3235.117644mp,ex,12,3.45e10mp,nuxy,12,0.2mp,dens,12,3221.585664mp,ex,13,3.45e10mp,nuxy,13,0.2mp,dens,13,3208.826871mp,ex,14,3.45e10mp,nuxy,14,0.2mp,dens,14,3194.279207mp,ex,15,3.45e10mp,nuxy,15,0.2mp,dens,15,3179.924673mp,ex,16,3.45e10mp,nuxy,16,0.2mp,dens,16,3166.445716mp,ex,17,3.45e10mp,nuxy,17,0.2mp,dens,17,3152.555731mp,ex,18,3.45e10mp,nuxy,18,0.2mp,dens,18,3138.312105mp,ex,19,3.45e10mp,nuxy,19,0.2mp,dens,19,3124.795334mp,ex,20,3.45e10mp,nuxy,20,0.2mp,dens,20,3110.7135mp,ex,21,3.45e10mp,nuxy,21,0.2mp,dens,21,3097.080875mp,ex,22,3.45e10mp,nuxy,22,0.2mp,dens,22,3083.186268mp,ex,23,3.45e10mp,nuxy,23,0.2mp,dens,23,3068.968824mp,ex,24,3.45e10mp,nuxy,24,0.2mp,dens,24,3055.612436mp,ex,25,3.45e10mp,nuxy,25,0.2mp,dens,25,3045.857147mp,ex,26,3.45e10mp,nuxy,26,0.2mp,dens,26,3035.174287mp,ex,27,3.45e10mp,nuxy,27,0.2mp,dens,27,3026.696551mp,ex,28,3.45e10mp,nuxy,28,0.2mp,dens,28,3015.795365mp,ex,29,3.45e10mp,nuxy,29,0.2mp,dens,29,3007.710181mp,ex,30,3.45e10mp,nuxy,30,0.2mp,dens,30,3000.513837mp,ex,31,3.45e10mp,nuxy,31,0.2mp,dens,31,2978.611375 mp,ex,32,3.45e10mp,nuxy,32,0.2mp,dens,32,2958.618861mp,ex,33,3.45e10mp,nuxy,33,0.2mp,dens,33,2937.888072mp,ex,34,3.45e10mp,nuxy,34,0.2mp,dens,34,2919.475751mp,ex,35,3.45e10mp,nuxy,35,0.2mp,dens,35,2903.359983 !!6700处mp,ex,36,3.45e10mp,nuxy,36,0.2mp,dens,36,3302.153125 !!可以不用mp,ex,37,3.45e10mp,nuxy,37,0.2mp,dens,37,3302.153125mp,ex,38,3.45e10 !!-700处梁mp,nuxy,38,0.2mp,dens,38,3180.578901mp,ex,39,3.45e10 !!边支点横隔板mp,nuxy,39,0.2mp,dens,39,2868.674818mp,ex,40,3.45e10 !!合拢段横隔板mp,nuxy,40,0.2mp,dens,40,2868.674818mp,ex,41,3.45e10 !!中支点横隔板mp,nuxy,41,0.2mp,dens,41,2757.470588mp,ex,1,3.25e10 !!!!定义墩材料属性mp,nuxy,1,0.2mp,dens,1,2650!!定义实常数编号，面积，IYY，IZZ，宽，高，，RMORE，，抗扭惯距!!主梁截面r,1,11.851,28.52,215.151,16.95,4,, !!直线段rmore,,56.8972r,2,11.862,28.57,215.458,16.95,4.001,, !!截面100处rmore,,56.905r,3,11.871,28.91,215.622,16.95,4.009,,rmore,,57.345r,4,11.903,29.03,215.78,16.95,4.025,,rmore,,57.938r,5,11.932,29.626,215.88,16.95,4.049,,rmore,,58.502r,6,12.043,30.402,217.093,16.95,4.081,,rmore,,59.768r,7,12.199,31.407,218.899,16.95,4.121,,rmore,,61.445r,8,12.441,32.716,221.806,16.95,4.169,,rmore,,63.652r,9,12.673,34.21,224.418,16.95,4.226,,rmore,,66.058r,10,12.932,35.94,227.279,16.95,4.290,,rmore,,68.801r,11,13.221,37.941,230.512,16.95,4.362,,rmore,,71.948r,12,13.534,40.233,233.872,16.95,4.442,,rmore,,75.432r,13,13.843,42.747,237.133,16.95,4.530,,rmore,,79.173r,14,14.213,45.688,241.054,16.95,4.627,,rmore,,83.484r,15,14.598,48.964,245.004,16.95,4.731,,rmore,,88.155r,16,14.979,52.547,248.817,16.95,4.843,,rmore,,93.111r,17,15.393,56.562,252.927,16.95,4.963,,rmore,,98.528r,18,15.842,61.069,257.347,16.95,5.092,,rmore,,104.47r,19,16.293,65.992,261.673,16.95,5.228,,rmore,,110.742r,20,16.791,71.484,266.426,16.95,5.372,,rmore,,117.648r,21,17.303,77.568,271.225,16.95,5.525,, rmore,,124.861r,22,17.858,84.282,276.363,16.95,5.685,, rmore,,132.631r,23,18.464,91.768,281.923,16.95,5.854,, rmore,,141.068r,24,19.072,99.936,287.364,16.95,6.03,, rmore,,149.875r,25,19.524,108.482,291.159,16.95,6.214,, rmore,,158.46r,26,20.084,118.046,295.597,16.95,6.407,, rmore,,167.857r,27,20.536,128.202,298.9,16.95,6.607,, rmore,,177.137r,28,21.148,139.669,303.86,16.95,6.816,, rmore,,187.495r,29,21.626,151.737,307.186,16.95,7.032,, rmore,,197.507r,30,22.07,164.734,310.057,16.95,7.257,, rmore,,207.681r,31,23.541,182.225,323.433,16.95,7.490,, rmore,,223.527r,32,25.066,201.932,336.467,16.95,7.730,, rmore,,241.046r,33,26.871,223.700,351.448,16.95,7.979,, rmore,,257.765r,34,28.707,247.928,365.663,16.95,8.235,, rmore,,275.922r,35,30.533,274.36,378.934,16.95,8.500,, rmore,,294.461r,46,14.58,35.599,238.702,16.95,4,, rmore,,238.702!!横隔板截面!!边支点r,36,35.376,54.393,339.84,16.95,4.00,, rmore,,128.747!!合拢段r,37,35.376,53.405,339.84,16.95,4.00,, rmore,,127.015!!中支点r,38,71.981,454.871,543.871,16.95,8.5,, rmore,,684.565!!主墩截面!!截面1（实心）r,39,39.932,82.876,212.138,8,5,, rmore,,203.709!!截面（1/2空心）r,40,23.863,72.571,168.024,8,5,, rmore,,168.491!!截面3(空心）r,41,14.236,52.366,114.529,8,5,, rmore,,116.391!!边墩截面!!截面1（实心1）r,42,36.113,66.217,177.703,7.7,4.7,, rmore,,165.119!!截面2（空心1/2）r,43,20.753,35.743,99.711,7,4,, rmore,,88.486!!截面3（空心）r,44,9.103,16.082,45.911,6.5,3.5,, rmore,,39.39!!截面4（实心2）r,45,22.673,23.005,79.361,6.5,3.5,, rmore,,61.527!!截面5（1/2）空心2r,49,16,21.958,68.708,6.5,3.5,, rmore,,55.682!!建立第一跨梁节点xl=atan(1.029/80) n,1,0,n,2,2.25,2.25*xln,3,6,6*xln,4,9,9*xln,5,11,11*xln,6,13,13*xl-0.0045 n,7,15,15*xl-0.0135 n,8,17,17*xl-0.023 n,9,19,19*xl-0.038 n,10,21,21*xl-0.062 n,11,23,23*xl-0.092 n,12,25,25*xl-0.1275 n,13,27,27*xl-0.167 n,14,29,29*xl-0.2115 n,15,31,31*xl-0.261 n,16,33,33*xl-0.314 n,17,35,35*xl-0.3715 n,18,37,37*xl-0.435 n,19,39,39*xl-0.502 n,20,41,41*xl-0.5725 n,21,43,43*xl-0.648 n,22,45,45*xl-0.728 n,23,47,47*xl-0.8115 n,24,49,49*xl-0.8995 n,25,51,51*xl-0.9925 n,26,53,53*xl-1.0895 n,27,55,55*xl-1.1905 n,28,57,57*xl-1.325 n,29,59,59*xl-1.433 n,30,61,61*xl-1.5155 n,31,63,63*xl-1.6325 n,32,65,65*xl-1.7535 n,33,67,67*xl-1.879 n,34,69,69*xl-2.007 n,35,71,71*xl-2.1325 n,36,73,73*xl-2.2635 n,37,75,75*xl-2.399 n,38,77,77*xl-2.533 n,39,79,79*xl-2.631xl2=atan(0.7/140)local,11,0,80,80*xl-2.631,0,xl2nsym,x,35,5,39,1 !!复制粱结点从5到39结点编号增加35 csys,0n,75,149.6,0.15852local,12,0,150,0.15852,0,xl2nsym,x,36,40,75,1local,13,0,220,-1.6027,0nsym,x,111,1,111,1csys,0n,500,80,80*xl-2.631n,501,220,-1.6027n,502,360,-1.6148csys,0 !!返回普通坐标!!建立墩结点!!PM112边墩截面n,300,0,-2.165 !!(实心）2米长n,301,0,-4.365n,302,0,-9.365 !!（1/2空心）5M长n,303,0,-13.365n,304,0,-17.365n,305,0,-21.365n,306,0,-24.93n,307,0,-28.93n,308,0,-30.63!!PM113主敦截面n,309,80,0.5-6.665-0.8n,310,80,0.5-3.5-6.665 !!(实心)3.5mn,311,80,0.5-8.5-6.665 !!(1/2空心）5mn,312,80,0.5-12.5-6.665 !!(空心）4Mn,313,80,0.5-16.5-6.665 !!n,314,80,0.5-18.594-6.665n,315,80,0.5-22.594-6.665n,316,80,0.5-24.594-6.665+0.4!(实心）2m!!pm114主墩截面n,350,220,0.5-6.665-0.8 !n,351,220,0.5-3.5-6.665 !!(1/2空心)5m n,352,220,0.5-8.5-6.665 !!(空心）5mn,353,220,0.5-12.5-6.665 !!(空心）6.064M n,354,220,0.5-16.5-6.665 !!(1/2空心）4M n,355,220,0.5-19.294-6.665 !!(实心）2M n,356,220,0.5-23.294-6.665n,357,220,0.5-25.294-6.665+0.4local,13,0,220,-1.6027,0nsym,x,17,300,316,1csys,0!!生成边跨1mat,39real,36e,1,2mat,38real,46e,2,3mat,2real,1e,3,4e,4,5mat,2real,2e,5,6*do,i,3,35,1mat,ireal,ie,3+i,4+i*enddo*do,i,3,35,1mat,ireal,ie,38+i,39+i*enddo*do,i,3,35,1 mat,ireal,ie,74+i,75+i *enddo*do,i,3,35,1 mat,ireal,ie,114+i,115+i *enddo*do,i,3,35,1 mat,ireal,ie,149+i,150+i *enddo*do,i,3,35,1 mat,ireal,ie,185+i,186+i *enddomat,41 real,38e,39,500mat,41 real,38e,500,74mat,2real,2e,41,40mat,2real,1e,40,75mat,40 real,37 e,75,111mat,2 real,1 e,111,76mat,2 real,2 e,76,77mat,41 real,38 e,110,501mat,41 real,38 e,501,221mat,2 real,2 e,188,187mat,2 real,1 e,187,222mat,41 real,37 e,222,186mat,2 real,1 e,186,151mat,2 real,2 e,151,152mat,41 real,38 e,185,502mat,41 real,38 e,502,150mat,39 real,36 e,112,113mat,38 real,46 e,113,114mat,2 real,1e,114,115 e,115,116mat,2 real,2e,116,117!!建立桥墩!!边墩1 mat,1 real,42 e,300,301mat,1 real,43 e,301,302mat,1 real,44 e,302,303mat,1 real,44e,303,304mat,1 real,44 e,304,305mat,1 real,44 e,305,306mat,1 real,49 e,306,307mat,1 real,45 e,307,308!!主墩1 mat,1 real,39 e,309,310mat,1 real,40 e,310,311mat,1 real,41 e,311,312mat,1 real,41 e,312,313mat,1 real,41 e,313,314mat,1 real,40 e,314,315mat,1 real,39 e,315,316!!主墩2 mat,1 real,39 e,350,351mat,1 real,40 e,351,352mat,1 real,41 e,352,353mat,1 real,41 e,353,354mat,1 real,41 e,354,355mat,1 real,40 e,355,356mat,1 real,39 e,356,357!!主墩3 mat,1 real,39 e,326,327mat,1 real,40 e,327,328mat,1 real,41 e,328,329mat,1 real,41 e,329,330mat,1 real,41 e,330,331mat,1 real,40 e,331,332mat,1 real,39 e,332,333!!边墩2 mat,1 real,42 e,317,318mat,1 real,43 e,318,319mat,1 real,44 e,319,320mat,1 real,44 e,320,321mat,1 real,44 e,321,322mat,1real,44e,322,323mat,1real,49e,323,324mat,1real,45e,324,325!!添加质量惯距!!两米段质量惯距r,100,,,,1293348.4 !!100米处和直线段的质量惯距r,101,,,,1296019.6r,102,,,,1297493r,103,,,,1301181.8r,104,,,,1311723r,105,,,,1326622r,106,,,,1348967r,107,,,,1370728r,108,,,,1395061r,109,,,,1422801r,110,,,,1452757r,111,,,,1483258r,112,,,,1519733r,113,,,,1558030r,114,,,,1597229r,115,,,,1640292r,116,,,,1687605r,117,,,,1736625r,118,,,,1790923r,119,,,,1848602r,120,,,,1911418r,121,,,,1980562r,122,,,,2052690r,123,,,,2118097r,124,,,,2192307r,125,,,,2263641r,126,,,,2350704r,127,,,,2432292r,128,,,,2516392r,129,,,,2679987r,130,,,,2853514r,131,,,,3048284r,132,,,,3252032r,133,,,,3462506!!边跨直线3米段r,148,,,,1293348.4*3/2!!-700处3.75米段的质量惯距r,134,,,,2725866.2!!边支点横隔r,135,,,,2246142.5!!合龙段横隔r,136,,,,833679/2!!中支点横隔r,137,,,,5292760!!添加边墩质量惯距从上到下!!变截面实心2米段r,138,,,,1163498!!1/2空心5米段r,139,,,,1794766!!空心段4米r,140,,,,657211!!1/2空心4米段r,141,,,,1085079!!实心2米段r,142,,,,1064567!!添加主墩质量惯性距从上到下!!实心3米段r,143,,,,2345361!!1/2空心5米段r,144,,,,2869095!!空心4米段r,145,,,,1769087!!1/2空心4米段r,146,,,,2869095*4/5!!实心2米段r,147,,,,2345361*2/3!!在结点处加入质量惯性距!!边跨合龙段!左边边跨type,2real,135e,2type,2real,134e,3type,2real,148e,4type,2real,100e,5!!右边边跨type,2real,135e,113type,2real,134e,114type,2real,148e,115type,2real,100e,116!!对第一个主跨的循环*do,i,6,39,1type,2real,94+ie,i*enddo!!对第二个主跨的循环*do,i,41,74,1type,2real,59+ie,i*enddo!!对第3个循环*do,i,77,110,1type,2real,33+ie,i*enddo!!对第4个循环*do,i,188,221,1 type,2real,(-88)+ie,i*enddo!!对第5个循环*do,i,152,185,1 type,2real,i-52e,i*enddo!!对第6个循环*do,i,117,150,1 type,2real,(-17)+ie,i*enddo!!添加质量惯距到桥墩!!左边墩type,2real,138e,301type,2real,139e,302type,2real,140e,303type,2real,140e,304type,2real,140e,305type,2real,140e,306type,2 real,141 e,307type,2 real,142 e,308!!右边墩type,2 real,138 e,318type,2 real,139 e,319type,2 real,140 e,320type,2 real,140 e,321type,2 real,140 e,322type,2 real,140 e,323type,2 real,141 e,324type,2real,142e,325!!左边主墩type,2real,143e,310type,2real,144e,311type,2real,145e,312type,2real,145e,313type,2r,400,,,,1769087*3/4 e,314type,2real,146e,315type,2real,147e,316!!中间主墩type,2real,143e,351type,2real,144e,352type,2real,145e,353type,2real,145e,354type,2r,400,,,,1769087*3/4 e,355type,2real,146e,356type,2real,147e,357!!右边主墩type,2real,143e,327type,2real,144e,328type,2real,145e,329type,2real,145e,330type,2r,400,,,,1769087*3/4e,331type,2real,146e,332type,2real,147e,333!!添加横隔板处质量惯性距!!中支点type,2real,137e,500type,2real,137e,501type,2real,137e,502!!合龙段1type,2real,136e,75type,2real,136e,111type,2r,411,,,,1293348.4*3/10 e,40type,2r,412,,,,1293348.4*3/10 e,76!!合龙段2type,2real,136e,222type,2real,136e,186type,2r,411,,,,1293348.4*3/10 e,187type,2r,412,,,,1293348.4*3/10 e,151!!偶合结点cp,1,uy,1,300cp,2,uz,1,300cp,3,rotx,1,300cp,4,uy,309,500cp,5,uz,309,500cp,6,rotx,309,500cp,7,ux,350,501cp,8,uy,350,501cp,9,uz,350,501cp,10,rotx,350,501 cp,11,roty,350,501cp,12,uy,326,502cp,13,uz,326,502cp,14,rotx,326,502cp,15,uy,317,112cp,16,uz,317,112cp,17,rotx,317,112!!!!加约束d,308,all,all,0d,316,all,all,0d,357,all,all,0d,333,all,all,0d,325,all,all,0allselfinish!!!!求解/soluantype,2 !模态分析acel,,9.8 !加载重力modopt,subsp,20 !提取前30阶模态mxpand,20,,,0LUMPM,1solvefinish/post1set,list。
Ansys经典命令流钢管混凝土拱桥算例/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.586RECTNG,0.35,0.364,0.8,2.2,RECTNG,0.85,0.836,0.8,2.2,asel,allaptn,allADELE,5,8,1 !删除多余的面积allselaadd,allsmrtsize,5amesh,allsecwrite,gg,sect,,1SECTYPE,1,BEAM,MESH, !截面代号1SECOFFSET,CENT,,,SECREAD,'gg','SECT',,MESH !读入空心钢管截面asel,allaclear,alladele,all,,,1 !将图形及面、线、点全部删除（只保留截面文件）/replotCYL4,0.6,0.6,0.586 !以下命令是用于建立一个空心钢管内的混凝土截面CYL4,0.6,2.4,0.586RECTNG,0.364,0.836,0.8,2.2,allselaadd,allsmrtsize,5amesh,allsecwrite,hnt,sect,,1SECTYPE,2,BEAM,MESH, !截面代号2SECOFFSET,CENT,,,SECREAD,'hnt','SECT',,MESH !读入混凝土（空心钢管内的）截面asel,allaclear,alladele,all,,,1 !将图形及面、线、点全部删除（只保留截面文件）/replot !重新生成图形k,1,-0.5,0 !以下命令是用于建立一个横梁标准截面k,2,-0.75,0.25k,3,-0.75,0.5k,4,-0.55,0.5k,5,-0.55,0.85k,6,0.55,0.85k,7,0.55,0.5k,8,0.75,0.5k,9,0.75,0.25k,10,0.5,0a,1,2,3,4,5,6,7,8,9,10 !生成横梁标准截面1 A-Asmrtsize,5amesh,allsecwrite,hl1,sect,,1SECTYPE,4,BEAM,MESH, !截面代号4SECOFFSET,CENT,,,SECREAD,'hl1','SECT',,MESH !读入横梁标准截面1 A-A asel,allaclear,alladele,all,,,1 !将图形及面、线、点全部删除（只保留截面文件）/replot !重新生成图形k,1,-0.5,0 !以下命令是用于建立一个横梁标准截面b-bk,2,-0.5,1k,3,-0.75,1.25k,4,-0.75,1.5k,5,0.75,1.5k,6,0.75,1.25k,7,0.5,1k,8,0.5a,1,2,3,4,5,6,7,8 !生成横梁标准截面2 b-bsmrtsize,5amesh,allsecwrite,hl2,sect,,1SECTYPE,5,BEAM,MESH, !截面代号5SECOFFSET,CENT,,,SECREAD,'hl2','SECT',,MESH !读入横梁标准截面2asel,allaclear,alladele,all,,,1 !将图形及面、线、点全部删除（只保留截面文件）/replot !重新生成图形k,1,-0.5,0 !以下命令是用于建立一个横梁标准截面c-ck,2,-0.5,1.04k,3,-0.75,1.04+0.25k,4,-0.75,1.04+0.5k,5,-0.55,1.04+0.5k,6,-0.55,1.04+0.5+0.35k,7,0.55,1.04+0.5+0.35k,8,0.55,1.04+0.5k,9,0.75,1.04+0.5k,10,0.75,1.04+0.25k,11,0.5,1.04k,12,0.5,0k,13,-0.2,0.58k,14,-0.2,1.04+0.5-0.3k,15,0.2,1.04+0.5-0.3k,16,0.2,0.58a,1,2,3,4,5,6,7,8,9,10,11,12a,13,14,15,16asba,1,2 !横梁标准截面3 C-Csmrtsize,5amesh,allsecwrite,hl3,sect,,1SECTYPE,6,BEAM,MESH, ! 截面代号6 SECOFFSET,CENT,,,SECREAD,'hl3','SECT',,MESH !横梁标准截面3asel,allaclear,alladele,all,,,1 !将图形及面、线、点全部删除（只保留截面文件）/replot !重新生成图形ldele,allkdele,allk,1,-0.5,0 !以下命令是用于建立一个横梁标准截面d-d k,2,-0.5,1.23k,3,-0.75,1.23+0.25k,4,-0.75,1.23+0.5k,5,-0.55,1.23+0.5k,6,-0.55,1.23+0.5+0.35k,7,0.55,1.23+0.5+0.35k,8,0.55,1.23+0.5k,9,0.75,1.23+0.5k,10,0.75,1.23+0.25k,11,0.5,1.23k,12,0.5,0k,13,-0.2,0.38k,14,-0.3,0.38+0.2k,15,-0.3,1.23+0.5-0.1-0.2k,16,-0.2,1.23+0.5-0.1k,17,0.2,1.23+0.5-0.1k,18,0.3,1.23+0.5-0.1-0.2k,19,0.3,0.38+0.2k,20,0.2,0.38a,1,2,3,4,5,6,7,8,9,10,11,12a,13,14,15,16,17,18,19,20asba,1,2 !横梁标准截面4 D-Dsmrtsize,5amesh,allsecwrite,hl4,sect,,1SECTYPE,7,BEAM,MESH, !截面代号7 SECOFFSET,CENT,,,SECREAD,'hl4','SECT',,MESH !读入横梁标准截面4 asel,allaclear,alladele,all,,,1 !将图形及面、线、点全部删除（只保留截面文件）/replot ! 重新生成图形k,1,-0.5,0 !以下命令是用于建立一个横梁标准截面e-ek,2,-0.5,1.23k,3,-0.75,1.23+0.25k,4,-0.75,1.23+0.5k,5,-0.55,1.23+0.5k,6,-0.55,1.23+0.5+0.35k,7,0.55,1.23+0.5+0.35k,8,0.55,1.23+0.5k,9,0.75,1.23+0.5k,10,0.75,1.23+0.25k,11,0.5,1.23k,12,0.5,0a,1,2,3,4,5,6,7,8,9,10,11,12 !横梁标准截面5 E-E smrtsize,5amesh,allsecwrite,hl5,sect,,1SECTYPE,8,BEAM,MESH, ! 截面代号8SECOFFSET,CENT,,,SECREAD,'hl5','SECT',,MESH !横梁标准截面5asel,allaclear,alladele,all,,,1/replotetdele,1 !将辅助单元类型删除/replot!*********************************************************************!定义所有材料特性et,1,beam44 !钢管特性mp,ex,1,2.1e11mp,dens,1,7800mp,prxy,1,0.3n,90000,0,0,30 !参考点et,2,beam44 !拱肋内混凝土特性 50#混凝土mp,ex,2,3.5e10mp,dens,2,2600mp,prxy,2,0.1667et,3,beam44 ! 纵梁特性 30#混凝土mp,ex,3,3.0e10mp,dens,3,2600mp,prxy,3,0.1667et,4,beam44 !横梁特性 30#混凝土mp,ex,4,3.0e10mp,dens,4,2600mp,prxy,4,0.1667et,5,beam44 !风撑特性mp,ex,5,2.1e11mp,dens,5,7800mp,prxy,5,0.3et,6,link10 !吊杆特性（钢绞线）mp,dens,6,7800mp,prxy,6,0.3keyopt,6,3,0 !只拉吊杆et,7,beam44 !盖梁特性 30# 混凝土mp,ex,7,3.0e10mp,dens,7,2600mp,prxy,7,0.1667et,8,beam44 !墩柱特性 40#混凝土mp,ex,8,3.3e10mp,dens,8,2600mp,prxy,8,0.1667et,9,beam44 !横系梁特性 30#混凝土mp,ex,9,3.0e10mp,dens,9,2600mp,prxy,9,0.1667et,10,link10 !系杆特性（钢绞线）mp,ex,10,2.0e11mp,dens,10,7800mp,prxy,10,0.3keyopt,10,3,0 !只拉系杆et,11,beam44 !承台桩基础特性 25#混凝土mp,ex,11,2.85e10mp,dens,11,2500n,90000,0,0,30!************************************************************** ********!开始建立模型!0# 横梁(端横梁)n,1,0,2.025n,2,6.25,2.025n,3,6.25+5.1,2.025n,4,6.25+5.1+0.2+1.9/2,2.025 !12.5米位置n,5,6.25+5.1+2.3,2.025n,6,6.25+5.1+2.3+2.23,2.025nsym,x,10,2,6,1type,4 !横梁r,1(定义一组单元实常数)real,1 (划分单元时指定的实常书)mat,4secnum,7! 单元1到10为0# 横梁e,1,2 !截面7 D-Dsecnum,6e,2,3 !截面6 D-Dsecnum,5e,3,4 !截面5 D-De,4,5 !截面5 D-Dsecnum,4e,5,6 !截面4 a-asecnum,7e,1,12 !截面7 D-Dsecnum,6e,12,13 !截面6 D-Dsecnum,5e,13,14 !截面5 D-De,14,15 !截面5 D-Dsecnum,4e,15,16 !截面4 a-a!1# 端横梁egen,2,100,1,10,1,,,,,,0,4.7+1.25/2+1.5/2!其余横梁2#~17#egen,17,100,11,20,1,,,,,,0,5.3+1.5! 横梁18#(端横梁)egen,2,100,171,180,1,,,,,,0,4.7+1.5/2+1.25/2cm,hl,elemesel,none!*********************************************************** ! 纵梁单元type,3r,2real,2mat,3SECTYPE,10,beam,rect,zl !纵梁截面代号10secdata,0.4,1.2secnum,10!0# 纵梁(端纵梁）e,1,101e,2,102e,3,103e,5,105e,12,112e,13,113e,15,115!中间的纵梁1#~16#*do,i,0,15,1e,101+i*100,201+i*100e,102+i*100,202+i*100e,104+i*100,204+i*100e,112+i*100,212+i*100e,114+i*100,214+i*100*enddo!17#纵梁(端纵梁）e,1701,1801e,1702,1802e,1703,1803e,1705,1805e,1712,1812e,1713,1813e,1715,1815cm,zl,elemesel,none!************************************************************** *********!定义拱轴线 !请参照《桥梁工程》p299页相关公式*dim,y,array,16 !将纵桥向距离定义为数组*dim,cs,array,16 !**值*dim,z,array,16 !将拱轴高度定义为数组!吊杆上点!首先定义吊杆的y向（顺桥向）坐标y(1)=8.1y(2)=14.9y(3)=21.7y(4)=28.5y(5)=35.3y(6)=42.1y(7)=48.9y(8)=55.7y(9)=62.5!风撑上点!定义风撑的y向（顺桥向）坐标y(11)=21.7y(12)=28.5y(13)=42.1y(14)=48.9y(15)=59.1!0#横梁的y向（顺桥向）坐标y(16)=2.025f1=25 !矢高m1=1.1 !悬链线拱轴系数k=log(m1+sqrt(m1*m1-1))*do,i,1,16,1cs(i)=1-y(i)*2/125z(i)=25-f1/(m1-1)*(cosh(k*cs(i))-1) !拱轴高度*enddo*dim,ygd,array,23 !定义预拱度（是在后续工况中反复调算的结果）初始运行时可将预拱度全部设置为0ygd(1) = 0.40004E-02+1.29E-02ygd(2) = 0.10685E-01+2.22E-02ygd(3) = 0.19483E-01+3.03E-02ygd(4) = 0.29219E-01+3.69E-02ygd(5) = 0.39022E-01+4.23E-02ygd(6) = 0.48243E-01+4.65E-02ygd(7) = 0.56390E-01+4.95E-02ygd(8) = 0.63070E-01+5.15E-02ygd(9) = 0.68014E-01+5.21E-02ygd(10)= 0.71003E-01+5.15E-02ygd(11)= 0.71739E-01+4.95E-02ygd(12)= 0.69962E-01+4.65E-02ygd(13)= 0.65514E-01+4.23E-02ygd(14)= 0.58289E-01+3.69E-02ygd(15)= 0.48248E-01+3.04E-02ygd(16)= 0.35507E-01+2.24E-02ygd(17)= 0.20415E-01+1.23E-02ygd(18)=0.65770E-01 +0.40234E-01ygd(19)=0.69769E-01 +0.41215E-01!生成各个吊杆的Z向节点（竖直方向）ngen,2,10000,104,114,10,0,0,z(1) +ygd(1) !1#吊杆节点ngen,2,10000,204,214,10,0,0,z(2) +ygd(2) !2#吊杆节点ngen,2,10000,304,314,10,0,0,z(3) +ygd(3) !3#吊杆节点ngen,2,10000,404,414,10,0,0,z(4) +ygd(4) !4#吊杆节点ngen,2,10000,504,514,10,0,0,z(5) +ygd(5) !5#吊杆节点ngen,2,10000,604,614,10,0,0,z(6) +ygd(6) !6#吊杆节点ngen,2,10000,704,714,10,0,0,z(7) +ygd(7) !7#吊杆节点ngen,2,10000,804,814,10,0,0,z(8) +ygd(8) !8#吊杆节点ngen,2,10000,904,914,10,0,0,z(9) +ygd(9) !9#吊杆节点ngen,2,10000,1004,1014,10,0,0,z(8) +ygd(10) !10#吊杆节点ngen,2,10000,1104,1114,10,0,0,z(7) +ygd(11) !11#吊杆节点ngen,2,10000,1204,1214,10,0,0,z(6) +ygd(12) !12#吊杆节点ngen,2,10000,1304,1314,10,0,0,z(5) +ygd(13) !13#吊杆节点ngen,2,10000,1404,1414,10,0,0,z(4) +ygd(14) !14#吊杆节点ngen,2,10000,1504,1514,10,0,0,z(3) +ygd(15) !15#吊杆节点ngen,2,10000,1604,1614,10,0,0,z(2) +ygd(16) !16#吊杆节点ngen,2,10000,1704,1714,10,0,0,z(1) +ygd(17) !17#吊杆节点!说明：之所以重复定义吊杆节点是因为以后吊杆转动要进行放松ngen,2,3000,104,114,10,0,0,0 !1#吊杆下节点ngen,2,3000,204,214,10,0,0,0 !2#吊杆下节点ngen,2,3000,304,314,10,0,0,0 !3#吊杆下节点ngen,2,3000,404,414,10,0,0,0 !4#吊杆下节点ngen,2,3000,504,514,10,0,0,0 !5#吊杆下节点ngen,2,3000,604,614,10,0,0,0 !6#吊杆下节点ngen,2,3000,704,714,10,0,0,0 !7#吊杆下节点ngen,2,3000,804,814,10,0,0,0 !8#吊杆下节点ngen,2,3000,904,914,10,0,0,0 !9#吊杆下节点ngen,2,3000,1004,1014,10,0,0,0 !10#吊杆下节点ngen,2,3000,1104,1114,10,0,0,0 !11#吊杆下节点ngen,2,3000,1204,1214,10,0,0,0 !12#吊杆下节点ngen,2,3000,1304,1314,10,0,0,0 !13#吊杆下节点ngen,2,3000,1404,1414,10,0,0,0 !14#吊杆下节点ngen,2,3000,1504,1514,10,0,0,0 !15#吊杆下节点ngen,2,3000,1604,1614,10,0,0,0 !16#吊杆下节点ngen,2,3000,1704,1714,10,0,0,0 !17#吊杆下节点ngen,2,6000,104,114,10,0,0,z(1) +ygd(1) !1#吊杆节点ngen,2,6000,204,214,10,0,0,z(2) +ygd(2) !2#吊杆节点ngen,2,6000,304,314,10,0,0,z(3) +ygd(3) !3#吊杆节点ngen,2,6000,404,414,10,0,0,z(4) +ygd(4) !4#吊杆节点ngen,2,6000,504,514,10,0,0,z(5) +ygd(5) !5#吊杆节点ngen,2,6000,604,614,10,0,0,z(6) +ygd(6) !6#吊杆节点ngen,2,6000,704,714,10,0,0,z(7) +ygd(7) !7#吊杆节点ngen,2,6000,804,814,10,0,0,z(8) +ygd(8) !8#吊杆节点ngen,2,6000,904,914,10,0,0,z(9) +ygd(9) !9#吊杆节点ngen,2,6000,1004,1014,10,0,0,z(8)+ygd(10) !10#吊杆节点ngen,2,6000,1104,1114,10,0,0,z(7)+ygd(11) !11#吊杆节点ngen,2,6000,1204,1214,10,0,0,z(6)+ygd(12) !12#吊杆节点ngen,2,6000,1304,1314,10,0,0,z(5)+ygd(13) !13#吊杆节点ngen,2,6000,1404,1414,10,0,0,z(4)+ygd(14) !14#吊杆节点ngen,2,6000,1504,1514,10,0,0,z(3)+ygd(15) !15#吊杆节点ngen,2,6000,1604,1614,10,0,0,z(2)+ygd(16) !16#吊杆节点ngen,2,6000,1704,1714,10,0,0,z(1)+ygd(17) !17#吊杆节点!定义吊杆的参数（每根吊杆都采用109根直径为7毫米的高强钢丝）dgarea=109*3.1415926*0.25*0.007*0.007 !吊杆面积(m2)dgyl=0.3*1670*1e6 !吊杆初应力取用0.3倍的抗拉强度(N/m2) dgl=dgyl*dgarea !吊杆力(N)edg=1.95e11 !吊杆的弹性模量(N/m2)*dim,L0,array,17 !定义吊杆的初始0应力长度数组*dim,L,array,17 !定义吊杆的受力后长度数组*dim,deltL,array,17 !定义吊杆的受力后长度增长量数组*dim,istrn,array,17 !定义吊杆的应变数组*dim,fdg,array,17 !定义在吊杆上施加的力的数组*do,i,1,9,1L0(i)=z(i) !0应力时单元长度L0(18-i)=z(i) !吊杆是沿中轴对称的*enddofdg(1)=dgl !1#吊杆施加的预应力fdg(2)=dgl !2#吊杆施加的预应力fdg(3)=dgl !3#吊杆施加的预应力fdg(4)=dgl !4#吊杆施加的预应力fdg(5)=dgl !5#吊杆施加的预应力fdg(6)=dgl !6#吊杆施加的预应力fdg(7)=dgl !7#吊杆施加的预应力fdg(8)=dgl !8#吊杆施加的预应力fdg(9)=dgl !9#吊杆施加的预应力fdg(10)=dgl !10#吊杆施加的预应力fdg(11)=dgl !11#吊杆施加的预应力fdg(12)=dgl !12#吊杆施加的预应力fdg(13)=dgl !13#吊杆施加的预应力fdg(14)=dgl !14#吊杆施加的预应力fdg(15)=dgl !15#吊杆施加的预应力fdg(16)=dgl !16#吊杆施加的预应力fdg(17)=dgl !17#吊杆施加的预应力type,6mat,6*do,i,1,17,1deltL(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,1real,i+1000 !赋予吊杆的实常数特性e,i*100+4+3000,i*100+4+6000 !生成吊杆单元e,i*100+10+4+3000,i*100+10+4+6000*enddocm,dg,elemesel,noneallsel!************************************************************** ********* r,1018real,1018 !定义一个全为0的实常数便于利用截面特性时，不产生冲突!************************************************************** *********!生成风撑上的节点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) !(中撑中点)!风撑1（倒K撑）type,5mat,5SECTYPE,500, beam,ctube,windc !定义风撑的形状为空心圆secdata,0.486,0.50 !赋予截面形状参数,!直径按1.0米计算厚度为1.4厘米secnum,500 !调用截面!生成风撑单元!风撑1（倒K撑）e,20002,10404,90000e,20002,10414,90000e,20002,10304,90000e,20002,10314,90000!风撑2（倒K撑）e,20004,10704,90000e,20004,10714,90000e,20004,10604,90000e,20004,10614,90000!风撑3（中撑X撑）e,20000,20015,90000e,20000,20016,90000e,20000,20025,90000e,20000,20026,90000!风撑4（k撑）e,20007,11104,90000e,20007,11114,90000e,20007,11204,90000e,20007,11214,90000!风撑5（k撑）e,20009,11404,90000e,20009,11414,90000e,20009,11504,90000e,20009,11514,90000cm,fc,elemesel,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)= 30002zgzx(3)= 10104zgzx(4)= 10204zgzx(5)= 10304zgzx(6)= 10404zgzx(7)= 10504zgzx(8)= 10604zgzx(9)= 10704zgzx(10)=10804zgzx(11)=20015zgzx(12)=10904zgzx(13)=20016zgzx(14)=11004zgzx(15)=11104zgzx(16)=11204zgzx(17)=11304zgzx(18)=11404zgzx(19)=11504zgzx(20)=11604zgzx(21)=11704zgzx(22)=30003zgzx(23)=30001ygzx(1)= 40000ygzx(2)= 40002ygzx(3)= 10114ygzx(4)= 10214ygzx(5)= 10314ygzx(6)= 10414ygzx(7)= 10514ygzx(8)= 10614ygzx(9)= 10714ygzx(10)=10814ygzx(11)=20025ygzx(12)=10914ygzx(13)=20026ygzx(14)=11014ygzx(15)=11114ygzx(16)=11214ygzx(17)=11314ygzx(18)=11414ygzx(19)=11514ygzx(20)=11614ygzx(21)=11714ygzx(22)=40003ygzx(23)=40001n,90001,12.5,0,30 !钢管左轴线参考点n,90002,-12.5,0,30 !钢管左轴线参考点type,1mat,1r,2000secnum,1*do,i,1,22 !生成拱轴单元（钢管）e,zgzx(i),zgzx(i+1),90001e,ygzx(i),ygzx(i+1),90002*enddocm,gzgg,elem !生成组esel,nonetype,2r,2001secnum,2*do,i,1,22 !生成拱轴单元（混凝土）e,zgzx(i),zgzx(i+1),90001e,ygzx(i),ygzx(i+1),90002*enddocm,gzhnt,elem !生成组esel,none!************************************************************** ***************!建立拱轴底座type,9 !30#混凝土r,6real,6mat,9SECTYPE,70, beam, rect,dizuo !定义截面形状secdata,1.4,4 !赋予截面形状参数secnum,70 !调用截面e,30001,30003,90001 !拱轴底座e,40001,40003,90002e,30000,30002,90001e,40000,40002,90002cm,gzdz,elem !生成组esel,none!************************************************************** ***************!建立盖梁ngen,2,2000,1,16,1,0,-2.025,0 !复制节点到拱轴线起点位置ngen,2,100,1801,1816,1,0,2.025,0 !复制节点到拱轴线终点位置mat,7r,7real,7SECTYPE,60, beam, rect,gail !定义截面形状secdata,4,1.6 !赋予截面形状参数secnum,60 !调用截面n,90005,0,0,20 !定义参考点n,90006,0,125,20e,1901,1902,90006 !盖梁单元125米处e,1902,1903,90006e,1903,1904,90006e,1904,1905,90006e,1905,1906,90006e,1901,1912,90006e,1912,1913,90006e,1913,1914,90006e,1914,1915,90006e,1915,1916,90006e,2001,2002,90005 !盖梁单元0米处e,2002,2003,90005e,2003,2004,90005e,2004,2005,90005e,2005,2006,90005e,2001,2012,90005e,2012,2013,90005e,2013,2014,90005e,2014,2015,90005e,2015,2016,90005cm,gailiang,elem !生成组esel,none!************************************************************** *********!建立墩柱n,90003,14.5,0,0 !参考点n,90004,-14.5,0,0type,8mat,8r,8real,8SECTYPE,50, beam,hrec,dun !定义截面形状secdata,3,4,0.8,0.8,0.8,0.8 !赋予截面形状参数ngen,2,10000,2004,2014,10,0,0,-10 !墩节点距离桥面-10米ngen,2,10000,1904,1914,10,0,0,-10secnum,50e,2004,12004,90003e,2014,12014,90004e,1904,11904,90003e,1914,11914,90004cm,dunzhu,elem !生成组esel,none!************************************************************** ************!承台单元type,11mat,11r,1SECTYPE,40, beam,rect,chengtaisecdata,12.7,9r,11real, 11secnum,40ngen,2,20000,2004,2014,10,0,0,-11.5ngen,2,20000,1904,1914,10,0,0,-11.5ngen,2,30000,2004,2014,10,0,0,-13ngen,2,30000,1904,1914,10,0,0,-13e,12004,22004,90003 !生成承台单元e,11904,21904,90003e,12014,22014,90004e,11914,21914,90004e,22004,32004,90003e,21904,31904,90003e,22014,32014,90004e,21914,31914,90004cm,chtai,elem !生成组esel,none!************************************************************** ******!横系梁type,9mat,9r,9real,9n,90007,12.5,0,10 !参考点n,90008,-12.5,0,10SECTYPE,90,beam, rect,hxl !定义截面形状secdata,3,3 !赋予截面形状参数secnum,90 !调用截面e,22004,22014,90007e,21904,21914,90008cm,hxl,elem !生成组esel,none!!************************************************************** *************!桩n,50000,6.15+1.85,2.7,-13 !定义桩节点n,50001,6.15+1.85,-2.7,-13n,50002,6.15+1.85+4.5,2.7,-13n,50003,6.15+1.85+4.5,-2.7,-13n,50004,6.15+12.7-1.85,2.7,-13n,50005,6.15+12.7-1.85,-2.7,-13nsel,s,,,50000,50005,1nsym,x,10,50000,50005,1allselngen,2,100,50000,50015,1,0,125,0 !首先生成平面内的24个节点! 每3米生成一层节点(共8层)ngen,10,1000,50000,51000,1,0,0,-2.5!生成桩单元esel,nonetype,11mat,11r,11real,11SECTYPE,20, beam,csolid,piersecdata,0.9secnum,20*do,i,0,8,1*do,k,0,5,1e,50000+k+i*1000,50000+k+(i+1)*1000 ,90003e,50000+k+10+i*1000,50000+k+10+(i+1)*1000 ,90003e,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*enddocm,zhuang,elem !生成组esel,none!************************************************************** ***!在桩上施加弹簧单元!第一层距离承台-3米 317-340!计算弹簧的弹性系数 m的具体计算方法参照本书斜拉桥建模部分d=1.8kf=0.9k0=1+1/db1=0.6 !和每排桩数有关的系数,当为2排时,为0.6L1=3.6 !桩间净距h1=3*(d+1) !桩在地面或冲刷线下的计算深k=b1+(1-b1)/0.6*L1/h1B=Kf*k0*k*d !换算宽度!第1层土粉质量粘土 m=8 厚度5!第2层土粘土（细砂） m=15 厚度5 !第3层土中砂 m=20 厚度5!第4层土粘土 m=25 厚度5!第5层基岩 c0=300 厚度10米*dim,m,array,5m(1)=8m(2)=15m(3)=20m(4)=25m(5)=4000 !基岩系数*dim,h,array,5h(1)=5 !第一层土的厚度h(2)=5h(3)=5h(4)=5h(5)=10*dim,c,array,4c(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,1e(i)=0.5*(c(i-1)+c(i))*h(i)*B*1e6*enddoet,20,combin14 !弹簧单元r,20,e(1) !第1层土的弹性系数r,30,e(2) !第2层土的弹性系数r,40,e(3) !第3层土的弹性系数r,50,e(4) !第4层土的弹性系数r,60,m(5)*3.1415926*B*B/4*1e6 !基础岩石弹性系数!第1层底部距离承台-5米ngen,2,10000,51000,51115,1,1,0,0 !x 方向ngen,2,20000,51000,51115,1,0,1,0 !y 方向!ngen,2,30000,51000,51115,1,0,0,1 !z 方向type,20real,20*do,i,0,5,1e,51000+i,51000+i+10000e,51000+i+10,51000+i+10+10000e,51100+i,51100+i+10000 ! xe,51100+i+10,51100+i+10+10000e,51000+i,51000+i+20000e,51000+i+10,51000+i+10+20000e,51100+i,51100+i+20000 !ye,51100+i+10,51100+i+10+20000*enddocm,spring1,elemesel,none!第2层底部距离承台-10米ngen,2,10000,53000,53115,1,1,0,0!x 方向ngen,2,20000,53000,53115,1,0,1,0!y 方向!ngen,2,30000,53000,53115,1,0,0,1!z 方向type,20real,30*do,i,0,5,1e,53000+i+10,53000+i+10+10000e,53000+i,53000+i+10000e,53100+i,53100+i+10000 ! xe,53100+i+10,53100+i+10+10000e,53000+i,53000+i+20000e,53000+i+10,53000+i+10+20000e,53100+i,53100+i+20000e,53100+i+10,53100+i+10+20000 !y *enddocm,spring2,elemesel,none!第3层距离承台-15米ngen,2,10000,55000,55115,1,1,0,0 !x 方向ngen,2,20000,55000,55115,1,0,1,0 !y 方向type,20real,40*do,i,0,5,1e,55000+i,55000+i+10000e,55000+i+10,55000+i+10+10000e,55100+i,55100+i+10000 ! xe,55100+i+10,55100+i+10+10000e,55000+i,55000+i+20000e,55000+i+10,55000+i+10+20000e,55100+i,55100+i+20000 !ye,55100+i+10,55100+i+10+20000*enddocm,spring3,elemesel,none!!第4层距离承台-20米ngen,2,10000,57000,57115,1,1,0,0 !x 方向ngen,2,20000,57000,57115,1,0,1,0 !y 方向type,20real,50*do,i,0,5,1e,57000+i+10,57000+i+10+10000e,57100+i,57100+i+10000 ! xe,57100+i+10,57100+i+10+10000e,57000+i,57000+i+10000e,57000+i,57000+i+20000e,57000+i+10,57000+i+10+20000e,57100+i,57100+i+20000 !ye,57100+i+10,57100+i+10+20000*enddocm,spring4,elemesel,none!!第5层距离承台-22.5米ngen,2,10000,58000,58115,1,1,0,0 !x 方向ngen,2,20000,58000,58115,1,0,1,0 !y 方向type,20real,60*do,i,0,5,1e,58000+i+10,58000+i+10+10000e,58000+i,58000+i+10000e,58100+i,58100+i+10000 ! xe,58100+i+10,58100+i+10+10000e,58000+i,58000+i+20000e,58000+i+10,58000+i+10+20000e,58100+i,58100+i+20000 !ye,58100+i+10,58100+i+10+20000*enddocm,spring5,elemesel,none!!第6层距离承台-25米ngen,2,10000,59000,59115,1,1,0,0 ngen,2,10000,59000,59115,1,1,0,0 !x 方向ngen,2,20000,59000,59115,1,0,1,0 ngen,2,20000,59000,59115,1,0,1,0 !y 方向ngen,2,30000,59000,59115,1,0,0,-1 ngen,2,30000,59000,59115,1,0,0,-1 !z 方向type,20real,60*do,i,0,5,1e,59000+i,59000+i+10000e,59000+i+10,59000+i+10+10000e,59100+i,59100+i+10000 ! xe,59100+i+10,59100+i+10+10000e,59000+i,59000+i+20000e,59000+i+10,59000+i+10+20000e,59100+i,59100+i+20000 !ye,59100+i+10,59100+i+10+20000e,59000+i,59000+i+30000e,59000+i+10,59000+i+10+30000e,59100+i,59100+i+30000 !ze,59100+i+10,59100+i+10+30000*enddocm,spring6,elemallsel!************************************************************** *********!剔除不需要进行对称复制的单元allselnsel,s,loc,z,-50,-0.01nsel,r,loc,y,-10,10cm,zdnode,nodeesln,s,,allcm,zdelem,elemallselnsel,u,,,zdnodensym,y,100000,all !节点对称allselESEL,U,,,zdelemESYM,,100000,ALL !见图11-5ALLSEL!图11-5 全桥单元模型图!************************************************************** ********* !定义系杆单元（全桥通长）!系杆单元!定义系杆的参数（每根系杆都采用8束27*7*d5高强钢丝）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,2xgyl=0.3*1860*1e6xgl=xgyl*xgareafxg(1)=xgl !施加的预应力fxg(2)=xgl !施加的预应力Lzero=250 !0应力时单元长度*do,i,1,2,1detL(i)=Lzero*Fxg(i)/(exg*xgarea)Llast(i)=Lzero+detL(i) !单元最终长度Listrn(i)=detl(i)/Llast(i)r,i+100,xgarea,Listrn(i)*enddotype,10mat,10real,101e,30001,130001 !系杆1 250米real,102e,40001,140001 !系杆2 250米.cm,xigan,elemesel,noneallsel!************************************************************** *******!将盖梁和拱轴交界点进行节点耦合,放松uy,和rotxcp,1,ux,2004,30000cp,1,uy,2004,30000cp,next,uz,2004,30000cp,next,roty,2004,30000 cp,next,rotz,2004,30000 cp,next,ux,1904,30001cp,next,uy,1904,30001cp,next,uz,1904,30001cp,next,roty,1904,30001 cp,next,rotz,1904,30001 cp,next,ux,2014,40000cp,next,uy,2014,40000cp,next,uz,2014,40000cp,next,roty,2014,40000 cp,next,rotz,2014,40000 cp,next,ux,1914,40001cp,next,uy,1914,40001cp,next,uz,1914,40001cp,next,roty,1914,40001 cp,next,rotz,1914,40001cp,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 allselnsel,s,,,30000,130000,100000nsel,a,,,40000,140000,100000nummrg,nodeallsel!************************************************************* ********do,i,1,3,1cp,next,all,1800+i,1900+i !盖梁和0#横梁耦合*enddoallsel*do,i,5,6,1cp,next,all,1800+i,1900+i*enddoallsel*do,i,1,2,1cp,next,all,1800+i+11,1900+i+11cp,next,all,i+11,2000+i+11,i+100000+11*enddoallsel*do,i,4,5,1cp,next,all,1800+i+11,1900+i+11cp,next,all,i+11,2000+i+11,i+100000+11*enddoallsel*do,i,1,3,1cp,next,all,i,2000+i,100000+i*enddoallsel*do,i,5,6,1cp,next,all,i,2000+i,100000+i*enddoallsel*do,i,1,3,1cp,next,all,101800+i,101900+i !盖梁和0#横梁耦合*enddoallsel*do,i,5,6,1cp,next,all,101800+i,101900+iallsel*enddo*do,i,1,2cp,next,all,101800+i+11,101900+i+11*enddoallsel*do,i,4,5,1cp,next,all,101800+i+11,101900+i+11*enddoallsel!************************************************************** ********!!桩节点和承台节点耦合nsel,s ,loc,z,-13nsel,r ,loc,y,-3,3nsel,r ,loc,x,8,20cp,next,all,allallselnsel,s ,loc,z,-13nsel,r ,loc,y,-3,3nsel,r ,loc,x,-8,-20cp,next,all,allallselnsel,s ,loc,z,-13nsel,r ,loc,y,122,128nsel,r ,loc,x,-8,-20cp,next,all,allallselnsel,s ,loc,z,-13nsel,r ,loc,y,122,128nsel,r ,loc,x,8,20cp,next,all,allallselnsel,s ,loc,z,-13nsel,r ,loc,y,-122,-128nsel,r ,loc,x,8,20cp,next,all,allallselnsel,s ,loc,z,-13nsel,r ,loc,y,-122,-128nsel,r ,loc,x,-8,-20cp,next,all,allallsel!************************************************************** ********!将吊杆转动进行放松*do,i,1,17,1cp,next,ux,i*100+4,i*100+4+3000cp,next,ux,i*100+4+6000,i*100+4+10000cp,next,uy,i*100+4,i*100+4+3000cp,next,uy,i*100+4+6000,i*100+4+10000cp,next,uz,i*100+4,i*100+4+3000cp,next,uz,i*100+4+6000,i*100+4+10000cp,next,ux,i*100+4+10,i*100+4+3000+10cp,next,ux,i*100+4+6000+10,i*100+4+10000+10cp,next,uy,i*100+4+10,i*100+4+3000+10cp,next,uy,i*100+4+6000+10,i*100+4+10000+10cp,next,uz,i*100+4+10,i*100+4+3000+10cp,next,uz,i*100+4+6000+10,i*100+4+10000+10cp,next,ux,i*100+4+100000,i*100+4+3000+100000cp,next,ux,i*100+4+6000+100000,i*100+4+10000+100000 cp,next,uy,i*100+4+100000,i*100+4+3000+100000cp,next,uy,i*100+4+6000+100000,i*100+4+10000+100000 cp,next,uz,i*100+4+100000,i*100+4+3000+100000cp,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+10,i*100+4+10000+100 000+10cp,next,uy,i*100+4+100000+10,i*100+4+3000+100000+10cp,next,uy,i*100+4+6000+100000+10,i*100+4+10000+100 000+10cp,next,uz,i*100+4+100000+10,i*100+4+3000+100000+10 cp,next,uz,i*100+4+6000+100000+10,i*100+4+10000+100 000+10*enddoallsel!************************************************************** ********!将土弹簧的一端约束allsel!nsel,u,,,erasensel,s,loc,x,13.5nsel,a,loc,x,9nsel,a,loc,x,18nsel,a,loc,x,-11.5nsel,a,loc,x,-7nsel,a,loc,x,-16nsel,a,loc,y,3.7nsel,a,loc,y,-1.7nsel,a,loc,y,128.7nsel,a,loc,y,123.3nsel,a,loc,y,-128.7nsel,a,loc,y,-123.3nsel,r,loc,z,-40,-13d,all,allallselnsel,s,loc,z,-36.5d,all,all/color,pbak,off/replotallselnsel,s,loc,y,0nsel,r,loc,z,0nummrg,node,allesel,noneallsel!************************************************************** ********!工况1 不包括桥面板时，求横梁和纵梁的总重量/soluallselmpdele,dens,allmp,dens,3,2600mp,dens,4,2600acel,,,9.8pstres,onsolvesave! 总重量：0.65064E+08 N!************************************************************** *********!工况2 恒载下调节预拱度首先只有拱肋自重作用下的拱的变位/soluallselmpdele,dens,3,11,1mp,dens,1,7800mp,dens,2,2600*do,i,1,17,1r,1000+i,dgarea,0*enddor,101,xgarea,0r,102,xgarea,0acel,,,9.8solve!************************************************************** *********!工况3一期恒载（包括桥面板）作用下的无预应力状态/soluallsel*do,i,1,17,1。
/TITLE, BRIDGE XX/PREP7! 定义单元类型ET,1,SHELL63ET,2,SOLID45ET,3,LINK8!定义1,钢，2，混凝土，3，无限刚度MP,EX,1,MP,DENS,1,MP,PRXY,1,mp,alpx,1,MP,EX,2,MP,DENS,2,MP,PRXY,2,MP,EX,3,MP,DENS,3,MP,PRXY,3,*set,nz,154!定义箱梁N1R,1,7,7,7,7!定义箱梁N2R,2,8,8,8,8!定义箱梁N3R,3,4,4,4,4!定义混凝土solid65的实常数R, 4!定义钢绞线的实常数R, 5, 140!降温=预加力/线膨胀系数*钢筋截面积*弹性模量tem1=*tem2=*!定义纵向梯度参数表*dim, z, array, nz*SET,z(1),130,260,371,482,593,702,811,920,1029,1138,1247,1356,146 5,1574,1683,,,*SET,z(19),,,,,,,,,,,,,*SET,z(31),,,,,,*SET,z(37),,,,,*SET,z(41),,,,,,,,,,,,,,*SET,z(55),,,,,,,,,,,*SET,z(65),,,,,,,,*SET,z(73),,,,,,,,*SET,z(80),,,,,,,,,,,*SET,z(91),,,,,,,,,,*SET,z(100),,,,,,,,,*SET,z(109),,,,,,,*SET,z(115),,,,,,,,,,,,*SET,z(127),,,,,,,,,*SET,z(135),,,,,,,,,,*SET,z(145),,,,,,,,,,!!定义横隔梁的梯度参数表*set,nd,18*dim, del, array, nd*SET,del(1), 2,5,15,31,47,63,76,82,92,95,98,108,114,125,132,142, *SET,del(17),150,153!!定义初始面座标参数表*set,ni,36*dim, initial, array,ni, 3*SET,initial(1,1),-633,-533,,,-133,0,,,400,,633,-633,-533,,,-133,*SET,initial(17,1),0,,,400,,633,-633,-633,633,633,-766,,-633,,-50 0,500,,633,*SET,initial(35,1),,766*SET,initial(1,2),0,0,0,0,0,0,0,0,0,0,0,0*SET,initial(12,2),60,60,60,60,60,60,60,60,60,60,60,*SET,Initial(23,2),184,276,184,276,368,368,368,368,368,368,368,36 8,368,368,*SET,initial(1,3), , , , , , , , , , , ,*SET,initial(12,3), , , , , , , , , , , , ,*SET,initial(24,3), , , , , , , , , , , , ,!定义横隔梁参数表*set,nv,25*dim,vertical,array,nv,3*SET,vertical(1 ,3), , , , , , , , , , , , ,*SET,vertical(14,3), , , , , , , , , , , ,*SET,vertical(1,2),184,184,184,184,184,184,184,184,184,*SET,vertical(10,2),276,276,276,276,276,276,276,276,276,*SET,vertical(19,2),368,368,368,368,368,368,368*SET,vertical(1,1),-533, , ,-133 ,0 , , , , , , , ,-133*SET,vertical(14,1),0, , , , , , , , , , ,!!定义混凝土的节点参数表*set,nc,135*dim,concrete,array,nc,3*SET,concrete(1,1), , , , , , , , , ,*SET,concrete(11,1), , , , , , , , , ,*SET,concrete(21,1), , , , , , , , , ,*SET,concrete(31,1), , , , , , , , , ,*SET,concrete(41,1), , , , , , , , , ,*SET,concrete(51,1), , , , , , , , , ,*SET,concrete(61,1), , , , , , , , , ,*SET,concrete(71,1), , , , , , , , , ,*SET,concrete(81,1), , , , , , , , , ,*SET,concrete(91,1), , , , , , , , , ,*SET,concrete(101,1), , , , , , , , , ,*SET,concrete(111,1), , , , , , , , , ,*SET,concrete(121,1), , , , , , , , , ,*SET,concrete(131,1), , , , ,*SET,concrete(1, 2), , , , , , , , , ,*SET,concrete(11, 2), , , , , , , , , ,*SET,concrete(21, 2), , , , , , , , , ,*SET,concrete(31, 2), , , , , , , , , ,*SET,concrete(41, 2), , , , , , , , , ,*SET,concrete(51, 2), , , , , , , , , ,*SET,concrete(61, 2), , , , , , , , , ,*SET,concrete(71, 2), , , , , , , , , ,*SET,concrete(81, 2), , , , , , , , , ,*SET,concrete(91, 2), , , ,,, ,,,368,*SET,concrete(101, 2), ,368,,,368, ,,,,*SET,concrete(111, 2), ,,,,, ,,,,*SET,concrete(121, 2), ,,,,,,,,,*SET,concrete(131, 2), ,,,,,*SET,concrete(1, 3), , , , , , , , , ,*SET,concrete(11, 3), , , , , , , , , ,*SET,concrete(21, 3), , , , , , , , , ,*SET,concrete(31, 3), , , , , , , , , ,*SET,concrete(41, 3), , , , , , , , , ,*SET,concrete(51, 3), , , , , , , , , ,*SET,concrete(61, 3), , , , , , , , , ,*SET,concrete(71, 3), , , , , , , , , ,*SET,concrete(81, 3), , , , , , , , , ,*SET,concrete(101, 3), ,,,,, ,,,,*SET,concrete(111, 3), ,,,,, ,,,,*SET,concrete(121, 3), ,,,,,,,,,*SET,concrete(131, 3), ,,,,,!节点编号，生成箱梁的节点*do, i, 1,1,1*do, j, 1, nin, j+ni*i, initial(j,1),initial(j,2),initial(j,3)+z(i) *enddo*enddo*do, j, 1, nin, j+ni*i, initial(j,1),initial(j,2),initial(j,3)+z(i)*enddo*enddo*do, i, 2,nz-1,1*do, j, 1, nin, j+ni*i, initial(j,1),initial(j,2),initial(j,3)+z(i)*enddo*enddo!生成横隔板节点*set, infl, ni*(nz+1)*set, nfl, 39*do, i, 1, 18, 1*SET,k,del(i)*do, j, 1, 25, 1n, infl+i*nfl+j, vertical(j,1), vertical(j,2), vertical(j,3)+z(k)*enddo*do, j, 26, 32n, infl+i*nfl+j, vertical(j-7,1), vertical(j-7,2), vertical(j-7,3)+z(k-1)*enddon, infl+i*nfl+j, vertical(j-14,1), vertical(j-14,2), vertical(j-14,3)+z(k+1)*enddo*enddo!节点编号，生成混凝土的节点*set,incr,infl+18*nfl+1*do, i, 1, 1,*do, j, 1, ncn, incr+j+nc*i, concrete(j,1),concrete(j,2),concrete(j,3)+z(i)*enddo*enddo*do, i, 2, nz*do, j, 1, ncn, incr+j+nc*i, concrete(j,1),concrete(j,2),concrete(j,3)+z(i)*enddo*enddo!生成钢绞线的节点*set,insty1,incr+nc*(nz+1)!Y1*do, i, 37, 127, 1n, insty1+(i-36),concrete(53,1),concrete(53,2),concrete(53,3)+z(i)n, insty1+(i-36)+127-37+1, concrete(68,1),concrete(68,2),concrete(68,3)+z(i)n, insty1+(i-36)+2*(127-37)+2, concrete(83,1),concrete(83,2),concrete(83,3)+z(i)*enddo!Y2*set,insty2,insty1+3*(127-37+1)*do, i, 66,114,1n, insty2+(i-65), concrete(56,1),concrete(56,2),concrete(56,3)+z(i)n, insty2+(i-65)+141-66+1, concrete(80,1),concrete(80,2),concrete(80,3)+z(i)*enddo!生成钢绞线锚固区节点*set,instb1,insty2+2*(141-66)*set,instb2,instb1+40!Y1*do, i, 1, 4*SET,k,27+in, instb1+i, , +12*i, +z(k)n, instb2+i, , +12*i, +z(k)*enddo*do, i, 1, 4*SET,k,27+in, instb1+i+4, 0 , +12*i, z(k) n, instb2+i+4, 0 , +12*i, z(k) *enddo*do, i, 1, 4*SET,k,27+in, instb1+i+8, , +12*i, +z(k) n, instb2+i+8, , +12*i, +z(k) *enddo*do, i, 1, 4*SET,k,128+i+6n, instb1+i+12, , *(i+6), +z(k) n, instb2+i+12, , *(i+6), +z(k) *enddo*do, i, 1, 4*SET,k,128+i+6n, instb1+i+16, 0 , *(i+6), z(k) n, instb2+i+16, 0 , *(i+6), z(k) *enddo*do, i, 1, 4*SET,k,128+i+6n, instb1+i+20, , *(i+6), +z(k) n, instb2+i+20, , *(i+6), +z(k) *enddo!Y2*do, i, 1, 4*SET,k,56+in, instb1+i+24, , +12*i, +z(k) n, instb2+i+24, , +12*i, +z(k) *enddo*do, i, 1, 4*SET,k,56+in, instb1+i+28, ,+12*i, +z(k) n, instb2+i+28, ,+12*i, +z(k) *enddo*do, i, 1, 4*SET,k,115+i+6n, instb1+i+32, ,*(i+6), +z(k) n, instb2+i+32, ,*(i+6), +z(k) *enddo*do, i, 1, 4*SET,k,115+i+6n, instb1+i+ni, ,*(i+6), +z(k)n, instb2+i+ni, ,*(i+6), +z(k)*enddo!生成箱梁的单元*do, i, 1,nz-1,1!底板，shell63，钢，厚度8type,1mat,1real,2*do, j, 1, 10, 1e, i*ni+j, i*ni+j+1, (i+1)*ni+j+1, (i+1)*ni+j, *enddo!顶板，shell63，钢，厚度7type,1mat,1real,1*do, j, 27, 30, 1e, i*ni+j, i*ni+j+1, (i+1)*ni+j+1, (i+1)*ni+j, *enddo*do, j, 32, 35, 1e, i*ni+j, i*ni+j+1, (i+1)*ni+j+1, (i+1)*ni+j, *enddo!腹板、底板加劲肋，shell63，钢，厚度4 type,1mat,1real,3*do, j, 1, 11, 1e, i*ni+j, (i+1)*ni+j, (i+1)*ni+j+11, i*ni+j+11*enddoe, i*ni+12, (i+1)*ni+12, (i+1)*ni+23, i*ni+23e, i*ni+23, (i+1)*ni+23, (i+1)*ni+24, i*ni+24e, i*ni+24, (i+1)*ni+24, (i+1)*ni+29, i*ni+29e, i*ni+22, (i+1)*ni+22, (i+1)*ni+25, i*ni+25e, i*ni+25, (i+1)*ni+25, (i+1)*ni+26, i*ni+26e, i*ni+26, (i+1)*ni+26, (i+1)*ni+34, i*ni+34*enddo!生成横隔板单元type,1mat,1real,3*do, i, 1, 18, 1!生成横隔板单元*do, l, 0, 9, 1e, del(i)*ni+1+l, del(i)*ni+2+l, del(i)*ni+13+l, del(i)*ni+12+l,*enddo*do, l, 0, 7,1e, infl+i*nfl+l+1, infl+i*nfl+l+2, infl+i*nfl+l+11, infl+i*nfl+l+10,e, del(i)*ni+l+13,del(i)*ni+l+14, infl+i*nfl+l+2, infl+i*nfl+l+1,*enddo*do, l, 0, 5, 1e, infl+i*nfl+l+11, infl+i*nfl+l+12, infl+i*nfl+l+20, infl+i*nfl+l+19,*enddo*do, l, 0, 5, 1e, infl+i*nfl+l+19, infl+i*nfl+l+20, infl+i*nfl+l+27, infl+i*nfl+l+26e, infl+i*nfl+l+33, infl+i*nfl+l+34, infl+i*nfl+l+20, infl+i*nfl+l+19*enddoe, del(i)*ni+31, infl+i*nfl+19, infl+i*nfl+26, (del(i)-1)*ni+31,e, (del(i)+1)*ni+31, infl+i*nfl+33, infl+i*nfl+19, del(i)*ni+31,e, infl+i*nfl+25, del(i)*ni+32, (del(i)-1)*ni+32, infl+i*nfl+32,e, infl+i*nfl+nfl, (del(i)+1)*ni+32, del(i)*ni+32, infl+i*nfl+25,e, del(i)*ni+12, del(i)*ni+13, infl+i*nfl+1,del(i)*ni+23e, del(i)*ni+23, infl+i*nfl+1, infl+i*nfl+10, del(i)*ni+24,e, del(i)*ni+24, infl+i*nfl+10, del(i)*ni+31, del(i)*ni+29,e, infl+i*nfl+10, infl+i*nfl+11, infl+i*nfl+19, del(i)*ni+31,e, del(i)*ni+21, del(i)*ni+22, del(i)*ni+25, infl+i*nfl+9,e, infl+i*nfl+9, del(i)*ni+25, del(i)*ni+26, infl+i*nfl+18,e, infl+i*nfl+18, del(i)*ni+26, del(i)*ni+34, del(i)*ni+32,e, infl+i*nfl+17, infl+i*nfl+18, del(i)*ni+32, infl+i*nfl+25,*enddo!生成混凝土单元type,2mat,2real, 4*do, i, 1, 153, 1*do, j, 0, 130, 3*do, k, 2, 3, 1e,incr+nc*i+j+k,incr+nc*(i+1)+j+k,incr+nc*(i+1)+j+k+3,incr+nc*i+j +k+3,incr+nc*i+j+k-1,incr+nc*(i+1)+j+k-1,incr+nc*(i+1)+j+k+2, incr+nc*i+j+k+2,*enddo*enddo*enddo!生成Y1钢绞线单元type,3real, 5*do, i, 37, 126, 1e, insty1+i-36,insty1+i-ni+1e, insty1+i-36+127-37+1, insty1+i-ni+127-37+2e, insty1+(i-36)+2*(127-37)+2,insty1+(i-ni)+2*(127-37)+3 *enddoe, instb1+4, insty1+1,e, instb1+8, insty1+1+127-37+1e, instb1+12, insty1+1+2*(127-37)+2e, instb1+13, insty1+(127-36)e, instb1+17, insty1+(127-36)+127-37+1,e, instb1+21, insty1+(127-36)+2*(127-37)+2*do, i, 1, 3*do, j ,0, 5e, instb1+i+j*4, instb1+i+j*4+1*enddo*enddofinish/SOL!进入求解器施加荷载ESEL,S,TYPE,,3 !选择现有单元类型为LINK8的单元BFE,all,TEMP,1,tem1 !对目前显示的所有单元施加温度荷载，大小为-tem1 ESEL,none !全不选择/PREP7!进入前处理器,生成Y2钢绞线的单元*set,insty2,insty1+3*(127-37+1)*do, i, 66,114,1n, insty2+(i-65), concrete(56,1),concrete(56,2),concrete(56,3)+z(i)n, insty2+(i-65)+141-66+1, concrete(80,1),concrete(80,2),concrete(80,3)+z(i)*enddo*do, i, 66, 113, 1e, insty2+i-65,insty2+i-65+1e, insty2+(i-65)+141-66+1, insty2+(i-65)+141-66+2*enddoe, instb1+28, insty2+1e, instb1+32, insty2+1+141-66+1e, instb1+33, insty2+(113-65)e, instb1+37,insty2+(113-65)+141-66+1*do, i, 1, 3*do, j ,0, 3e, instb1+i+j*4+24, instb1+i+j*4+1+24*enddo*enddoFINISH/SOL!进入求解器施加荷载BFE,all,TEMP,1,tem2 !对目前选择的单元即Y2施加温度荷载，大小为-tem2 ESEL,all !选择所有的单元FINISH/PREP7!生成钢筋锚固区type,2mat,2real, 4*do, i, 1, 3e, instb2+i, instb2+4+i, instb2+5+i, instb2+1+i, incr+54+nc*(i+27),incr+69+nc*(i+27),incr+69+nc*(i+28),incr+54+nc*(i+28),e, instb2+4+i,instb2+8+i, instb2+9+i,instb2+5+i, incr+nc*(i+27)+69,incr+nc*(i+27)+84, incr+nc*(i+28)+84, incr+nc*(i+28)+69e, instb2+12+i, instb2+16+i, instb2+17+i, instb2+13+i,incr+54+nc*(i+134),incr+69+nc*(i+134),incr+69+nc*(i+135),incr+54+nc*(i+135)e, instb2+16+i, instb2+20+i, instb2+21+i, instb2+17+i, incr+69+nc*(i+134),incr+84+nc*(i+134),incr+84+nc*(i+135),incr+69+nc*(i+135)*enddo*do, i,1, 3e, instb2+24+i, instb2+28+i, instb2+29+i, instb2+25+i,incr+57+nc*(i+56),incr+81+nc*(i+56),incr+81+nc*(i+57),incr+57+nc*(i+57)e, instb2+32+i, instb2+36+i, instb2+37+i, instb2+33+i,incr+57+nc*(i+121),incr+81+nc*(i+121),incr+81+nc*(i+122),incr+57+nc*(i+122)*enddo!加载FINISH/SOL*do, i, 1, 2*do, j, 4,8d, ni*i+j, all,*enddo*enddo*do, i, 94, 96 *do, j, 4,8d, ni*i+j, all,*enddo*enddo*do, i, 153, 154 *do, j, 4,8d, ni*i+j, ux, *enddo*enddo*do, i, 153, 154 *do, j, 4,8d, ni*i+j, uz, *enddo*enddo FINISH/PREP7 cpintf, ux, cpintf, uy, cpintf, uz, nsel, alleplot/solu solve。
大桥全长2996.8m，其中主桥采用跨度为101.5+188.5+580+217.5+159.5+116m的钢桁梁斜拉桥；非通航孔正桥采用6孔跨径64m预应力混凝土简支箱梁；东引桥采用16孔梁长32.6m预应力混凝土简支箱梁；跨大堤桥采用48.9+86+48.8m预应力混凝土连续箱梁；西引桥采用15孔梁长32.6m预应力混凝土简支梁及2孔梁长24.6m预应力混凝土简支梁，其中宁安线采用箱梁，阜景线采用T梁。
主桥采用103+188.5+580+217.5+159.5+117.5m两塔钢桁斜拉桥方案，全长1366m。
主梁为三片主桁钢桁梁，桁间距2x14m，节间长14.5m，桁高15m。
主塔为钢筋混凝土结构，塔顶高程+204.00m，塔底高程-6.00m，斜拉索为空间三索面，立面上每塔两侧共18对索，全桥216根斜拉索。
所有桥墩上均设竖向和横向约束，4#塔与主梁之间设纵向水平约束，3#塔与梁间使用带限位功能的粘滞阻尼器。
主梁为”N”字型桁式，横向采用三片桁结构，主桁的横向中心距各为14m，桁高15m，节间距14.5m[2]。
结构构造主桥采用两塔钢桁斜拉桥方案，主梁为三片主桁钢桁梁，主桁上下弦杆均为箱型截面，上弦杆内高1000mm，内宽1200mm，板厚20～48mm。
下弦杆内高1400mm，宽1200mm，板厚20～56mm。
下弦杆顶板向桁内侧加宽700mm与整体桥面板焊接。
腹杆主要采用H型截面。
H型杆件宽1200mm，高720和760mm，板厚20～48mm。
根据不同的受力区段选用不同的杆件截面，在辅助墩附近的压重区梁段，腹杆采用箱型截面杆件。
主桁采用焊接杆件，整体节点。
在节点外以高强度螺栓拼接的结构形式，上下弦杆四面等强对接拼装。
H型腹杆采用插入式连接。
箱型腹杆采用四面与主桁节点对拼的连接形式。
主桁拼接采用M30高强螺栓。
表1.2 安庆长江大桥主要工程数量表表 2.1 设计构件节点和单元数表 2.2安庆长江铁路大桥模型主要参数列表表2.3 模型实常数/nerr,0!将参数读入数据库finish/clear/PREP7!单元类型编号et,1,beam44et,2,shell63et,3,link8!材料类型编号Q420/Q370/拉索-1，混凝土C40/C50-2，刚臂-3 mp,ex,1,2.1e11mp,prxy,1,0.3mp,dens,1,7850mp,alpx,1,1e-5mp,ex,2,3.5e7mp,prxy,2,0.16mp,dens,2,2500mp,alpx,2,1e-5mp,ex,3,1e15mp,prxy,3,1mp,dens,3,0mp,alpx,3,0!读取自定义截面并编号*dim,name,char,24name(1)='Sec_1a'name(2)='Sec_2a'name(3)='Sec_3a'name(4)='Sec_4a'name(5)='Sec_5a'name(6)='Sec_6a'name(7)='Sec_7a'name(8)='Sec_1'name(9)='Sec_2'name(10)='Sec_3'name(11)='Sec_4'name(12)='Sec_5'name(13)='Sec_6'name(14)='Sec_7'name(15)='Sec_8'name(16)='Sec_9'name(17)='Sec_10'name(18)='Sec_11'name(19)='Sec_12'name(20)='Sec_13'name(21)='Sec_14'name(22)='Sec_15'name(23)='Sec_16'name(24)='Sec_x2'*do,i,1,24,1sectype,i,beam,mesh,secread,name(i),,,mesh*enddo!将ta_sec参数读入数据库!读取自定义截面并编号*dim,name2,char,15name2(1)='ta_1'name2(2)='ta_2'name2(3)='ta_3'name2(4)='ta_4'name2(5)='ta_5'name2(6)='ta_6'name2(7)='ta_7'name2(8)='ta_8'name2(9)='ta_9'name2(10)='ta_10'name2(11)='ta_11'name2(12)='ta_12'name2(13)='ta_13'name2(14)='ta_14'*do,i,1,14,1sectype,40+i,beam,mesh,secread,name2(i),,,mesh*enddoSECTYPE, 26, BEAM, HREC, , 0 SECOFFSET, USER, 0.1, 0.266SECDATA,0.2,0.266,0.006,0.006,0.006,0.006,0,0,0,0 SECTYPE, 27, BEAM, I, , 0SECOFFSET, USER, 0, 0.584SECDATA,0.240,0.240,0.6,0.016,0.016,0.012,0,0,0,0 SECTYPE, 28, BEAM, T, T_nin, 0 SECOFFSET, USER, 0, 1.4SECDATA,0.6,1.4,0.032,0.018,0,0,0,0,0,0SECTYPE, 29, BEAM, T, T_max, 0SECOFFSET, USER, 0, 2.1SECDATA,0.6,2.1,0.032,0.016,0,0,0,0,0,0 SECTYPE, 31, BEAM, HREC, , 0SECOFFSET, CENTSECDATA,0.548,0.500,0.016,0.016,0.016,0.016,0,0,0,0 SECTYPE, 32, BEAM, I, , 0SECOFFSET, CENTSECDATA,0.560,0.560,0.500,0.028,0.028,0.016,0,0,0,0 !刚臂的截面SECTYPE, 100, BEAM, CSOLID, , 0 SECOFFSET, CENTSECDATA,0.06,0,0,0,0,0,0,0,0,0r,1,r,2,0.016,!拉索截面控制!7-349r,3,134.31e-4,!7-301r,4,115.84e-4,!7-283r,5,108.91e-4,!7-241r,6,100.000e-4,finish!开始建立模型/prep7!********************************************* !!!xiaXianGan!!!********************************************* *SET,j,-14*do,i,1,(94*25+1)*2,94*25+1n,i,0,j,0n,94*25+i,1363.000,j,0FILL,i,94*25+i,94*25-1,,1,,,!中桁*if,i,eq,94*25+2,then!7截面type,1mat,1secnum,14*do,a,1,3*25,1E,i-1+a,i+a*enddo*do,b,1,11*25,1E,35*25-1+i+b,35*25+i+b *enddo*do,c,1,8*25,1E,77*25-1+i+c,77*25+i+c *enddo!6截面type,1mat,1secnum,13*do,a,1,3*25,1E,i+3*25-1+a,i+3*25+a*enddo*do,b,1,3*25,1E,32*25-1+i+b,32*25+i+b *enddo*do,c,1,3*25,1E,46*25-1+i+c,46*25+i+c *enddo*do,d,1,1*25,1E,i+76*25-1+d,i+76*25+d *enddo*do,e,1,7*25,1E,87*25-1+i+e,87*25+i+e *enddo!5截面type,1mat,1secnum,12*do,a,1,25,1E,8*25-1+i+a,8*25+i+a*enddo*do,b,1,2*25,1E,30*25-1+i+b,30*25+i+b *enddo*do,c,1,2*25,1*enddo*do,d,1,2*25,1E,72*25-1+i+d,72*25+i+d *enddo!4截面type,1mat,1secnum,11*do,a,1,3*25,1E,9*25-1+i+a,9*25+i+a*enddo*do,b,1,4*25,1E,26*25-1+i+b,26*25+i+b *enddo*do,c,1,3*25,1E,51*25-1+i+c,51*25+i+c *enddo*do,d,1,4*25,1E,68*25-1+i+d,68*25+i+d *enddo!3截面type,1mat,1secnum,10*do,a,1,1*25,1E,7*25-1+i+a,7*25+i+a*enddo*do,b,1,5*25,1E,12*25-1+i+b,12*25+i+b *enddo*do,c,1,4*25,1E,22*25-1+i+c,22*25+i+c *enddo*do,d,1,4*25,1E,54*25-1+i+d,54*25+i+d *enddo*do,e,1,5*25,1E,63*25-1+i+e,63*25+i+e *enddo!2截面type,1mat,1secnum,9*enddo*do,b,1,2*25,1E,17*25-1+i+b,17*25+i+b*enddo*do,c,1,1*25,1E,21*25-1+i+c,21*25+i+c*enddo*do,d,1,1*25,1E,58*25-1+i+d,58*25+i+d*enddo*do,e,1,2*25,1E,61*25-1+i+e,61*25+i+e*enddo*do,f,1,2*25,1E,74*25-1+i+f,74*25+i+f*enddo*do,g,1,1*25,1E,85*25-1+i+g,85*25+i+g*enddo!1截面type,1mat,1secnum,8*do,a,1,2*25,1E,19*25-1+i+a,19*25+i+a*enddo*do,b,1,2*25,1E,59*25-1+i+b,59*25+i+b*enddo*do,c,1,1*25,1E,86*25-1+i+c,86*25+i+c*enddo!左边桁*elseif,i,eq,(94*25+1)*2+1,then !7截面type,1mat,1secnum,7*do,a,1,3*25,1E,i+a,i-1+a*enddo*do,b,1,11*25,1*do,c,1,8*25,1E,77*25+i+c,77*25-1+i+c *enddo!6截面type,1mat,1secnum,6*do,a,1,3*25,1E,i+3*25+a,i+3*25-1+a*enddo*do,b,1,3*25,1E,32*25+i+b,32*25-1+i+b *enddo*do,c,1,3*25,1E,46*25+i+c,46*25-1+i+c *enddo*do,d,1,1*25,1E,i+76*25+d,i+76*25-1+d *enddo*do,e,1,7*25,1E,87*25+i+e,87*25-1+i+e *enddo!5截面type,1mat,1secnum,5*do,a,1,25,1E,8*25+i+a,8*25-1+i+a*enddo*do,b,1,2*25,1E,30*25+i+b,30*25-1+i+b *enddo*do,c,1,2*25,1E,49*25+i+c,49*25-1+i+c *enddo*do,d,1,2*25,1E,72*25+i+d,72*25-1+i+d *enddo!4截面type,1mat,1secnum,4*do,a,1,3*25,1*do,b,1,4*25,1E,26*25+i+b,26*25-1+i+b *enddo*do,c,1,3*25,1E,51*25+i+c,51*25-1+i+c *enddo*do,d,1,4*25,1E,68*25+i+d,68*25-1+i+d *enddo!3截面type,1mat,1secnum,3*do,a,1,1*25,1E,7*25+i+a,7*25-1+i+a*enddo*do,b,1,5*25,1E,12*25+i+b,12*25-1+i+b *enddo*do,c,1,4*25,1E,22*25+i+c,22*25-1+i+c *enddo*do,d,1,4*25,1E,54*25+i+d,54*25-1+i+d *enddo*do,e,1,5*25,1E,63*25+i+e,63*25-1+i+e *enddo!2截面type,1mat,1secnum,2*do,a,1,1*25,1E,6*25+i+a,6*25-1+i+a*enddo*do,b,1,2*25,1E,17*25+i+b,17*25-1+i+b *enddo*do,c,1,1*25,1E,21*25+i+c,21*25-1+i+c *enddo*do,d,1,1*25,1E,58*25+i+d,58*25-1+i+dE,61*25+i+e,61*25-1+i+e *enddo*do,f,1,2*25,1E,74*25+i+f,74*25-1+i+f *enddo*do,g,1,1*25,1E,85*25+i+g,85*25-1+i+g *enddo!1截面type,1mat,1secnum,1*do,a,1,2*25,1E,19*25+i+a,19*25-1+i+a *enddo*do,b,1,2*25,1E,59*25+i+b,59*25-1+i+b *enddo*do,c,1,1*25,1E,86*25+i+c,86*25-1+i+c *enddo!右边桁*else!7截面type,1mat,1secnum,7*do,a,1,3*25,1E,i-1+a,i+a*enddo*do,b,1,11*25,1E,35*25-1+i+b,35*25+i+b *enddo*do,c,1,8*25,1E,77*25-1+i+c,77*25+i+c *enddo!6截面type,1mat,1secnum,6*do,a,1,3*25,1E,i+3*25-1+a,i+3*25+aE,32*25-1+i+b,32*25+i+b *enddo*do,c,1,3*25,1E,46*25-1+i+c,46*25+i+c *enddo*do,d,1,1*25,1E,i+76*25-1+d,i+76*25+d *enddo*do,e,1,7*25,1E,87*25-1+i+e,87*25+i+e *enddo!5截面type,1mat,1secnum,5*do,a,1,25,1E,8*25-1+i+a,8*25+i+a*enddo*do,b,1,2*25,1E,30*25-1+i+b,30*25+i+b *enddo*do,c,1,2*25,1E,49*25-1+i+c,49*25+i+c *enddo*do,d,1,2*25,1E,72*25-1+i+d,72*25+i+d *enddo!4截面type,1mat,1secnum,4*do,a,1,3*25,1E,9*25-1+i+a,9*25+i+a*enddo*do,b,1,4*25,1E,26*25-1+i+b,26*25+i+b *enddo*do,c,1,3*25,1E,51*25-1+i+c,51*25+i+c *enddo*do,d,1,4*25,1E,68*25-1+i+d,68*25+i+d *enddomat,1secnum,3*do,a,1,1*25,1E,7*25-1+i+a,7*25+i+a*enddo*do,b,1,5*25,1E,12*25-1+i+b,12*25+i+b *enddo*do,c,1,4*25,1E,22*25-1+i+c,22*25+i+c *enddo*do,d,1,4*25,1E,54*25-1+i+d,54*25+i+d *enddo*do,e,1,5*25,1E,63*25-1+i+e,63*25+i+e *enddo!2截面type,1mat,1secnum,2*do,a,1,1*25,1E,6*25-1+i+a,6*25+i+a*enddo*do,b,1,2*25,1E,17*25-1+i+b,17*25+i+b *enddo*do,c,1,1*25,1E,21*25-1+i+c,21*25+i+c *enddo*do,d,1,1*25,1E,58*25-1+i+d,58*25+i+d *enddo*do,e,1,2*25,1E,61*25-1+i+e,61*25+i+e *enddo*do,f,1,2*25,1E,74*25-1+i+f,74*25+i+f *enddo*do,g,1,1*25,1E,85*25-1+i+g,85*25+i+g *enddo!1截面secnum,1*do,a,1,2*25,1E,19*25-1+i+a,19*25+i+a*enddo*do,b,1,2*25,1E,59*25-1+i+b,59*25+i+b*enddo*do,c,1,1*25,1E,86*25-1+i+c,86*25+i+c*enddo*endif*SET,j,j+14*enddo!*********************************************!!!T_beam AND shell63's nodes!!!*********************************************/prep7ngen,2,7053,1,2351,1,0,0,0.7,ngen,4,2351,7054,7054+2350,1,0,0.57,0,ngen,4,2351,7054+3*2351,7054+3*2351+2350,1,0,0.60,0, ngen,3,2351,7054+6*2351,7054+6*2351+2350,1,0,0.45,0, ngen,2,2351,7054+8*2351,7054+8*2351+2350,1,0,0.60,0, ngen,3,2351,7054+9*2351,7054+9*2351+2350,1,0,0.45,0, ngen,2,2351,7054+11*2351,7054+11*2351+2350,1,0,0.70,0,local,11,0,0,-7,0nsym,y,2351*12,7054+2351,32914+2*2351,1csys,0nsym,y,58775,7054,65828,1ngen,2,122252,2352,2352+2350,1,0,0,0.7,!生成边上4条U肋type,1*do,h,6,9,1*do,k,1,2350,1e,2351*h+k,2351*h+1+k *enddo*enddo!生成T梁1type,1mat,1secnum,27*do,k,1,2350,1e,2351*10+k,2351*10+1+k *enddo!生成U肋2条type,1mat,1secnum,26*do,h,11,12,1*do,k,1,2350,1e,2351*h+k,2351*h+1+k *enddo*enddo!生成T梁2type,1mat,1secnum,27*do,k,1,2350,1e,2351*13+k,2351*13+1+k *enddo!生成U肋另外2条type,1mat,1secnum,26*do,h,14,15,1*do,k,1,2350,1e,2351*h+k,2351*h+1+k *enddo*enddo!镜像四分之一纵向的单元!（最好用*get 最大节点号单元号）esym,,2351*12,7051,30550,1!生成横梁type,1mat,1secnum,29*do,h,7054,9404,25*do,k,h,h+11*2351,2351e,k,k+2351*enddo*enddo*do,h,37617,39967,25*do,k,h,h+10*2351,2351e,k,k+2351*enddo*enddo*do,h,35266,37616,25e,h,h+2351*12*enddo!生成横肋type,1mat,1secnum,28*do,ah,7059,9384,25*do,ac,0,15,5*do,ak,ah+ac,ah+ac+11*2351,2351 e,ak,ak+2351*enddo*enddo*enddo*do,ah,37622,39947,25*do,ac,0,15,5*do,ak,ah+ac,ah+ac+10*2351,2351 e,ak,ak+2351*enddo*enddo*enddo*do,ac,0,15,5e,ah+ac,ah+ac+2351*12 *enddo*enddo!镜像一半单元csys,0esym,,58775,7051,65354,1!补充横梁type,1mat,1secnum,29*do,h,37617,39967,25 e,h,h+86987*enddo*do,h,96392,98742,25 e,h,h+28212*enddo!补充横肋type,1mat,1secnum,28*do,ah,37622,39947,25*do,ac,0,15,5e,ah+ac,ah+ac+86987*enddo*enddo*do,ah,96397,98722,25*do,ac,0,15,5e,ah+ac,ah+ac+28212*enddo*enddo!刚臂type,1mat,3secnum,100*do,h,1,2351,1e,h,h+7053*enddo*do,h,4703,7053,1*enddo*do,h,2352,4702,1e,h,h+124604-2352*enddo!********************************************* !!!shell!!!********************************************* type,2mat,1real,2*do,i,0,11,1*do,h,7054,7054+2349,1e,h+i*2351,h+1+i*2351,h+2352+i*2351,h+2351+i*2351 *enddo*enddo*do,h,35266,35266+2349,1e,h,h+1,h+63479-35266,h+63478-35266*enddo*do,i,0,-10,-1*do,h,63478,63478+2349,1e,h+i*2351,h+1+i*2351,h+i*2351-2350,h+i*2351-2351 *enddo*enddo*do,h,37617,37617+2349,1e,h,h+1,h+1+124604-37617,h+124604-37617*enddocsys,0esym,,58775,131654,190403,1*do,h,124604,124604+2349,1e,h,h+1,h+1-(124604-96392),h-(124604-96392)*enddo!!!ShangXianGan!!!********************************************* /prep7real,1*SET,Nmax1,300000*SET,j,-14*do,i,Nmax1+1,Nmax1+(92*25+1)*2,92*25+1n,i,14.50,j,15.00n,92*25+i,1348.50,j,15.00FILL,i,92*25+i,,,,,,!中桁*if,i,eq,Nmax1+92*25+2,then!13截面type,1mat,1secnum,20*do,a,1,2*25,1E,i-1+a,i+a*enddo*do,b,1,6*25,1E,3*25-1+i+b,3*25+i+b*enddo*do,c,1,3*25,1E,70*25-1+i+c,70*25+i+c*enddo*do,d,1,10*25,1E,74*25-1+i+d,74*25+i+d*enddo*do,e,1,3*25,1E,85*25-1+i+e,85*25+i+e*enddo*do,f,1,1*25,1E,91*25-1+i+f,91*25+i+f*enddo!12截面type,1mat,1secnum,19*do,a,1,1*25,1E,i+2*25-1+a,i+2*25+a*do,b,1,1*25,1E,9*25-1+i+b,9*25+i+b *enddo*do,c,1,8*25,1E,35*25-1+i+c,35*25+i+c *enddo*do,d,1,1*25,1E,i+69*25-1+d,i+69*25+d *enddo*do,e,1,1*25,1E,73*25-1+i+e,73*25+i+e *enddo*do,f,1,3*25,1E,88*25-1+i+f,88*25+i+f *enddo!11截面type,1mat,1secnum,18*do,a,1,25,1E,10*25-1+i+a,10*25+i+a *enddo*do,b,1,1*25,1E,18*25-1+i+b,18*25+i+b *enddo*do,c,1,1*25,1E,26*25-1+i+c,26*25+i+c *enddo*do,d,1,4*25,1E,31*25-1+i+d,31*25+i+d *enddo*do,e,1,1*25,1E,43*25-1+i+e,43*25+i+e *enddo*do,f,1,5*25,1E,45*25-1+i+f,45*25+i+f *enddo*do,g,1,1*25,1E,59*25-1+i+g,59*25+i+g *enddo*do,h,1,1*25,1E,68*25-1+i+h,68*25+i+h *enddo!10截面mat,1secnum,17*do,a,1,1*25,1E,11*25-1+i+a,11*25+i+a *enddo*do,b,1,1*25,1E,17*25-1+i+b,17*25+i+b *enddo*do,c,1,1*25,1E,20*25-1+i+c,20*25+i+c *enddo*do,d,1,1*25,1E,25*25-1+i+d,25*25+i+d *enddo*do,e,1,4*25,1E,27*25-1+i+e,27*25+i+e *enddo*do,f,1,3*25,1E,50*25-1+i+f,50*25+i+f*enddo*do,g,1,1*25,1E,57*25-1+i+g,57*25+i+g *enddo*do,h,1,1*25,1E,60*25-1+i+h,60*25+i+h *enddo*do,hh,1,1*25,1E,67*25-1+i+hh,67*25+i+hh *enddo!9截面type,1mat,1secnum,16*do,a,1,5*25,1E,12*25-1+i+a,12*25+i+a *enddo*do,b,1,1*25,1E,19*25-1+i+b,19*25+i+b *enddo*do,c,1,4*25,1E,21*25-1+i+c,21*25+i+c *enddo*do,d,1,4*25,1E,53*25-1+i+d,53*25+i+dE,58*25-1+i+e,58*25+i+e*enddo*do,g,1,6*25,1E,61*25-1+i+g,61*25+i+g*enddo*do,h,1,1*25,1E,84*25-1+i+h,84*25+i+h*enddo!8截面type,1mat,1secnum,15*do,a,1,1*25,1E,44*25-1+i+a,44*25+i+a*enddo!左边桁*elseif,i,eq,Nmax1+(92*25+1)*2+1,then !13截面type,1mat,1secnum,20*do,a,1,2*25,1E,i-1+a,i+a*enddo*do,b,1,6*25,1E,3*25-1+i+b,3*25+i+b*enddo*do,c,1,3*25,1E,70*25-1+i+c,70*25+i+c*enddo*do,d,1,10*25,1E,74*25-1+i+d,74*25+i+d*enddo*do,e,1,3*25,1E,85*25-1+i+e,85*25+i+e*enddo*do,f,1,1*25,1E,91*25-1+i+f,91*25+i+f*enddo!12截面type,1mat,1E,i+2*25-1+a,i+2*25+a*enddo*do,b,1,1*25,1E,9*25-1+i+b,9*25+i+b *enddo*do,c,1,8*25,1E,35*25-1+i+c,35*25+i+c *enddo*do,d,1,1*25,1E,i+69*25-1+d,i+69*25+d *enddo*do,e,1,1*25,1E,73*25-1+i+e,73*25+i+e *enddo*do,f,1,3*25,1E,88*25-1+i+f,88*25+i+f *enddo!11截面type,1mat,1secnum,18*do,a,1,25,1E,10*25-1+i+a,10*25+i+a *enddo*do,b,1,1*25,1E,18*25-1+i+b,18*25+i+b *enddo*do,c,1,1*25,1E,26*25-1+i+c,26*25+i+c *enddo*do,d,1,4*25,1E,31*25-1+i+d,31*25+i+d *enddo*do,e,1,1*25,1E,43*25-1+i+e,43*25+i+e *enddo*do,f,1,5*25,1E,45*25-1+i+f,45*25+i+f *enddo*do,g,1,1*25,1E,59*25-1+i+g,59*25+i+g *enddo*do,h,1,1*25,1*enddo!10截面type,1mat,1secnum,17*do,a,1,1*25,1E,11*25-1+i+a,11*25+i+a *enddo*do,b,1,1*25,1E,17*25-1+i+b,17*25+i+b *enddo*do,c,1,1*25,1E,20*25-1+i+c,20*25+i+c *enddo*do,d,1,1*25,1E,25*25-1+i+d,25*25+i+d *enddo*do,e,1,4*25,1E,27*25-1+i+e,27*25+i+e *enddo*do,f,1,3*25,1E,50*25-1+i+f,50*25+i+f*enddo*do,g,1,1*25,1E,57*25-1+i+g,57*25+i+g *enddo*do,h,1,1*25,1E,60*25-1+i+h,60*25+i+h *enddo*do,hh,1,1*25,1E,67*25-1+i+hh,67*25+i+hh *enddo!9截面type,1mat,1secnum,16*do,a,1,5*25,1E,12*25-1+i+a,12*25+i+a *enddo*do,b,1,1*25,1E,19*25-1+i+b,19*25+i+b *enddo*do,c,1,4*25,1E,21*25-1+i+c,21*25+i+c*do,d,1,4*25,1E,53*25-1+i+d,53*25+i+d *enddo*do,e,1,1*25,1E,58*25-1+i+e,58*25+i+e *enddo*do,g,1,6*25,1E,61*25-1+i+g,61*25+i+g *enddo*do,h,1,1*25,1E,84*25-1+i+h,84*25+i+h *enddo!8截面type,1mat,1secnum,15*do,a,1,1*25,1E,44*25-1+i+a,44*25+i+a *enddo!右边桁*else!13截面type,1mat,1secnum,20*do,a,1,2*25,1E,i-1+a,i+a*enddo*do,b,1,6*25,1E,3*25-1+i+b,3*25+i+b *enddo*do,c,1,3*25,1E,70*25-1+i+c,70*25+i+c *enddo*do,d,1,10*25,1E,74*25-1+i+d,74*25+i+d *enddo*do,e,1,3*25,1E,85*25-1+i+e,85*25+i+e *enddo*do,f,1,1*25,1E,91*25-1+i+f,91*25+i+f *enddotype,1mat,1secnum,19*do,a,1,1*25,1E,i+2*25-1+a,i+2*25+a*enddo*do,b,1,1*25,1E,9*25-1+i+b,9*25+i+b *enddo*do,c,1,8*25,1E,35*25-1+i+c,35*25+i+c *enddo*do,d,1,1*25,1E,i+69*25-1+d,i+69*25+d *enddo*do,e,1,1*25,1E,73*25-1+i+e,73*25+i+e *enddo*do,f,1,3*25,1E,88*25-1+i+f,88*25+i+f *enddo!11截面type,1mat,1secnum,18*do,a,1,25,1E,10*25-1+i+a,10*25+i+a *enddo*do,b,1,1*25,1E,18*25-1+i+b,18*25+i+b *enddo*do,c,1,1*25,1E,26*25-1+i+c,26*25+i+c *enddo*do,d,1,4*25,1E,31*25-1+i+d,31*25+i+d *enddo*do,e,1,1*25,1E,43*25-1+i+e,43*25+i+e *enddo*do,f,1,5*25,1E,45*25-1+i+f,45*25+i+f *enddo*do,g,1,1*25,1*enddo*do,h,1,1*25,1E,68*25-1+i+h,68*25+i+h *enddo!10截面type,1mat,1secnum,17*do,a,1,1*25,1E,11*25-1+i+a,11*25+i+a *enddo*do,b,1,1*25,1E,17*25-1+i+b,17*25+i+b *enddo*do,c,1,1*25,1E,20*25-1+i+c,20*25+i+c *enddo*do,d,1,1*25,1E,25*25-1+i+d,25*25+i+d *enddo*do,e,1,4*25,1E,27*25-1+i+e,27*25+i+e *enddo*do,f,1,3*25,1E,50*25-1+i+f,50*25+i+f*enddo*do,g,1,1*25,1E,57*25-1+i+g,57*25+i+g *enddo*do,h,1,1*25,1E,60*25-1+i+h,60*25+i+h *enddo*do,hh,1,1*25,1E,67*25-1+i+hh,67*25+i+hh *enddo!9截面type,1mat,1secnum,16*do,a,1,5*25,1E,12*25-1+i+a,12*25+i+a *enddo*do,b,1,1*25,1E,19*25-1+i+b,19*25+i+b*do,c,1,4*25,1E,21*25-1+i+c,21*25+i+c*enddo*do,d,1,4*25,1E,53*25-1+i+d,53*25+i+d*enddo*do,e,1,1*25,1E,58*25-1+i+e,58*25+i+e*enddo*do,g,1,6*25,1E,61*25-1+i+g,61*25+i+g*enddo*do,h,1,1*25,1E,84*25-1+i+h,84*25+i+h*enddo!8截面type,1mat,1secnum,15*do,a,1,1*25,1E,44*25-1+i+a,44*25+i+a*enddo*endif*SET,j,j+14*enddo!********************************************* !!!FuGan!!!********************************************* /prep7*SET,j,-14*do,i,1,(94*25+1)*2,94*25+1k,i,0,j,0k,94*25+i,1363.000,j,0kFILL,i,94*25+i,93,,,,,j=j+14*enddo*SET,Nmax1,10000*SET,j,-14k,i,14.50,j,15.00k,92*25+i,1348.50,j,15.00 kFILL,i,92*25+i,91,,,,,j=j+14*enddo!连接右斜线1j=Nmax1+1*do,n,1,(94*25+1)*2,94*25+1 kk=0*do,i,n,n-1+20*25,25lstr,i,jkk=kk+1*if,kk,eq,1,thenlatt,1,1,1,,,,21esize,,25lmesh,all*elseif,kk,ge,4,and,kk,le,6,then latt,1,1,1,,,,21esize,,25lmesh,all*elseif,kk,ge,8,and,kk,le,9,then latt,1,1,1,,,,21esize,,25lmesh,all*elseif,kk,eq,7,thenlatt,1,1,1,,,,23esize,,25lmesh,all*elselatt,1,1,1,,,,23esize,,25lmesh,all*endifj=j+25*enddoj=j+2301-(20-0)*25*enddo!连接右斜线2j=Nmax1+1+40*25*do,n,1,(94*25+1)*2,94*25+1 kk=0*do,i,n+40*25,n-1+60*25,25 lstr,i,jkk=kk+1latt,1,1,1,,,,21esize,,25lmesh,all*elseif,kk,ge,3,and,kk,le,9,then latt,1,1,1,,,,21esize,,25lmesh,all*elselatt,1,1,1,,,,22esize,,25lmesh,all*endifj=j+25*enddoj=j+2301-(60-40)*25*enddo!连接左斜线1j=Nmax1+19*25+1*do,n,1,(94*25+1)*2,94*25+1 kk=0*do,i,n+21*25,n-1+41*25,25 lstr,i,jkk=kk+1*if,kk,eq,1,thenlatt,1,1,1,,,,21esize,,25lmesh,all*elseif,kk,ge,13,and,kk,le,18,then latt,1,1,1,,,,21esize,,25lmesh,all*elselatt,1,1,1,,,,22esize,,25lmesh,all*endifj=j+25*enddoj=j+2301-(41-21)*25*enddo!连接左斜线2j=Nmax1+59*25+1*do,n,1,(94*25+1)*2,94*25+1 kk=0kk=kk+1*if,kk,ge,1,and,kk,le,12,then latt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,ge,21,and,kk,le,22,then latt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,ge,15,and,kk,le,16,then latt,1,1,1,,,,23esize,,25lmesh,all*elseif,kk,ge,24,and,kk,le,28,then latt,1,1,1,,,,23esize,,25lmesh,all*elseif,kk,ge,33,and,kk,le,34,then latt,1,1,1,,,,23esize,,25lmesh,all*elselatt,1,1,1,,,,21esize,,25lmesh,all*endifj=j+25*enddoj=j+2301-(95-61)*25*enddo!连接shu线j=Nmax1+1*do,n,1,(94*25+1)*2,94*25+1 kk=0*do,i,n+1*25,n-1+94*25,25 lstr,i,jkk=kk+1*if,kk,ge,1,and,kk,le,3,then latt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,ge,22,and,kk,le,58,then latt,1,1,1,,,,22*elseif,kk,ge,62,and,kk,le,74,then latt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,ge,78,and,kk,le,83,then latt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,ge,89,and,kk,le,93,then latt,1,1,1,,,,22esize,,25lmesh,all*elseif,kk,eq,7,thenlatt,1,,1,,,,23esize,,25lmesh,all*elseif,kk,eq,60,thenlatt,1,1,1,,,,23esize,,25lmesh,all*elseif,kk,eq,75,thenlatt,1,1,1,,,,23esize,,25lmesh,all*elseif,kk,eq,86,thenlatt,1,1,1,,,,23esize,,25lmesh,all*elseif,kk,eq,87,thenlatt,1,1,1,,,,23esize,,25lmesh,all*elselatt,1,1,1,,,,21esize,,25lmesh,all*endifj=j+25*enddoj=j+2301-93*25*enddomodmsh,detackldele,all!********************************************* !!!2X+上平联!!!********************************************* /prep7k,1,14.5,-14,15.00k,2,14.5,-7,15.00k,3,14.5,-14,10.2k,4,14.5,-7,10.2lstr,1,4lstr,2,3lstr,3,4lgen,2,all,,,,7,nummrg,kpnumcmp,kplstr,2,4lgen,93,all,,,14.5,,lgen,2,all,,,,14,,latt,1,,1,,,,24esize,,10lmesh,allmodmsh,detackldele,allkdele,all*SET,j,-14*do,i,1,(92*25+1)*2,92*25+1k,i,14.50,j,15.00k,92*25+i,1348.50,j,15.00kFILL,i,92*25+i,91,,,,,j=j+14*enddo*do,i,1,92*25+1,25lstr,i,i+92*25+1lstr,i+92*25+1,i+(92*25+1)*2*enddolatt,1,,1,,,,31esize,,20lmesh,allmodmsh,detackldele,all*do,i,1,92*25,25lstr,i+92*25+1,i+(92*25+1)*2+1*25*enddo*do,i,1+1*25,92*25+1,25lstr,i,i+92*25+1-25lstr,i+92*25+1,i+(92*25+1)*2-25*enddonummrg,kpnumcmp,kplatt,1,,1,,,,32esize,,20lmesh,allmodmsh,detackldele,allkdele,all!************************************************! !!1#ta!!!************************************************ /prep7nn=400000n,1+nn,290,0,158.9n,13+nn,290,0,132.5n,20+nn,290,0,115n,22+nn,290,0,110.5n,63+nn,290,0,110.5n,52+nn,290,-19,-6n,62+nn,290,-11.5,-34.5n,93+nn,290,19,-6n,103+nn,290,11.5,-34.5n,104+nn,290,-19,-6n,107+nn,290,19,-6fill,1+nn,13+nn,11,fill,13+nn,20+nn,6,fill,20+nn,22+nn,1,fill,22+nn,52+nn,29,fill,52+nn,62+nn,9,fill,63+nn,93+nn,29,fill,93+nn,103+nn,9,fill,104+nn,107+nn,2,real,1type,1mat,2*do,a,1+nn,12+nn,1 E,a,a+1*enddotype,1mat,2secnum,52*do,b,13+nn,19+nn,1 E,b,b+1*enddotype,1mat,2secnum,51*do,b,20+nn,21+nn,1 E,b,b+1*enddotype,1mat,2secnum,46*do,b,22+nn,51+nn,1 E,b,b+1*enddotype,1mat,2secnum,42*do,b,52+nn,61+nn,1 E,b,b+1*enddotype,1mat,2secnum,46*do,b,63+nn,92+nn,1 E,b,b+1*enddotype,1mat,2secnum,42*do,b,93+nn,102+nn,1 E,b,b+1*enddotype,1mat,2secnum,54*do,b,104+nn,106+nn,1 E,b,b+1!************************************************ !!2#ta!!************************************************* nn=450000n,1+nn,870,0,158.9n,13+nn,870,0,132.5n,20+nn,870,0,115n,22+nn,870,0,110.5n,63+nn,870,0,110.5n,52+nn,870,-19,-6n,62+nn,870,-11.5,-34.5n,93+nn,870,19,-6n,103+nn,870,11.5,-34.5n,104+nn,870,-19,-6n,107+nn,870,19,-6fill,1+nn,13+nn,11,fill,13+nn,20+nn,6,fill,20+nn,22+nn,1,fill,22+nn,52+nn,29,fill,52+nn,62+nn,9,fill,63+nn,93+nn,29,fill,93+nn,103+nn,9,fill,104+nn,107+nn,2,real,1type,1mat,2secnum,53*do,a,1+nn,12+nn,1E,a,a+1*enddotype,1mat,2secnum,52*do,b,13+nn,19+nn,1E,b,b+1*enddotype,1mat,2secnum,51*do,b,20+nn,21+nn,1E,b,b+1*enddosecnum,46*do,b,22+nn,51+nn,1E,b,b+1*enddotype,1mat,2secnum,42*do,b,52+nn,61+nn,1E,b,b+1*enddotype,1mat,2secnum,46*do,b,63+nn,92+nn,1E,b,b+1*enddotype,1mat,2secnum,42*do,b,93+nn,102+nn,1E,b,b+1*enddotype,1mat,2secnum,54*do,b,104+nn,106+nn,1E,b,b+1*enddo!************************************************ !!link!!************************************************ /prep7nn=300000mm=400000!*****************************ES*do,j,nn+1,nn+2301*2+1,2301*do,a,j,10*25+j,25type,3mat,1real,3E,a,mm+3mm=400000*enddo*do,j,nn+1+11*25,nn+2301*3,2301 *do,a,j,2*25+j,25type,3mat,1real,4E,a,mm+14mm=mm+1*enddomm=400000*enddo*do,j,nn+1+14*25,nn+2301*3,2301 *do,a,j,2*25+j,25type,3mat,1real,5E,a,mm+17mm=mm+1*enddomm=400000*enddo*do,j,nn+1+17*25,nn+2301*3,2301 type,3mat,1real,6E,j,mm+20mm=400000*enddo!****************************EM mm=400000*do,j,nn+1+21*25,nn+2301*3,2301 type,3mat,1real,6E,j,mm+20mm=400000*enddo*do,j,nn+1+22*25,nn+2301*3,2301 *do,a,j,2*25+j,25type,3mat,1real,5mm=mm-1*enddomm=400000*enddo*do,j,nn+1+25*25,nn+2301*3,2301*do,a,j,2*25+j,25type,3mat,1real,4E,a,mm+16mm=mm-1*enddomm=400000*enddo*do,j,nn+1+28*25,nn+2301*3,2301*do,a,j,10*25+j,25type,3mat,1real,3E,a,mm+13mm=mm-1*enddomm=400000*enddo!***************************************WM mm=450000*do,j,nn+1+40*25,nn+2301*3,2301*do,a,j,10*25+j,25type,3mat,1real,3E,a,mm+3mm=mm+1*enddomm=450000*enddo*do,j,nn+1+51*25,nn+2301*3,2301*do,a,j,2*25+j,25type,3mat,1real,4E,a,mm+14mm=mm+1*enddo*enddo*do,j,nn+1+54*25,nn+2301*3,2301*do,a,j,2*25+j,25type,3mat,1real,5E,a,mm+17mm=mm+1*enddomm=450000*enddo*do,j,nn+1+57*25,nn+2301*3,2301type,3mat,1real,6E,j,mm+20mm=450000*enddo!***************************************WS mm=450000*do,j,nn+1+61*25,nn+2301*3,2301type,3mat,1real,6E,j,mm+20mm=450000*enddo*do,j,nn+1+62*25,nn+2301*3,2301*do,a,j,2*25+j,25type,3mat,1real,5E,a,mm+19mm=mm-1*enddomm=450000*enddo*do,j,nn+1+65*25,nn+2301*3,2301*do,a,j,2*25+j,25type,3mat,1real,4E,a,mm+16mm=mm-1mm=450000*enddo*do,j,nn+1+68*25,nn+2301*3,2301 *do,a,j,10*25+j,25type,3mat,1real,3E,a,mm+13mm=mm-1*enddomm=450000*enddonummrg,nodenumcmp,nodeFINISH。
ansys命令流中文说明展开全文KB、KE: 待划分线的定向关键点起始、终止号SECNUM: 截面类型号u SECPLOT,SECID,MESHKEY 画梁截面的几何形状及网格划分SECID:由SECTYPE命令分配的截面编号MESHKEY：0：不显示网格划分1：显示网格划分u /ESHAPE, SCALE 按看似固体化分的形式显示线、面单元SCALE: 0:简单显示线、面单元1：使用实常数显示单元形状u esurf, xnode, tlab, shape 在已存在的选中单元的自由表面覆盖产生单元xnode: 仅为产生surf151 或surf152单元时使用tlab: 仅用来生成接触元或目标元top 产生单元且法线方向与所覆盖的单元相同，仅对梁或壳有效，对实体单元无效Bottom产生单元且法线方向与所覆盖的单元相反，仅对梁或壳有效，对实体单元无效Reverse 将已产生单元反向Shape: 空与所覆盖单元形状相同Tri 产生三角形表面的目标元注意：选中的单元是由所选节点决定的，而不是选单元，如同将压力加在节点上而不是单元上u Nummrg,label,toler, Gtoler,action,switch 合并相同位置的itemlabel: 要合并的项目node: 节点， Elem,单元，kp: 关键点（也合并线，面及点）mat: 材料，type: 单元类型，Real: 实常数cp：耦合项，CE：约束项，CE: 约束方程，All：所有项toler: 公差Gtoler：实体公差Action: sele 仅选择不合并空合并switch: 较低号还是较高号被保留（low, high）注意：可以先选择一部分项目，再执行合并。
如果多次发生合并命令，一定要先合并节点，再合并关键点。
合并节点后，实体荷载不能转化到单元，此时可合并关键点解决问题。
u Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线type: s 从全部线中选一组线r 从当前选中线中选一组线a 再选一部线附加给当前选中组aunoneu(unselect)inve: 反向选择item: line 线号loc 坐标length 线长comp: x,y,zkswp: 0 只选线1 选择线及相关关键点、节点和单元u Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0” 使用正负号“1”仅用绝对值u NSLL,type, nkey 选择与所选线相联系的节点u nsla, type, nkey: 选择与选中面相关的节点type：s 选一套新节点r 从已选节点中再选a 附加一部分节点到已选节点u 从已选节点中去除一部分nkey: 0 仅选面内的节点1 选所有和面相联系的节点（如面内线，关键点处的节点）u esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元Type: S: 选择一组单元（缺省）R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组（？）Stat: 显示当前选择状态Item： Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号u ALLSEL, LABT, ENTITY 选中所有项目LABT: ALL: 选所有项目及其低级项目BELOW: 选指定项目的直接下属及更低级项目ENTITY: ALL: 所有项目（缺省）VOLU:体高级AREA:面LINE :线KP:关键点ELEM:单元NODE:节点低级u Tshap,shape 定义接触目标面为2D、3D的简单图形Shape: line:直线Arc:顺时针弧Tria:3点三角形Quad:4点四边形………….2.6 根据需要耦合某些节点自由度u cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotznode1-node17: 待耦合的节点号。
