装配式混凝土柱钢梁节点ANSYS命令流

合集下载

钢管混凝土ansys命令流

碳纤维加固桥梁结构的健康监测摘要：提出了一种针对老化桥梁的健康监测的可靠性方法。

首先, 根据桥梁的类型提出了临界失效标准。

然后基于这些标准定义安全界限。

假定大量的失效标准是正态分布的随机变量。

基于这样的假设,基本可靠性指标和基本失效概率为每个关键失效模式估计。

已经发现基本失效概率值,桥梁系统失效概率值按这段时间考虑计算。

最后,用系统失效概率来获得桥梁系统的可靠性指标。

系统的可靠性指数用来表达的桥梁在使用期间的状况较为片面。

选择斯里兰卡国家铁路桥网络的一座铁路桥梁作为一个案例研究来说明该可靠度程序。

对于这座桥,失效的最重要的形式是疲劳和腐蚀。

应用程序的建议方法表明,桥梁在现有负载的情况下的现状是令人满意的。

对可靠性的预测,不是为了抵御破坏,而是一种对基础设施管理进行有效的积极的健康监测方式。

关键字：碳纤维结构健康监测状态评估维护可靠性指标失效概率1、介绍在当今迅速变化的世界中,人类的生活质量和国家经济发展依赖于基础设施的数量、质量、和有效性]1[。

在世界各国，民用基础设施通常是需要维护的,康复,或更换的]2[。

特别是对基础设施构成一个相当大的投资的高速公路和相关的地面运输系统，直接影响到一个国家的工业产值。

目前大多数桥梁养护策略基于由桥梁检查员在不同时间间隔进行的视觉检查]3[。

由于人类检测取决于个人检查员,其结果存在一定程度的不确定性。

因此,绝对的依赖检验报告是不可靠的，而且可能导致错误的决策和桥梁拥有者的高维护成本。

基于桥梁的剩余使用寿命维护成本的优化。

另一方面,桥梁的检查员发现结构和非结构化缺陷,破坏和恶化过程的能力不应该被忽视。

事实上,人类经验应该被纳入维护决策,但要有足够的主观的知识。

只在那时,任何可接受的战略才能成为可能。

当前桥梁结构健康监测方法存在许多缺点和普遍的不确定性。

对结构健康监测和桥梁维护的知识不足造成由于频繁的维护活动,不可预见的桥梁损害,和偶然的桥梁失效引起的重大的经济损失。

ansys命令流中文说明

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: 待耦合的节点号。

ANSYS工程结构数值分析---命令流实例

!------------------------------------------------------!EX8.26 钢筋混凝土简支梁数值分析!分离式模型,关闭压碎,keyopt(1)=0,keyopt(7)=1 !力加载,位移收敛准则，误差1.5%,1/4模型分析!--------------------------------------------finish/clear/config,nres,2000/prep7!1.定义单元与材料性质--------------------et,1,solid65,,,,,,,1 !K1=0,k7=1et,2,link8mp,ex,1,13585mp,prxy,1,0.2fc=14.3ft=1.43tb,concr,1tbdata,,0.5,0.95,ft,-1tb,miso,1,,11tbpt,,0.0002,fc*0.19tbpt,,0.0004,fc*0.36tbpt,,0.0006,fc*0.51tbpt,,0.0008,fc*0.64tbpt,,0.001,fc*0.75tbpt,,0.0012,fc*0.84tbpt,,0.0014,fc*0.91tbpt,,0.0016,fc*0.96tbpt,,0.0018,fc*0.99tbpt,,0.002,fctbpt,,0.0033,fc*0.85mp,ex,2,2.0e5mp,prxy,2,0.3tb,biso,2tbdata,,300,0pi=acos(-1)r,1,0.25*pi*22*22r,2,0.25*pi*22*22/2r,3,0.25*pi*10*10r,4,0.25*pi*10*10/2!2.创建几何模型blc4,,,150/2,300,2000/2*do,i,1,9wpoff,,,100vsbw,all*enddo wpcsys,-1 wpoff,,,50 vsbw,all wpcsys,-1 wprota,,-90 wpoff,,,30 vsbw,all wpoff,,,240 vsbw,all wpcsys,-1 wpoff,30 wprota,,,90 vsbw,all wpcsys,-1!3.划分钢筋网格elemsiz=50 lsel,s,loc,x,30 lsel,r,loc,y,30 cm,zj,linelatt,2,1,2 lesize,all,elemsiz lsel,s,loc,x,75 lsel,r,loc,y,30 cm,zjb,line latt,2,2,2 lesize,all,elemsiz lsel,s,loc,x,30 lsel,r,loc,y,270 cm,jlj,linelatt,2,3,2 lesize,all,elemsiz lsel,s,tan1,z lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 lsel,u,loc,z,50 cm,gj,linelatt,2,3,2 lesize,all,elemsiz lsel,s,loc,z,0 lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 cm,gjb,line latt,2,4,2lesize,all,elemsizlsel,allcmsel,s,zjcmsel,a,zjbcmsel,a,jljcmsel,a,gjcmsel,a,gjbcm,gj,linelmesh,alllsel,all!4.划分混凝土网格vatt,1,,1mshkey,1esize,elemsizvmesh,allallsel,all!5.施加荷载和约束lsel,s,loc,y,0lsel,r,loc,z,900dl,all,,uyasel,s,loc,z,0da,all,symmasel,s,loc,x,75da,all,symm!-----------------p0=180000q0=p0/150/100asel,s,loc,z,0,50asel,r,loc,y,300sfa,all,1,pres,q0allsel,all!6.求解控制设置/soluantype,0nsubst,80outres,all,allautos,onneqit,50cnvtol,u,,0.015solve!7.进入POST1查看结果/post1set,lastesel,s,type,,2etable,saxl,ls,1plls,saxl,saxlesel,s,type,,1/device,vector,onplcrackeall!8.进入时程后处理/post26nsol,2,205,u,yprod,3,1,,,,,,p0/1000prod,4,2,,,,,,-1xvar,4plvar,3!EX2.6D 60kg/m钢轨实体建模FINISH/CLEAR/PREP7!1.定义参数H1=30.5H2=48.5H3=48.5B1=45.75B2=29.25B3=16.5B4=36.5B5=10XD1=1/3XD2=1/9XD3=1/3XD4=1/20R0=400R1=2R2=4R3=40R4=20R5=25R7=13R8=80R9=300R10=2FAI=43H=H1+2*H2+H3!2.创建轨底主要组成部分的线K,1K,2,B1+B2K,3,B1+B2,H1-B1*XD1-B2*XD2K,4,B1,H1-B1*XD1K,5,0,H1*DO,I,1,4L,I,I+1*ENDDO!3.计算确定轨头下三个圆弧的位置CTA=A TAN(XD3)BTA=(90-FAI)*ACOS(-1)/180-CTAA=1/TAN(CTA)C1=B3/2+R0-R6*SIN(CTA)-(R5-R6)*COS(BTA)C2=R5*SIN(BTA)+R6*COS(CTA)-R6*SIN(BTA)-H2C=(C1-C2/TAN(CTA))/(R0-R5)C3=A*C-SQRT(A*A-C*C+1)C3=C3/(A*A+1)REFA=ASIN(C3)C1=R6*SIN(CTA)+(R5-R6)*COS(BTA)+(R0-R5)*COS(REFA) !4.创建轨腰和轨头结构线ROUX=B3/2+R0-C1ROUY=ROUX*TAN(CTA)ANGI=ASIN((H1+H2)/R0)K,6,B3/2+R0-R0*COS(ANGI)K,7,B3/2+R0-R0*COS(REFA),H1+H2+R0*SIN(REFA) LARC,6,7,2,R0DX1=R6*(COS(BTA)-SIN(CTA))DY1=R6*(COS(CTA)-SIN(BTA))K,8,ROUX-DX1,H1+2*H2+ROUY-DY1K,9,ROUX,H1+2*H2+ROUYLARC,7,8,2,R5LARC,8,9,2,R6K,10,B4,H1+2*H2+B4*XD3K,11,B4-(H3-B4*XD3)*XD4,HL,9,10L,10,11ANGI=ASIN(B5/R9)K,12,0,HK,13,B5,H-(R9-R9*COS(ANGI)) LARC,12,13,1,R9DX1=B4-B5+R8*SIN(ANGI)ANGI1=ASIN(DX1/R8)Y1=R8*COS(ANGI)-R8*COS(ANGI1) K,14,B4,KY(13)-Y1LARC,13,14,1,R8!5.倒角或弧线连接LFILLT,1,2,R1LFILLT,2,3,R2LFILLT,3,4,R3LPTN,4,5LFILLT,17,18,R4LFILLT,8,9,R10LPTN,9,11LFILLT,20,22,R7LDELE,15,16,1,1LDELE,19,21,2,1LSYMM,X,ALLNUMMRG,ALLNUMCMP,ALLAL,ALLASUMVOFFST,1,300!EX2.9B 斜向圆台建模finish/clear/prep7r0=3r1=2h=4cyl4,,,r0wpoff,r1,,hcyl4,,,r1*do,i,1,4askin,i,i+4*enddova,all!斜向延伸创建圆台vgen,1,1,,,,3*r0,,,,1wpcsyscyl4,,,r0vext,7,,,r1,,h,r1/r0,r1/r0!EX2.9B 已知函数方程时的建模!玫瑰线1==================================== !极坐标方程为ρ=asin(4Φ/3)finish/clear/prep7a=20csys,1*do,i,0,360*3fei=i*acos(-1)/180r0=a*sin(fei*4/3)k,,r0,i*enddo*do,i,1,360*3l,i,i+1*enddo!玫瑰线2=================================== !极坐标方程为ρ=asin(4Φ)finish/clear/prep7a=20csys,1*do,i,0,360fei=i*acos(-1)/180r0=a*sin(fei*4)k,,r0,i*enddo*do,i,1,360l,i,i+1*enddo!长辐圆内旋轮线============================ !参数方程x=(a-b)cost+λcos(a-b)/bt!参数方程y=(a-b)sint-λsin(a-b)/btfinish/clear/prep7a=50b=5lmda=15*afun,deg*do,i,0,360x1=(a-b)*cos(i)+lmda*cos((a-b)/b*i)y1=(a-b)*sin(i)-lmda*sin((a-b)/b*i)k,,x1,y1*enddo*do,i,1,360l,i,i+1*enddo!椭圆抛物面================================ !方程为z=x*x/a/a+y*y/b/b!使用蒙皮方法finish/clear/prep7a=20b=10n=20*do,i,1,niz=i/10lsel,nonewpoff,,,izcyl4,,,1adele,alllsscale,all,,,a*sqrt(iz),b*sqrt(iz),,,,1*enddoallsel*do,i,1,4*(n-1),4askin,i,i+4askin,i+1,i+5askin,i+2,i+6askin,i+3,i+7*enddo!EX2.9C 旋转图素建模!类花瓶建模--------------------finish/clear/prep7K,1K,2,100K,3,300,400K,4,150,550k,5,130,880K,6,300,1150k,7,400,1100k,8,500,1150bsplin,allk,1001,,1000arotat,all,,,,,,1,1001,360,5!类轮建模----------------------finish/clear/prep7!先创建类轮结构的一个断面，然后旋转之blc4,,,4,6cyl4,2,6,1.5asba,1,2k,100,-3k,101,-3,10vrotat,all,,,,,,100,101!EX3.1 自定义箱形截面finish/clear/prep7k,1k,2,2k,3,2,2.2k,4,3,2.3k,5,3,2.5k,6,0,2.5a,1,2,3,4,5,6blc4,,0.2,1.7,2asba,1,2wpoff,1.7wprota,,,90asbw,allwpoff,,,0.3asbw,allwpoff,,0.2wprota,,90asbw,allwpoff,,,-2asbw,allarsym,x,allwpcsysaglue,allet,1,plane82lesize,all,,,1amesh,allsecwrite,myboxfinish/clear/prep7et,1,beam189mp,ex,1,3.0e10mp,prxy,1,0.167sectype,1,beam,mesh secread,mybox,,,meshk,1k,2,10k,100,5,5l,1,2lesize,all,,,20latt,1,,1,,100,,1lmesh,all/eshape,1Eplot!EX3.2 自定义多种材料截面finish/clear/prep7Ro=1.5Ri=1.0csys,1cyl4,,,ricyl4,,,roaptn,allwprota,,90asbwa,allwprota,,,90asbw,allwpcsyset,1,plane82mymat1=4mymat2=7mp,ex,mymat1,1.0mp,ex,mymat2,2.0asel,s,loc,x,0,riaatt,mymat1,,1asel,s,loc,x,ri,roaatt,mymat2,,1allselesize,0.25mshape,0,2dmshkey,1amesh,allsecwrite,mycsolid,sectfinish/clear/prep7et,1,beam189mym1=4mym2=7mp,ex,mym1,3.0e10mp,prxy,mym1,0.167 mp,ex,mym2,2.1e11mp,prxy,mym2,0.3 sectype,1,beam,mesh secread,mycsolid,sect,,mesh k,1k,2,,,10l,1,2lesize,all,,,20latt,,,1,,,,1lmesh,all/eshape,1/pnum,mat,1eplot!EX3.3A 单个方位关键点示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,ISECDATA,100,40,160,10,10,8K,1K,2,,,1000L,1,2K,100,,500,500LATT,1,,1,,100,,1LGEN,4,1,,,500LESIZE,ALL,,,10LMESH,ALL/ESHAPE,1EPLOT!EX3.3B 同时使用KB和KE方位关键点示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,ISECDATA,100,40,160,10,10,8L0=1000DL=500DXC=400K,1K,2,,,L0L,1,2K,100,,DLK,200,DXC,-DLK,300,2*DXC,DLK,301,2*DXC+DLK,400K,500,8*DXCLGEN,5,1,,,DXCLATT,1,,1,,100,,1LSEL,S,,,2LATT,1,,1,,200,,1LSEL,S,,,3LATT,1,,1,,300,301,1LSEL,S,,,4LATT,1,,1,,400,,1LSEL,S,,,5LATT,1,,1,,500,,1LSEL,ALLLESIZE,ALL,,,50LMESH,ALL/ESHAPE,1EPLOT!EX3.3C 双方位关键点翻转示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,I SECDATA,100,40,160,10,10,8L0=1000DL=500DXC=400K,1K,2,,,L0K,3,DXCK,4,DXC,,L0/2K,5,DXC,,L0L,1,2L,3,4L,4,5K,100,,DLk,101,,-DLK,200,DXC,DLk,201,2*DXCK,202,DXC,-DLLATT,1,,1,,100,101,1lesize,all,,,100LSEL,S,,,2LATT,1,,1,,200,201,1lesize,all,,,50LSEL,S,,,3LATT,1,,1,,201,202,1lesize,all,,,50LSEL,ALLLMESH,ALL/ESHAPE,1EPLOT!EX3.4A 两种单元形状和两种网格划分比较finish/clear/prep7et,1,plane82k,1k,2,8k,3,7,6k,4,1,6a,1,2,3,4esize,1mshape,0mshkey,0!mshape,0!mshkey,1!mshape,1!mshkey,0!mshape,1!mshkey,1amesh,all!EX3.4B 中间节点位置控制网格划分比较finish/clear/prep7cyl4,,,4,,8,60lesize,all,,,2mshape,0mshkey,1mshmid,0!mshmid,1!mshmid,2amesh,all!EX3.5 线上单元尺寸设置示例!下边密上边稀finish/clear/prep7et,1,plane82blc4,,,10,10lsel,s,tan1,ylesize,all,,,10lsel,s,loc,x,0lesize,all,,,9,1/8lsel,s,loc,x,10lesize,all,,,9,8lsel,allmshape,0mshkey,1amesh,all!中间密外边稀finish/clear/prep7et,1,plane82blc4,,,10,10lsel,s,tan1,ylesize,all,,,10,-1/5lsel,s,tan1,xlesize,all,,,9,-1/8lsel,allmshape,0mshkey,1amesh,all!EX3.6 合并线和连接线以进行映射网格划分finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,3mshape,0mshkey,1lccat,1,2lcomb,4,5amesh,all!EX3.7 合并线或连接线的网格划分设置finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,,10lesize,6,,,8lesize,1,,,4lesize,2,,,3lesize,4,,,5lesize,5,,,2mshape,0mshkey,1lcomb,1,2!lccat,1,2!如采用lcatt,1,2则采用7个划分数lcomb,4,5amesh,all!EX3.7B 简化映射网格划分finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,,6mshape,0amap,1,2,5,3,4!EX3.8A 过渡四边形映射网格finish/clear/prep7k,1k,2,10,-1k,3,8,6k,4,1,3a,1,2,3,4lesize,1,,,8lesize,3,,,3lesize,4,,,7lesize,2,,,2mshape,0,2dmshkey,1amesh,all!EX3.8B 过渡四边形映射网格finish/clear/prep7et,1,plane42k,1k,2,10,-1k,3,8,6k,4,1,3a,1,2,3,4lesize,1,,,11lesize,3,,,3lesize,4,,,2lesize,2,,,2mshape,0,2dmshkey,1amesh,all!EX3.9A 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,all,,,3lesize,5,,,15allselmshape,0,3dvmesh,all!EX3.9B 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,all,,,4lesize,7,,,12mshape,0,3dmshkey,1vmesh,all!EX3.9C 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lsel,s,,,1,8,7lesize,all,,,6lsel,s,,,2,7,5lesize,all,,,3lsel,s,,,4,5lesize,all,,,7lsel,s,,,3,6,3lesize,all,,,2lsel,s,,,9,12,1lesize,all,,,5allselmshape,0,3dmshkey,1vmesh,all!EX3.9D 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,2,,,2lesize,4,,,2lesize,5,,,2lesize,1,,,3lesize,3,,,3lesize,7,,,6lesize,6,,,7lesize,11,,,7lesize,9,,,3lesize,10,,,3lesize,12,,,3allselmshape,0,3dmshkey,1vmesh,all!EX3.10 带两孔的长方体的扫掠网格划分finish/clear/prep7a=10r=2et,1,mesh200,6et,2,solid45blc4,,,2*a,a,acyl4,a/2,a/2,r,,,,awprota,,90cyl4,1.5*a,a/2,r,,,,-avsbv,1,2vsbv,4,3wprota,,,90wpoff,,,a/2vsbw,allwpoff,,,a/2vsbw,allwpoff,,,a/2vsbw,allwpcsyswpoff,,a/2,a/2vsbw,allwprota,,90vsbw,allwpcsysesize,1amap,105,15,16,26,63amap,107,16,13,60,26amap,108,13,14,28,60amap,103,15,14,28,63lesize,94,,,4lesize,79,,,5vsel,s,loc,x,0,avsweep,allasel,s,loc,y,aasel,r,loc,x,a,2*alccat,2,45lccat,57,71lccat,68,78lccat,65,72mshape,0,2dmshkey,1amesh,allvsel,s,loc,x,a,2*avsweep,allallsel/view,1,1,2,3Eplot!EX3.11 拉伸类命令生成体单元网格finish/clear/prep7et,1,82et,2,95blc4,,,4,4blc4,6,,4,4blc4,12,,4,4blc4,18,,4,4esize,1amesh,allesize,,8vrotat,1,,,,,,1,4,90vext,2,,,,,10,0.5,0.5voffst,3,10vdrag,4,,,,,,35!EX3.12 圆的网格划分finish/prep7et,1,plane82r0=10cyl4,,,r0cyl4,3*r0,,,,r0 wprota,,90asbw,allwprota,,,90asbw,allwpoff,,,3*r0asbw,allwpcsys,-1asel,s,loc,x,-r0,r0lsla,slesize,all,,,8mshape,0,2dmshkey,1amesh,allallselasel,s,loc,x,2*r0,4*r0 lesize,all,,,8lsel,r,length,,r0lesize,all,,,8,0.1,1 amesh,allallsel!EX3.13 圆环的网格划分finish/clear/prep7et,1,plane82r0=10cyl4,,,r0/3,,r0,90cyl4,2*r0,,r0/10,,r0,90 asel,s,loc,x,-r0,r0lsla,slesize,all,,,8lsel,r,length,,r0*2/3 lesize,all,,,3,,1 mshape,0,2dmshkey,1amesh,allasel,s,loc,x,2*r0,4*r0lesize,5,,,12lesize,7,,,6lsel,s,length,,r0*9/10lesize,all,,,7amesh,all!EX3.14 圆柱面和圆柱体的网格划分!圆柱面finish/clear/prep7r0=10h0=50et,1,shell63cyl4,,,r0adele,1cm,l1cm,linek,50k,51,,,h0l,50,51adrag,l1cm,,,,,,5lsel,s,loc,z,0lesize,all,,,6lsel,s,length,,h0lesize,all,,,8mshape,0,2dmshkey,1amesh,all!圆柱体finish/clear/prep7r0=10h0=50et,1,solid95cyl4,,,r0,,,,h0wprota,,90vsbw,allwprota,,,90vsbw,allmshape,0,3dmshkey,1lsel,s,loc,z,0lsel,s,length,,h0lesize,all,,,8vmesh,all!EX3.15 圆锥的网格划分!扫掠网格划分finish/clear/prep7cone,10,,,15,,90et,1,200,7et,2,95lsel,s,,,5,6lesize,all,,,12,0.5lsel,alllesize,3,,,6mshape,0,2dmshkey,1amesh,3vsweep,1,3,4vsymm,x,allvsymm,y,allvglue,all!六面体映射网格划分finish/clear/prep7cone,10,,,15,,90et,2,95lesize,all,,,8mshape,0,3dmshkey,1vmesh,allvsymm,x,allvsymm,y,allvglue,all!EX3.16 正多边形面的通用网格划分finish/clear/prep7ns=11rpr4,ns,,,10kp0=100+nsk,kp0*if,mod(ns,2),eq,0,then*do,i,1,ns/2l,kp0,2*i-1*enddo*else*do,i,1,nsl,kp0,i*enddo*endiflsel,s,,,ns+1,2*nscm,l1cm,lineallselasbl,1,l1cm*if,mod(ns,2),eq,0,then cmsel,s,l1cmlesize,all,,,8lsel,invelesize,all,,,4*elselesize,all,,,8*endifallselmshape,0mshkey,1amesh,all!EX3.17 球体及球面网格划分finish/clear/prep7r0=10sphere,,r0,,90vsbw,allvdele,2,,,1esize,2et,1,solid95mshape,0,3dmshkey,1vsymm,x,all vsymm,y,all vsymm,z,allvglue,allfinish/clear/prep7r0=10sphere,,r0,,90 vsbw,allvdele,2,,,1vdele,allasel,s,loc,x,0asel,a,loc,y,0asel,a,loc,z,0adele,all,,,1allselesize,2et,1,shell63 mshape,0,2d mshkey,1amesh,allarsym,x,allarsym,y,allarsym,z,all nummrg,all!EX3.18 孔板网格划分finish/clear/prep7a0=300b0=800r0=15blc4,,,a0,b0cyl4,a0/4,b0/8,r0 agen,2,2,,,a0/2 agen,2,2,3,1,,b0/8 agen,2,2,5,1,,b0*5/8 asel,s,,,2,9,1cm,a2cm,areaallselwprota,,-90*do,i,1,5wpoff,,,b0/16asbw,all*enddowpoff,,,b0*5/16*do,i,1,5wpoff,,,b0/16asbw,all*enddowprota,,,90*do,i,1,3wpoff,,,a0/4asbw,all*enddowpcsys,-1numcmp,alllsel,s,radius,,r0lesize,all,,,8lsel,invelesize,all,,,4lsel,allet,1,82mshape,0,2dmshkey,1asel,u,loc,y,b0/16,b0*5/16 asel,u,loc,y,b0*11/16,b0*15/16 lsla,slsel,r,tan1,xlesize,all,50,,,,1amesh,allallselamap,21,30,31,54,62 amap,22,31,32,54,64 amap,32,29,30,62,76 amap,33,29,32,64,76amap,42,34,35,76,82 amap,43,35,36,76,84 amap,41,33,34,53,82 amap,44,33,36,53,84amap,19,22,23,49,64 amap,20,23,24,49,60 amap,9,21,22,64,74 amap,31,21,24,60,74amap,38,26,27,74,84 amap,39,27,28,74,80 amap,37,25,26,50,84 amap,40,25,28,50,80amap,17,14,15,43,59 amap,18,15,16,43,63 amap,7,13,14,59,73 amap,8,13,16,63,73amap,34,18,19,73,79 amap,30,17,18,79,44 amap,35,19,20,73,83 amap,36,17,20,44,83amap,15,6,7,39,63amap,16,7,8,39,58amap,4,6,5,63,71amap,5,5,8,58,71amap,27,10,11,71,83 amap,28,11,12,71,78 amap,26,9,10,83,40 amap,29,9,12,40,78!EX3.19 角支架的网格划分finish/clear/prep7blc4,,,150,50blc4,100,,50,-50cyl4,,25,25cyl4,125,-50,25aadd,allnumcmp,allcyl4,,25,10cyl4,125,-50,10 asel,s,,,2,3cm,a1cm,area asel,allasba,1,a1cm lcomb,1,6lfillt,1,2,20 asbl,4,6adele,1,,,1lfillt,3,4,20al,18,19,20 aadd,all numcmp,all!---------------------- wprota,,90 wpoff,,,-25 asbw,all wpoff,,,75 asbw,all wprota,,,90 asbw,all wpoff,,,125 asbw,all wpcsys,-1 wpoff,25 wprota,,,90 asbw,all kwpave,18 asbw,all kwpave,3 wprota,,90 asbw,all kwpave,21 wprota,,-45 asbw,8 wpcsys,-1!-------------------- et,1,plane82 mshape,0,2d mshkey,1 esize,6lesize,33,,,6 lesize,37,,,6 lesize,42,,,6amesh,3,5,2amap,6,9,10,4,23amap,7,9,12,1,23amesh,11,13,2amesh,1,9,8amesh,2,4,2amap,16,3,18,26,28amap,12,14,15,5,28amap,15,13,14,6,28lcomb,19,22lcomb,27,46amesh,10,14,4!========================= !EX3.20 大板小孔的网格划分A finish/clear/prep7a0=100blc4,,,a0,a0cyl4,,,a0/100asba,1,2csys,1k,50,a0/20k,51,a0/20,90l,50,51asbl,all,1wprota,,90,45asbw,allet,1,plane82mshape,0,2dmshkey,1lsel,s,loc,x,0,a0/20lesize,all,,,8lsel,alllesize,7,,,12,0.1lesize,8,,,12,0.1lesize,9,,,12,10amesh,allcsys,0arsym,x,allarsym,y,allnummrg,all!============================= !EX3.20 大板小孔的网格划分Bfinish/clear/prep7a0=100blc4,,,a0,a0cyl4,,,a0/100asba,1,2csys,1k,50,a0/20k,51,a0/20,90l,50,51asbl,all,1et,1,plane82mshape,0,2dmshkey,1lesize,5,,,8lesize,1,,,10lesize,4,,,8,5lesize,6,,,8,5amesh,1lesize,7,,,20,0.1lesize,8,,,20,0.1amap,2,50,51,2,4csys,0arsym,x,allarsym,y,allnummrg,all!EX3.21 Z形面网格划分及局部细分finish/clear/prep7blc4,,,15,10blc4,10,6,14,12aadd,allwprota,,-90wpoff,,,6asbw,allwpoff,,,4wprota,,,90wpoff,,,10asbw,allwpoff,,,5asbw,allwpcsys,-1et,1,82esize,2mshape,0,2dmshkey,1amesh,allkrefine,9,10,1,1,,off !不光滑处理（节点位置不变）aclear,allamesh,allkrefine,9,10,1,1,,smooth !光滑处理，节点位置可改变，但既有单元不删除aclear,allamesh,allkrefine,9,10,1,1 !光滑处理，可删除既有单元重细分，节点位置改变!EX3.22 分布小孔环形面的网格划分finish/clear/prep7r1=15r2=25r3=3n=8cyl4,,,r1,,r2,180/ncyl4,0.5*(r1+r2),,r3asba,1,2kl,2larc,6,8,4,0.5*(r1+r2)asbl,all,4lesize,4,,,4lesize,10,,,4lesize,6,,,8lesize,3,,,6lesize,8,,,4lesize,5,,,8lesize,9,,,4lesize,7,,,4et,1,plane82mshape,0,2damap,2,6,7,4,8amap,1,5,6,8,1!lrefine,5,6,1,1,smootharsym,y,allcsys,1agen,n,all,,,,360/nnummrg,all!EX3.23 同心圆轴体的网格划分finish/clear/prep7cylind,50,0,100,150cylind,40,0,200,250cylind,30,0,100,250cylind,20,0,50,250cylind,10,0,0,250vptn,allwprota,,,90vsbw,allwprota,,90vsbw,allet,1,95csys,1lsel,s,length,,50lesize,all,,,5lsel,invelsel,u,radius,,10,50lesize,all,,,4n1=6n2=6n3=8n4=10n5=12lsel,r,loc,x,0,10lesize,all,,,n1,,1lsel,s,radius,,10lesize,all,,,n1lsel,s,radius,,20lsel,s,radius,,30lesize,all,,,n3lsel,s,radius,,40lesize,all,,,n4lsel,s,radius,,50lesize,all,,,n5allselmshape,0,3dmshkey,1vmesh,all!EX3.24A 圆柱与长方体组合模型的网格划分!采用六面体映射网格划分finish/clear/prep7a0=30h1=15h2=25r0=7blc4,,,a0/2,a0/2,h1cyl4,,,r0,,,90,h1+h2vptn,allaccat,4,6et,1,95esize,3mshape,0,3dmshkey,1vmesh,allasel,s,accaadele,alllsel,s,lccaldele,allallselvsymm,x,allvsymm,y,allnummrg,all!EX3.24B 圆柱与长方体组合模型的网格划分!采用过渡六面体映射网格划分finish/cleara0=30h1=15h2=25r0=7blc4,,,a0/2,a0/2,h1cyl4,,,r0,,,90,h1+h2vptn,alllsel,s,radius,,r0lsel,a,length,,r0lesize,all,,,6lsel,s,loc,z,0lsel,a,loc,z,h1lesize,all,,,7asel,s,loc,x,a0/2asel,a,loc,y,a0/2accat,allallselet,1,95esize,3mshape,0,3dmshkey,1vmesh,allasel,s,accaadele,alllsel,s,lccaldele,allallselvsymm,x,allvsymm,y,allnummrg,all!EX3.25-任意圈数圆柱形螺旋弹簧的扫略网格划分FINISH/CLEAR/PREP7!1.定义弹簧参数（同上）D=4C=8DZ=C*DT=DZ/6N=4.7*IF,T,LT,D,THEN*ENDIFTKPD=90!2.创建全部螺旋线CSYS,1TDEG=N*360TDEG1=MOD(TDEG,TKPD) N0=(TDEG-TDEG1)/TKPD+1 *AFUN,DEG*DO,I,1,N0CTA=(I-1)*TKPDZ=T/360*CTAK,I,DZ/2,CTA,Z*ENDDO*IF,TDEG1,LT,1.0E-2,THEN *ELSEN0=N0+1CTA=CTA+TDEG1Z=T/360*CTAK,N0,DZ/2,CTA,Z*ENDIF*DO,I,1,N0-1L,I,I+1*ENDDOCM,L1,LINE!3.在螺旋线端部创建簧丝截面CSYS,0WPOFF,DZ/2 WPROTA,,90CYL4,,,D/2wprota,,90asbw,allwprota,,,90asbw,allcm,a1cm,areaVDRAG,a1cm,,,,,,L1et,1,mesh200,7et,2,solid95cmsel,s,a1cmlsla,slesize,all,,,6amesh,allallselesize,d/2vsweep,all!EX3.26 具坑缺陷圆柱的网格划分finish/clear/prep7l0=50r0=10r1=5r2=1h0=5et,1,solid45cyl4,,,r0,,,,l0wpave,0,r0,l0/2wprota,,90con4,,,r1,r2,h0vsbv,1,2wprota,,90vsbw,allwpoff,,,r1+3vsbw,allwpoff,,,-2*(r1+3)vsbw,allwprota,,,90vsbw,allwpoff,,,r1+1vsbw,allwpoff,,,-2*(r1+1)vsbw,allwpcsys,-1wprota,,90vsbw,allwpoff,,,-r1vsbw,allwpoff,,,2*r1vsbw,allwpcsys,-1numcmp,allksel,s,loc,z,l0/2ksel,r,loc,x,0ksel,u,loc,y,0*get,kp1,kp,0,num,min kp2=kpnext(kp1) ksel,s,loc,z,l0/2+r2*get,kp3,kp,0,num,min allselk,1000,,-r0,l0/2+5l,kp3,1000asel,none*get,l1,line,0,num,max arotat,l1,,,,,,kp1,kp2 cm,a1cm,areaallselvsba,all,a1cm!vsel,s,loc,z,17,33!vsel,r,loc,x,-6,6!aslv,s!aplotaccat,128,61accat,160,138 accat,224,202 accat,183,35accat,31,134accat,151,164 accat,215,231 accat,82,195accat,130,5accat,162,145 accat,229,209 accat,192,89accat,12,123accat,117,155 accat,109,219 accat,14,178mshape,0,3d mshkey,1esize,1vmesh,all!EX3.27长方体开圆柱槽的网格划分finish/clear/prep7blc5,,,20,10,30wpoff,,5,5cyl4,,,8,,,,10wpoff,,,10cyl4,,,4,,,,12vsel,s,,,2,3cm,v1cm,voluallselwpcsys,-1vsbv,1,v1cmvoffst,15,5voffst,16,3voffst,17,15vptn,allwprota,,,90vsbw,allwpcsys,-1wpoff,,,5vsbw,allwpoff,,,10vsbw,allwpoff,,,12vsbw,allet,1,mesh200,7esize,2mshape,0,2dmshkey,1asel,s,loc,z,30lccat,6,74lccat,8,68amesh,allasel,s,loc,z,0lccat,3,73lccat,1,67amesh,allallselet,2,solid95mshape,0,3dmshkey,1vsweep,all!ex4.1累加节点自由度约束值finish/clear/prep7et,1,beam3k,1k,2,10l,1,2esize,1lmesh,alld,1,ux,1e-2,,,,uyd,1,rotzd,2,alldlistnsel,s,d,u,0,0.1dofsel,s,uydscale,2.5dlistdofsel,alldcum,addd,1,ux,2e-2dlistdcum,ignod,2,uy,1.0dlist!ex4.2 对线施加约束并转换finish/clear/prep7et,1,95blc4,,,10,10,10dl,7,,ux,0.1dl,5,,alldl,11,6,symmdl,10,6,asymdl,6,,symmdllistesize,2vmesh,alldtrandlist!ex4.3 施加集中荷载与节点坐标系finish/clear/prep7et,1,beam4k,1k,2,5k,3,10l,1,2l,2,3local,12,0,,,,90nrotat,alldk,1,allfk,2,fy,-1000esize,1lmesh,allnrotat,alllplotfk,3,fy,1000f,6,fx,-1000sbctraneplot!ex4.4a 3d单元sf加载示例finish/clear/prep7et,1,95blc4,,,10,10,20esize,,4vmesh,allasel,s,loc,y,10sf,all,pres,1000asel,s,loc,z,20nsla,s,1sf,all,pres,1000!ex4.4b 2d单元sf加载示例finish/clear/prep7et,1,82blc4,,,100,200blc4,30,60,40,80asba,1,2wprota,,-90wpoff,,,60asbw,allwpoff,,,80asbw,allwprota,,,90wpoff,,,30asbw,allwpoff,,,40asbw,allwpcsys,-1esize,5amesh,all/psf,pres,norm,2sf,all,pres,100sfdele,all,presnsel,s,loc,x,0sf,all,pres,100nsel,s,loc,x,15,20esln,s,1nsel,r,loc,x,15sf,all,pres,110nsel,s,loc,x,40,60nsel,r,loc,y,10,30esln,s,1sf,all,pres,100lsel,s,loc,x,100nsll,s,1esln,snsel,s,loc,x,95sf,all,pres,-100eplot!ex4.5 节点号及其荷载函数finish/clear/prep7et,1,45blc4,,,10,10,20esize,5vmesh,all*dim,mypres,,100*do,i,1,100mypres(i)=i*10.0*enddosffun,pres,mypres(1)nsel,s,loc,y,10sf,all,pres,10sflist*do,i,1,100mypres(i)=i*50.0*enddonsel,s,loc,z,20sf,all,pres,0allselsflist!ex4.6a 2d平面单元plane82 finish/clear/prep7et,1,82blc4,,,10,50esize,2amesh,allsfe,1,4,pres,,100,50sfe,6,4,pres,,200,40nsel,s,loc,x,10esln,ssfe,all,2,pres,,-100/psf,pres,norm,2,0,1eplot!ex4.6b 3d体单元solid95 finish/clear/prep7et,1,95blc4,,,10,10,30esize,5vmesh,all/psf,pres,norm,2,0,1sfe,22,3,pres,,100sfe,22,6,pres,,100,110,120,130 sfelist!ex4.6c 3d壳单元shell63 finish/clear/prep7et,1,63wprota,,90blc4,,,10,10esize,5amesh,all/psf,pres,norm,2,0,1sfe,1,1,pres,,100sfe,1,3,pres,,100sfe,1,6,pres,,100sfe,4,1,pres,,100,110,120,130!ex4.7 在梁单元上施加荷载finish/clear/prep7et,1,beam3k,1k,2,10l,1,2esize,,10lmesh,all/pnum,elem,1sfbeam,3,1,pres,50,100 sfbeam,5,1,pres,100 sfbeam,7,1,pres,50,100,,,0.2,0.1sfbeam,9,1,pres,100,,,,0.4,-1sfbeam,3,2,pres,50,100,,,0.2,0.1!ex4.8a 在线上施加面荷载finish/clear/prep7et,1,82blc4,,,10,30esize,5amesh,all/psf,pres,norm,2sfl,4,pres,10,60sfl,2,pres,60sftraneplot!ex4.8b 在线上施加面荷载finish/clear/prep7et,1,shell63wprota,,90blc4,,,10,30esize,5amesh,allsfl,3,pres,100sftraneplot!ex4.9 利用荷载梯度在直角坐标系下的施加方法finish/clear/prep7et,1,82blc4,,,10,60esize,2amesh,all/psf,pres,norm,2sfgrad,pres,,y,0,-5nsel,s,loc,x,0nsel,r,loc,y,0,40sf,all,pres,600sfgrad,pres,,y,30,-20nsel,s,loc,x,10sf,all,pres,0allseleplot!ex4.10利用荷载梯度在柱坐标系下的施加方法finish/clear/prep7csys,1k,1,10,-90k,2,10,90k,3,10,90,30l,1,2l,2,3adrag,1,,,,,,2ldele,2,,,1numcmp,allet,1,63esize,2amesh,allsfgrad,pres,1,y,-90,1sfa,all,2,pres,400sftran!可考察下述命令及其施加方法!①!sfgrad,pres,1,y,270,1sfa,all,2,pres,400sftran!②!local,12,1cscir,12,1sfgrad,pres,12,y,270,1sfa,all,2,pres,400sftran!ex4.11 ldread命令的示例finish/clear/filname,ldtest1/prep7et,1,plane82mp,ex,1,2.1e5mp,prxy,1,0.3blc4,,,10,40esize,2amesh,alld,all,allsfl,3,pres,100/solusolve/post1prrsolfinish/filename,ldtest2/solulsclear,allnsel,s,loc,y,0d,all,allallselldread,reac,,,,,ldtest1,rst solve/post1plnsol,s,y!ex4.12 初应力荷载finish/clear/filname,colu1/prep7et,1,plane82mp,ex,1,2e5mp,nuxy,1,0.3blc4,,,1,10esize,2amesh,allnsel,s,loc,y,0d,all,uyd,1,uxnsel,s,loc,y,10sf,all,pres,-10。

Ansys经典命令流钢管混凝土拱桥算例

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。

Ansys命令流大全（整理）

Ansys命令流大全（整理）1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9此命令用已知的一组关键点点（P1~P9）来定义面（Area），最少使用三个点才能围成面，同时产生转围绕些面的线。

点要依次序输入，输入的顺序会决定面的法线方向。

如果超过四个点，则这些点必须在同一个平面上。

Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>Through KPs2、*ABBR，Abbr,String－－定义一个缩略语．Abbr:用来表示字符串＂String＂的缩略语，长度不超过８个字符．String：将由＂Abbr＂表示的字符串，长度不超过６０个字符．3、ABBRES，Lab，Fname，Ext－从一个编码文件中读出缩略语．Lab：指定读操作的标题，NEW：用这些读出的缩略语重新取代当前的缩略语（默认）CHANGE：将读出的缩略语添加到当前缩略语阵列，并替代现存同名的缩略语．Ext：如果＂Fname＂是空的，则缺省的扩展命是＂ABBR＂．4、ABBSA V，Lab，Fname，Ext－将当前的缩略语写入一个文本文件里Lab：指定写操作的标题，若为ALL，表示将所有的缩略语都写入文件（默认）5、add, ir, ia,ib,ic,name,--,--,facta, factb, factc将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic：变量号name: 变量的名称6、Adele,na1,na2,ninc,kswp ！kswp=0时只删除掉面积本身，＝1时低单元点一并删除。

7、Adrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 ！面积的建立，沿某组线段路径，拉伸而成。

8、Afillt,na1,na2,rad ！建立圆角面积，在两相交平面间产生曲面，rad为半径。

ansys混凝土命令流

vsel,s,loc,x,-0.002,0 !选择左边建筑结构胶
vatt,2,2,1
esize,0.01
mshape,0,3d$mshkey,1
vmesh,all
vsel,s,loc,x,0.102,0.112 !选择右边钢板
vatt,3,3,1
wpoff,0.002,, !工作平面平移
blc4,,,0.01,0.05 ,0.16 !钢板
wpcsys,-1
wpoff,-0.002,0.025,0.05 !工作平面平移
blc4,,,0.002,0.05 ,0.15 !建筑结构胶
blc4,,,-0.01,0.05 ,0.16 !钢板
fini
/cle
/FILNAME, hntd, 1
/TITLE, tiegang
/prep7
!----------定义单元及材料等---------------------
Et,1,solid45 !定义混凝土单元类型
Et,2,solid45 !定义建筑结构胶单元类型
mp,prxy,1,0.2 !定义混凝土泊松比
mp,prxy,2,0.317 !定义建筑结构胶泊松比
mp,prxy,3,0.3 !定义钢板泊松比
mp,dens,1,2550 !定义混凝土密度
mp,dens,2,1700 !定义建筑结构胶密度
sfa,all,,pres,-1200e3/(2*6)
save
!-----------求解---------------------------
/solu
/nerr,10000,100000
solve

钢筋混凝土梁ansys分析附命令流

钢筋混凝⼟梁ansys分析附命令流钢筋混凝⼟⾮线性分析2015⼤作业上海交通⼤学陈明1、参数选择梁的截⾯宽度为200mm，上部配置2Φ8受压筋，混凝⼟的净保护层厚度为25 mm（从纵向钢筋外边缘算起），箍筋两端区采⽤8@100的双肢箍，中间区取8@200 双肢箍1）梁的截⾯⾼度选300mm；2）两加载间的距离选1000mm；3）混凝⼟选C30；4）纵向受拉钢筋配筋选218；2、描述选⽤的有限元模型及单元的特点采⽤ansys软件进⾏模拟计算，钢筋混凝⼟模型采⽤分离式模型，不考虑钢筋与混凝⼟之间的相对滑移。

混凝⼟采⽤solid65单元模拟，solid65⽤于模拟三维有钢筋或⽆钢筋的混凝⼟模型。

该单元能够计算拉裂和压碎。

在混凝⼟应⽤中，该单元的实体功能可以⽤于建⽴混凝⼟模型，同时，还可⽤加筋功能建⽴钢筋混凝⼟模型。

另外，该单元还可以应⽤于加强复合物和地质材料。

该单元由⼋个节点定义，每个节点有三个⾃由度:节点坐标系的x，y，z⽅向的平动。

⾄多可以定义三种不同规格的钢筋。

钢筋单元采⽤link180单元模拟，link180是⼀个适⽤于各类⼯程应⽤的三维杆单元。

根据具体情况，该单元可以被看作桁架单元、索单元、链杆单元或弹簧单元等等。

本单元是⼀个轴向拉伸⼀压缩单元，每个节点有三个⾃由度:节点坐标系的x，y，z⽅向的平动。

本单元是⼀种顶端铰接结构，不考虑单元弯曲。

本单元具有塑性、蠕变、旋转、⼤变形和⼤应变功能。

缺省时，当考虑⼤变形时任何分析中LINK180单元都包括应⼒刚化选项。

3、描述选⽤的混凝⼟与钢筋粘结滑移本构关系的具体形式、参数等。

钢筋的应⼒应变关系曲线考虑到极限塑性应变最⼤值为0.01,钢筋本构模型采⽤多线性模型kinh，初始弹性模量为Es=200000Mpa，强化系数为0.001。

混凝⼟的应⼒应变关系曲线混凝⼟选⽤各向同性的miso模型，当计⼊下降端时，程序报错，所以只取了前⾯的上升段，⽤5段折线模拟混凝⼟应⼒应变曲线。

(完整版)ANSYS命令流总结(全)

ANSYS结构分析单元功能与特性/可以组成一一些命令，一般是一种总体命令（session),三十也有特殊，比如是处理/POST1! 是注释说明符号，，与其他软件的说明是一样的，ansys不作为命令读取，* 此符号一般是APDL的标识符，也就是ansys的参数化语言，如*do ,,,*enddo等等NSEL的意思是node select，即选择节点。

