自编ANSYS培训讲义

主讲人简介：谢宝来，男，汉族，1975年生，1999年7月毕业于西南交通大学交通土建结构工程专业，现任分院总工程师兼副院长，高级工程师。2003年开始研究ANSYS程序，编写了《弯桥计算程序》和《直桥计算程序》，为本单位ANSYS软件技术带头人，在发布的多篇论文中应用到了ANSYS分析技术，该讲义为多年研究总结的结果，为ANSYS初学者不可多得的好资料。第一讲ANSYS初探（4）关于本培训的目的：-使用图形用户界面和APDL语言来操作ANSYS软件-策划及进行基本的分析（线性、静态的结构分析）ANSYS功能概览：•结构分析•热分析•电磁分析•流体分析（CFD）•耦合场分析-多物理场ANSYS结构分析概览：结构分析用于确定结构的变形、应变、应力及反作用力等。结构分析的类型：•静力分析-用于静态载荷，可以考虑结构的线性及非线性行为，例如：大变形、大应变、应力刚化、接触、塑性、超弹及蠕变等。•模态分析-计算线性结构的自振频率及振形。谱分析是模态分析的扩展，用于计算由于随机载荷引起的结构应力和应变（也叫作响应谱或PSD）。•谐响应分析-确定线性结构对随时间按正弦曲线变化的载荷的响应。•瞬态动力学分析-确定结构对随时间任意变化的载荷的响应，可以考虑与静力分析相同的结构非线性行为。•特征屈曲分析-用于计算线性屈曲临界载荷并确定屈曲形状（结合瞬态动力学分析可以实现非线性屈曲分析）。•专项分析-断裂分析，复合材料分析，疲劳分析。ANSYS除了提供标准的隐式动力学以外，还提供了显示动力学分析模块ANSYS/LS-DYNA。•用于模拟高度非线性，惯性力占支配地位的问题，并可以考虑所有的非线性行为。•它的显式方程是求解冲击、碰撞、快速成型等问题，是目前求解这类问题最有效的方法。有限元分析：有限元分析（FEA）是对物理现象（几何及载荷工况）的模拟，是对真实情况的数值近似。通过划分单元，求解有限个数值来近似模拟真实环境的无限个未知数。ANSYS的分析方法：创建有限元模型-创建或读入几何模型.-定义材料属性.-划分单元（节点及单元）.施加载荷进行求解-施加载荷及载荷选项.-求解.查看结果-查看分析结果.-检验结果（分析是否正确）.1.建立有限元模型Preprocessor/Prep72.施加载荷求解Solution/Solu3.查看结果GeneralPostproc/Post1例1.1一个简支梁的例子：（仅由老师演示使用）

⑴、分析目的利用梁的理论，求跨径为L的简支梁的跨中弯矩(2)、定义问题(理论解)材料扬氏系数E=210x109Pa截面参数H=0.5m，W=0.2m几何参数L=10m边界及外力Ux=Uy=0在x=0Uy=0在x=Lq=1000N/m(3)、系统组合及求解M=qL2/8⑷、检查结果M=1000x102/8=12500Nm/prep7!建立有限兀模型k,1!定义关键点1k,2,10!定义关键点2lstr,1,2!定义直线et,1,beam3!定义单元类型为平面梁单元mp,ex,1,210e9!定义材料的弹性模量mp,prxy,1,1/3!定义材料的泊松比r,1,0.2*0.5,0.5**3*.2/12,.5!定义截面尺寸(A,I,H)latt,1,1,1!线的属性(材料，实常数，单元类型)lesize,1,1,,1!定义单元长度为1米