s就是select，选择。

DIM是定义数组的意思。

array 数组。

MP命令用来定义材料参数。

K是建立关键点命令。

K,关键点编号,x坐标,y坐标，z坐标。

K, NPT, X, Y, Z是定义关键点，K是命令，NPT是关键点编号，XYZ是坐标。

NUMMRG, keypoint 用这个命令，要保证关键点的位置完全一样，只是关键点号不一样的才行。

这个命令对于重复的线面都可以用。

这个很简单，压缩关键。

Ngen 复制节点e,节点号码：这个命令式通过节点来形成单元NUMCMP,ALL：压缩所有编号，这样你所有的线都会按次序重新编号~你要是需要固定的线固定的标号NSUBST,100,500,50：通过指定子步数来设置载荷步的子步LNSRCH线性搜索是求解非线性代数方程组的一种技巧，此法会在一段区间内，以一定的步长逐步搜索根，相比常用的牛顿迭代法所要耗费的计算量大得多，但它可以避免在一些情况下牛顿迭代法出现的跳跃现象。

LNSRCH激活线性搜索PRED 激活自由度求解预测NEQIT指定一个荷载步中的最大子步数AUTOTS 自动求解控制打开自动时间步长.KBC -指定阶段状或者用跳板装载里面一个负荷步骤。

SPLINE:P1，P2，P3，P4，P5，P6，XV1，YV1，ZV1，XV6，YV6，ZV6（生成分段样条曲线）*DIM，Par，Type，IMAX，JMAX，KMAX，Var1，Var2，Var3（定义载荷数组的名称）【注】Par: 数组名Type：array 数组，如同fortran,下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8个字符）tableIMAX，JMAX，KMAX各维的最大下标号Var1，Var2，Var3 各维变量名，缺省为row,column,plane(当type为table时)/config是设置ansys配置参数的命令格式为/CONFIG, Lab, V ALUELab为参数名称value为参数值例如：/config，MXEL，10000的意思是最大单元数为10000杆单元:LINK1、8、10、11、180梁单元：BEAM3、4、23、24，44，54，188，189管单元:PIPE16，17，18，20，59，602D实体元:PLANE2，25，42，82，83，145，146，182，1833D实体元:SOLID45，46，64，65，72，73，92，95，147，148，185，186，187，191壳单元:SHELL28，41，43，51，61，63，91，93，99，143，150，181，208，209弹簧单元:COMBIN7，14，37，39，40质量单元:MASS21接触单元:CONTAC12，52，TARGE169，170，CONTA171，172，173，174，175，178矩阵单元:MATRIX27，50表面效应元:SURF153，154粘弹实体元:VISCO88，89，106，107，108，超弹实体元:HYPER56，58，74，84，86，158耦合场单元:SOLID5，PLANE13，FLUID29，30，38，SOLID62，FLUID79，FLUID80，81， SOLID98，FLUID129，INFIN110，111，FLUID116，130 界面单元:INTER192，193，194，195 显式动力分析单元:LINK160，BEAM161，PLANE162，SHELL163，SOLID164，COMBI16杆单元(Large deflection)，F-大应变(Large strain)或有限应变(Finite strain)，B-单元生死(Birth and dead),G-应力刚化(Stress stiffness)或几何刚度(Geometric stiffening)，A-自适应下降(Adaptive descent )等。