lmesh,1!分单元/eshape,1eplot!定义单元显示比例!显示单元/solu!施加载荷求解dk,1,ux,0,,,uy!边界条件dk,2,uy,0!边界条件sfbeam,all,1,pres,1000,1000!施加线载荷1000N/mSolve!求解/post1!查看结果etable,mi,smisc,6!定义I端弯矩etable,mj,smisc,12!定义J端弯矩plls,mi,mj,-1!画弯矩图〖重要说明〗程序注解：可采用“！”加注解内容，也可以采用“/com”和“C***”加注解内容，但必须自成一行。例题1.2悬臂梁的挠度(在以后的章节中进行练习)：、分析目的利用梁的理论，求解x=L变形量、定义问题(理论解)材料扬氏系数截面参数几何参数边界及外力E=200x10歹a、定义问题(理论解)材料扬氏系数截面参数几何参数边界及外力E=200x10歹aH=0.01m，W=0.03mL=4m，a=2m，b=2mUx=Uy=Rotz=0在x=0在x=a在a^x^L、系统的组合及求解y=Fa2(3L-a)/6EI+q(3L4-4a3L+aD/24EI⑷、检查结果y(L)=0.02667+0.02733=0.05400m第二讲建立节点(8)有限元分析的第一步便是如何将土木结构系统转化成节点(Node)及单元(Element)所组合的有限元模型(FiniteElementModeling)，该有限元模型与土建结构系统的几何外形一致。有限元模型的建立基本上可以分直接法及间接法。直接法-按照土木结构系统的几何外形建立节点及单元。间接法-通过点、线、面、体，先建立系统的几何外形，然后进行网格(Meshing)分割。坐标系统：空间中任何一点通常可用卡式坐标(Cartesian)、圆柱坐标(Cylinder)或球面坐标(Sphericity)来表示该点的坐标位置，不管是哪种坐标系统都需要3个参数表示该点的正确位置。坐标形式参数确定代码卡式x,y,z0圆柱R,0,z1球面R,e逐2【命令介绍】LOCAL,KCN,KCS,XC,YC,ZC,THXY,THYZ,THZX,PAR1,PAR2定义区域(LOCAL)坐标系统，该命令执行后，ANSYS坐标系统自动更改为新建立的坐标系统。KCN:系统确定代码，大于10的整数KCS:系统属性代码，0(卡式)、1(圆柱)、2(球面)及其他XC,YC,ZC:该区域坐标与整体坐标系统原点的关系THXY,THYZ,THZX:该区域坐标与整体坐标系统Z、X、Y轴的关系MenuPaths:UtilityMenu>WorkPlane>LocalCoordinateSystems>CreateLocalCS>AtSpecifiedLocCSYS,KCN声明坐标系统(CoordinateSystem)oMenuPaths:UtilityMenu>WorkPlane>ChangeActiveCSto>SpecifiedCoordSys相关参考命令：CSLIST,DSYS建立节点：有限元模型的建立是将土木结构转换成很多节点和单元相连接，所以节点即为土木结构中的一个坐标点，指定一个号码和坐标位置。【命令介绍】N,NODE,X,Y,Z,THXY,THYZ,THZX定义节点(Nodedefinition)，号码编排顺序不影响分析结果，节点的建立也不一定要连号，但强烈建议连号。在圆柱坐标系统下X，Y，Z相对应于R，0，Z，在球面坐标系统下X，Y，Z对应于R，0，中。节点重复定义，位置以后定义的为准，即使它已经隶属于单元。NODE:欲建立节点的号码，缺省系统将自动编号。X,Y,Z:节点在目前坐标系统下的坐标位置。MenuPaths:MainMenu>Preprocessor>Modeling>Create>Nodes>InActiveCSNDELE,NODE1,NODE2,NINC删除已建立的节点(NodeDelete)，如果此节点已经隶属于单元，则不能删除，除非先删除隶属的单元。NODE1,NODE2,NINC:预删除节点的范围，起始码NODE1，终止码NODE2，隔号码NINC。(见本讲的重要说明)MenuPaths:MainMenu>Preprocessor>Modeling>Delete>NodesNPLOT,KNUM节点显示(NodePlot)，KNUM=0表示不显示号码，KNUM=1表示显示号码。MenuPaths:UtilityMenu>Plot>NodesNLIST,NODE1,NODE2,NINC,Lcoord,SORT1,SORT2,SORT3节点列表(NodeList)。MenuPaths:UtilityMenu>List>Nodes相关参考命令：NDIST，NKPT

〖重要说明〗范围控制码：(有很多命令均采用此文法)NODE1NODE2NINC说明12011到20所有节点、120同上12021到20所有奇数点22022到20所有偶数点12041,5,9,13,1733号节点例题2.1：第一讲的例1.1/prep7!前处理N,1,0,0,0!在(0,0,0)位置建立节点1N,2,1!同“N,2,1,0,0”N,3,2$N,4,3$N,5,4!“$”可以把几个命令写在一行N,6,5$N,7,6$N,8,7$N,9,8N,,9$N,,10!编号缺省Nlis!节点列表，命令仅写了前4个字母以上参数位置为零时可以空白表示，例如“N,1,0,0,0”可以表示为“N,1”，也可以表示为“N”，再例如“N,5,0,5,0”可以表示为“N,5,,5”。例题2.2:第一讲的例1.2/prep7!前处理N,,.5$N,,1$N,,1.5$N,,2!建立节点1-4N,,2.5$N,,3$N,,3.5$N,,4!建立节点5-8Nplo!显示节点Nplo,1!显示节点及编号Nlis!节点列表

例题2.3：平面桁架1000N/prep7!前处理N$N,,2$N,,4N,,6$N,,8$N,,2,1N,,4,2$N,,6,1例题2.4：平面梁/prep7!前处理Local,11,1,15,-20!定义当前坐标系为柱坐标系N,1,25,126.87$N,2,25,119.50$N,3,25,112.12$N,4,25,104.75N,5,25,97.37$N,6,25,90$N,7,25,82.63$N,8,25,75.25N,9,25,67.88$N,10,25,60.50$N,11,25,53.13Local,12,1,50,-23!定义当前坐标系为柱坐标系

N,12,30.5,124.18N,14,30.5,110.51N,16,30.5,96.84N,18,30.5,83.17N,20,30.5,69.50N,22,30.5,55.83Local,13,1,85,-20$N,13,30.5,117.34$N,15,30.5,103.67$N,17,30.5,90$N,19,30.5,76.33$N,21,30.5,62.66$N,23,30.5,48.99N,24,25,119.50N,12,30.5,124.18N,14,30.5,110.51N,16,30.5,96.84N,18,30.5,83.17N,20,30.5,69.50N,22,30.5,55.83Local,13,1,85,-20N,24,25,119.50$N,25,25,112.12$N,26,25,104.75$N,27,25,97.37N,28,25,90N,32,25,60.50Csys,0N,34,30,-2.5N,38,30,-7.5$N,29,25,82.63$N,33,25,53.13$N,35,70,-2.5$N,39,70,-7.5$N,30,25,75.25$N,31,25,67.88!转换当前坐标为直角坐标系$N,36,30,-5$N,37,70,-5$N,40,30,-10$N,41,70,-10例题2.5：平面挖孔板/prep7!前处理N,1,,-2$N,2,,-1$N,3$N,4,,1N,5,,2$N,6,4,-2$N,7,4,-1$N,8,4N,9,4,1$N,10,4,2$N,11,8,-2$N,12,8,-1N,13,8$N,14,8,1$N,15,8,2$N,16,12,-2N,17,12,-1$N,18,12$N,19,12,1$N,20,12,2N,21,2,-2$N,22,6,-2$N,23,10,-2$N,24,2,2N,25,6,2$N,26,10,2$Local,11,1,2N,27,1$N,28,1,30$N,29,1,60$N,30,1,90N,31,1,120$N,32,1,150$N,33,1,180$N,34,1,210N,35,1,240$N,36,1,270$N,37,1,300$N,38,1,330Local,12,1,6,,,,,,1.6!定义柱坐标，并把Y轴扩大1.6倍N,39,1$N,40,1,30$N,41,1,60$N,42,1,90N,43,1,120$N,44,1,150$N,45,1,180$N,46,1,210N,47,1,240$N,48,1,270$N,49,1,300$N,50,1,330Local,13,1,10$N,51,1$N,52,1,30$N,53,1,60N,54,1,90$N,55,1,120$N,56,1,150$N,57,1,180N,58,1,210$N,59,1,240$N,60,1,270$N,61,1,300N,62,1,330【命令介绍】FILL,NODE1,NODE2,NFILL,NSTRT,NINC,ITIME,INC,SPACE节点填充(Fill)命令是自动将两节点间在现有坐标系统下填充许多节点，两节点间填充的节点个数及分布状态视其参数而定。MenuPaths：MainMenu>Preprocessor>Modeling>Create>Nodes>FillbetweenNdsNGEN,ITIME,INC,NODE1,NODE2,NINC,DX,DY,DZ,SPACE节点复制(NodeGeneration)命令是将一组节点在现有坐标系统下复制到其他位置。ITIME:复制次数，包含自己本身。INC:每次复制时节点号码的增加量。NODE1,NODE2,NINC:选取要复制的节点。DX,DY,DZ:每次复制时在现有坐标系统下，几何位置的改变量。MenuPaths：MainMenu>Preprocessor>Modeling>Copy>Nodes>Copy例题2.6第一讲的例1.1

/prep7!前处理N,1$N,11,10Fill,1,11!填充节点Nlis例题2.7本讲的例2.4/prep7!前处理N,1$N,11,30$N,23,70$N,33,100Local,11,1,15,-20!定义当前坐标系为柱坐标系Fill,1,11!按当前坐标系进行填充Local,12,1,50,-23!定义当前坐标系为柱坐标系Fill,11,23!按当前坐标系进行填充Local,13,1,85,-20!定义当前坐标系为柱坐标系Fill,23,33!按当前坐标系进行填充Csys,0!转换当前坐标为直角坐标系N,34,30,-2.5$N,35,70,-2.5Ngen,4,2,34,35,1,,-2.5!复制34,35两节点例题2.8本讲的例2.5/prep7!前处理N,1,,-2$N,5,,2$Fill,1,5Ngen,4,5,1,5,1,4!复制节点Local,11,1,2$N,27,1Ngen,12,1,27,,,,30!复制节点Local,12,1,6,,,,,,1.6$N,39,1Ngen,12,1,39,,,,30!复制节点Local,13,1,10$N,51,1Ngen,12,1,51,,,,30!复制节点Csys,0$N,21,2,-2$Ngen,3,1,21,,,4Ngen,2,3,21,23,1,,4第三讲建立单元(8)当节点建立完成后，必须使用适当元素，将土建结构按照节点连接成单元，并完成其有限元模型。单元选择正确与否直接影响分析结果。每一个单元都有不同编号，例如Linkl是第一号单元、Beam3是第三号单元，单元编号也可以代表此单元。每个单元号码前的名称可判断适用范围或其形状，基本上单元形状可分为线单元、平面单元和实体单元。建立单元前必须先行定义使用者欲选择的单元形式代码、单元材料特性、单元几何特性等。熟记下表中最常用的10种单元:编号单元名称实常数用途1Beam3AREA,IZZ,HEIGHT平面梁2Link8AREA,ISTRN索3Pipe16OD,TKWALL管4Mass21质量5Beam44AREA1,IZ1,IY1,TKZB1,TKYB1,IX1空间梁6Shell63TK(I),TK(J),TK(K),TK(L)薄板4节点7Solid65MAT1,VR1,THETA1,PHI1,…混凝土8节点8Plane82THK当(KEYOPT(3)=3)8节点9Solid186无6面体20节点10Beam188无变截面梁【命令介绍】ET,ITYPE,Ename,KOP1,KOP2,KOP3,KOP4,KOP5,KOP6,INOPR单元类型(ElementType)为构成土建结构系统所含的单元种类。ITYPE:单元类型号码。Ename：单元库名称，例如：Beam3，也可以用“3”代替。KOP1-KOP6:单元特性编码。MenuPaths:MainMenu>Preprocessor>ElementType>Add/Edit/DeleteMP,Lab,MAT,C0,C1,C2,C3,C4定义材料特性(MaterialProperty)oMAT：对应前面所定义单元类型号码(ITYPE),表示该组材料隶属于此单元。Lab：材料特性的类别，常用的类别见下表：Lab物理意义EX弹性模量ALPX热膨胀伸缩系数PRXY泊松比DAMP阻尼系数DENS密度C0：材料对应的值。MenuPaths:MainMenu>Preprocessor>MaterialProps>MaterialModelsR,NSET,R1,R2,R3,R4,R5,R6R命令(RealConstant)为所定义单元类型的几何特性。NSET:几何特性编码。R1-R6:几何特性值，要参考相应的单元说明。MenuPaths:MainMenu>Preprocessor>RealConstants>Add/Edit/DeleteE,I,J,K,L,虬N,O,P定义单元(Element)的连接方式，梁单元也可能用到参数K，但此时节点K与I、J不在一条线上，通过K节点来确定单元主惯性轴的方向，与I、J、K三点确定的平面的方向相同。MenuPaths:MainMenu>Preprocessor>Modeling>Create>Elements>AutoNumbered>ThruNodesEPLOT单元显示(ElementPlot)。MenuPaths:UtilityMenu>Plot>ElementsEGEN,ITIME,NINC,IEL1,IEL2,IEINC,MINC,TINC,RINC,CINC,SINC,DX,DY,DZ单元复制(ElementGeneration)，是将一组单元在现有坐标系统下复制到其他位置。如果事先建立了节点，位置的确定根据节点走，则忽略DX，DY，DZ参数，如果事先没建立节点，则根据DX，DY，DZ参数自动建立节点。ITIME:复制的次数，包含自己本身。NINC:每次复制单元时，相对应节点号码的增加量。IEL1,IEL2,IEINC:选取复制单元。DX,DY,DZ:新建节点几何位置的改变量，如果节点已经存在，则按现有节点考虑。MenuPaths:MainMenu>Preprocessor>Modeling>Copy>Elements>AutoNumbered/PNUM,Label,KEY在有限元模块图形中显示号码(PlotNumber)oLabel:欲显示对象的名称Label意义Label意义NODE节点SEC截面ELEM单元MAT材料KP关键点TYPE单元类型LINE线REAL实常数AREA面LOC坐标系VOLU体ESYS单元坐标系KEY:控制码，0不显示号码(系统默认)，1显示号码。MenuPaths:UtilityMenu>PlotCtrls>Numbering相关参考命令：ELIST，EDEL，/Eshape，EMODIF例题3.1：第一讲例1.1/prep7!前处理Et,1,beam3!定义平面梁单元Mp,ex,1,210e9!定义材料的弹性模量Mp,prxy,1,1/3!定义材料的泊松比r,1,0.2*0.5,0.5**3*.2/12,.5!定义截面尺寸(A,I,H)N$N,2,1$E,1,2Egen,10,1,1,,,,,,,,1!复制单元/Eshape,1!定义单元显示比例/Pnum,elem,1$Eplot!显示单元编号

例题3.2：第二讲例2.3/prep7!前处理Et,1,Link1!杆单元Mp,ex,1,80e9!定义材料的弹性模量Mp,prxy,1,1/3!定义材料的泊松比R,1,10e-4!实常数(A)N$N,5,8$Fill,1,5$N,6,2,1N,7,4,2$N,8,6,1$E,1,2$E,1,6E,8,5$E,7,2$E,7,3$E,7,4Egen,7,1,1$/Pnum,elem,1$Edele,10$E,2,6E,4,8$/eshape,1$Eplot例题3.3：第二讲例2.4接例题2.4Et,1,44!空间梁单元Beam44Mp,ex,1,30e9!定义材料的弹性模量Mp,prxy,1,1/6!定义材料的泊松比R,,1*4,1**3*4/12,1*4**3/12,2,.5!实常数(A,IZ,IY,Z/2,Y/2)E,1,2$egen,32,1,1$E,11,34$E,23,35E,34,36$E,35,37$Egen,3,2,35,36,1/eshape,1$eplot例题3.4：第二讲例2.5接例题2.5Et,1,Shell63!空间壳单元Shell63Mp,ex,1,210e9!定义材料的弹性模量Mp,prxy,1,1/3!定义材料的泊松比R,1,.1!厚度0.1m

E,1,2,34,35E,5,24,30,31E,8,7,38,27E,6,7,46,47E,10,25,42,43E,13,12,50,39E,11,12,58,59E,15,26,54,55E,18,17,62,51/eshape,1$E,2,3,33,34$E,24,10,29,30$E,7,6,37,38$E,7,8,45,46E,1,2,34,35E,5,24,30,31E,8,7,38,27E,6,7,46,47E,10,25,42,43E,13,12,50,39E,11,12,58,59E,15,26,54,55E,18,17,62,51/eshape,1$E,2,3,33,34$E,24,10,29,30$E,7,6,37,38$E,7,8,45,46$E,25,15,41,42$E,12,11,49,50$E,12,13,57,58$E,26,20,53,54$E,17,16,61,62$Eplot$E,3,4,32,33$E,10,9,28,29$E,6,21,36,37$E,8,9,44,45$E,15,14,40,41$E,11,22,48,49$E,13,14,56,57$E,20,19,52,53$E,16,23,60,61$E,4,5,31,32$E,9,8,27,28$E,21,1,35,36$E,9,10,43,44$E,14,13,39,40$E,22,6,47,48$E,14,15,55,56$E,19,18,51,52$E,23,11,59,60铰接处两个梁的端点有两个不同的Node，但坐标一样，意义在于它们的平动位移一致，则把这两个Node的平动自由度耦合即可。ANSYS可采用两种方法来实现铰接：在同一位置用2个Node，然后CP。只用一个Node，然后根据需要用一些可以释放某些自由度的单元（通过设置keyopt来实现），如Beam44。推荐采用第一种，因直观易懂，不容易出错。【命令介绍】CP,NSET,Lab,NODE1,NODE2,NODE3,NODE4,NODE5,NODE6,NODE7,NODE8,NODE9,NODE1。，NODE11,NODE12,NODE13,NODE14,NODE15,NODE16,NODE17把相关节点进行自由度藕合。NSET:藕合编号，当Lab=All时可缺省，当等于NEXT可以自动增加。Lab:要藕合的自由度，Ux,Uy,Uz,Rotx,Roty,Rotz,All等。NODE1-NODE17:要耦合的节点。

MenuPaths:MainMenu>Preprocessor>Coupling/Ceqn>CoupleDOFs相关参考命令：CPLIST,CPDELE例题3.5做下列连续梁的弯矩图4000N10000N6000N/mQ1n/h;21;21b2Jwb2>44**/Prep7Et,1,Beam3!平面梁单元Beam3Mp,ex,1,210e9!定义材料的弹性模量Mp,prxy,1,1/3!定义材料的泊松比R,1,1,1/12,1!定义实常数，静定，不影响结果N,1$N,3,2$N,4,2$N,8,6N,9,6$N,17,14$Fill,1,3$Fill,4,8Fill,9,17$E,1,2$Egen,16,1,1Edele,3,8,5!删除0长度单元Cp,1,ux,3,4$Cp,2,uy,3,4$Cp,3,ux,8,9$Cp,4,uy,8,9〖重要说明〗定义截面类型对于Beam44单元，可以定义截面类型，但是如果定义了实常数，则定义的截面类型将被忽略。但对于Beam188/Beam199单元，就必须定义截面类型。截面可以是组合的，例如叠合梁和钢管混凝土，也可以是变截面的，这些结构形式均能在ANSYS中实现，后面章节将作为详细学习。ANSYS提供了11种常见的截面类型，我们可以选择与实际结构截面相近的类型来使用，下表中列出了Beam单元的11种常见截面：

编号Subtype意义编号Subtype意义1RECT矩形7ZZ形2QUAD4边形8LL形3CSOLID圆形9TT形4CTUBE圆环10HATS几形5CHAN槽形11HRES箱形6I工形【命令介绍】SECTYPE,SECID,Type,Subtype,Name,REFINEKEY设置截面类型(SectionType)信息。SECID：截面号码。Type：单元类型，等截面梁(BEAM),变截面梁(TAPER)等。Subtype：截面类型，参见本讲11种常见截面表。Name:截面名字，可以缺省，建议填写截面说明。MenuPaths:MainMenu>Preprocessor>Sections>Beam>CommonSectionsSECOFFSET,Location,OFFSET1,OFFSET2,CG-Y,CG-Z,SH-Y,SH-Z设置截面形心位置。Location:见下表Location意义CENT质心SHRC剪切中心ORIGIN边缘USER用户定义OFFSET1,OFFSET2:当Location为USER时，设定偏心值。MenuPaths:MainMenu>Preprocessor>Sections>Beam>CommonSectionsSECDATA,VAL1,VAL2,VAL3,VAL4,VAL5,VAL6,VAL7,VAL8,VAL9,VAL10设置截面数据，必须根据截面类型确定参数，使参数与截面类型匹配。MenuPaths:MainMenu>Preprocessor>Sections>Beam>CommonSections

相关参考命令：SECWRITE,SECREAD,SDELETE,SLIST,SECPLOT,SECNUM例题3.6第一讲例1.1定义两种截面形式，先采用第二种截面，显示单元，然后更改单元属性，把单元截面改为第一种，然后显示单元。截面1：工形，顶底板宽0.2m，厚度0.01m，高度0.5m，腹板厚度0.01m截面2：环形，外径0.5m，内径0.4m/Prep7Et,1,beam44/Prep7Et,1,beam44Mp,ex,1,210e9Mp,prxy,1,1/3$NNgen,2,11,1,10,1,,1Sectype,1,beam,ISecdata,.2,.2,.5,.01,.01,.01Sectype,2,beam,CTUBESecdata,.4,.5Secnum,2E,1,2,12$Egen,10,1,1/Pnum,sec,1Emodif,all,secnum,1!空间梁单元，用他代替beam3!定义材料的弹性模量$N,11,10$Fill,1,11!建立控制单元主惯方向的节点$Secoffset,cent!上3条命令定义“工”形截面$Secoffset,cent!上3条命令定义“圆”形截面!指定截面号为2$/Eshape,1$Eplot!显示截面编号$Eplot第四讲加载求解(10)有限元最主要目的，在于了解结构系统组件受外力载荷后的反映，顾明确定义适当、正确的载荷，这对于设计分析工作是非常重要的。载荷可分为边界条件和实际外力两大类。不同领域的载荷：结构力学位移、集中力、压力、温度、重力热八、、学温度、热流率、热源、对流、无限表面磁学磁场、磁通量、磁源密度、无限表面电学电位、电流、电荷、电荷密度流体力学速度、压力【命令介绍】荷载定义ANTYPE,Antype,Status,LDSTEP,SUBSTEP,Action声明分析类型(AnalysisType)，进入求解状态(/Solu)，必须先确定分析类型系统默认为静态分析。Antype:分析类型，也可以用代码代替：Antype代码意义STATIC0静态分析(系统默认)BUCKLE1屈曲分析MODAL2模态分析HARMIC3谐响应分析TRANS4瞬态动力学分析SUBSTR7子结构分析SPECTR8响应谱分析MenuPaths:MainMenu>Solution>AnalysisType>NewAnalysisF,NODE,Lab,VALUE,VALUE2,NEND,NINC定义节点上的集中力(Force)。NODE:节点号，也可以为ALL(当前激活的所以节点)。Lab:外力形式：FX,FY,FZ,MX,MY,MZ。VALUE:外力的大小。NODE,NEND,NINC:施加节点的范围，在建立节点时应先规划节点号码，方便选择。MenuPaths:MainMenu>Solution>DefineLoads>Apply>Structural>Force/Moment>OnNodesNODE,Lab,VALUE,VALUE2,NEND,NINC,Lab2,Lab3,Lab4,Lab5,Lab6定义节点自由度。NODE,NEND,NINC:节点的范围。Lab:约束形式：DX,DY,DZ线位移,ROTX,ROTY,ROTZ角位移，也可以是ALL。VALUE:位移大小。Lab2-Lab6:每个节点依单元不同最多6个自由度，Lab只约束了一个，其他自由度约束由此定义。MenuPaths:MainMenu>Solution>DefineLoads>Apply>Structural>Displacement>OnNodesSFBEAM,ELEM,LKEY,Lab,VALI,VALJ,VAL2I,VAL2J,IOFFST,JOFFST定义作用于梁单元上的均布力。ELEM:元素编号，也可以是ALL。LKEY:依节点顺序梁单元有4个面，该参数为定义分布力所施加面的号码。Lab:PRES(表示压力)。VALI,VALJ:在I点及J点均布力的大小。MenuPaths:MainMenu>Solution>DefineLoads>Apply>Structural>Pressure>OnBeamsOMEGA,OMEGX,OMEGY,OMEGZ,KSPIN定义结构等速旋转(OMEGA)负载效应，OMEGX,OMEGY,OMEGZ为结构相对于整体坐标系统X,Y,Z方向的旋转角速度的值。MenuPaths:MainMenu>Solution>DefineLoads>Apply>Structural>Inertia>AngularVeloc>GlobalACEL,ACELX,ACELY,ACELZ定义结构线性加速度负载效应，ACEX,ACEY,ACEZ为结构相对应整体坐标系统X,Y,Z方向的加速度的值。注意方向关系，Z轴负方向的加速度，则ACEL,,,9.81。MenuPaths:MainMenu>Solution>DefineLoads>Apply>Structural>Inertia>

Gravity>Global相关参考命令：SFE,SF【命令介绍】后处理器PLESOL,Item,Comp,KUND,Fact显示单元的解答，以轮廓线方式表达，通常2-D单Item,Comp:见下表：元及3-D单元才适用。ItemComp意义SX,Y,Z,XY,YZ,XZ应力S1,2,3主应力SEQV,INT等效应力FX,Y,Z结构力MX,Y,Z结构力矩KUND:KUND=0,显示变形后的结构形式KUND=1,同时显示变形前和变形后的结构形式KUND=2,同时显示变形前和变形后的结构形式，但仅显示结构外观MenuPaths:MainMenu>GeneralPostproc>PlotResults>ContourPlot>ElementSoluPLNSOL,Item,Comp,KUND,Fact,FileID显示节点的解答，以轮廓线方式表达。Item,Comp:见下表：ItemComp意义SX,Y,Z,XY,YZ,XZ应力S1,2,3主应力SEQV,INT等效应力UX,Y,Z,SUM位移分量及向量位移ROTX,Y,Z,SUM旋转位移分量及向量旋转位移TEMP温度MenuPaths:MainMenu>GeneralPostproc>PlotResults>ContourPlot>NodalSoluPRESOL,Item,Comp

列表显示单元的解答，以单元为单位，将Item所声明项次的结果显示在窗口中，使用者可以保存该资料。Item,Comp:见下表：ItemComp意义SX,Y,Z,XY,YZ,XZ应力FX,Y,Z结构力MX,Y,Z结构力矩MenuPaths:MainMenu>GeneralPostproc>ListResults>ElementSolutionPRNSOL,Item,Comp列表显示节点的解答，以节点为单位，将Item所声明项次的结果显示在窗口中，使用者可以保存该资料。Item,Comp:见下表：ItemComp意义UX,Y,ZX,Y,Z位移UCOMPX,Y,Z位移及向量位移SCOMP应力SPRIN主应力，等效应力MenuPaths:MainMenu>GeneralPostproc>ListResults>NodalSolutionPRRSOL,Lab列表显示支反力计算结果。Lab:显示的FX,FY,FZ,MX,MY,MZ分量，缺省显示全部MenuPaths:MainMenu>GeneralPostproc>ListResults>ReactionSoluETABLE,Lab,Item,Comp将单元(Element)某项结果制作成表格(Table)形式。Lab:字段名称，最多8个字符。Item,Comp:见下表：ItemComp意义UX,Y,ZX,Y,Z位移ROTX,Y,ZX,Y,Z转角TEMP温度SX,Y,Z,XY,YZ,XZ应力S1,2,3主应力SEQV,INT等效应力FX,Y,Z结构力MX,Y,Z结构力矩SMISCsnum参见单元帮助NMISCsnum参见单元帮助MenuPaths:MainMenu>GeneralPostproc>ElementTable>DefineTablePLLS,LabI,LabJ,Fact画出节点表数据，LabI为I节点表，LabJ为J节点表，Fact为图形放大系数，不放大结果，当为-1时则改变图形的显示方向。MenuPathsMainMenu>GeneralPostproc>PlotResults>ContourPlot>LineElemResPRETAB,Lab1,Lab2,Lab3,Lab4,Lab5,Lab6,Lab7,Lab8,Lab9列表显示表格资料，Lab1-Lab9为前面定义的表格字段名称，缺省显示全部的表。MenuPaths:MainMenu>GeneralPostproc>ElementTable>ListElemTablePLETAB,Itlab,Avglab绘出结果表格。Itlab：为前面所定义的表格字段名称。MenuPaths:MainMenu>GeneralPostproc>ElementTable>PlotElemTable〖重要说明〗隐循环-建模的利器当我们学习到APDL语言的时候，会讲到Do循环，其实ANSYS还提供了一种隐循环，其操作速度更快。例如：N,(1:11:1),(0:10:1)，一条命令生成了11个节点！数据表达式(IVAL:FVAL:INC)，IVAL为开始变量，FVAL为结束变量，INC为步长，循环次数为(FVAL-IVAL)/INC+1，注意用冒号分隔，容易误写逗号。

N,(1:11),(0:10)!生成11个节点，步长1可缺省E,(1:10),(2:11)!生成10个单元F,(1:11:2),fy,-1000!在奇数节点上施加力Sfbeam,(4:7),2,pres,1000!在4,5,6,7梁单元上加均布力例题4.1第一讲（例1.1），画弯矩和剪力图，并显示支反力，最后列表显示/Prep7!前处理Et,1,beam3$R,1,.2*.5,.2*.5**3/12,.5Mp,ex,1,210e9$Mp,prxy,1,1/3N,(1:11),(0:10)$E,(1:10),(2:11)/Eshape,1$Eplo/Solu!求解D,1,ux,0,,,,uy!约束1节点的ux,uyD,11,uy,0!约束11节点的uySfbeam,all,1,pres,1000!施加均布力1000N/mSolve!求解/Post1!后处理Etable,Mi,smisc,6!定义I端弯矩Etable,Mj,smisc,12!定义J端弯矩Etable,Qi,smisc,2!定义I端剪力Etable,Qj,smisc,8!定义J端剪力Plls,Mi,Mj,-1!绘制弯矩图Plls,Qi,Qj,-1!绘制剪力图Prrsol!列表显示支反力Pretab!列表显示：Elem,Mi,Mj,Qi,Qj例题4.2第一讲（例1.2），求悬臂梁梁端位移/Prep7!前处理Et,1,beam44$mp,ex,1,200e9

mp,prxy,1,1/3R,1,.03*.01,.03*.01**3/12,.03**3*.01/12,.03/2,.01/2N,(1:9),(0:4:.5)$E,(1:8),(2:9)/solu!求解器$Sfbeam,(5:8),2,pres,.5!定义均布力D,1,all$F,5,fy,-2Solve/post1!后处理Plnsol,u,y!画节点位移图Prnsol,u,y!列表显示节点位移例题4.3第三讲（例3.2），求桁架各杆内力、应力及应变/Prep7!前处理Et,1,Link1$mp,ex,1,80e9mp,prxy,1,1/3$r,1,10e-4N,(1:5),(0:8:2)$N,(6:8:2),(2:6:4),1N,7,4,2$E,(1:4),(2:5)E,1,6$E,8,5$E,(6:7),(7:8)E,2,6$E,4,8$E,2,7$E,4,7E,3,7/solu!求解器D,1,all,,,5,4$F,all,fy,-1000Fdele,(2:4),fy!删除多余的荷载Solve/post1!后处理Etable,nn,smisc,1!轴力列表Etable,sx,ls,1!应力列表Etable,dt,lepel,1!应变列表Pretab!列表显示：Elem,nn,sx,dt

例题4.4第三讲（例3.3），列出在自重下的弯矩、剪力、轴力表接第三讲（例3.3）!前处理Mp,dens,1,2500!定义密度，混凝土标号C30/solu!求解器Acel,,9.81!Y轴负方向的重力加速度，正的9.81D,(1:33:32),all$D,(40:41),allSolve/post1Etable,mi,smisc,6!定义I端弯矩Etable,mj,smisc,12!定义J端弯矩Etable,qi,smisc,2!定义I端剪力Etable,qj,smisc,8!定义J端剪力Etable,ni,smisc,1!定义I端轴力Etable,nj,smisc,7!定义J端轴力Pretab!列出弯矩、剪力、轴力表例题4.5第三讲（例3.4），求在q=2000kN/m下的应力分布图接第三讲（例3.3）!前处理D,(1:5),all$F,(16:20:4),fx,2000000*.5F,(17:19),fx,2000000*1/solu$solve/post1!后处理Plnsol,s,x!x方向正应力例题4.6第三讲（例3.5），画连续梁的弯矩剪力图/Prep7Et,1,Beam3!平面梁单兀Beam3Mp,ex,1,210e9!定义材料的弹性模量Mp,prxy,1,1/3R,1,1,1/12,1N,(1:3),(0:2)N,(9:17),(6:14)(4:7),(5:8)Cp,1,ux,3,4$Cp,2,uy,3,4/SoluD,1,all3,fy,-4000Mp,prxy,1,1/3R,1,1,1/12,1N,(1:3),(0:2)N,(9:17),(6:14)(4:7),(5:8)Cp,1,ux,3,4$Cp,2,uy,3,4/SoluD,1,all3,fy,-4000Sfbeam,(11:14),1,pres,6000/post1Etable,mi,smisc,6Etable,mj,smisc,12Etable,qi,smisc,2Etable,qj,smisc,8Plls,mi,mj,-1Plls,qi,qj,-1!定义材料的泊松比!定义实常数，静定，不影响结果$N,(4:8),(2:6)$E,(1:2),(2:3)$E,(9:16),(10:17)$Cp,3,ux,8,9$Cp,4,uy,8,9!求解器$D,(11:15:4),uy$F,6,fy,-10000$solve!后处理!定义I端弯矩!定义J端弯矩!定义I端剪力!定义J端剪力!绘制弯矩图!绘制剪力图6000N/mM/Mooog/Prep7Et,1,beam3Mp,ex,1,210e9/Prep7Et,1,beam3Mp,ex,1,210e9N,(1:3),(0:2),4N,(6:10),(2:6),4E,(1:4),(2:5)/soluSfbeam,(1:2),1,pres,6000!前处理$r,1,1,1/12,1$Mp,prxy,1,1/3$N,(4:5),2,(3:2:-1)$N,(11:14),6,(3:0:-1)$E,(6:13),(7:14)!求解器$Sfbeam,(5:8),1,pres,6000$D,14,ux,,,,,uy$Cp,2,uy,3,6!后处理!定义I端弯矩!定义J端弯矩!定义I端剪力!定义J端剪力!定义I端轴力!定义J端轴力!弯矩图!剪力图!轴力图Sfbeam,(9:12),1,pres,(0:4500:1500),(1500:6000:1500)D,5,uy,,,,,rotzCp,1,ux,3,6Solve/post1Etable,mi,smisc,6Etable,mj,smisc,12Etable,qi,smisc,2Etable,qj,smisc,8Etable,ni,smisc,1Etable,nj,smisc,7Plls,mi,mj,-1Plls,qi,qj,-1Plls,ni,nj,-1第五讲基本命令(6)有限元模型的直接建立法中主要的对象有节点及单元，我们已经学习了许多相关命令，实体模型包含关键点、线、面和体等对象，从对象的角度来说许多命令非常相似。通过前面几讲的学习，我们能通过直接法建立一些简单的模型，接下来我们将把直接法建模的技术提高到最高层次，建立复杂的桥梁模型，然后将学到建立实体模型，这样我们将接触更多的命令，为了更好的记忆这些命令及为下一步建立复杂的模型做准备，特在此讲来介绍群组的基本命令。1.用“*”来代表下述字符：K,L,A,V,N,E,CM,CP,CS,ET,MP,R,ICKKeypoint关键点LLine线AArea面VVolume体NNode节点EElement单元CMCoMponent组CPCouPled耦合CSCoordinateSystem坐标系ETElementType单元类型MPMaterialProperty材料特性RRealconstant实常数ICInitialConditions节点初始条件2.用“$”代表下述字符：D,F,SF,BFDDOFconstraint节点约束FForceload节点施加集中力SFSurfaceForce节点施加面何载BFBodyForce节点施加体荷载3.许多命令由“$*”字符构成DKDOFconstraintsatKp关键点上施加约束

DLDOFconstraintsonLines线上施加约束DADOFconstraintsonAreas面上施加约束FKForceloadsatKp关键点上施加集中力SFLSurfaceForceonLines线上施加面荷载SFASurfaceForceonAreas面上施加面荷载SFESurfaceForceonElement单元上施加面荷载BFKBodyForceasKp关键点上施加体荷载BFLBodyForceonLines线上施加体荷载BFABodyForceonAreas面上施加体荷载BFVBodyForceonVolumes体上施加体荷载BFEBodyForceonElements单元上施加体荷载4.列表显示命令由“*Lis”、"$Lis”和“$*Lis”组成5.图形显示命令由“*Plo”组成KPlo绘制关键点NPlo绘制节点LPlo绘制线EPlo绘制单元APlo绘制面CMPlo绘制组VPlo绘制体6.删除命令由“*DeT'、"$Del”和“$*Del”组成Kdel,P删除关键点SFDel,P删除节点面荷载LDel,P删除线BFDel,P删除节点体荷载ADel,P删除面DKDel,P删除关键点约束VDel,P删除体DLDel,P删除线约束NDel,P删除节点DADel,P删除面约束EDel,P删除单元FKDel,P删除关键点集中力CMDel,P删除组SFLDel,P删除线上均布力CPDel,P删除耦合SFADel,P删除面上均布力CSDel,P删除坐标系SFEDel,P删除单元上均布力ETDel,P删除单元类型BFKDel,P删除关键点体荷载MPDel,P删除材料特性BFLDel,P删除线上体荷载

RDel,P删除实常数BFADel,P删除面上体荷载ICDel,P删除节点初始条件BFVDel,P删除体上体荷载DDel,P删除节点约束BFEDel,P删除单元上体荷载FDel,P删除节点集中力7.两个实体的距离由“*dis,P”组成Ndis,PDistancebetweentwonodes节点距离Kdis,PDistancebetweentwoKPs关键点距离8.几何分网由“*MES,P”组成KMES,P关键点分网LMES,P线分网AMES,P面积分网VMES,P体积分网分网前划分线段和面积大小指定：“*size,P”Lesize,P直线划分Aesize,P面积划分清除网格：“*CLE,P”KCLE,P清除关键点网格LCLE,P清除线网格ACLE,P清除面网格VCLE,P清除体网格布尔运算：相交“*IN*”AINA,P面面相交VINV,P体体相交LINL,P线线相交AINV,P面体相交LINV,P线体相交LINA,P线面相交布尔运算：相加“*ADD”AADD,P面相加VADD,P体相加布尔运算：粘合"*GLUE”VGLUE,P体粘合AGLUE,P面粘合LGLUE,P线粘合布尔运算：减、分开“*sb*,p”AsbA,PSubtractAreafromArea面积相减AsbL,PDivideAreabyLine面被线分开VsbA,PDivideVolumebyArea体被面分开LsbW,PDivideLinebyWorkplane线被工作平面分开

VsbW,PDivideVolumebyWorkplane体被工作平面分开AsbW,PDivideAreabyWorkplane面被工作平面分开VsbV,PSubtractVolumebyVolume体体相减15.布尔运算：交叠"*OVLAP,P”AOVLAP,P面积交叠VOVLAP,P体积交叠16.拖拉操作："*drag,P”Vdrag,PDragareasalongalineto面沿直线拖拉生成体Adrag,PcreateanewvolumeDraglinealongalineto直线沿直线拖拉生成面Ldrag,PcreateanewareaDragKPalongalinetocreate关键点沿直线拖拉生成线anewline17.复制对象："*GEN,,P”NGEN,,P节点复制EGEN,,P单元复制KGEN,,P关键点复制LGEN,,P线复制AGEN,,P面复制VGEN,,P体复制18.从下到上建模：“*,P&**,P”L,PCreatelinesfromKP通过关键点连线AK,PCreateareafromKP通过关键点建面AL,PCreateareafromlines通过线建面V,PCreatevolumefromKP通过关键点建体VA,PCreatevolumefromareas通过面建体E,PCreateelementfromnodes通过节点建立单元EN,PCreateelementfromnodes通过节点建立单元（可编号）19.重新分网：“*ref,P”Kref,P关键点重新分网Lref,P线重新分网Aref,P面重新分网Vref,P体重新分网Nref,P节点重新分网Eref,p单元重新分网20.将荷载从实体模型转换到有限元模型：“*tran”Dtran约束转换Ftran集中力转换sftran表面荷载转换bftran体荷载转换〖重要说明〗选择命令，用“*”代表K,L,A,V,N,E对像，选择命令由*Sel组成，谈到选择，就不得不提到编组命令CM,然后用CMsel命令选择这些编好的组，编组命令CM其语法为：CM,Cname,EntityCname：本编组的名称，可以用单引号引上也可以不用Entity：VOLU、AREA、LINE、KP、ELEM和NODE，分别代表给当前激活的体、面、线、关键点、单元和节点编组*SEL,Type,Item,Comp,VMIN,VMAX,VINC,KSWPType:选择的方式Type意义S选择一组对象为激活对像（缺省）R在现有激活对象中，重新选择某些激活对像A再在没有选择的对像中，选择某些对像加入现在的激活对像中U在现在的激活对像中排除一些对像ALL选择所有对像为激活对像NONE与ALL相反INVE逆激活STAT显示选择信息Item：选择依据Item意义P点选组名CM命令定义的组名对像VOLU、AREA、LINE、KP、ELEM、NODE坐标LOC（ElEM不能通过坐标选）属性MAT、TYPE、REAL、SECN等Comp:当Item为LOC时，Comp为x,y,z；其他将缺省本参数VMIN,VMAX,VINC：范围

ItemVMIN,VMAX,VINC意义P缺省组名组名，VMAX和VINC缺省对像对像编号集坐标表示Comp(x,y,z)的范围，VINC缺省属性属性编号集KSWP:是否选择隶属于本对像的低级元素MenuPaths:UtilityMenu>Select>Entities*SL$,Type,ARKEY通过对像$来选择对像*，参见下表：*$NE，K，L，A，VEN，L，A，VKN，LLK，AAL,VVAType:选择的方式，仅能是S、R、A和U,缺省为S。KSWP:是否选择隶属于本对像的低级元素MenuPaths:UtilityMenu>Select>EntitiesALLSEL,LabT,Entity选择和激活所有对像，LabT和Entity可以缺省MenuPaths:UtilityMenu>Select>Everything第六讲高级使用（12）目前为止，我们所建立的模型均狠简单，单一的单元类型，单一的材料属性，单一的几何参数，例子中也没有使用截面命令，复杂的例题会浪费大量的建模时间，这样会降低我们的学习兴趣，但现在我敢说，我们已经喜欢上了此软件，而且通过前几讲的学习，现在我们有能力建立高级的模型了。复杂的模型中经常包含多种单元类型、多种材料、多种几何参数和多种截面类型，这时我们在建立单元时，需要经常更改单元类型、材料类型、几何参数和截面类型，完成更改操作用以下命令：更改单元类型TYPE,ITYPE建立单元类型ET,ITYPE,Ename,…更改几何参数REAL,NSET建立几何参数R,NSET,…更改截面类型SECNUM,SECID建立截面类型SECTYPE,SECID,…更改材料类型MAT,MAT建立材料类型MP,Lab,MAT,…对于同时定义了几何参数和截面类型的单元，其几何参数优先，例如BEAM44，除非此单元不支持几何参数，例如BEAM188。如果实常数编号从2开始定义就可以避免几何参数和截面类型相混淆的问题，因为缺省的实常数1不存在，也就无法优先了，只能找定义的截面了。复杂的模型建模过程：ET,1,…!定义三种单元类型ET,2,…ET,3,…R,2,…!定义两种实常数R,3,…SECTYPE,1,…!定义两种截面SECTYPE,2,…MP,•••,1,…!定义两种材料MP,-„,2,…N,…!建立节点TYPE,2!指定单元类型2,实常数2,材料1REAL,2

MAT,1!可缺省，系统初始为1E,…MAT,1TYPE,1!指定单元类型1,实常数2,材料2REAL,2MAT,2!可缺省，因为前面指定为2了E,…TYPE,3!指定单元类型3,截面2,材料2REAL!指定材料类型为空，或REAL,0SECNUM,2MAT,2!可缺省，因为前面指定为2了E,…下面以某特大桥为例，用来掌握高级建模的方法。此桥受力模式为桁架拱，计算跨径156米，上下弦为二次抛物线。以桥梁中心为坐标原点的方程如下：下弦拱的抛物线方程:y=y=23（1+三）（1-三）6464xeL78,78〕xeLxeL56,56〕y=27（1+—）（1-—）8484—（x土112）2+3—xeL79,-56〕u156,79〕26433下弦拱承受全部静载及部分活载；上弦和下弦构成桁架，起加劲的作用，以承受局部载重产生的力矩。Finish/clear!清除当前数据库/prep7!前处理

Et,1,Beam188!定义截面梁单元Et,2,Link8!定义索单元R,2,0.005!定义吊杆实常数R,3,0.015,3.74e-3!定义拉杆实常数（初拉力734t）Mp,ex,1,210e9!定义材料弹性模量（钢材Q345qC）Mp,prxy,1,1/3!定义材料泊松比Mp,dens,1,7850*1.1!定义材料密度（考虑节点板附加重）Mp,alpx,1,1.2e-5!定义材料温度伸长率Mp,ex,2,195e9!定义材料弹性模量（钢绞线）Mp,prxy,2,1/3!定义材料泊松比Mp,dens,2,7850!定义材料密度Mp,alpx,2,1.2e-5!定义材料温度伸长率/VIEW,1,1,1,1!改变视角/ANG,1/REP,FASTN,1!建立桥面系节点N,2,3N,3,9N,21,153N,22,159N,23,162Fill,3,21Ngen,2,23,1,23,1,,,0.925Ngen,8,23,24,46,1,,,22.65/7Ngen,2,23,185,207,1,,,0.925N,231,3,-11.164!建立下弦拱节点N,232,9,-6.11N,233,17,0N,234,25,5.39N,235,33,10.062N,236,41,14.015N,237,49,17.25N,238,57,19.765N,239,65,21.562N,240,73,22.64N,241,81,23N,242,89,22.64N,243,97,21.562N,244,105,19.765N,245,113,17.25N,246,121,14.015N,247,129,10.062N,248,137,5.39N,249,145,0N,250,153,-6.11N,251,159,-11.164Ngen,3,21,231,251,1,,,24.5/2!建立上弦拱节点!建立上弦拱节点N,295,9,9.181N,296,17,11.848N,297,25,15N,298,33,18.183N,299,41,20.877N,300,49,23.081N,301,57,24.795N,302,65,26.02N,303,73,26.755N,304,81,27

N,305,89,26.755N,306,97,26.02N,307,105,24.795N,308,113,23.081N,309,121,20.877N,310,129,18.183N,311,137,15N,312,145,11.848N,313,153,9.181N,314,159,7.5Ngen,3,21,294,314,1,,,24.5/2Ngen,2,333,24,185,23,,1!次纵梁方向控制节点SECTYPE,1,BEAM,I,'ZZL',0SECDATA,2,2,1,.01,.01,0.01!主纵梁截面SECTYPE,2,BEAM,T,'CZL',0!次纵梁截面SECOFFSET,USER,0,0.295SECDATA,.4,1.3,.01,.01,.016,.01!OFFSETZ=1.3-0.005-1=0.295SECTYPE,3,BEAM,I,'HL',0SECDATA,.6,.6,2,.04,.04,0.02!横梁截面SECTYPE,4,BEAM,HREC,'G',0SECDATA,1,1,.025,.025,.025,.025!拱截面SECTYPE,5,BEAM,I,'FG',0SECDATA,.6,.6,1,.01,.01,0.01!腹杆截面SECTYPE,6,BEAM,HREC,'JQFG',0SECDATA,.6,1,.02,.02,.02,.02!加强腹杆截面SECTYPE,7,BEAM,I,'GPL',0SECDATA,.4,.4,.4,.01,.01,0.01!拱平联截面/eshape,1!单元外形显示比例为1E,(1:22),(2:23)!主纵梁，单元特性缺省E,(208:229),(209:230)Secnum,2!次纵梁E,(24:45),(25:46),357Egen,8,23,45,66,1Secnum,3!横梁E,(2:186:23),(25:209:23),294Egen,21,1,221,229,1Secnum,4!拱E,(231:250),(232:251)E,(273:292),(274:293)E,(294:313),(295:314)E,(336:355),(337:356)Secnum,5!腹杆E,(296:312),(233:249)E,(296:303),(234:241)E,(305:312),(241:248)Egen,2,42,490,522,1Secnum,6!加强腹杆E,(231:232),(2:3)E,(2:3),(294:295)E,(232:233),(294:295)E,(250:251),(21:22)E,(21:22),(313:314)E,(313:314),(249:250)E,(273:274),(209:210)E,(209:210),(336:337)E,(274:275),(336:337)E,(292:293),(228:229)E,(228:229),(355:356)E,(291:292),(355:356)Secnum,7!拱平联E,(294:315:21),(315:336:21),2!上平连Egen,21,1,580,581,1E,(235:256:21),(256:277:21),6!下平连Egen,13,1,622,623,1E,294,316,253!上K撑E,316,336,253Egen,10,1,648,649,1E,314,334,271E,334,356,271Egen,10,-1,668,669,1E,234,256,319!下K撑E,256,276,319Egen,7,1,688,689,1E,248,268,331E,268,290,331Egen,7,-1,702,703,1Type,2!吊杆Real,2Mat,2E,5,234E,212,276Egen,15,1,716,717,1Real,3!拉杆E,(4:211:207),(20:227:207)Cp,,all,(4:20:16),(233:249:16)!自由度偶合Cp,,all,(211:227:16),(275:291:16)D,(231:273:42),all!边界条件D,(251:293:42),allD,(294:336:42)