ansys混凝土浇筑命令流

ansys混凝土浇筑命令流/CLEAR/FILNAME,CT,1/TITLE,CT01*SET,T1,1 !第一层开始浇筑时间*SET,T2,8 !第二层开始浇筑时间*SET,TZZ,30 !终止时间*SET,cdiam,0.028 !水管径*SET,carea,3.14*cdiam*cdiam/4 !水管面积*SET,massflow,28800 !管内流体质量流率（kg/d）*SET,KX1,234 !混凝土导热系数*SET,C1,0.963 !混凝土比热*SET,DENSE1,2400 !砼密度*SET,KX2,50.1 !水导热系数*SET,C2,4.2 !水比热*SET,DENSE2,1000 !水密度*SET,DLD,588 !上层表面放热系数$*SET,DLB,994 !底层放热系数$*SET,DLC,1983 !侧面放热系数*SET,Q_SUM,326*377 !单位水化热总量*SET,TEMP_q,19 !大气温度*SET,TEMP_c,22 !砼温度*SET,TEMP_w,17 !水流温度/PREP7ET,1,70 !定义单元solid70et,2,116,1,2 !单元fluid116MP,KXX,1,KX1 !单元70 材料MP,C,1,C1MP,DENS,1,DENSE1MP,KXX,2,KX2 !单元116材料MP,C,2,C2MP,DENS,2,DENSE2MP,HF,2,3574r,1,cdiam,carea,1 !单元116实常数!建立实体模型K,1,0,0,0$K,2,12.19,0,0$K,3,12.19,0,8.25$K,4,0,0,8.25K,5,12.65,0,8.25K,6,15.15,0,0$A,1,2,3,4A,2,3,5,6$ASEL,ALL$VEXT,ALL,,,,5 !通过偏移量由面生成体ALLSEL,ALL !选择所有的实体VGLUE,ALL !体粘接NUMMRG,ALL !对实体进行合并$NUMCMP,ALL !压缩所定义项的编号TYPE, 1 !划分模型 !指定单元类型属性指示器FLST,5,4,4,ORDE,4FITEM,5,1FITEM,5,3FITEM,5,8FITEM,5,10CM,_Y,LINE !由所选择的实体生成一个元件LSEL, , , ,P51X !p51x说包含你用鼠标选取的所有线CM,_Y1,LINECMSEL,,_Y !选择一组元件或部件子集LESIZE,_Y1, , ,23, , , , ,1 !对选择的线设置网格单元大小23为分段数1修改所有已选择的线FLST,5,6,4,ORDE,6 !为某一命令选择操作实体时产生的FITEM,5,2FITEM,5,4FITEM,5,6FITEM,5,11FITEM,5,17CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINE CMSEL,,_YLESIZE,_Y1, , ,33, , , , ,1 FLST,5,6,4,ORDE,4 FITEM,5,12FITEM,5,-15FITEM,5,19FITEM,5,-20CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINE CMSEL,,_YLESIZE,_Y1, , ,10, , , , ,1 FLST,5,4,4,ORDE,4 FITEM,5,5FITEM,5,7FITEM,5,16FITEM,5,18CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINE CMSEL,,_YLESIZE,_Y1, , ,5, , , , ,1 FLST,5,2,6,ORDE,2 !6为实体FITEM,5,-2CM,_Y,VOLUVSEL, , , ,P51XCM,_Y1,VOLUCHKMSH,'VOLU' !用来检查网格质量的，可查看你画的网格是否有警告单元，或err单元CMSEL,S,_YVSWEEP,_Y1 !利用与体邻近的面单元采用扫琼方式对体进行网格划分CMDELE,_Y !删除一个现存的元件或部件CMDELE,_Y1CMDELE,_Y2NSEL,S,LOC,Y,1 !选择第一层水管节点FLST,5,242,1,ORDE,212FITEM,5,1916FITEM,5,2174FITEM,5,2184FITEM,5,2200FITEM,5,2209FITEM,5,2218FITEM,5,2227FITEM,5,2236FITEM,5,2245FITEM,5,2254FITEM,5,2263FITEM,5,2272FITEM,5,2281FITEM,5,2299 FITEM,5,2308 FITEM,5,2317 FITEM,5,2326 FITEM,5,2335 FITEM,5,2344 FITEM,5,2353 FITEM,5,2362 FITEM,5,2371 FITEM,5,2380 FITEM,5,2389 FITEM,5,2414 FITEM,5,2429 FITEM,5,2438 FITEM,5,2447 FITEM,5,2456 FITEM,5,2465 FITEM,5,2474 FITEM,5,2483 FITEM,5,2492 FITEM,5,2501 FITEM,5,2510 FITEM,5,2519 FITEM,5,2528 FITEM,5,2537 FITEM,5,2546 FITEM,5,2555 FITEM,5,2564 FITEM,5,2573FITEM,5,2591 FITEM,5,2600 FITEM,5,2609 FITEM,5,2618 FITEM,5,2643 FITEM,5,2658 FITEM,5,2946 FITEM,5,2955 FITEM,5,2964 FITEM,5,2973 FITEM,5,2981 FITEM,5,3021 FITEM,5,-3052 FITEM,5,3311 FITEM,5,3326 FITEM,5,3335 FITEM,5,3344 FITEM,5,3353 FITEM,5,3360 FITEM,5,3369 FITEM,5,3378 FITEM,5,3387 FITEM,5,3396 FITEM,5,3405 FITEM,5,3414 FITEM,5,3423 FITEM,5,3432 FITEM,5,3441 FITEM,5,3450FITEM,5,3468 FITEM,5,3477 FITEM,5,3486 FITEM,5,3495 FITEM,5,3504 FITEM,5,3513 FITEM,5,3522 FITEM,5,3531 FITEM,5,3540 FITEM,5,3549 FITEM,5,3558 FITEM,5,3567 FITEM,5,3576 FITEM,5,3585 FITEM,5,3594 FITEM,5,3603 FITEM,5,3612 FITEM,5,3621 FITEM,5,3630 FITEM,5,3639 FITEM,5,3648 FITEM,5,3657 FITEM,5,3666 FITEM,5,3675 FITEM,5,3684 FITEM,5,3693 FITEM,5,3702 FITEM,5,3711 FITEM,5,3720FITEM,5,3738 FITEM,5,3747 FITEM,5,3756 FITEM,5,3765 FITEM,5,3774 FITEM,5,3783 FITEM,5,3792 FITEM,5,3801 FITEM,5,3810 FITEM,5,3819 FITEM,5,3828 FITEM,5,3837 FITEM,5,3846 FITEM,5,3855 FITEM,5,3864 FITEM,5,3873 FITEM,5,3882 FITEM,5,3891 FITEM,5,3900 FITEM,5,3909 FITEM,5,3918 FITEM,5,3927 FITEM,5,3936 FITEM,5,3945 FITEM,5,3954 FITEM,5,3963 FITEM,5,3972 FITEM,5,3981 FITEM,5,3990FITEM,5,4008 FITEM,5,4017 FITEM,5,4026 FITEM,5,4035 FITEM,5,4044 FITEM,5,4053 FITEM,5,4062 FITEM,5,4071 FITEM,5,4080 FITEM,5,4089 FITEM,5,4098 FITEM,5,4107 FITEM,5,4116 FITEM,5,4125 FITEM,5,4134 FITEM,5,4143 FITEM,5,4152 FITEM,5,4161 FITEM,5,4170 FITEM,5,4179 FITEM,5,4188 FITEM,5,4197 FITEM,5,4206 FITEM,5,4215 FITEM,5,4224 FITEM,5,4233 FITEM,5,4242 FITEM,5,4251 FITEM,5,4260FITEM,5,4278 FITEM,5,4287 FITEM,5,4296 FITEM,5,4305 FITEM,5,4314 FITEM,5,4323 FITEM,5,4332 FITEM,5,4341 FITEM,5,4350 FITEM,5,4359 FITEM,5,4368 FITEM,5,4377 FITEM,5,4386 FITEM,5,4395 FITEM,5,4404 FITEM,5,4413 FITEM,5,4422 FITEM,5,4431 FITEM,5,4440 FITEM,5,4449 FITEM,5,4458 FITEM,5,4467 FITEM,5,4476 FITEM,5,4485 FITEM,5,4494 FITEM,5,4503 FITEM,5,4518 FITEM,5,4806 FITEM,5,5094FITEM,5,5670 FITEM,5,5958 FITEM,5,6246 FITEM,5,6534 FITEM,5,6822 FITEM,5,7110 FITEM,5,7398 FITEM,5,7686 FITEM,5,7974 FITEM,5,8262FITEM,5,8550 FITEM,5,8838 FITEM,5,9126 FITEM,5,9414 FITEM,5,9702 FITEM,5,9990 FITEM,5,10278 FITEM,5,10566 NSEL,U, , ,P51X NPLOT !显示节点编号FLST,5,362,1,ORDE,333 FITEM,5,1885 FITEM,5,-1915 FITEM,5,4527 FITEM,5,4536 FITEM,5,4545 FITEM,5,4554 FITEM,5,4563FITEM,5,4581 FITEM,5,4590 FITEM,5,4599 FITEM,5,4608 FITEM,5,4617 FITEM,5,4626 FITEM,5,4635 FITEM,5,4644 FITEM,5,4653 FITEM,5,4662 FITEM,5,4671 FITEM,5,4680 FITEM,5,4689 FITEM,5,4698 FITEM,5,4707 FITEM,5,4716 FITEM,5,4725 FITEM,5,4734 FITEM,5,4743 FITEM,5,4752 FITEM,5,4761 FITEM,5,4770 FITEM,5,4779 FITEM,5,4788 FITEM,5,5112 FITEM,5,5121 FITEM,5,5130 FITEM,5,5139 FITEM,5,5148FITEM,5,5166 FITEM,5,5175 FITEM,5,5184 FITEM,5,5193 FITEM,5,5202 FITEM,5,5211 FITEM,5,5220 FITEM,5,5229 FITEM,5,5238 FITEM,5,5247 FITEM,5,5256 FITEM,5,5265 FITEM,5,5274 FITEM,5,5283 FITEM,5,5292 FITEM,5,5301 FITEM,5,5310 FITEM,5,5319 FITEM,5,5328 FITEM,5,5337 FITEM,5,5346 FITEM,5,5355 FITEM,5,5364 FITEM,5,5373 FITEM,5,5679 FITEM,5,5688 FITEM,5,5697 FITEM,5,5706 FITEM,5,5715FITEM,5,5733 FITEM,5,5742 FITEM,5,5751 FITEM,5,5760 FITEM,5,5769 FITEM,5,5778 FITEM,5,5787 FITEM,5,5796 FITEM,5,5805 FITEM,5,5814 FITEM,5,5823 FITEM,5,5832 FITEM,5,5841 FITEM,5,5850 FITEM,5,5859 FITEM,5,5868 FITEM,5,5877 FITEM,5,5886 FITEM,5,5895 FITEM,5,5904 FITEM,5,5913 FITEM,5,5922 FITEM,5,5931 FITEM,5,5940 FITEM,5,6264 FITEM,5,6273 FITEM,5,6282 FITEM,5,6291 FITEM,5,6300FITEM,5,6318 FITEM,5,6327 FITEM,5,6336 FITEM,5,6345 FITEM,5,6354 FITEM,5,6363 FITEM,5,6372 FITEM,5,6381 FITEM,5,6390 FITEM,5,6399 FITEM,5,6408 FITEM,5,6417 FITEM,5,6426 FITEM,5,6435 FITEM,5,6444 FITEM,5,6453 FITEM,5,6462 FITEM,5,6471 FITEM,5,6480 FITEM,5,6489 FITEM,5,6498 FITEM,5,6507 FITEM,5,6516 FITEM,5,6525 FITEM,5,6831 FITEM,5,6840 FITEM,5,6849 FITEM,5,6858 FITEM,5,6867FITEM,5,6885 FITEM,5,6894 FITEM,5,6903 FITEM,5,6912 FITEM,5,6921 FITEM,5,6930 FITEM,5,6939 FITEM,5,6948 FITEM,5,6957 FITEM,5,6966 FITEM,5,6975 FITEM,5,6984 FITEM,5,6993 FITEM,5,7002 FITEM,5,7011 FITEM,5,7020 FITEM,5,7029 FITEM,5,7038 FITEM,5,7047 FITEM,5,7056 FITEM,5,7065 FITEM,5,7074 FITEM,5,7083 FITEM,5,7092 FITEM,5,7407 FITEM,5,7416 FITEM,5,7425 FITEM,5,7434 FITEM,5,7443FITEM,5,7461 FITEM,5,7470 FITEM,5,7479 FITEM,5,7488 FITEM,5,7497 FITEM,5,7506 FITEM,5,7515 FITEM,5,7524 FITEM,5,7533 FITEM,5,7542 FITEM,5,7551 FITEM,5,7560 FITEM,5,7569 FITEM,5,7578 FITEM,5,7587 FITEM,5,7596 FITEM,5,7605 FITEM,5,7614 FITEM,5,7623 FITEM,5,7632 FITEM,5,7641 FITEM,5,7650 FITEM,5,7659 FITEM,5,7668 FITEM,5,7677 FITEM,5,7983 FITEM,5,7992 FITEM,5,8001 FITEM,5,8010FITEM,5,8028 FITEM,5,8037 FITEM,5,8046 FITEM,5,8055 FITEM,5,8064 FITEM,5,8073 FITEM,5,8082 FITEM,5,8091 FITEM,5,8100 FITEM,5,8109 FITEM,5, 8118 FITEM,5,8127 FITEM,5,8136 FITEM,5,8145 FITEM,5,8154 FITEM,5,8163 FITEM,5,8172 FITEM,5,8181 FITEM,5,8190 FITEM,5,8199 FITEM,5,8208 FITEM,5,8217 FITEM,5,8226 FITEM,5,8235 FITEM,5,8244 FITEM,5,8568 FITEM,5,8577 FITEM,5,8586FITEM,5,8604 FITEM,5,8613 FITEM,5,8622 FITEM,5,8631 FITEM,5,8640 FITEM,5,8649 FITEM,5,8658 FITEM,5,8667 FITEM,5,8676 FITEM,5,8685 FITEM,5,8694 FITEM,5,8703 FITEM,5,8712 FITEM,5,8721 FITEM,5,8730 FITEM,5,8739 FITEM,5,8748 FITEM,5,8757 FITEM,5,8766 FITEM,5,8775 FITEM,5,8784 FITEM,5,8793 FITEM,5,8802 FITEM,5,8811 FITEM,5,8820 FITEM,5,8829 FITEM,5,9135 FITEM,5,9144 FITEM,5,9153FITEM,5,9171 FITEM,5,9180 FITEM,5,9189 FITEM,5,9198 FITEM,5,9207 FITEM,5,9216 FITEM,5,9225 FITEM,5,9234 FITEM,5,9243 FITEM,5,9252 FITEM,5,9261 FITEM,5,9270 FITEM,5,9279 FITEM,5,9288 FITEM,5,9297 FITEM,5,9306 FITEM,5,9315 FITEM,5,9324 FITEM,5,9333 FITEM,5,9342 FITEM,5,9351 FITEM,5,9360 FITEM,5,9369 FITEM,5,9378 FITEM,5,9387 FITEM,5,9396 FITEM,5,9720 FITEM,5,9729 FITEM,5,9738FITEM,5,9756 FITEM,5,9765 FITEM,5,9774 FITEM,5,9783 FITEM,5,9792 FITEM,5,9801 FITEM,5,9810 FITEM,5,9819 FITEM,5,9828 FITEM,5,9837 FITEM,5,9846 FITEM,5,9855 FITEM,5,9864 FITEM,5,9873 FITEM,5,9882 FITEM,5,9891 FITEM,5,9900 FITEM,5,9909 FITEM,5,9918 FITEM,5,9927 FITEM,5,9936 FITEM,5,9945 FITEM,5,9954 FITEM,5,9963 FITEM,5,9972 FITEM,5,9981 FITEM,5,10287 FITEM,5,10296 FITEM,5,10305FITEM,5,10323 FITEM,5,10332 FITEM,5,10341 FITEM,5,10350 FITEM,5,10359 FITEM,5,10368 FITEM,5,10377 FITEM,5,10386 FITEM,5,10395 FITEM,5,10404 FITEM,5,10413 FITEM,5,10422 FITEM,5,10431 FITEM,5,10440 FITEM,5,10449 FITEM,5,10458 FITEM,5,10467 FITEM,5,10476 FITEM,5,10485 FITEM,5,10494 FITEM,5,10503 FITEM,5,10512 FITEM,5,10521 FITEM,5,10530 FITEM,5,10539 FITEM,5,10548 NSEL,U, , ,P51X NPLOTFLST,5,20,1,ORDE,20FITEM,5,4815FITEM,5,5085FITEM,5,5391FITEM,5,5661FITEM,5,5967FITEM,5,6237FITEM,5,6543FITEM,5,6813FITEM,5,7389FITEM,5,7965FITEM,5,8271FITEM,5,8541FITEM,5,8847FITEM,5,9117FITEM,5,9423FITEM,5,9693FITEM,5,9999FITEM,5,10269FITEM,5,10845NSEL,U, , ,P51XNSEL,A,NODE,,10573,10574,1NSEL,A,NODE,,7693,7694,1NSEL,A,NODE,,7117,7118,1NSEL,A,NODE,,2665,2666,1CM,NODE1,NODENGEN,2,100000,ALL !创建第一层水管单元ngen复制节点TYPE,2MAT,2REAL,1 !单元实常数属性指示器E,110573,110574,10573,10574 !通过节点相连生成一个单元E,110574,110575,10574,10575*do,i,1,29,1e,110575+(i-1)*9,110584+(i-1)*9,10575+(i-1)*9,10584+(i-1)*9*enddoe,110836,110557,10836,10557e,110557,110260,10557,10260*do,i,1,28,1e,110260-(i-1)*9,110251-(i-1)*9,10260-(i-1)*9,10251-(i-1)*9 *enddoe,110008,109711,10008,9711e,109711,109432,9711,9432*do,i,1,28,1e,109432+(i-1)*9,109441+(i-1)*9,9432+(i-1)*9,9441+(i-1)*9 *enddoe,109684,109405,9684,9405e,109405,109108,9405,9108*do,i,1,28,1e,109108-(i-1)*9,109099-(i-1)*9,9108-(i-1)*9,9099-(i-1)*9*enddoe,108856,108559,8856,8559e,108559,108280,8559,8280*do,i,1,28,1e,108280+(i-1)*9,108289+(i-1)*9,8280+(i-1)*9,8289+(i-1)*9 *enddoe,108532,108253,8532,8253e,108253,107956,8253,7956*do,i,1,29,1e,107956-(i-1)*9,107947-(i-1)*9,7956-(i-1)*9,7947-(i-1)*9 *enddoE,107695,107694,7695,7694E,107694,107693,7694,7693E,102665,102666,2665,2666E,102666,102667,2666,2667*do,i,1,29,1e,102667+(i-1)*9,102676+(i-1)*9,2667+(i-1)*9,2676+(i-1)*9 *enddoe,102928,104797,2928,4797e,104797,105076,4797,5076*do,i,1,28,1e,105076-(i-1)*9,105067-(i-1)*9,5076-(i-1)*9,5067-(i-1)*9 *enddoe,104824,105103,4824,5103e,105103,105400,5103,5400*do,i,1,28,1e,105400+(i-1)*9,105409+(i-1)*9,5400+(i-1)*9,5409+(i-1)*9 *enddoe,105652,105949,5652,5949e,105949,106228,5949,6228*do,i,1,28,1e,106228-(i-1)*9,106219-(i-1)*9,6228-(i-1)*9,6219-(i-1)*9 *enddoe,105976,106255,5976,6255e,106255,106552,6255,6552*do,i,1,28,1e,106552+(i-1)*9,106561+(i-1)*9,6552+(i-1)*9,6561+(i-1)*9 *enddoe,106804,107101,6804,7101e,107101,107380,7101,7380*do,i,1,29,1e,107380-(i-1)*9,107371-(i-1)*9,7380-(i-1)*9,7371-(i-1)*9 *enddoE,107119,107118,7119,7118E,107118,107117,7118,7117 ! 第一层水管单元创建完毕! 选择第二层水管节点allsel,allNSEL,S,LOC,Y,3NSEL,U,LOC,X,0NSEL,U,LOC,z,0NSEL,U,LOC,Z,7.9,8.25NSEL,R,LOC,X,0.2,12.20FLST,5,351,1,ORDE,351FITEM,5,2933FITEM,5,4523FITEM,5,4532FITEM,5,4541FITEM,5,4550FITEM,5,4559FITEM,5,4568FITEM,5,4577FITEM,5,4586FITEM,5,4595FITEM,5,4604FITEM,5,4613FITEM,5,4622FITEM,5,4640 FITEM,5,4649 FITEM,5,4658 FITEM,5,4667 FITEM,5,4676 FITEM,5,4685 FITEM,5,4694 FITEM,5,4703 FITEM,5,4712 FITEM,5,4721 FITEM,5,4730 FITEM,5,4739 FITEM,5,4748 FITEM,5,4757 FITEM,5,4766 FITEM,5,4775 FITEM,5,4784 FITEM,5,4811 FITEM,5,5081 FITEM,5,5108 FITEM,5,5117 FITEM,5,5126 FITEM,5,5135 FITEM,5,5144 FITEM,5,5153 FITEM,5,5162 FITEM,5,5171 FITEM,5,5180 FITEM,5,5189FITEM,5,5207 FITEM,5,5216 FITEM,5,5225 FITEM,5,5234 FITEM,5,5243 FITEM,5,5252 FITEM,5,5261 FITEM,5,5270 FITEM,5,5279 FITEM,5,5288 FITEM,5,5297 FITEM,5,5306 FITEM,5,5315 FITEM,5,5324 FITEM,5,5333 FITEM,5,5342 FITEM,5,5351 FITEM,5,5360 FITEM,5,5369 FITEM,5,5387 FITEM,5,5657 FITEM,5,5675 FITEM,5,5684 FITEM,5,5693 FITEM,5,5702 FITEM,5,5711 FITEM,5,5720 FITEM,5,5729 FITEM,5,5738FITEM,5,5756 FITEM,5,5765 FITEM,5,5774 FITEM,5,5783 FITEM,5,5792 FITEM,5,5801 FITEM,5,5810 FITEM,5,5819 FITEM,5,5828 FITEM,5,5837 FITEM,5,5846 FITEM,5,5855 FITEM,5,5864 FITEM,5,5873 FITEM,5,5882 FITEM,5,5891 FITEM,5,5900 FITEM,5,5909 FITEM,5,5918 FITEM,5,5927 FITEM,5,5936 FITEM,5,5963 FITEM,5,6233 FITEM,5,6260 FITEM,5,6269 FITEM,5,6278 FITEM,5,6287 FITEM,5,6296FITEM,5,6314 FITEM,5,6323 FITEM,5,6332 FITEM,5,6341 FITEM,5,6350 FITEM,5,6359 FITEM,5,6368 FITEM,5,6377 FITEM,5,6386 FITEM,5,6395 FITEM,5,6404 FITEM,5,6413 FITEM,5,6422 FITEM,5,6431 FITEM,5,6440 FITEM,5,6449 FITEM,5,6458 FITEM,5,6467 FITEM,5,6476 FITEM,5,6485 FITEM,5,6494 FITEM,5,6503 FITEM,5,6512 FITEM,5,6521 FITEM,5,6539 FITEM,5,6809 FITEM,5,6827 FITEM,5,6836 FITEM,5,6845FITEM,5,6863 FITEM,5,6872 FITEM,5,6881 FITEM,5,6890 FITEM,5,6899 FITEM,5,6908 FITEM,5,6917 FITEM,5,6926 FITEM,5,6935 FITEM,5,6944 FITEM,5,6953 FITEM,5,6962 FITEM,5,6971 FITEM,5,6980 FITEM,5,6989 FITEM,5,6998 FITEM,5,7007 FITEM,5,7016 FITEM,5,7025 FITEM,5,7034 FITEM,5,7043 FITEM,5,7052 FITEM,5,7061 FITEM,5,7070 FITEM,5,7079 FITEM,5,7088 FITEM,5,7385 FITEM,5,7403 FITEM,5,7412FITEM,5,7430 FITEM,5,7439 FITEM,5,7448 FITEM,5,7457 FITEM,5,7466 FITEM,5,7475 FITEM,5,7484 FITEM,5,7493 FITEM,5,7502 FITEM,5,7511 FITEM,5,7520 FITEM,5,7529 FITEM,5,7538 FITEM,5,7547 FITEM,5,7556 FITEM,5,7565 FITEM,5,7574 FITEM,5,7583 FITEM,5,7592 FITEM,5,7601 FITEM,5,7610 FITEM,5,7619 FITEM,5,7628 FITEM,5,7637 FITEM,5,7646 FITEM,5,7655 FITEM,5,7664 FITEM,5,7673 FITEM,5,7961FITEM,5,7988 FITEM,5,7997 FITEM,5,8006 FITEM,5,8015 FITEM,5,8024 FITEM,5,8033 FITEM,5,8042 FITEM,5,8051 FITEM,5,8060 FITEM,5,8069 FITEM,5,8078 FITEM,5,8087 FITEM,5,8096 FITEM,5,8105 FITEM,5,8114 FITEM,5,8123 FITEM,5,8132 FITEM,5,8141 FITEM,5,8150 FITEM,5,8159 FITEM,5,8168 FITEM,5,8177 FITEM,5,8186 FITEM,5,8195 FITEM,5,8204 FITEM,5,8213 FITEM,5,8222 FITEM,5,8231 FITEM,5,8240FITEM,5,8537 FITEM,5,8564 FITEM,5,8573 FITEM,5,8582 FITEM,5,8591 FITEM,5,8600 FITEM,5,8609 FITEM,5,8618 FITEM,5,8627 FITEM,5,8636 FITEM,5,8645 FITEM,5,8654 FITEM,5,8663 FITEM,5,8672 FITEM,5,8681 FITEM,5,8690 FITEM,5,8699 FITEM,5,8708 FITEM,5,8717 FITEM,5,8726 FITEM,5,8735 FITEM,5,8744 FITEM,5,8753 FITEM,5,8762 FITEM,5,8771 FITEM,5,8780 FITEM,5,8789 FITEM,5,8798 FITEM,5,8807FITEM,5,8825 FITEM,5,8843 FITEM,5,9113 FITEM,5,9131 FITEM,5,9140 FITEM,5,9149 FITEM,5,9158 FITEM,5,9167 FITEM,5,9176 FITEM,5,9185 FITEM,5,9194 FITEM,5,9203 FITEM,5,9212 FITEM,5,9221 FITEM,5,9230 FITEM,5,9239 FITEM,5,9248 FITEM,5,9257 FITEM,5,9266 FITEM,5,9275 FITEM,5,9284 FITEM,5,9293 FITEM,5,9302 FITEM,5,9311 FITEM,5,9320 FITEM,5,9329 FITEM,5,9338 FITEM,5,9347 FITEM,5,9356FITEM,5,9374 FITEM,5,9383 FITEM,5,9392 FITEM,5,9419 FITEM,5,9689 FITEM,5,9716 FITEM,5,9725 FITEM,5,9734 FITEM,5,9743 FITEM,5,9752 FITEM,5,9761 FITEM,5,9770 FITEM,5,9779 FITEM,5,9788 FITEM,5,9797 FITEM,5,9806 FITEM,5,9815 FITEM,5,9824 FITEM,5,9833 FITEM,5,9842 FITEM,5,9851 FITEM,5,9860 FITEM,5,9869 FITEM,5,9878 FITEM,5,9887 FITEM,5,9896 FITEM,5,9905 FITEM,5,9914FITEM,5,9932 FITEM,5,9941 FITEM,5,9950 FITEM,5,9959 FITEM,5,9968 FITEM,5,9977 FITEM,5,9995 FITEM,5,10265 FITEM,5,10283 FITEM,5,10292 FITEM,5,10301 FITEM,5,10310 FITEM,5,10319 FITEM,5,10328 FITEM,5,10337 FITEM,5,10346 FITEM,5,10355 FITEM,5,10364 FITEM,5,10373 FITEM,5,10382 FITEM,5,10391 FITEM,5,10400 FITEM,5,10409 FITEM,5,10418 FITEM,5,10427 FITEM,5,10436 FITEM,5,10445 FITEM,5,10454 FITEM,5,10463FITEM,5,10481 FITEM,5,10490 FITEM,5,10499 FITEM,5,10508 FITEM,5,10517 FITEM,5,10526 FITEM,5,10535 FITEM,5,10544 FITEM,5,10841 NSEL,U, , ,P51X NSEL,A,NODE,,1506 NSEL,A,NODE,,10568 NSEL,A,NODE,,10569 NSEL,A,NODE,,10570NSEL,A,NODE,,7690 NSEL,A,NODE,,7689 NSEL,A,NODE,,7688 NSEL,A,NODE,,1186 NSEL,A,NODE,,7114 NSEL,A,NODE,,7113 NSEL,A,NODE,,7112 NSEL,A,NODE,,1122 NSEL,A,NODE,,2662 NSEL,A,NODE,,2661 NSEL,A,NODE,,2660 NSEL,A,NODE,,1567 cm,node2,nodengen,2,200000,allcm,node2_,node!创建第二层水管单元TYPE,2MAT,2REAL,1E,201506,210568,1506,10568E,210568,210569,10568,10569E,210569,210570,210569,210570E,210570,210571,10570,10571*do,i,1,29,1$e,200000+10571+(i-1)*9,200000+10580+(i-1)*9,10571+(i-1)*9,10580+(i-1)*9$*enddo$e,210832,210553,10832,10553$e,210553,210256,10553,10256*do,i,1,28,1$e,200000+10256-(i-1)*9,200000+10247-(i-1)*9,10256-(i-1)*9,10247-(i-1)*9$*enddo$e,200000+10004,200000+9707,10004,9707$e,200000+9707,200000+9428,9707,9428*do,i,1,28,1$e,200000+9428+(i-1)*9,200000+9437+(i-1)*9,9428+(i-1)*9,9437+(i-1)*9$*enddo$e,200000+9680,200000+9401,9680,9401$e,200000+9401,200000+9104,9401,9104*do,i,1,28,1$e,200000+9104-(i-1)*9,200000+9095-(i-1)*9,9104-(i-1)*9,9095-(i-1)*9$*enddo$e,200000+8852,200000+8555,8852,8555$e,200000+8555,200000+8276,8555,8276*do,i,1,28,1$e,200000+8276+(i-1)*9,200000+8285+(i-1)*9,8276+(i-1)*9,8285+(i-1)*9$*enddo$e,200000+8528,200000+8249,8528,8249$e,200000+8249,200000+7952,8249,7952*do,i,1,29,1$e,200000+7952-(i-1)*9,200000+7943-(i-1)*9,7952-(i-1)*9,7943-(i-1)*9$*enddoE,207691,207690,7691,7690E,207690,207689,7690,7689E,207689,207688,7689,7688E,207688,201186,7688,1186E,201567,202660,1567,2660E,202660,202661,2660,2661E,202661,202662,2661,2662E,202662,202663,2662,2663*do,i,1,29,1$e,200000+2663+(i-1)*9,200000+2672+(i-1)*9,2663+(i-1)*9,2672+(i-1)*9$*enddo$e,200000+2924,200000+4793,2924,4793$e,200000+4793,200000+5072,4793,5072$*do,i,1,28,1$e,200000+5072-(i-1)*9,200000+5063-(i-1)*9,5072-(i-1)*9,5063-(i-1)*9$*enddo$e,200000+4820,200000+5099,4820,5099$e,200000+5099,200000+5396,5099,5396$*do,i,1,28,1$e,200000+5396+(i-1)*9,200000+5405+(i-1)*9,5396+(i-1)*9,5405+(i-1)*9$*enddo$e,200000+5648,200000+5945,5648,5945$e,200000+5945,200000+6224,5945,6224$*do,i,1,28,1$e,200000+6224-(i-1)*9,200000+6215-(i-1)*9,6224-(i-1)*9,6215-(i-1)*9$*enddo$e,200000+5972,200000+6251,5972,6251$e,200000+6251,200000+6548,6251,6548$*do,i,1,28,1$e,200000+6548+(i-1)*9,200000+6557+(i-1)*9,6548+(i-1)*9,6557+(i-1)*9$*enddo$e,200000+6800,200000+7097,6800,7097$e,200000+7097,200000+7376,7097,7376$*do,i,1,29,1$e,200000+7376-(i-1)*9,200000+7367-(i-1)*9,7376-(i-1)*9,7367-(i-1)*9$*enddoE,207115,207114,7115,7114E,207114,207113,7114,7113E,207113,207112,7113,7112E,207112,201122,7112,1122!第二层水管创建完毕!选择第三层水管节点ALLSEL,ALLNSEL,S,LOC,Y,4NSEL,U,LOC,Z,0NSEL,U,LOC,Z,7.9,8.25NSEL,R,LOC,X,0.2,12.2FLST,5,351,1,ORDE,351FITEM,5,2931FITEM,5,4521FITEM,5,4530FITEM,5,4539FITEM,5,4548FITEM,5,4557FITEM,5,4566FITEM,5,4575FITEM,5,4584FITEM,5,4593FITEM,5,4602FITEM,5,4611FITEM,5,4620FITEM,5,4629FITEM,5,4647 FITEM,5,4656 FITEM,5,4665 FITEM,5,4674 FITEM,5,4683 FITEM,5,4692 FITEM,5,4701 FITEM,5,4710 FITEM,5,4719 FITEM,5,4728 FITEM,5,4737 FITEM,5,4746 FITEM,5,4755 FITEM,5,4764 FITEM,5,4773 FITEM,5,4782 FITEM,5,4809 FITEM,5,5079 FITEM,5,5106 FITEM,5,5115 FITEM,5,5124 FITEM,5,5133 FITEM,5,5142 FITEM,5,5151 FITEM,5,5160 FITEM,5,5169 FITEM,5,5178 FITEM,5,5187 FITEM,5,5196FITEM,5,5214 FITEM,5,5223 FITEM,5,5232 FITEM,5,5241 FITEM,5,5250 FITEM,5,5259 FITEM,5,5268 FITEM,5,5277 FITEM,5,5286 FITEM,5,5295 FITEM,5,5304 FITEM,5,5313 FITEM,5,5322 FITEM,5,5331 FITEM,5,5340 FITEM,5,5349 FITEM,5,5358 FITEM,5,5367 FITEM,5,5385 FITEM,5,5655 FITEM,5,5673 FITEM,5,5682 FITEM,5,5691 FITEM,5,5700 FITEM,5,5709 FITEM,5,5718 FITEM,5,5727 FITEM,5,5736 FITEM,5,5745FITEM,5,5763 FITEM,5,5772 FITEM,5,5781 FITEM,5,5790 FITEM,5,5799 FITEM,5,5808 FITEM,5,5817 FITEM,5,5826 FITEM,5,5835 FITEM,5,5844 FITEM,5,5853 FITEM,5,5862 FITEM,5,5871 FITEM,5,5880 FITEM,5,5889 FITEM,5,5898 FITEM,5,5907 FITEM,5,5916 FITEM,5,5925 FITEM,5,5934 FITEM,5,5961 FITEM,5,6231 FITEM,5,6258 FITEM,5,6267 FITEM,5,6276 FITEM,5,6285 FITEM,5,6294 FITEM,5,6303 FITEM,5,6312FITEM,5,6330 FITEM,5,6339 FITEM,5,6348 FITEM,5,6357 FITEM,5,6366 FITEM,5,6375 FITEM,5,6384 FITEM,5,6393 FITEM,5,6402 FITEM,5,6411 FITEM,5,6420 FITEM,5,6429 FITEM,5,6438 FITEM,5,6447 FITEM,5,6456 FITEM,5,6465 FITEM,5,6474 FITEM,5,6483 FITEM,5,6492 FITEM,5,6501 FITEM,5,6510 FITEM,5,6519 FITEM,5,6537 FITEM,5,6807 FITEM,5,6825 FITEM,5,6834 FITEM,5,6843 FITEM,5,6852 FITEM,5,6861FITEM,5,6879 FITEM,5,6888 FITEM,5,6897 FITEM,5,6906 FITEM,5,6915 FITEM,5,6924 FITEM,5,6933 FITEM,5,6942 FITEM,5,6951 FITEM,5,6960 FITEM,5,6969 FITEM,5,6978 FITEM,5,6987 FITEM,5,6996 FITEM,5,7005 FITEM,5,7014 FITEM,5,7023 FITEM,5,7032 FITEM,5,7041 FITEM,5,7050 FITEM,5,7059 FITEM,5,7068 FITEM,5,7077 FITEM,5,7086 FITEM,5,7383 FITEM,5,7401 FITEM,5,7410 FITEM,5,7419FITEM,5,7437 FITEM,5,7446 FITEM,5,7455 FITEM,5,7464 FITEM,5,7473 FITEM,5,7482 FITEM,5,7491 FITEM,5,7500 FITEM,5,7509 FITEM,5,7518 FITEM,5,7527 FITEM,5,7536 FITEM,5,7545 FITEM,5,7554 FITEM,5,7563 FITEM,5,7572 FITEM,5,7581 FITEM,5,7590 FITEM,5,7599 FITEM,5,7608 FITEM,5,7617 FITEM,5,7626 FITEM,5,7635 FITEM,5,7644 FITEM,5,7653 FITEM,5,7662 FITEM,5,7671 FITEM,5,7959 FITEM,5,7977FITEM,5,7995 FITEM,5,8004 FITEM,5,8013 FITEM,5,8022 FITEM,5,8031 FITEM,5,8040 FITEM,5,8049 FITEM,5,8058 FITEM,5,8067 FITEM,5,8076 FITEM,5,8085 FITEM,5,8094 FITEM,5,8103 FITEM,5,8112 FITEM,5,8121 FITEM,5,8130 FITEM,5,8139 FITEM,5,8148 FITEM,5,8157 FITEM,5,8166 FITEM,5,8175 FITEM,5,8184 FITEM,5,8193 FITEM,5,8202 FITEM,5,8211 FITEM,5,8220 FITEM,5,8229 FITEM,5,8238 FITEM,5,8265。

1.梁柱焊接ANSYS命令流(梁焊接过程)

et,1,solid70 mptemp,1,20,100,200,300,400,500 mptemp,7,600,700,800,900,1000,1500 mpdata,ex,1,1,2.06e11,2.06e11,1.85e11,1.65e11,1.45e11,1.24e11 mpdata,ex,1,7,6.39e10,2.68e10,1.86e10,1.39e10,9.27e9,1.32e9 tb,bkin,1,6 tbtemp,20,1 tbdata,1,235e6,0.206e11 tbtemp,200,2 tbdata,1,235e6,0.185e11 tbtemp,500,3 tbdata,1,183e6,0.124e11 tbtemp,800,4 tbdata,1,26e6,0.186e10 tbtemp,1000,5 tbdata,1,9.4e6,0.927e9 tbtemp,1500,6 tbdata,1,1.5e6,0.132e9 mp,dens,1,7850 mp,alpx,1,1.4e-5 600
l=1.2 wweld=15e-3 lnum=60 !******************************************焊接参数
I=200 U=30 v=1.2/600 pi=3.1415926 eta=0.7 k=1.2e4 them0=25 !**********************************************定义高斯表面移动热源参数
/solu /VIEW,1,1,1,1 /ANG,1 /REP,FAST /DIST,1,0.924021086472,1 /REP,FAST /DIST,1,0.924021086472,1 /REP,FAST /DIST,1,0.924021086472,1 /REP,FAST /DIST,1,0.924021086472,1 /REP,FAST

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

finish$/clear$/config,fsplit,750$/prep7!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义材料!!!!!!!!!下面定义材料 1，Q235A型钢钢材!!!!!!!!!mp,ex,1,2.1e11$mp,prxy,1,0.3$mp,dens,1,7850tb,bkin,1tbdata,,268.85e6,0!!下面定义材料 9，Q235A型钢钢材!!!!!!!!!mp,ex,3,2.1e11$mp,prxy,3,0.3$mp,dens,3,7850tb,bkin,3$tbdata,,235e6,0!!!下面定义材料 2,混凝土C40!!!!!!!!!fck=26*0.85*0.8e6$ftk=2.3e6mp,ex,2,fck*0.44/0.0005$mp,prxy,2,0.2$mp,dens,2,2600tb,concr,2$tbdata,,0.125,0.7,FTK,-1tb,mkin,2$tbtemp,,strain ! 混凝土屈服准则Tbdata,1,0.0005,0.001,0.0015,0.002,0.0038$tbtemp,0 $Tbdata,1,0.44*fck,0.75*fck,0.94*fck,fck,0.85*fck !混凝土应力应变曲线!!!下面定义材料 4,纵钢筋直径20!!!!!!!!! mp,ex,4,2.1e11$mp,prxy,4,0.3$mp,dens,4,7850 tb,bkin,4tbdata,,365e6,0!!!下面定义材料7,箍筋直径8!!!!!!!!!mp,ex,5,2.1e11$mp,prxy,5,0.3$mp,dens,5,7850 tb,bkin,5tbdata,,243e6,0!!!!定义实常数!!!!!!Pi=acos(-1)$ r,1$ r,2$r,3,0.25*pi*0.020*0.020r,4,0.25*pi*0.008*0.008R,5!!!!定义单元!!!!!!!et,1,solid45$et,2,solid65$et,3,LINK8keyopt,2,7,1$keyopt,2,1,1$keyopt,2,3,2!!!!建模型!!!!!!!!柱子wpstyl,defa$Block,-0.2,0.2,-1.5,1.5,0,0.2cm,zhuzi,volu!工字钢vsel,none$wpstyl,defa$$block,-2,2,-0.125,0.125,0,0.003$bloc k,-2,2,0.125-0.015,0.125,0,0.078$$block,-2,2,-0.125,-0.125+0.015,0,0.078$$wpstyl,defa$wpoff,-0.2$$block,0,0.010,-0.125,0.125,0,0.078 $wpstyl,defa$wpoff,0.2$$block,-0.010,0,-0.125,0.125,0,0.078 $wpstyl,defa$wpoff,,0.125-0.015$block,-0.075,0.075,0,1.5-0.125+0.015,0,0.003$$block,-0.075,-0.075+0.008,0,1.39,0,0.078$$block,0.075-0.008,0.075,0,1.39,0,0.078$$wpstyl,defa$wpoff,,-0.125+0.015$block,-0.075,0.075,-1.39,0,0,0.003$$block,-0.075,-0.075+0.008,-1.39,0,0,0.078$$block,0.075-0.008,0.075,-1.39,0,0,0.078$$wpstyl,defa$wpoff,,-0.125+0.015$block,-0.2,0.2,-0.06-0.015,0,0,0.003$$wpstyl,defa$wpoff,,0.125-0.015$block,-0.2,0.2,0,0.060+0.015,0,0.003$$wpstyl,defa$wpoff,-0.2$$block,0,0.010,0.125,0.125+0.06,0,0 .2$wpstyl,defa$wpoff,0.2$$block,-0.010,0,0.125,0.125+0.06,0,0 .2$wpstyl,defa$wpoff,-0.2$$block,0,0.010,-0.125-0.06,-0.125,0 ,0.2$wpstyl,defa$wpoff,0.2$$block,-0.010,0,-0.125-0.06,-0.125,0 ,0.2$wpstyl,defa$wpoff,,-0.125$wpoff,,,0.2$block,-0.2,0.2,-0.06,0,-0.010,0$$wpstyl,defa$wpoff,,0.125$wpoff,,,0.2$block,-0.2,0.2,0,0.06,-0.010,0$$wpstyl,defa$wpoff,-0.2$$block,0,0.010,-0.11,0.11,0,0.2$wps tyl,defa$wpoff,0.2$$block,-0.010,0,-0.11,0.11,0,0.2 $wpstyl,defa$wpoff,,,0.2$block,-0.2,0.2,-0.11,0.11,-0.010,0$Nummrg,kp$vptn,all$cm,gzg,volu$numcmp,all$allsel,all$vptn,all$$allsel,allVsel,u,,,gzgCm,zhuzi,volu!切分gzgCmsel,s,gzg,volu$Nummrg,kp$!画纵筋$wpstyl,defa$wpoff,0.2$wpoff,,,0.160$vsbw,all$ $wpstyl,defa$wpoff,0.160$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,-0.16$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,0.110$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,-0.11$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,0.2$wpoff,,,0.078$vsbw,all$$numcmp,all!画箍筋$wpstyl,defa$wprota,,90 $wpoff,,,0.04$vsbw,all $wpoff,,,0.070$vsbw,all$wpoff,,,0.015$vsbw,all$wpoff,,,0.030$vsbw,all$wpoff,,,0.045$vsbw,all$wpoff,,,0.05$vsbw,all$*do,i,1,12,1$wpoff,,,0.10$vsbw,all$*enddo$wpstyl,defa$wprota,,-90 $$wpoff,,,0.04$vsbw,all $wpoff,,,0.070$vsbw,all$wpoff,,,0.015$vsbw,all$wpoff,,,0.030$vsbw,all$wpoff,,,0.045$vsbw,all$wpoff,,,0.05$vsbw,all$*do,i,1,12,1$wpoff,,,0.10$vsbw,all$*enddo$$wpstyl,defa$wprota,,,90$wpoff,,,0.075-0.008$vsbw,all $wpoff,,,0.008$vsbw,all$wpstyl,defa$wprota,,,-90$wpoff,,,0.075-0.008$vsbw,all $wpoff,,,0.008$vsbw,all$wpstyl,defa$wpoff,,,0.003$vsbw,all$wpstyl,defa$wpoff,,,0.078$vsbw,all$wpstyl,defa$wprota,,,90$vsbw,all$wpstyl,defa$wprota,,90 $wpoff,,,0.125+0.06$vsbw,all $wpstyl,defa$wprota,,-90 $wpoff,,,0.125+0.06$vsbw,all$numcmp,allcm,gzg,volu!切分柱cmsel,s,zhuzi,volu$!画纵筋$wpstyl,defa$wpoff,0.2$wpoff,,,0.160$vsbw,all$ $wpstyl,defa$wpoff,0.160$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,-0.16$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,0.110$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,-0.11$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,,,0.078$vsbw,all$$numcmp,all!画箍筋$wpstyl,defa$wprota,,90 $wpoff,,,0.04$vsbw,all $wpoff,,,0.070$vsbw,all$wpoff,,,0.015$vsbw,all$wpoff,,,0.030$vsbw,all$wpoff,,,0.045$vsbw,all$wpoff,,,0.05$vsbw,all$*do,i,1,12,1$wpoff,,,0.10$vsbw,all$*enddo$wpstyl,defa$wprota,,-90 $$wpoff,,,0.04$vsbw,all $wpoff,,,0.070$vsbw,all$wpoff,,,0.015$vsbw,all$wpoff,,,0.030$vsbw,all$wpoff,,,0.045$vsbw,all$wpoff,,,0.05$vsbw,all$*do,i,1,12,1$wpoff,,,0.10$vsbw,all$*enddo$$wpstyl,defa$wprota,,,90$wpoff,,,0.075-0.008$vsbw,all$wpoff,,,0.008$vsbw,all$wpstyl,defa$wprota,,,-90$wpoff,,,0.075-0.008$vsbw,all $wpoff,,,0.008$vsbw,all$wpstyl,defa$wpoff,,,0.003$vsbw,all$wpstyl,defa$wprota,,,90$vsbw,all$wpstyl,defa$wpoff,0.2-0.01$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,-0.2+0.01$wprota,,,90$ vsbw,all $wpstyl,defa$wpoff,0.2$wpoff,,,0.2-0.01$vsbw,all$ $wpstyl,defa$wprota,,90 $wpoff,,,0.125+0.06$vsbw,all $wpstyl,defa$wprota,,-90 $wpoff,,,0.125+0.06$vsbw,all$numcmp,allcm,zhuzi,volu!!!!!!!!!!赋钢筋单元属性并划分网格!!!!!!!!!elemsiz=50e-3lsel,s,loc,x,-0.16$lsel,r,loc,z,0.078$cm,waizj1,line$latt,4 ,3,3,0$lesize,all,elemsizlsel,s,loc,x,-0.16$lsel,r,loc,z,0.16$cm,waizj2,line$latt,4,3,3,0$lesize,all,elemsizlsel,s,loc,x,-0.11$lsel,r,loc,z,0.16$cm,waizj3,line $latt,4,3,3,0$lesize,all,elemsizlsel,s,loc,x,0.16$lsel,r,loc,z,0.078$cm,ywaizj1,line$latt,4 ,3,3,0$lesize,all,elemsizlsel,s,loc,x,0.16$lsel,r,loc,z,0.16$cm,ywaizj2,line $latt,4,3,3,0$lesize,all,elemsizlsel,s,loc,x,0.11$lsel,r,loc,z,0.16$cm,ywaizj3,line $latt,4,3,3,0$lesize,all,elemsiz$cmsel,s,waizj1,line$$cmsel,a,waizj2,line$$cmsel,a,waizj3,l ine$$cmsel,a,ywaizj1,line$$cmsel,a,ywaizj2,line$$cmsel,a,ywaizj 3,line$$cm,zj,line!节点外域箍筋lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x,-2.25,-0.161 $lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015 $lsel,u,loc,y,-0.295,0.295 $lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,u,loc,y,1.5$lsel,u,loc,y,-1.5$lsel,r,loc,z,0,0.16$lsel,r,loc,x,-0.16$cm,sgj1,line$lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x,-2.25,-0.161 $lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015 $lsel,u,loc,y,-0.295,0.295$lsel,u,loc,y,1.5$lsel,u,loc,y,-1.5$lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,r,loc,z,0,0.16$lsel,r,loc,x,0.16$cm,sgj2,line$lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x,-2.25,-0.161 $lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015 $lsel,u,loc,y,-0.295,0.295 $lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,u,loc,y,1.5$lsel,u,loc,y,-1.5$lsel,r,loc,z,0.16$cm,sgj3,line$lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x,-2.25,-0.161$lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015$lsel,u,loc,y,-0.295,0.295$lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,u,loc,y,1.5$lsel,u,loc,y,-1.5$Lsel,r,loc,z,0.078$cm,sgj4,line$$cmsel,s,sgj1,line$$cmsel,a,sgj2,line$ cmsel,a,sgj3,line$$cmsel,a,sgj4,line$cm,guj1,line$!节点区域箍筋!节点区域上下四肢箍箍筋lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015 $lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,u,loc,y,-0.095,0.095 $lsel,u,loc,y,0.251,1.5$lsel,u,loc,y,-1.5,-0.251$lsel,u,loc,y, 0.125$lsel,u,loc,y,-0.125$Lsel,r,loc,z,0.16$cm,ssgj1,line$lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x,-2.25,-0.161 $lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015 $lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,u,loc,y,0.251,1.5$lsel,u,loc,y,-1.5,-0.251$lsel,u,loc,y, 0.125$lsel,u,loc,y,-0.125$Lsel,r,loc,z,0.078$cm,ssgj2,line$lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x,-2.25,-0.161 $lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015$lsel,u,loc,y,-0.095,0.095 $lsel,u,loc,y,0.251,1.5$lsel,u,loc,y,-1.5,-0.251$lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,u,loc,y, 0.125$lsel,u,loc,y,-0.125$Lsel,r,loc,x,-0.16$Lsel,r,loc,z,0,0.16$cm,ssgj3,line$lsel,s,tan1,y !对称外定义箍筋$lsel,u,loc,x,-2.25,-0.161$lsel,u,loc,x, 0.161,2.25$lsel,u,loc,y,0.125-0.015$lsel,u,loc,y,-0.125+0.015$lsel,u,loc,y,-0.125-0.06$lsel,u,loc,y,0.125+0.06$lsel,u,loc,y,-0.095,0.095$lsel,u,loc,y,0.251,1.5$lsel,u,loc,y,-1.5,-0.251$lsel,u,loc,y, 0.125$lsel,u,loc,y,-0.125$Lsel,r,loc,x,0.16$Lsel,r,loc,z,0,0.16$cm,ssgj4,line$$cmsel,s,ssgj1,line$$cmsel,a,ssgj2,line$cmsel,a,ssgj3,line$cmsel,a,ssgj4,line$$cm,guj3,line$$cmsel,s,guj1,line$$cmsel,a,guj3,line$$cm,guj,line$latt,5,4 ,3,0$lesize,all,elemsiz$cmsel,s,guj,line$cmsel,a,zj,line$cm,gj,lineCmsel,s,gzg,volu$VSEL,r,LOC,Y,-0.125+0.015,0.125-0.015 $VSEL,r,LOC,z,0,0.003$cm,fb,VOLUCmsel,s,gzg,volu$VSEL,u,LOC,Y,-0.125,0.125$cm,GZhu,VOLU Cmsel,s,gzg,volu$$cmSEL,U,fb,VOLU$VSEL,r,LOC,Y,-0.125+0.015 ,0.125-0.015$cm,mianban,VOLUCmsel,s,gzg,volu$$cmSEL,U,fb,VOLU$$cmSEL,U,gzhu,VOLU$$cmSEL ,U,mianban,VOLU$$cm,yy,VOLU!!!!!!!!!!!!!!!!!!!!!!!! Mesh钢筋!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!cmsel,s,gj,line$lmesh,all $lsel,all$/eshape,1$eplot !钢筋划分完毕!!!!!!!!!!!!!!!!!!!! mesh 钢筋混凝土柱!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Cmsel,s,zhuzi,volu $esize,50e-3$vatt,2,1,2,0$MSHAPE,0,3D$MSHKEY,1$VMESH,allallsel,all!!!!!!!!!!!!!!!!!! Mesh 翼缘!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Cmsel,s,yy,volu$$LSEL,s,length,,1.8LESIZE,all, , ,12vatt,1,1,1,0$esize,50e-3$$MSHAPE,0,3D$MSHKEY,1$VMESH,all$allsel,all!!!!!!!!!!!!!!!!!! Mesh fb!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Cmsel,s,fb,volu$vatt,1,1,1,0$esize,50e-3$$MSHAPE,0,3D$MSHKEY,1$VMESH,all$allsel,all!!!!!!!!!!!!!!!!!! Mesh mianban gzhu!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Cmsel,s,mianban,volu$Cmsel,a,gzhu,volu$vatt,3,5,1,0$esize,50e-3$$MSHAPE,0,3D$MSHKEY,1 $VMESH,all$allsel,all$nummrg,node$numcmp,all!定义远距离加载点$ET,4,TARGE170ET,5,CONTA174R,6keyopt,4,2,1$keyopt,5,2,2$keyopt,5,4,2$keyopt,5,12,5 asel,s,loc,x,2$nsla,s,1$real,6$type,5$esurfnmax=ndinqr(0,14)+1$n,nmax,2,,type,4$tshap,pilo$e,nmaxR,7keyopt,4,2,1$keyopt,5,2,2$keyopt,5,4,2$keyopt,5,12,5 asel,s,loc,x,-2$nsla,s,1$real,7$type,5$esurfnmax=ndinqr(0,14)+1$n,nmax,-2,,type,4$tshap,pilo$e,nmaxR,8keyopt,4,2,1$keyopt,5,2,2$keyopt,5,4,2$keyopt,5,12,5 asel,s,loc,y,1.5$nsla,s,1$real,8$type,5$esurfnmax=ndinqr(0,14)+1$n,nmax,0,1.5,type,4$tshap,pilo$e,nmaxAllsel,all!施加约束和荷载Allsel,allAsel,s,loc,z,0Da,all,symm$D,8665,uy,,,,,uz$D,8666,uy,,,,,uz$Asel,s,loc,y,-1.5Da,all,allAllsel,all/solu$ANTYPE,0$OUTRES,ALL,ALL$$pred,off$neqit,30$lnsrch,on$nlgenom,off$Cnvtol,f,,0.05$kbc,0$eqslv,sparse,1.0E-5,2 $p=0.01$time,1$nsubst,50,,10$Asel,s,loc,y,1.5f,8667,fy,-582000$$allsel,all$OUTRES,ALL,ALL $ALLSEL,ALL $LSWRITE,1TIME,2$D,8667,Ux,19*p$$kbc,0$ NSUBST,100,1000000,1$AUTOTS,o n$OUTRES,ALL,ALL $ALLSEL,ALL$LSWRITE,2$Lssolve,1,2$save$finish/post1SET,LAST$SET,PREVIOUS*GET,RTIME,ACTIVE,0,SET,TIME !获取最终收敛时间Pldisp$plnsol,s,x !显示变形图与X方向应力分布图Plnsol,epel,x !X弹性应变分布图Plnsol,eppl,x !x塑性应变分布图Plnsol,epto,x !X方向总应变分布图Plnsol,nl,srat !显示塑性区域/device,vector,off$plcrack,Cmsel,s,gzgEslv,s/POST26nsel,s,loc,x,0.975NUMVAR,200$NSOL,2,8667,U,x$nsel,s,loc,y,-1.55$RFORCE,3,8667,F,xcmsel,s,gj,line$Esll,s$/output,Z1,csv$prod,4,2,,,,,,1000$Prod,5,3,,,,,,1/1000PRVAR,4,5!!! 钢梁应力的显示/post1$Cmsel,s,gzg$cmsel,a,FUBAN$vsel,r,loc,z,0,0.0045$vsel,r,l oc,y,-0.125+0.006,0.125-0.006$vsel,r,loc,x,-0.175,0.175$eslv,s$nsle,splnsol,s,eqvPlnsol,epel,x !X弹性应变分布图Plnsol,eppl,x !x塑性应变分布图Plesol,epto,1 !X方向总应变分布图/post1$Cmsel,s,gzg$$eslv,s$nsle,s$plnsol,s,eqv/post1$Cmsel,s,zhuzi$eslv,s$nsle,splnsol,s,eqv$Cmsel,s,zhuzi$cmsel,a,FUBAN$vsel,r,loc,z,0,0.0045$vsel,r ,loc,y,-0.125+0.006,0.125-0.006$vsel,r,loc,x,-0.175,0.175$e slv,s$nsle,splesol,s,eqv/post1$Cmsel,s,zhuzi$eslv,s$nsle,splnsol,s,eqv/post1SET,LAST$SET,PREVIOUS*GET,RTIME,ACTIVE,0,SET,TIME !获取最终收敛时间Pldisp$plnsol,s,x !显示变形图与X方向应力分布图Plnsol,epel,x !X弹性应变分布图Plnsol,eppl,x !x塑性应变分布图Plnsol,epto,eqv !X方向总应变分布图Plnsol,nl,srat !显示塑性区域/post1Cmsel,s,gj$esll,sETABLE,rebar,ls,1PLETAB,rebar/post1$wpstyl,defa$$path,23,2,30,120$ppath,1,,0.185,0.125,0$ppath ,2,,0.975,0.125,0pdef,myseqv,s,eqv$pdef,mysy,s,1plpath,myseqv,mysyPLPAGM,myseqv,1,blankPLPAGM,myseqv,1,nodeprpath,myseqv,,xg,yg,zg/device,vector,on$plcrack,Cmsel,s,gzgEslv,snsel,s,loc,x,0.975NUMVAR,200$NSOL,2,9989,U,Y$nsel,s,loc,y,-1.55$RFORCE,3,9989,F,Y/output,Z1,csv$prod,4,2,,,,,,-1000$Prod,5,3,,,,,,-1/1000 PRVAR,4,5/POST26nsel,s,loc,x,0.975NUMVAR,200$NSOL,6,10678,U,Y$nsel,s,loc,y,-1.55$RFORCE,7,10678,F,Y/output,Z2,csv$prod,8,6,,,,,,1000$Prod,9,7,,,,,,1/1000!后处理阶段/post1SET,LAST$SET,PREVIOUS*GET,RTIME,ACTIVE,0,SET,TIME !获取最终收敛时间Pldisp$plnsol,s,x !显示变形图与X方向应力分布图Plnsol,epel,x !X弹性应变分布图Plnsol,eppl,x !x塑性应变分布图Plnsol,epto,eqv !X方向总应变分布图Plnsol,nl,srat !显示塑性区域esel,s,type,,2cmsel,s,gj,line$esll,sETABLE,rebar,ls,1PLETAB,rebar!混凝土裂缝的显示!外混凝土$Cmsel,s,zhuzi,$vsel,r,loc,z,0.0625+0.001,0.175$vsel,r,lo c,y,-0.125+0.006,0.125-0.006$eslv,s$$/dscale,,off$/device,vector,on$plcrack,1,1 !第一条裂纹!内混凝土$Cmsel,s,zhuzi$vsel,r,loc,z,0,0.0625$vsel,r,loc,y,-0.125+ 0.006,0.125-0.006$eslv,s$/dscale,,off$$/device,vector,on$pl crack/post1!混凝土应力的显示$/gformat,defa !恢复缺省格式$Cmsel,s,zhuzi$eslv,s$plnsol,s,1 !主应力$Cmsel,s,zhuzi$eslv,s$plnsol,s,3 !最大压应力$Cmsel,s,zhuzi$eslv,s$plnsol,s,eqv !等效应力$Cmsel,s,gzg$Cmsel,a,fb$eslv,s$plnsol,s,eqv/post1!内混凝土矢量图$Cmsel,s,zhuzi,$vsel,r,loc,z,0,0.0625$vsel,r,loc,y,-0.125 +0.006,0.125-0.006$eslv,s$$/device,vector,off$/dscale,,off$$plnsol,s,eqv $Plvect,u$$Plvect,s$/vscale,,1$plvect,s,,,,vect !调整箭头大小$/gline,,-1$plvect,s,,,,vect ! 取消单元边界$/dscale,,off$$plnsol,s,1 !关闭结构变形$/dscale,defas$/gline,1 !恢复变形比例缺省，显示单元边界$/Device,vector,1 !设置矢量显示模式$/clabel,,5 !设置等高线上的文字间隔$/contour,,18,-16,,500 !设置等高线数目$Plnsol,s,x！外混凝土矢量图$Cmsel,s,zhuzi,$vsel,r,loc,z,0.0625+0.001,0.175$vsel,r,lo c,y,-0.125+0.006,0.125-0.006$eslv,s$/device,vector,off$/dscale,,off$$plesol,s,3$Plvect,u$$Plvect,s$/vscale,,1$plvect,s,,,,vect !调整箭头大小$/gline,,-1$plvect,s,,,,vect ! 取消单元边界$/dscale,,off$$plnsol,s,1 !关闭结构变形$/dscale,defas$/gline,1 !恢复变形比例缺省，显示单元边界$/Device,vector,1 !设置矢量显示模式$/clabel,,5 !设置等高线上的文字间隔$/contour,,18,-16,,500 !设置等高线数目$Plnsol,s,x!混凝土主压应变与等效应变曲线$Cmsel,s,zhuzi,$vsel,r,loc,z,0,0.175$vsel,r,loc,y,-0.125+ 0.006,0.125-0.006$vsel,r,loc,x,-0.135+0.03,0.175$eslv,s$nsle,s/POST26NUMVAR,200$$Esol,3,1080,26,s,eqv$esol,4,1080,26,epel,eqv$esol,5,1080,26,eppl,eqv$add,4,4,5$prod,5,3,,,,,,1/1000000$/axlab,x,strain-eqv$/axlab,y,sigx-eqv(MPa)$PRVAR,4,5/POST26NUMVAR,200$$Esol,3,1080,45,s,eqv$esol,4,1080,45,epel,eqv$esol,5,1080,45,eppl,eqv$add,4,4,5$prod,5,3,,,,,,1/1000000$/axlab,x,strain-eqv$/axlab,y,sigx-eqv(MPa) $PRVAR,4,5/POST26NUMVAR,200$$Esol,3,2706,46,s,eqv$esol,4,2706,46,epel,eqv$esol,5,2706,46,eppl,eqv$add,4,4,5$prod,5,3,,,,,,1/1000000$/axlab,x,strain-eqv$/axlab,y,sigx-eqv(MPa) $PRVAR,4,5/POST26NUMVAR,200$$Esol,3,3157,296,s,eqv$esol,4,3157,296,epel,eqv$esol,5,3157,296,eppl,eqv$add,4,4,5$prod,5,3,,,,,,1/1000000$/axlab,x,strain-eqv$/axlab,y,sigx-eqv(MPa) $PRVAR,4,5/POST26NUMVAR,200$$Esol,3,3157,291,s,eqv$esol,4,3157,291,epel,eqv$esol,5,3157,291,eppl,eqv$add,4,4,5$prod,5,3,,,,,,1/1000000$/axlab,x,strain-eqv$/axlab,y,sigx-eqv(MPa)$PRVAR,4,5/POST26NUMVAR,200$$Esol,3,3160,1169,s,eqv$esol,4,3160,1169,epel,eqv$esol,5,3160,1169,eppl,eqv$add,4,4,5$prod,5,3,,,,,,1/1000000$/axlab,x,strain-eqv$/axlab,y,sigx-eqv(MPa)$PRVAR,4,5! 剪应变与梁端荷载曲线$Cmsel,s,gzg$vsel,r,loc,z,0,0.0045$vsel,r,loc,y,-0.125+0. 006,0.125-0.006$eslv,s$nsle,s/POST26NUMVAR,200$$RFORCE,3,7150,F,Y$esol,4,6351,2661,epel,xy$esol,5,6351,2661,eppl,xy$add,4,4,5$prod,5,3,,,,,,-1/1000$PRVAR,4,5/post1$Cmsel,s,gzg$vsel,r,loc,z,0,0.0045$vsel,r,loc,y,-0.125+0. 006,0.125-0.006$$vsel,r,loc,x,-0.175,0.175$eslv,s$nsle,sSET,LAST$SET,PREVIOUS*GET,RTIME,ACTIVE,0,SET,TIME !获取最终收敛时间Pldisp$plnsol,s,x !显示变形图与X方向应力分布图Plnsol,epel,xy !Xy弹性应变分布图Plnsol,eppl,xy !xy塑性应变分布图Plnsol,epto,xy !Xy方向总应变分布图/post1Allsel,all!腹板剪力$Cmsel,s,gzg$Nslv,s,1$Nsel,r,loc,x,-0.175,0.175$Nsel,r,loc,z,-0.0045,0.0045 $Nsel,r,loc,y,0,0.05$Esln,s,1$Nsel,r,loc,y,0$fsum!内混凝土剪力Allsel,all$Cmsel,s,zhuzi$Nsel,r,loc,x,-0.175,0.175$Nsel,r,loc,z,-0.0625,0.0625$Nsel,r,loc,y,0,0.05$Esln,s,1$Nsel,r,loc,y,0$fsum/post1Allsel,all!外混凝土剪力$Cmsel,s,zhuzi$Nslv,s,1$Nsel,r,loc,x,-0.175,0.175$Nsel,u,loc,z,-0.0625+0.001,0.0625-0.001 $Nsel,r,loc,y,0,0.05$Esln,s,1$Nsel,r,loc,y,0$fsum!混凝土总剪力Allsel,all$Cmsel,s,zhuzi$Nsel,r,loc,x,-0.175,0.175 $Nsel,r,loc,y,0,0.05$Esln,s,1$Nsel,r,loc,y,0$fsum!总剪力/post1Allsel,all$Nsel,s,loc,x,-0.175,0.175 $Nsel,r,loc,y,0,0.05$Esln,s,1$Nsel,r,loc,y,0$fsum。