某钢桥的ansys命令流
- 格式:doc
- 大小:211.00 KB
- 文档页数:20
碳纤维加固桥梁结构的健康监测摘要:提出了一种针对老化桥梁的健康监测的可靠性方法。
首先, 根据桥梁的类型提出了临界失效标准。
然后基于这些标准定义安全界限。
假定大量的失效标准是正态分布的随机变量。
基于这样的假设,基本可靠性指标和基本失效概率为每个关键失效模式估计。
已经发现基本失效概率值,桥梁系统失效概率值按这段时间考虑计算。
最后,用系统失效概率来获得桥梁系统的可靠性指标。
系统的可靠性指数用来表达的桥梁在使用期间的状况较为片面。
选择斯里兰卡国家铁路桥网络的一座铁路桥梁作为一个案例研究来说明该可靠度程序。
对于这座桥,失效的最重要的形式是疲劳和腐蚀。
应用程序的建议方法表明,桥梁在现有负载的情况下的现状是令人满意的。
对可靠性的预测,不是为了抵御破坏,而是一种对基础设施管理进行有效的积极的健康监测方式。
关键字:碳纤维结构健康监测状态评估维护可靠性指标失效概率1、介绍在当今迅速变化的世界中,人类的生活质量和国家经济发展依赖于基础设施的数量、质量、和有效性]1[。
在世界各国,民用基础设施通常是需要维护的,康复,或更换的]2[。
特别是对基础设施构成一个相当大的投资的高速公路和相关的地面运输系统,直接影响到一个国家的工业产值。
目前大多数桥梁养护策略基于由桥梁检查员在不同时间间隔进行的视觉检查]3[。
由于人类检测取决于个人检查员,其结果存在一定程度的不确定性。
因此,绝对的依赖检验报告是不可靠的,而且可能导致错误的决策和桥梁拥有者的高维护成本。
基于桥梁的剩余使用寿命维护成本的优化。
另一方面,桥梁的检查员发现结构和非结构化缺陷,破坏和恶化过程的能力不应该被忽视。
事实上,人类经验应该被纳入维护决策,但要有足够的主观的知识。
只在那时,任何可接受的战略才能成为可能。
当前桥梁结构健康监测方法存在许多缺点和普遍的不确定性。
对结构健康监测和桥梁维护的知识不足造成由于频繁的维护活动,不可预见的桥梁损害,和偶然的桥梁失效引起的重大的经济损失。
八天学会 Ansys 命令流为方便大家的交流和学习, 特推出 " 跟我学命令流 " 课程本课程分为三部分: 前处理 , 加载求解 , 后处理每部分的学习时间:10 天 , 共计 30 天每天学习大约10 个命令希望本课程对大家能有所帮助第一天目标 : 熟悉 ANSYS基本关键字的含义k --> Keypoints 关键点l --> Lines 线a --> Area 面v --> Volumes 体e --> Elements 单元n --> Nodes 节点cm --> component 组元et --> element type 单元类型mp --> material property 材料属性r --> real constant 实常数d --> DOF constraint约束f --> Force Load 集中力sf --> Surface load on nodes 表面载荷bf --> Body Force on Nodes 体载荷ic --> Initial Conditions 初始条件第二天目标 : 了解命令流的整体结构, 掌握每个模块的标识!文件说明段/BATCH/TITILE,test analysis !定义工作标题/FILENAME,test !定义工作文件名/PREP7 ! 进入前处理模块标识!定义单元 , 材料属性 , 实常数段ET,1,SHELL63 ! 指定单元类型ET,2,SOLID45 !指定体单元MP,EX,1,2E8 !指定弹性模量MP,PRXY,1,0.3 !入泊松比MP,DENS,1,7.8E3 ! 入材料密度R,1,0.001 !指定壳元常数- 厚度......!建立模型K,1,0,0,, !定关点K,2,50,0,,K,3,50,10,,K,4,10,10,,K,5,10,50,,K,6,0,50,,A,1,2,3,4,5,6, !由关点生成面......!划分网格ESIZE,1,0,AMESH,1......FINISH ! 前理束/SOLU ! 入求解模!施加束和荷DL,5,,ALLSFL,3,PRES,1000SFL,2,PRES,1000......SOLVE !求解FINISH ! 求解模束/POST1 ! 入通用后理器....../POST26 ! 入程后理器⋯⋯/EXIT,SAVE ! 退出并存以下是日志文件中常出的一些命令的明, 希望能大家在整理LOG文件有所帮助/ANGLE ! 指定绕轴旋转视图/DIST !说明对视图进行缩放/DEVICE ! 设置图例的显示, 如: 风格 , 字体等/REPLOT ! 重新显示当前图例/RESET ! 恢复缺省的图形设置/VIEW ! 设置观察方向/ZOOM !对图形显示窗口的某一区域进行缩放第三天生成关键点和线部分1.生成关键点K, 关键点编号 ,X 坐标 ,Y 坐标 ,Z 坐标例:K,1,0,0,02.在激活坐标系生成直线LSTR,关键点 P1, 关键点 P2例LSTR,1,23.在两个关键点之间连线L, 关键点 P1, 关键点 P2例L,1,2注: 此命令会随当前的激活坐标系不同而生成直线或弧线4.由三个关键点生成弧线LARC,关键点 P1, 关键点 P2, 关键点 PC,半径 RAD例LARC,1,3,2,0.05注: 关键点 PC是用来控制弧线的凹向5.通过圆心半径生成圆弧CIRCLE,关键点圆心 , 半径 RAD,,,, 圆弧段数NSEG例:CIRCLE,1,0.05,,,,46.通过关键点生成样条线BSPLIN,关键点 P1, 关键点 P2,关键点 P3, 关键点 P4, 关键点 P5, 关键点 P6 例:BSPLIN,1,2,3,4,5,67.生成倒角线LFILLT, 线 NL1, 线 NL2,倒角半径RAD例LFILLT,1,2,0.0058.通过关键点生成面A, 关键点 P1, 关键点 P2, 关键点 P3, 关键点 P4, 关键点 P5, 关键点 P6,P7,P8...例:A,1,2,3,49.通过线生成面AL, 线 L1, 线 L2, 线 L3, 线 L4, 线 L5, 线 L6, 线 L7, 线 L8, 线 L9, 线 L10例:AL,5,6,7,810.通过线的滑移生成面ASKIN, 线 NL1,线 NL2,线 NL3,线 NL4, 线 NL5, 线 NL6, 线 NL7,线 NL8, 线 NL9例:ASKIN,1,4,5,6,7,8注:线 1 为滑移的导向线第四天目标 : 掌握常用的实体 - 面的生成生成矩形面1.通过矩形角上定位点生成面BLC4,定位点 X 方向坐标 XCORNER,定位点 Y 方向坐标 YCORNER,矩形宽度 WIDTH,矩形高度HEIGHT,矩形深度 DEPTH例:BLC4,0,0,5,3,02.通过矩形中心定位点生成面BLC5,定位点 X 方向坐标 XCENTER,定位点 Y 方向坐标 YCENTER,矩形宽度 WIDTH,矩形高度HEIGHT,矩形深度 DEPTH注:与上条命令的不同就在于矩形的定位点不一样例:BLC5,2.5,1.5,5,3,03. 通过在工作平面定义矩形X.Y 坐标生成面RECTNG,矩形左边界 X 坐标 X1,矩形右边界X 坐标 X2, 矩形下边界Y坐标 Y1, 矩形上边界Y 坐标Y2例:RECTNG,0,5,0,3生成圆面4.通过中心定位点生成实心圆面CYL4,定位点 X 方向坐标 XCENTER,定位点 Y 方向坐标 YCENTER,圆面的内半径 RAD1,内圆面旋转角度 THETA1,圆面的外半径 RAD2,外圆面旋转角度THETA2,圆面的深度 DEPTH注:如要实心的圆面则不用 RAD2,THETA2,DEPTH例:CYL4,0,0,5,3605.生成扇形圆面命令介绍如上例 1 实心扇形 :CYL4,0,0,5,60例2 扇形圆环 :CYL4,0,0,5,60,10,60例3 整的圆环 :CYL4,0,0,5,360,10,360注: 同时可通过定义圆面的深度以生成柱体6.通过在工作平面定义起始点生成圆面CYL5,开始点 X坐标 XEDGE1,开始点 Y坐标 YEDGE1,结束点 X坐标 XEDGE2,结束点 Y坐标YEDGE2, 圆面深度 DEPTH例:CYL5,0,0,2,2,7.通过在工作平面定义内外半径和起始角度来生成圆面PCIRC,内半径 RAD1,外半径 RAD2,起始角度THETA1,结束角度THETA2例LCIRC,2,5,30,1808.生成面与面的倒角AFILLT, 面 1 的编号 NA1, 面 2 的编号 NA2,倒角半径RAD例:AFILLT,2,5,2下一讲 : 多边形面的生成第五天目标 : 掌握多边形面和体的生成1.生成多边形面命令 :RPR4, 多边形的边数NSIDES,中心定位点X 坐标 XCENTER,中心定位点Y 坐标 YCENTER, 中心定位点距各边顶点的距离RADIUS,多边形旋转角度THETA例:RPR4,4,0,0,0.15,30注: 这条命令可通过定义不同的NSIDES生成三边形 , 四边形 ,...,八边形2.生成多边形体命令 :RPR4, 多边形的边数NSIDES,中心定位点X 坐标 XCENTER,中心定位点Y 坐标 YCENTER, 中心定位点距各边顶点的距离RADIUS,多边形旋转角度THETA,多边形的深度DEPTH例:RPR4,4,0,0,0.15,30,0.1注: 多边形体和面命令唯一的不同就在于深度DEPTH的定义到此 , 关键点 , 线 , 面的生成讲解已结束, 下一讲 : 体的生成第六天目标 : 掌握体的生成命令1.通过关键点生成体命令 :V, 关键点 P1, 关键点 P2, P3, P4, P5, P6, P7, P8例:V,4,5,6,7,15,24,252.通过面生成体命令 :VA, 面 A1, 面 A2, A3, A4, A5, A6, A7, A8, A9, A10例:VA,3,4,5,8,103.通过长方形角上定位点生成体命令 :BLC4该命令前面在讲生成面的时候已作介绍, 唯一的不同在于深度DEPTH的定义 .4.通过长方形中心定位点生成面命令 :BLC55.通过定义长方体起始位置生成体命令 :BLOCK,开始点 X 坐标 X1,结束点 X 坐标 X2, Y1, Y2, Z1, Z2例:BLOCK,2,5,0,2,1,36.生成圆柱体基本命令通生成圆形面, 不同在于DEPTH的定义基本命令 :CYL4基本命令 :CYL5基本命令 :CYLIND7.生成棱柱基本命令通生成多边形, 不同在于DEPTH的定义基本命令 :RPR48.通过球心半径生成球体命令 :SPH4, 球心 X 坐标 XCENTER,球心 Y 坐标 YCENTER,半径 RAD1,半径 RAD2例:SPH4,1,1,2,59.通过直径上起始点坐标生成球体命令 :SPH5, 起点 X 坐标 XEDGE1,起点 Y坐标 YEDGE1,结束点 X 坐标 XEDGE2,结束点 Y 坐标YEDGE2例:SPH5,2,5,7,610.在工作平面起点通过半径和转动角度生成球体命令 :SPHERE,半径 RAD1,半径 RAD2,转动角度THETA1,转动角度THETA2例:SPHERE,2,5,0,6011.生成圆锥体命令 :CONE,底面半径 RBOT,顶面半径 RTOP,底面高 Z1, 顶面高 Z2, 转动角度 THETA1,转动角度THETA2例:CONE,10,20,0,50,0,180下一讲 : 布尔操作第七天目标 : 掌握常用的布尔操作命令1.沿法向延伸面生成体命令 :VOFFST,面的编号NAREA,面拉伸的长度DIST, 关键点增量KINC例:VOFFST,1,2,,2.通过坐标的增量延伸面生成体命令 :VEXT, 面 1 的编号 NA1,面 2 的编号 NA2,增量 NINC,X 方向的增量 DX,Y 方向的增量 DY,Z 方向的增量 DZ, RX, RY, RZ例:VEXT,1,5,1,1,2,2,3. 面绕轴旋转生成体命令 :VROTAT,面 1 的编号PAX1,定位轴关键点 2 编号NA1,面 2 的编号PAX2,旋转角度NA2,NA3, NA4, NA5, NA6,ARC,生成体的段数 NSEG定位轴关键点 1 编号例:VROTAT,1,2,,,,,4,5,360,44. 沿线延伸面生成体命令 :VDRAG,面 1 的编号 NA1,面 2 的编号 NA2, NA3, NA4, NA5, NA6,导引线 2 的编号 NLP2, NLP3, NLP4, NLP5, NLP6例:VDRAG,2,3,,,,,8,导引线 1 的编号NLP1, 5. 线绕轴旋转生成面命令 :AROTAT,线 1 的编号 NL1, NL2, NL3, NL4, NL5, NL6, 位轴关键点 2 的编号 PAX2,旋转角度ARC,生成面的段数例:AROTAT,3,4,,,,,6,8,360,4定位轴关键点NSEG1 的编号PAX1,定6.沿线延伸线生成面命令 :ADRAG,线 1 的编号 NL1,NL2, NL3, NL4, NL5, NL6, 导引线 1 的编号 NLP1, NLP2, NLP3, NLP4, NLP5, NLP6例:ADRAG,3,,,,,,87.同理可以延伸关键点 , 相应的命令如下 :LROTAT, NK1, NK2, NK3, NK4, NK5, NK6, PAX1, PAX2, ARC, NSEGLDRAG, NK1, NK2, NK3, NK4, NK5, NK6, NL1, NL2, NL3, NL4, NL5, NL6各选项的含义雷同于上 .8.延伸一条线命令 LEXTND,线的编号 NL1, 定位关键点编号 NK1,延伸的距离 DIST, 原有线是否保留控制项KEEP例LEXTND,5,2,1.5,09.布尔操作 : 加命令 LCOMB,线编号 NL1,线编号 NL2, 是否修改控制项KEEP例LCOMB,2,5注: 对面和体的相应为:VADD,AADD.选项的含义都类似10.布尔操作 : 粘接和搭接搭接的核心关键字为:OVLAP,随实体的不同略有不同, 如 :对体为 VOVLAP对面为 AOVLAP对线为 LOVLAP粘接的核心关键字为:GLUE, 随实体的不同略有不同, 如 :对体为 VGLUE对面为 AGLUE对线为 LGLUE但其他的选项的含义是类似的, 这里就不再累述.下一讲 : 移动 , 复制 , 映射 , 删除 ...第八天目标 : 掌握体素的移动, 复制 , 删除 , 映射一. 移动关键点命令 :KMODIF,关键点编号NPT,移动后的坐标X, 移动后的坐标Y, 移动后的坐标Z例:KMODIF,5,0,0,2二. 移动复制关键点命令 :KGEN,复制次数选项 ITIME, 起始关键点编号 NP1,结束关键点编号 NP2,增量 NINC,偏移DX,偏移 DY,偏移 DZ,关键点编号增量 KINC, 生成节点单元控制项 NOELEM,原关键点是否被修改选项 IMOVE例:KGEN,2,1,10,1,2,2,2,,,,注:IMOVE 选项说明 , 设置为 0 时, 不修改原关键点 , 即为复制 , 设置为 1 时 , 修改原关键点 , 即为移动 , 从而通过控制 IMOVE选项实现移动或复制 .三. 移动复制线命LGEN,ITIME,NL1,NL2,NINC,DX,DY,DZ,KINC,NOELEM,IMOVE各选项的含义同上四. 移动复制面命:AGEN,ITIME,NA1,NA2,NINC,DX,DY,DZ,KINC,NOELEM,IMOVE各选项的含义同上五. 移动复制体命令 :VGEN,ITIME,NV1,NV2,NINC,DX,DY,DZ,KINC,NOELEM,IMOVE各选项的含义同上六. 修改面的法向方向命令 :ANORM,面的编号 ANUM,单元的法向方向是否修改选项NOEFLIP例:ANORM,2七.体素的删除基本的命令为 :*DELE组合不同的关键字形成不同的命令如:KDELE,LDELE,ADELE,VDELE基本的命令格式为 :*DELE, 起始体素编号 N*1, 结束体素编号 N*2, 增量 NINC,是否删除体素下层的元素选项 KSWP如LDELE,2,5,1,1八.体素的映射基本的命令为 :*SYMM组合不同的关键字形成不同的命令如:KSYMM,LSYMM,ARSYM,VSYMM基本的命令格式为 :*SYMM,映射轴选项 NCOMP,起始体素编号 N*1, 结束体素编号 N*2, 增量NINC,关键点编号增量 KINC,NOELEM, IMOVE如:VSYMM,X,1,10,1,,,,。
ansys命令流使用方法
在ANSYS中,命令流是一种用于执行特定操作的自动化工具。
以下是ANSYS命令流使用的一般步骤:
1. 打开ANSYS软件并加载您要使用的工程文件。
2. 在ANSYS Graphical User Interface (GUI) 中,将鼠标指针放
在工具栏上。
在“Run”下拉菜单中选择“Command Line”。
3. 在命令行窗口中,输入和编辑您想执行的命令。
您可以使用ANSYS的命令语言以及相关命令进行模型操作、网格生成、
求解等。
4. 您可以通过多种方式输入命令:直接在命令行中输入、从脚本文件中读取、从ANSYS GUI中的日志文件中复制粘贴等。
5. 您可以使用命令流中的参数和变量来进行自动化操作。
使用“!VARIABLE”语句定义变量,并通过“!VARIABLE = value”语
句赋值。
6. 使用ANSYS的各种功能命令对模型进行操作。
例如,在预
处理阶段,您可以使用命令生成几何体、定义材料属性、设定网格、添加边界条件等。
7. 在求解阶段,使用命令启动求解器,设置求解器选项,运行求解器,并监视求解器的输出。
8. 在结果后处理阶段,使用命令读取并处理结果数据,生成图形、报告等。
9. 执行命令流,您可以一次性执行整个命令流,或者逐个执行命令。
10. 您还可以将命令流保存为脚本文件,以便将来再次使用。
以上是ANSYS命令流的一般用法,具体的命令和语法取决于您的特定需求和ANSYS的版本。
建议您参考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"的缩略语,长度不超过8个字符.String:将由"Abbr"表示的字符串,长度不超过60个字符.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最常用命令流+中文注释VSBV, NV1, NV2, SEPO, KEEP1, KEEP2 —Subtracts volumes from volumes,用于2个solid相减操作,最终目的是要nv1-nv2=?通过后面的参数设置,可以得到很多种情况:sepo项是2个体的边界情况,当缺省的时候,是表示2个体相减后,其边界是公用的,当为sepo的时候,表示相减后,2个体有各自的独立边界。
keep1与keep2是询问相减后,保留哪个体?当第一个为keep时,保留nv1,都缺省的时候,操作结果最终只有一个体,比如:vsbv,1,2,sepo,,keep,表示执行1-2的操作,结果是保留体2,体1被删除,还有一个1-2的结果体,现在一共是2个体(即1-2与2),且都各自有自己的边界。
如vsbv,1,2,,keep,,则为1-2后,剩下体1和体1-2,且2个体在边界处公用。
同理,将v换成a 及l是对面和线进行减操作!mp,lab, mat, co, c1,…….c4 定义材料号及特性lab: 待定义的特性项目(ex,alpx,reft,prxy,nuxy,gxy,mu,dens) ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号(缺省为当前材料号)co: 材料特性值,或材料之特性,温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项,2次项,3次项,4次项的系数定义DP材料:首先要定义EX和泊松比:MP,EX,MA T,……MP,NUXY,MAT,……定义DP材料单元表(这里不考虑温度):TB,DP,MA T进入单元表并编辑添加单元表:TBDATA,1,CTBDATA,2,ψTBDATA,3,……如定义:EX=1E8,NUXY=0.3,C=27,ψ=45的命令如下:MP,EX,1,1E8MP,NUXY,1,0.3TB,DP,1TBDATA,1,27TBDATA,2,45这里要注意的是,在前处理的最初,要将角度单位转化到“度”,即命令:*afun,degVSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Type,是选择的方式,有选择(s),补选(a),不选(u),全选(all)、反选(inv)等,其余方式不常用Item, Comp 是选取的原则以及下面的子项如volu 就是根据实体编号选择,loc 就是根据坐标选取,它的comp就可以是实体的某方向坐标!其余还有材料类型、实常数等MIN, VMAX, VINC,这个就不必说了吧!,例:vsel,s,volu,,14vsel,a,volu,,17,23,2上面的命令选中了实体编号为14,17,19,21,23的五个实体VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体nv1:初始体号nv2:最终的体号ninc:体号之间的间隔kswp=0:只删除体kswp=1:删除体及组成关键点,线面如果nv1=all,则nv2,ninc不起作用其后面常常跟着一条显示命令VPLO,或aplo,nplo,这个湿没有参数的命令,输入后直接回车,就可以显示刚刚选择了的体、面或节点,很实用的哦!Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点(缺省)R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0” 使用正负号“1”仅用绝对值下面是单元生死第一个载荷步中命令输入示例:!第一个载荷步TIME,... !设定时间值(静力分析选项)NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿-拉夫森选项ESTIF,... !设定非缺省缩减因子(可选)ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点(不与活动单元相连的结点)D,ALL,ALL,0 !约束所有不活动的结点自由度(可选)NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SA VESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D, F,SF和BF命令得到更详细的解释。
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"的缩略语,长度不超过8个字符.String:将由"Abbr"表示的字符串,长度不超过60个字符.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为半径。
9、*AFUN,Lab在参数表达式中,为角度函数指定单位.Lab:指定将要使用的角度单位.有3个选项.RAD:在角度函数的输入与输出中使用弧度单位(默认)DEG:在角度函数的输入与输出中使用度单位.STAT:显示该命令当前的设置(即是度还是弧度).10、Agen, itime,na1,na2,ninc,dx,dy,dz,kinc,noelem,imove !面积复制命令。
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, VALUELab为参数名称 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杆单元单元名称简称节点数节点自由度特性备注LINK1 2D杆 2 Ux,Uy EPCSDGB常用杆元LINK8 3D杆Ux,Uy,Uz EPCSDGBLINK103D仅受拉或仅受压杆EDGB模拟缆索的松弛及间隙LINK11 3D线性调节器EGB模拟液压缸和大转动LINK183D有限应EPCDFG另可考虑粘0 变杆 B 弹塑性E-弹性(Elasticity),P-塑性(Plasticity),C-蠕变(Creep),S-膨胀(Swelling),D-大变形或大挠度(Large deflection),F-大应变(Large strain)或有限应变(Finite strain),B-单元生死(Birth and dead),G-应力刚化(Stress stiffness)或几何刚度(Geometric stiffening),A-自适应下降(Adaptive descent)等。
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命令流大全(整理)1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9此命令用已知的一组关键点点(P1~P9)来定义面(Area), 最少使用三个点才能围成面,同时产生转围绕些面的线。
点要依次序输入,输入的顺序会决定面的法线方向。
如果超过四个点,则这些点必须在同一个平面上。
MenuPaths:MainMenu>Preprocessor>Create>ArbitrarP>Th roughK Ps 2、GABBR,Abbr,String ――定义一个缩略语.Abbr:用来表示字符串"String "的缩略语,长度不超过8个字符. String :将由"Abbr "表示的字符串,长度不超过6 0个字符.3、ABBRES,Lab,Fname,EGt —从一个编码文件中读出缩略语. Lab:指定读操作的标题,NEW :用这些读出的缩略语重新取代当前的缩略语(默认) CHANGE :将读岀的缩略语添加到当前缩略语阵列,并替代现存同名的缩略语.EGt:如果"Fname"是空的,则缺省的扩展命是" ABBR ".4、ABBSAV,Lab,Fname,EGt —将当前的缩略语写入一个文本文件里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=O 时只删除掉面积本身,= 1 时低单元点一并删除。
7、Adrag,nl1,nl2,nl3,nl4,nl5,nl6,nlp1,nlp2,nlp3,nlp4,nlp5,nl p6 !面积的建立,沿某组线段路径,拉伸而成。
Finish(退出四大模块,回到BEGIN层)/clear (清空内存,开始新的计算)1.定义参数、数组,并赋值.2./prep7(进入前处理)定义几何图形:关键点、线、面、体定义几个所关心的节点,以备后处理时调用节点号。
设材料线弹性、非线性特性设置单元类型及相应KEYOPT设置实常数设置网格划分,划分网格根据需要耦合某些节点自由度定义单元表3./solu加边界条件设置求解选项定义载荷步求解载荷步4./post1(通用后处理)5./post26 (时间历程后处理)6.PLOTCONTROL菜单命令7.参数化设计语言8.理论手册Finish(退出四大模块,回到BEGIN层)/clear (清空内存,开始新的计算)1.定义参数、数组,并赋值.dim, par, type, imax, jmax, kmax, var1, vae2, var3 定义数组par: 数组名type:array 数组,如同fortran,下标最小号为1,可以多达三维(缺省)char 字符串组(每个元素最多8个字符)tableimax,jmax, kmax 各维的最大下标号var1,var2,var3 各维变量名,缺省为row,column,plane(当type为table时) 2./prep7(进入前处理)2.1 设置单元类型及相应KEYOPTET, itype, ename, kop1……kop6, inopr 设定当前单元类型Itype:单元号Ename:单元名设置实常数Keyopt, itype, knum, valueitype: 已定义的单元类型号knum: 单元的关键字号value: 数值注意:如果,则必须使用keyopt命令,否则也可在ET命令中输入2.2 定义几个所关心的节点,以备后处理时调用节点号。
n,node,x,y,z,thxy, thyz, thzx 根据坐标定义节点号如果已有此节点,则原节点被重新定义,一般为最大节点号。
1模型简介模型以壳单元shell63为主,横梁采用beam44单元,吊杆采用link10单元,横梁、吊杆与拱肋的连接采用节点耦合的方式。
不考虑纵坡和横披。
111213141516X Y ZROTCPACEL图1 三维图111213141516X Y Z 图2 正立面图2建模命令流finish/clear/config,nproc,2/prep7!双核运行!材料!钢材EX,1,2.06e11PRXY ,1,0.3DENS,1,7800MPDA TA,ALPX,1,,1.2e-5!吊杆钢材EX,2,2.06e11PRXY ,2,0.3DENS,2,7800MPDA TA,ALPX,2,,1.2e-5!桥面钢材包括铺装层重量EX,3,2.06e11PRXY,3,0.3DENS,3,7800*3MPDA TA,ALPX,3,,1.2e-5!单元ET,1,plane82ET,2,shell63!设置shell63关键字KEYOPT,2,1,0KEYOPT,2,2,0KEYOPT,2,3,0KEYOPT,2,5,2KEYOPT,2,6,0KEYOPT,2,7,0KEYOPT,2,8,0KEYOPT,2,9,0KEYOPT,2,11,0et,3,beam189!设置beam189关键字KEYOPT,3,1,1KEYOPT,3,2,0KEYOPT,3,4,2KEYOPT,3,6,3KEYOPT,3,7,2KEYOPT,3,8,3KEYOPT,3,9,3KEYOPT,3,10,0KEYOPT,3,11,0KEYOPT,3,12,0ET,4,solid45et,5,link10r,1,0.010r,2,0.016r,3,0.02r,4,0.04!建立各种截面!钢横梁!HL1/HL3a=1.0b=0.680th=0.02ts=0.016RECTNG,-a/2,a/2,-b/2,b/2, RECTNG,-a/2+ts,a/2-ts,-b/2+th,b/2-th,RECTNG,-a/2+ts,-a/2+ts+0.16,-0.01,0.01, RECTNG,a/2-ts,a/2-ts-0.16,-0.01,0.01, RECTNG,-0.01,0.01,-b/2+0.02,-b/2+0.02+0.16, RECTNG,-0.01,0.01,b/2-0.02-0.16,b/2-0.02, AGLUE,allwprot,0,90,0wpoff,0,0,-b/2+thASBW,allwpoff,0,0,b/2-th-0.01ASBW,allwpoff,0,0,0.02ASBW,allwpoff,0,0,b/2-th-0.01ASBW,allwprot,0,0,90wpoff,0,0,-a/2+tsASBW,allwpoff,0,0,a-2*tsASBW,allwpoff,0,0,-a/2+ts+0.01ASBW,allwpoff,0,0,-0.02ASBW,allESIZE,0.3,0,amesh,allSECWRITE,sec3SECTYPE,3,BEAM,MESH, SECOFFSET,CENT,,,SECREAD,'sec3','SECT',' ',MESHaclear,allADELE,all, , ,1WPCSYS,-1,0!HL2!吊杆SECTYPE, 4, BEAM, CTUBE, , 0 SECOFFSET, CENTSECDA TA,0.032,0.04,,0,0,0,0,0,0,0!建立跨中截面所有线k,1,0,-2.058,0+5.45k,2,0,-2,0+5.45k,,0,-1.6,0+5.45k,,0,-0.8,0+5.45k,,0,-0.4,0+5.45k,,0,-1.6,0.108+5.45k,,0,-1.2,0.108+5.45k,,0,-0.8,0.108+5.45k,,0,-0.4,0.108+5.45k,,0,-2.541661,0.200339+5.45 k,,0,-2.37,0.576+5.45k,,0,-1.6,0.576+5.45k,,0,-1.2,0.576+5.45k,,0,-0.8,0.576+5.45k,,0,-0.4,0.576+5.45k,,0,-2.742,0.684+5.45k,,0,-2.37,0.684+5.45k,,0,-2,0.684+5.45k,,0,-1.6,0.684+5.45k,,0,-1.2,0.684+5.45k,,0,-0.8,0.684+5.45k,,0,-0.4,0.684+5.45k,,0,0,0.684+5.45k,,0,0,0+5.45*do,i,1,5,1l,i,i+1*enddol,6,25*do,i,17,24,1l,i,i+1*enddo*do,i,3,6,1l,i,i+4*enddo*do,i,20,23,1l,i,i-7*enddol,18,12larc,1,17,11l,2,19!建立拖拉的路径线k,26,-40,0,0+5.45l,25,26!拖拉成箱梁*do,i,1,25,1ADRAG,i, , , , , ,26*enddoallselNUMMRG,KP, , , ,LOWwprot,0,0,90*do,i,1,15,1wpoff,0,0,-2.5ASBW,all*enddo!通过循环生成横隔板,选择封闭的为轮廓线来生成!生成矩形部分*do,i,0,16,1lsel,s,loc,x,-2.5*i-0.001,-2.5*i+0.001lsel,r,loc,y,-2.001,0.001lsel,u,loc,z,0.001+5.45,0.2+5.45lsel,u,loc,z,0.683+5.45,0.6+5.45al,all*enddo!生成扇形部分*do,i,0,16,1lsel,s,loc,x,-2.5*i-0.001,-2.5*i+0.001lsel,r,loc,y,-1.999,-2.749lsel,u,loc,z,0.683+5.45,0.6+5.45al,all*enddoallsel!多次分割WPCSYS,-1,0asel,u,loc,z,0.001+5.45,0.2+5.45asel,u,loc,z,0.683+5.45,0.6+5.45wprot,0,90,0*do,i,1,6,1KWPA VE,iasbw,all*enddoKWPA VE,18asbw,allallselWPA VE,0,0,0wpoff,0,0,0.2asbw,allwpoff,0,0,1.58asbw,allallselWPCSYS,-1,0KWPA VE,10asbw,allKWPA VE,16asbw,allWPCSYS,-1,0wpoff,0,0,0.2+5.45asbw,allwpoff,0,0,0.268asbw,all!删除需要空心的地方*do,j,0,3*do,i,1,3,1asel,s,loc,x,-2.5*(i+4*j)-0.001,-2.5*(i+4*j)+0.001 asel,r,loc,z,0.2+5.45,0.2+0.268+5.45asel,r,loc,y,-0.2,-0.2-1.580ADELE,all, , ,1*enddo*enddoallselaglue,allNUMMRG,KP, , , ,LOWNUMCMP,KPNUMCMP,LINENUMCMP,AREAWPCSYS,-1,0!建立壳单元拱肋!上拱肋跨中断面关键点k,10100,0,-3.25-0.5,-0.35+10.1k,10101,0,-3.25-0.5,0+10.1k,10102,0,-3.25-0.5,0.35+10.1k,10103,0,-3.25-0.34,0+10.1k,10104,0,-3.25-0.25,-0.35+10.1k,10105,0,-3.25-0.25,-0.19+10.1k,10106,0,-3.25-0.25,0.19+10.1k,10107,0,-3.25-0.25,0.35+10.1k,10108,0,-3.25-0.16,0+10.1k,10109,0,-3.25+0,-0.35+10.1k,10110,0,-3.25+0,0+10.1k,10111,0,-3.25+0,0.35+10.1k,10112,0,-3.25+0.16,0+10.1k,10113,0,-3.25+0.25,-0.35+10.1k,10114,0,-3.25+0.25,-0.19+10.1k,10115,0,-3.25+0.25,0.19+10.1 k,10116,0,-3.25+0.25,0.35+10.1 k,10117,0,-3.25+0.34,0+10.1 k,10118,0,-3.25+0.5,-0.35+10.1 k,10119,0,-3.25+0.5,0+10.1 k,10120,0,-3.25+0.5,0.35+10.1 !下拱肋跨中断面关键点k,10200,0,-3.25-0.5,-0.35+5 k,10201,0,-3.25-0.5,0+5k,10202,0,-3.25-0.5,0.35+5k,10203,0,-3.25-0.34,0+5k,10204,0,-3.25+0,-0.35+5k,10205,0,-3.25+0,-0.19+5k,10206,0,-3.25+0,0.19+5k,10207,0,-3.25+0,0.35+5k,10208,0,-3.25+0.34,0+5k,10209,0,-3.25+0.5,-0.35+5 k,10210,0,-3.25+0.5,0+5k,10211,0,-3.25+0.5,0.35+5!连成线NUMCMP,LINELSTR, 10118, 10113 LSTR, 10113, 10109 LSTR, 10109, 10104 LSTR, 10104, 10100 LSTR, 10100, 10101 LSTR, 10101, 10102 LSTR, 10102, 10107 LSTR, 10107, 10111 LSTR, 10111, 10116 LSTR, 10116, 10120 LSTR, 10120, 10119 LSTR, 10119, 10118 LSTR, 10111, 10110 LSTR, 10110, 10109 LSTR, 10116, 10115 LSTR, 10113, 10114 LSTR, 10119, 10117 LSTR, 10110, 10112 LSTR, 10107, 10106 LSTR, 10104, 10105 LSTR, 10101, 10103 LSTR, 10110, 10108LSTR, 10209, 10204LSTR, 10204, 10200LSTR, 10200, 10201LSTR, 10201, 10202LSTR, 10202, 10207LSTR, 10207, 10211LSTR, 10211, 10210LSTR, 10210, 10209LSTR, 10210, 10208LSTR, 10201, 10203LSTR, 10204, 10205LSTR, 10207, 10206!通过旋转生成面!建立上拱肋圆心关键点k,10300,0,0,-27.141k,10301,0,-3,-27.141!建立下拱肋圆心关键点k,10302,0,0,-67.525+5k,10303,0,-3,-67.525+5 WPCSYS,-1,0!上拱肋*do,i,2499,2520,1AROTA T, i,,,,,,10300,10301,43.21447, *enddo!下拱肋*do,i,2521,2532,1AROTA T, i,,,,,,10302,10303,22.18745, *enddo!拱脚的处理ASEL,S,LOC,Y,-3+0.251,-10 ALLSEL,BELOW,AREA NUMMRG,KP, , , ,LOW NUMCMP,KPNUMCMP,LINENUMCMP,AREANUMMRG,KP, , , SELE,LOW ASEL,R,LOC,Z,0,4.5wprot,0,90,0wpoff,0,0,3.00asbw,allwpoff,0,0,0.50asbw,allASEL,S,LOC,Y,-3+0.251,-10WPCSYS,-1,0wprot,0,0,90*do,i,1,5,1wpoff,0,0,-3asbw,all*enddoWPA VE,0,0,0*do,i,1,2,1wpoff,0,0,-10asbw,all*enddoasel,r,loc,x,-25,-20AOVLAP,1593,1625KWPLAN,-1, 1441, 1440, 1112ASBW,allasel,r,loc,x,-25,-23.5ADELE,all, , ,1ASEL,S,LOC,Y,-3+0.251,-10ALLSEL,BELOW,AREANUMMRG,KP, , , SELE,LOWWPCSYS,-1,0!拱脚加劲板LSTR, 1409, 1431al,2710,2703,2533LSTR, 1432, 1405al,3177,2536,2992!上拱肋横隔板*do,i,0,5,1lSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5.001+0.35,10.101+0.35lsel,r,loc,x,-3*i-0.001,-3*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-3*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-3*i)**2)** 0.5)+0.05lsel,r,loc,y,-2.749,-3.251al,alllSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5.001+0.35,10.101+0.35lsel,r,loc,x,-3*i-0.001,-3*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-3*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-3*i)**2)** 0.5)+0.05lsel,r,loc,y,-3.249,-4al,all*enddoi=1lSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5.001,10.101+0.35lsel,r,loc,x,-10*i-0.001,-10*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-10*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-10*i)**2) **0.5)+0.05lsel,r,loc,y,-2.749,-3.251al,alllSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5.001,10.101+0.35lsel,r,loc,x,-10*i-0.001,-10*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-10*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-10*i)**2) **0.5)+0.05lsel,r,loc,y,-3.249,-4al,alli=2lSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,3.8,10.101+0.35lsel,r,loc,x,-10*i-0.001,-10*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-10*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-10*i)**2) **0.5)+0.05lsel,r,loc,y,-2.749,-3.251al,alllSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,3.8,10.101+0.35lsel,r,loc,x,-10*i-0.001,-10*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-10*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-10*i)**2) **0.5)+0.05lsel,r,loc,y,-3.249,-4al,all!下拱肋横隔板*do,i,0,5,1lSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5-(67.525-(67.525**2-(-3*i)**2)**0.5)-0.6,5-(67.525-(67.525**2-(-3*i)**2)**0.5)+0. 6lsel,r,loc,x,-3*i-0.001,-3*i+0.001lsel,u,loc,y,-3.251,-3.249lsel,u,loc,z,5-(67.525-(67.525**2-(-3*i)**2)**0.5)-0.05,5-(67.525-(67.525**2-(-3*i)**2)**0.5)+ 0.05al,all*enddo*do,i,10,20,10lSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5-(67.525-(67.525**2-(-i)**2)**0.5)-0.6,5-(67.525-(67.525**2-(-i)**2)**0.5)+0.6 lsel,r,loc,x,-i-0.001,-i+0.001lsel,u,loc,y,-3.251,-3.249lsel,u,loc,z,5-(67.525-(67.525**2-(-i)**2)**0.5)-0.05,5-(67.525-(67.525**2-(-i)**2)**0.5)+0.05 al,all*enddoallsel!删除上拱肋多余部分asel,s,loc,y,-3.251,-10asel,r,loc,x,-12,-50KWPA VE, 1411asbw,allasel,r,loc,z,kz(997),20ADELE,all, , ,1asel,s,loc,y,-3.251,-10asel,r,loc,z,5.5,20ADELE,all, , ,1asel,s,loc,y,-2.749,-10NUMMRG,KP, , , SELE,LOW!建立新吊杆线LSTR, 1174, 1152LSTR, 1215, 1190LSTR, 1252, 1229LSTR, 1289, 1266LSTR, 1326, 1303LSTR, 1363, 1340FLST,5,6,4,ORDE,2FITEM,5,2501FITEM,5,-2506LSEL,R, , ,P51Xcm,diaogan,lineallsel!建立横梁线allselk,10001,0,0,kz(1158)k,10002,0,-2.73,kz(1158) k,10003,-10,0,kz(1387) k,10004,-10,-2.73,kz(1387) k,10005,-20,0,kz(1410) k,10006,-20,-2.73,kz(1410)l,10001,10002l,10003,10004l,10005,10006lsel,s,loc,z,-10,5.01lsel,r,loc,y,0,-2.73 WPCSYS,-1,0KWPA VE, 10001 wprot,0,90,0wpoff,0,0,1.6lsbw,allcm,hengliang,line!划分!划分横梁type,3mat,1real,1SECNUM,3esize,0.5lmesh,hengliang!划分吊杆allselTYPE,5mat,2real,1SECNUM,4esize,20lmesh,diaoganallsel!划分横隔板type,2mat,1*do,i,0,40,2.5asel,s,loc,x,-i+0.001,-i-0.001 asel,r,loc,y,0,-2.75 ALLSEL,BELOW,AREA LESIZE,all,0.25, , , , , , ,1 aclear,allamesh,all*enddo!划分纵隔板allselLESIZE,all, , ,-1, ,1asel,s,loc,y,-0.001,0.001 asel,a,loc,y,-2.001,-1.999 ALLSEL,BELOW,AREA LESIZE,all,0.4, , , , , , ,1 type,2mat,1real,3amesh,all!划分加劲肋asel,s,loc,y,-2.370,-2.370*do,i,0.4,1.6,0.4asel,a,loc,y,-i-0.001,-i+0.001 *enddoALLSEL,BELOW,AREA LESIZE,all, , ,-1, ,1 LESIZE,all,0.4, , , , , , ,1 type,2mat,1real,1amesh,allallsel!划分箱体外包面ASEL,S,MA T,,0asel,r,loc,y,0.001,-2.75+0.001 LESIZE,all, , ,-1, ,1 LESIZE,all,0.4, , , , , , ,1 real,2amesh,allallsel!拱肋划分allselLESIZE,all, , ,-1, ,1mat,1real,3esize,0.5asel,s,loc,y,-2.749,-10 amesh,allallsel!镜像生成全桥LSEL,S,TYPE,,3LSEL,a,TYPE,,5 LSYMM,Y,all, , , ,0,0 asel,s,loc,y,-0.001,-10 ARSYM,Y,all, , , ,0,0LSEL,S,TYPE,,3LSEL,a,TYPE,,5lsel,r,loc,x,-50,-0.001 LSYMM,x,all, , , ,0,0 asel,s,loc,x,-0.001,-50 ARSYM,x,all, , , ,0,0allselNUMMRG,NODE, , , ,LOW NUMMRG,KP, , , ,LOW NUMMRG,NODE,,, ,LOW !统一所有单元坐标allselcsys,0EMODIF,P51X,ESYS, ,!耦合横梁和箱梁!跨中nsel,s,loc,x,-0.001,0.001 nsel,r,loc,y,-1.601,-1.599n1=node(0,-1.60,0)nsel,u,,,n1n2=node(0,-1.60,0)nsel,a,,,n1cp,1,ux,n1,n2cp,next,uy,n1,n2cp,next,uz,n1,n2nsel,s,loc,x,-0.001,0.001 nsel,r,loc,y,1.601,1.599n1=node(0,1.60,0)nsel,u,,,n1n2=node(0,1.60,0)nsel,a,,,n1cp,next,ux,n1,n2cp,next,uz,n1,n2!H2横梁处nsel,s,loc,x,-10.001,-9.999 nsel,r,loc,y,-1.601,-1.599n1=node(0,-1.60,0) nsel,u,,,n1n2=node(0,-1.60,0) nsel,a,,,n1cp,next,uy,n1,n2cp,next,uz,n1,n2nsel,s,loc,x,-10.001,-9.999 nsel,r,loc,y,1.601,1.599n1=node(0,1.60,0)nsel,u,,,n1n2=node(0,1.60,0)nsel,a,,,n1cp,next,uz,n1,n2!H2'横梁处nsel,s,loc,x,10.001,9.999 nsel,r,loc,y,-1.601,-1.599n1=node(0,-1.60,0) nsel,u,,,n1n2=node(0,-1.60,0) nsel,a,,,n1cp,next,uy,n1,n2cp,next,uz,n1,n2nsel,s,loc,x,10.001,9.999 nsel,r,loc,y,1.601,1.599n1=node(0,1.60,0)nsel,u,,,n1n2=node(0,1.60,0)nsel,a,,,n1cp,next,uz,n1,n2!H1横梁处nsel,s,loc,x,-20.001,-19.999 nsel,r,loc,y,-1.601,-1.599n1=node(0,-1.60,0) nsel,u,,,n1n2=node(0,-1.60,0) nsel,a,,,n1cp,next,uy,n1,n2cp,next,uz,n1,n2nsel,s,loc,x,-20.001,-19.999 nsel,r,loc,y,1.601,1.599n1=node(0,1.60,0)nsel,u,,,n1n2=node(0,1.60,0)nsel,a,,,n1cp,next,uz,n1,n2!H1'横梁处nsel,s,loc,x,20.001,19.999 nsel,r,loc,y,-1.601,-1.599n1=node(0,-1.60,0) nsel,u,,,n1n2=node(0,-1.60,0) nsel,a,,,n1cp,next,uy,n1,n2cp,next,uz,n1,n2nsel,s,loc,x,20.001,19.999 nsel,r,loc,y,1.601,1.599n1=node(0,1.60,0)nsel,u,,,n1n2=node(0,1.60,0)nsel,a,,,n1cp,next,uz,n1,n2allsel!耦合横梁与拱肋!跨中横梁n1=node(0,-2.72,5)n2=node(-0.5,-2.75,5)n3=node(0.5,-2.75,5)n4=node(0,-2.75,5)n5=node(0,-2.75,5+0.35) cp,next,all,n1,n4n1=node(0,2.72,5)n2=node(-0.5,2.75,5)n3=node(0.5,2.75,5)n4=node(0,2.75,5)n5=node(0,2.75,5+0.35) cp,next,all,n1,n4!H2n1=node(-10,-2.72,4.255) n2=node(-10.5,-2.75,4.255) n3=node(-9.5,-2.75,4.255)n4=node(-10,-2.75,4.255)n5=node(-10,-2.75,4.255+0.35) cp,next,all,n1,n4n1=node(-10,2.72,4.255)n2=node(-10.5,2.75,4.255)n3=node(-9.5,2.75,4.255)n4=node(-10,2.75,4.255)n5=node(-10,2.75,4.255+0.35) cp,next,all,n1,n4n1=node(10,-2.72,4.255)n2=node(10.5,-2.75,4.255)n3=node(9.5,-2.75,4.255)n4=node(10,-2.75,4.255)n5=node(10,-2.75,4.255+0.35) cp,next,all,n1,n4n1=node(10,2.72,4.255)n2=node(10.5,2.75,4.255)n3=node(9.5,2.75,4.255)n4=node(10,2.75,4.255)n5=node(10,2.75,4.255+0.35) cp,next,all,n1,n4!H1n1=node(-20,-2.72,4.274)n2=node(-20.5,-2.75,4.274-0.35) n3=node(-20.5,-2.75,4.274+0.35) n4=node(-19.5,-2.75,4.274-0.35) n5=node(-19.5,-2.75,4.274+0.35) n6=node(-20,-2.75,4.274)n7=node(-20,-2.75,4.274+0.35) cp,next,all,n1,n6n1=node(-20,2.72,4.274)n2=node(-20.5,2.75,4.274-0.35) n3=node(-20.5,2.75,4.274+0.35) n4=node(-19.5,2.75,4.274-0.35) n5=node(-19.5,2.75,4.274+0.35) n6=node(-20,2.75,4.274)n7=node(-20,2.75,4.274+0.35) cp,next,all,n1,n6n1=node(20,-2.72,4.274)n2=node(20.5,-2.75,4.274-0.35) n3=node(20.5,-2.75,4.274+0.35) n4=node(19.5,-2.75,4.274-0.35) n5=node(19.5,-2.75,4.274+0.35) n6=node(20,-2.75,4.274)n7=node(20,-2.75,4.274+0.35)cp,next,all,n1,n6n1=node(20,2.72,4.274)n2=node(20.5,2.75,4.274-0.35)n3=node(20.5,2.75,4.274+0.35)n4=node(19.5,2.75,4.274-0.35)n5=node(19.5,2.75,4.274+0.35)n6=node(20,2.75,4.274)n7=node(20,2.75,4.274+0.35)cp,next,all,n1,n6! 施加约束nsel,s,loc,z,5,50D,node(-40,-1.60,0), , , , , ,UY,UZ, , , ,D,node(-40,1.60,0), , , , , ,Uz,, , , ,D,node(-30,-1.60,0), , , , , ,UY,UZ, , , ,D,node(-30,1.60,0), , , , , ,Uz,, , , ,D,node(30,-1.60,0), , , , , ,UY,UZ, , , ,D,node(30,1.60,0), , , , , ,Uz,, , , ,D,node(40,-1.60,0), , , , , ,UY,UZ, , , ,D,node(40,1.60,0), , , , , ,Uz,, , , ,nsel,s,loc,x,-25.366,-50nsel,r,loc,z,-10,0.256D,all, , , , , ,all,, , , ,nsel,s,loc,x,25.366,50nsel,r,loc,z,-10,0.256D,all, , , , , ,all,, , , ,allsel!修改桥面单元,使之重量包括铺装层重量asel,s,loc,z,6.133,6.135asel,r,loc,y,2.73,-2.73ALLSEL,BELOW,AREAEMODIF,all,MA T,3,allsel3影响线求解(静力法)!在给一排单元上施加-1000kN/m^2的均布荷载,每排单元在纵向上的长度为2.5/7m。
ANSYS桥梁受力分析命令流ANSYS是一款广泛应用于工程领域的有限元分析软件,它可以帮助工程师、设计师、科学家等进行各种力学或结构分析。
在工程领域中,桥梁是一种比较常见的结构,因此本文将介绍在ANSYS中进行桥梁受力分析的命令流。
1. 建立桥梁模型首先,我们需要建立桥梁模型。
在ANSYS中,可以使用多种模型建立工具,例如现在比较流行的三维CAD软件,如SolidWorks、CATIA等,均可以建立出桥梁模型。
建立好桥梁模型后,需要导入到ANSYS中进行后续分析。
2. 分配材料属性桥梁模型建立完后,需要根据实际情况为它分配材料属性。
以混凝土桥梁为例,我们可以使用“MP(材料属性)”命令为桥梁定义混凝土的相应属性。
比如,定义混凝土的弹性模量、泊松比等。
3. 设定支座约束桥梁的受力分析需要考虑桥梁的支座约束情况。
在ANSYS中,我们可以通过使用“SUPPORT”命令设定不同类型的支座约束。
4. 设定荷载桥梁承受的荷载对于力学分析至关重要,因此在进行桥梁受力分析时,必须设定荷载的种类和大小。
ANSYS提供了多种荷载设定方式,如集中荷载、分布荷载、自重荷载等。
5. 载荷应用载荷应用是桥梁受力分析中的一个关键步骤,它能够准确模拟桥梁所受荷载。
在ANSYS中,可以使用“D”命令定义荷载的应用方式,如定义一个位移荷载、力荷载等。
6. 求解过程ANSYS中提供了多种求解器,可以快速准确地解析桥梁模型受力情况。
在进行桥梁受力分析时,需要选择合适的求解器并进行计算。
ANSYS中提供了两种类型求解器,一种是基于传统矩阵进行计算的Direct Solver求解器,另一种是基于有限元模型优化的Iterative Solver求解器。
7. 结果分析在桥梁受力分析计算完成后,需要对应力和位移等结果进行分析并产生可视化的输出。
在ANSYS中,可以使用“POST1”命令来输出结果,并对结果进行后处理。
8.经过以上步骤进行受力分析后,可以得到桥梁受力的详细情况,也可以对桥梁的结构进行合理的优化。
钢管混凝⼟拱桥ansys命令流——很实⽤,很经典fini/clear/title,concrete-filled tube arch bridge /prep7截⾯1et,1,82cyl4,0.6,0.6,0.6,,0.586cyl4,0.6,2.4,0.6,,0.586rectng,0.35,0.364,0.8,2.2rectng,0.85,0.836,0.8,2.2allselaptn,alladele,5,8,1allselaadd,allsmrtsize,5amesh,allsecwrite,gg,sect,,1sectype,1,beam,meshsecoffset,cent,,,secread,'gg','sect',,meshasel,allaclear,alladele,all,,,1截⾯2cyl4,0.6,0.6,0.586cyl4,0.6,2.4,0.586rectng,0.364,0.836,0.8,2.2allselaadd,allsmrtsize,5amesh,allsecwrite,hnt,sect,,1sectype,2,beam,meshsecoffset,cent,,,secread,'hnt','sect',,meshasel,allaclear,alladele,all,,,1建⽴横梁标准截⾯a-a k,1,-0.5k,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,10smrtsize,5amesh,allsecwrite,hl1,sect,,1 sectype,4,beam,mesh secoffset,cent,,, secread,'hl1','sect',,mesh asel,all aclear,alladele,all,,,1建⽴横梁标准截⾯b-bk,1,-0.5k,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,8smrtsize,5amesh,allsecwrite,hl2,sect,,1 sectype,5,beam,meshsecoffset,cent,,, secread,'hl2','sect',,mesh asel,allaclear,alladele,all,,,1建⽴横梁标准截⾯c-ck,1,-0.5k,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.5k,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,12 a,13,14,15,16asba,1,2smrtsize,5amesh,allsecwrite,hl3,sect,,1 sectype,6,beam,mesh secoffset,cent,,, secread,'hl3','sect',,mesh asel,all aclear,alladele,all,,,1建⽴横梁标准截⾯d-d ldele,allkdele,allk,1,-0.5k,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.35 k,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.5k,13,-0.2,0.38k,14,-0.3,0.38+0.2k,15,-0.3,1.23+0.5-0.1-0.2 k,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.2 k,19,0.3,0.38+0.2k,20,0.2,0.38a,1,2,3,4,5,6,7,8,9,10,11,12 a,13,14,15,16,17,18,19,20 asba,1,2smrtsize,5amesh,allsecwrite,hl4,sect,,1 sectype,7,beam,mesh secoffset,cent,,, secread,'hl4','sect',,mesh asel,all aclear,alladele,all,,,1建⽴横梁标准截⾯e-ek,1,-0.5k,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.5a,1,2,3,4,5,6,7,8,9,10,11,12smrtsize,5amesh,allsecwrite,hl5,sect,,1sectype,8,beam,meshsecoffset,cent,,,secread,'hl5','sect',,meshasel,allaclear,alladele,all,,,1etdele,1 !!将辅助单元删除***********************!定义所有材料特性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.5e10 mp,dens,2,2600mp,prxy,2,0.1667et,3,beam44 !!纵梁30#混凝⼟钢管特性mp,ex,3,3.0e10 mp,dens,3,2600mp,prxy,3,0.1667et,4,beam44 !!横梁30#混凝⼟钢管特性mp,ex,4,3.0e10 mp,dens,4,2600mp,prxy,4,0.1667et,5,beam44 !!风撑特性mp,ex,5,2.1e11mp,dens,5,7800mp,prxy,5,0.3et,6,link10 !!吊杆特性(钢绞线)mp,ex,6,1.9e11mp,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.85e10 mp,dens,11,2500mp,prxy,11,0.1667##############################!开始建⽴模型!0号横梁(端横梁)n,1,,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.5m位置n,5,6.25+5.1+2.3,2.025 n,6,6.25+5.1+2.3+2.23,2.025nsym,x,10,2,6,1type,4r,1real,1mat,4secnum,7!单元1到10号为0号横梁e,1,2secnum,6e,2,3secnum,5e,3,4e,4,5e,5,6secnum,7e,1,12secnum,6e,12,13secnum,5e,13,14e,14,15secnum,4e,15,16!!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.25/2+1.5/2 cm,hl,elemesel,none!!********************************!纵梁单元type,3r,2real,2mat,3sectype,10,beam,rect,zlsecdata,0.4,1.2secnum,10!0号纵梁e,1,101e,2,102e,3,103e,5,105e,12,112e,13,113!中间纵梁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****************************!定义拱轴线*dim,y,array,16 !将纵桥向距离定义为数组*dim,cs,array,16 !定义分值cs=(1-2x/l)*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(12)=28.5y(13)=42.1y(14)=48.9y(15)=59.1!0号横梁的y向坐标y(16)=2.025fl=25 !⽮⾼ml=1.1k=log(ml+sqrt(ml*ml-1)) !k=ch-1(m)*do,i,1,16,1cs(i)=1-y(i)*2/125z(i)=25-fl/(ml-1)*(cosh(k*cs(i))-1) !拱轴⾼度*enddo!上述定义的拱轴线可参照姚玲森《桥梁⼯程》p1299*dim,ygd,array,23 !定义预⼯拱度(是在后续⼯况中反复调算的结果)!初始运⾏中可将预拱度全部设为0ygd(1)=0.4000e-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) ngen,2,10000,204,214,10,0,0,z(2)+ygd(2) ngen,2,10000,304,314,10,0,0,z(3)+ygd(3) ngen,2,10000,404,414,10,0,0,z(4)+ygd(4) ngen,2,10000,504,514,10,0,0,z(5)+ygd(5) ngen,2,10000,604,614,10,0,0,z(6)+ygd(6) ngen,2,10000,704,714,10,0,0,z(7)+ygd(7) ngen,2,10000,804,814,10,0,0,z(8)+ygd(8) ngen,2,10000,904,914,10,0,0,z(9)+ygd(9) ngen,2,10000,1004,1014,10,0,0,z(8)+ygd(10) ngen,2,10000,1104,1114,10,0,0,z(7)+ygd(11) ngen,2,10000,1204,1214,10,0,0,z(6)+ygd(12) ngen,2,10000,1304,1314,10,0,0,z(5)+ygd(13) ngen,2,10000,1404,1414,10,0,0,z(4)+ygd(14) ngen,2,10000,1504,1514,10,0,0,z(3)+ygd(15) ngen,2,10000,1604,1614,10,0,0,z(2)+ygd(16) ngen,2,10000,1704,1714,10,0,0,z(1)+ygd(17)ngen,2,3000,104,114,10,0,0,0ngen,2,3000,204,214,10,0,0,0ngen,2,3000,304,314,10,0,0,0ngen,2,3000,404,414,10,0,0,0ngen,2,3000,504,514,10,0,0,0ngen,2,3000,604,614,10,0,0,0ngen,2,3000,704,714,10,0,0,0ngen,2,3000,804,814,10,0,0,0ngen,2,3000,904,914,10,0,0,0ngen,2,3000,1004,1014,10,0,0,0ngen,2,3000,1104,1114,10,0,0,0ngen,2,3000,1204,1214,10,0,0,0ngen,2,3000,1304,1314,10,0,0,0ngen,2,3000,1404,1414,10,0,0,0ngen,2,3000,1504,1514,10,0,0,0ngen,2,3000,1604,1614,10,0,0,0ngen,2,3000,1704,1714,10,0,0,0ngen,2,6000,104,114,10,0,0,z(1)+ygd(1)ngen,2,6000,204,214,10,0,0,z(2)+ygd(2)ngen,2,6000,304,314,10,0,0,z(3)+ygd(3)ngen,2,6000,404,414,10,0,0,z(4)+ygd(4)ngen,2,6000,504,514,10,0,0,z(5)+ygd(5)ngen,2,6000,604,614,10,0,0,z(6)+ygd(6)ngen,2,6000,704,714,10,0,0,z(7)+ygd(7)ngen,2,6000,804,814,10,0,0,z(8)+ygd(8)ngen,2,6000,904,914,10,0,0,z(9)+ygd(9)ngen,2,6000,1004,1014,10,0,0,z(8)+ygd(10)ngen,2,6000,1104,1114,10,0,0,z(7)+ygd(11)ngen,2,6000,1204,1214,10,0,0,z(6)+ygd(12)ngen,2,6000,1304,1314,10,0,0,z(5)+ygd(13)ngen,2,6000,1404,1414,10,0,0,z(4)+ygd(14)ngen,2,6000,1504,1514,10,0,0,z(3)+ygd(15)ngen,2,6000,1604,1614,10,0,0,z(2)+ygd(16)ngen,2,6000,1704,1714,10,0,0,z(1)+ygd(17)!*do,i,1,9,1!ngen,2,10000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i) !吊杆节点!*enddo!!jj=1!*do,i,10,17,1!ngen,2,10000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i) !吊杆节点!jj=jj+1!*enddo!!为了对以后吊杆转动进⾏放松,重复吊杆节点!*do,i,1,9,1!ngen,2,3000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i) !吊杆下节点!*enddo!!jj=1!*do,i,10,17,1!ngen,2,3000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i) !吊杆节点!jj=jj+1!*enddo!!*do,i,1,10,1!ngen,2,6000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i) !吊杆节点!*enddo!!jj=1!*do,i,10,17,1!ngen,2,6000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i) !吊杆节点!jj=jj+1!*enddo!定义吊杆的参数(每根吊杆均采⽤109根直径7mm⾼强钢丝)dgarea=109*3.141593*0.25*0.007*0.007 dgyl=0.3*1670*1e6 !吊杆除应⼒采⽤0.3倍抗拉强度(n/mm*mm)dgl=dgyl*dgareaedg=1.95e11 !吊杆弹模*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) !吊杆沿中轴对称*enddo*do,i,1,17,1fdg(i)=dgl !各吊杆施加的预应⼒*enddotype,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!!*********************************real,1018 !定义⼀个全为零的实常数便于利⽤截⾯特性时,不产⽣冲突!!********************************* !⽣成风撑上的节点n,20001,0,21.7,z(11)n,20002,0,28.5,z(12)n,20003,0,42.1,z(13)n,20004,0,48.9,z(14)n,20005,0,59.1,z(15)n,20006,0,65.9,z(15)n,20007,0,76.1,z(14)n,20008,0,82.9,z(13)n,20009,0,96.5,z(12)n,20010,0,103.3,z(11)n,20011,12.5,21.7,z(11)n,20012,12.5,28.5,z(12)n,20013,12.5,42.1,z(13)n,20014,12.5,48.9,z(14)n,20015,12.5,59.1,z(15)+ygd(18)n,20016,12.5,65.9,z(15)+ygd(19)n,20017,12.5,76.1,z(14)n,20018,12.5,82.9,z(13)n,20019,12.5,96.5,z(12)n,20020,12.5,103.3,z(11)n,20021,-12.5,21.7,z(11)n,20022,-12.5,28.5,z(12)n,20023,-12.5,42.1,z(13)n,20024,-12.5,48.9,z(14)n,20025,-12.5,59.1,z(15)+ygd(18)n,20026,-12.5,65.9,z(15)+ygd(19)n,20027,-12.5,76.1,z(14)n,20028,-12.5,82.9,z(13)n,20029,-12.5,96.5,z(12)n,20030,-12.5,103.3,z(11)n,20000,0,62.5,z(15) !中撑中点!风撑1(倒k撑)type,5mat,5sectype,500,beam,ctube,windc !定义风撑形状为空⼼圆secdata,0.486,0.50secnum,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,90000esel,s,real,,1018cm,fc,elemesel,none!!******************************************* !连接拱轴线!定义拱轴上的坐标点n,30000,12.5 !左拱轴起点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 !右拱轴起点n,40001,-12.5,125,0 !右拱轴端点n,40002,-12.5,2.025,z(16) !右横梁起点n,40003,-12.5,125-2.025,z(16) !右横梁端点*dim,zgzx,array,23 !定义拱轴线,左轴线*dim,ygzx,array,23 !定义拱轴线,右轴线zgzx(1)=30000 !以下为拱轴线上的节点号zgzx(2)=30002 *do,i,3,10,1zgzx(i)=10004+(i-2)*100*enddozgzx(11)=20015zgzx(12)=10904zgzx(13)=20016*do,i,14,21,1zgzx(i)=11004+(i-14)*100*enddozgzx(22)=30003zgzx(23)=30001ygzx(1)=40000ygzx(2)=40002*do,i,3,10,1ygzx(i)=10014+(i-2)*100*enddoygzx(11)=20025ygzx(12)=10914ygzx(13)=20026*do,i,14,21,1ygzx(i)=11014+(i-14)*100*enddoygzx(22)=40003ygzx(23)=40001n,90001,12.5,0,30 !钢管左轴线参考点n,90002,-12.5,0,30 !钢管右轴线参考点type,1mat,1secnum,1*do,i,1,22 !⽣成拱轴单元(钢管)e,zgzx(i),zgzx(i+1),90001e,ygzx(i),ygzx(i+1),90002*enddocm,gzgg,elemesel,nonetype,2mat,2secnum,2*do,i,1,22e,zgzx(i),zgzx(i+1),90001 !⽣成拱轴单元(混凝⼟)e,ygzx(i),ygzx(i+1),90002*enddocm,gzhnt,elemesel,none*****************************************!建⽴拱轴底座type,9r,6real,6mat,9sectype,70,beam,rect,dizuo !定义截⾯形状secdata,1.4,4secnum,70e,30001,30003,90001e,40001,40003,90002e,30000,30002,90001e,40000,40002,90002cm,gzdz,elemesel,none*****************************************!建⽴盖梁ngen,2,2000,1,16,1,0,-2.025,0 !复制节点到拱轴线起点位置ngen,2,100,1801,1816,1,0,2.025,0 !复制节点到拱轴线终点位置type,7mat,7r,7sectype,60,beam,rect,gail secdata,4,1.6secnum,60n,90005,0,0,20n,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,90005e,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,elemesel,none!*****************************!建⽴墩柱n,90003,14.5 ! 参考点n,90004,-14.5type,8mat,8r,8sectype,50,beam,hrec,dun !定义截⾯形状secdata,3,4,0.8,0.8,0.8,0.8 ngen,2,10000,2004,2014,10,0,0,-10ngen,2,10000,1904,1914,10,0,0,-10 !墩柱节点,距离桥⾯-10m secnum,50e,2004,12004,90003e,2014,12014,90004e,1904,11904,90003e,1914,11914,90004cm,dunzhu,elemesel,none*****************************************!承台单元type,11mat,11r,1sectype,40,beam,rect,chengtaisecdata,12.7,9r,11real,11secnum,40secdata,3,4,0.8,0.8,0.8,0.8ngen,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,elemesel,none!!**********************************!横系梁type,9mat,9r,9real,9n,90007,12.5,0,10 !参考点n,90008,12.5,0,10sectype,90,beam,rect,hxlsecdata,3,3secnum,90e,22004,22014,90007e,21904,21914,90008cm,hxl,elemesel,none!!**********************************!桩n,50000,6.15+1.85,2.7,-13 !定义桩节点n,50001,6.15+1.85,-2.7,-13 n,50002,6.15+1.85+4.5,2.7,-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个节点!每3m⽣成⼀层节点(共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,90003 e,50000+k+100+i*1000,50000+k+100+(i+1)*1000,90003 e,50000+k+110+i*1000,50000+k+110+(i+1)*1000,90003 *enddo*enddocm,zhuang,elemesel,none***********************************!在桩上施加弹簧单元!第⼀层距离承台-3m 317-340!计算弹簧的弹性系数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/h1。
ANSYS最常用命令流+中文注释VSBV, NV1, NV2, SEPO, KEEP1, KEEP2 —Subtracts volumes from volumes,用于2个solid相减操作,最终目的是要nv1-nv2=?通过后面的参数设置,可以得到很多种情况:sepo项是2个体的边界情况,当缺省的时候,是表示2个体相减后,其边界是公用的,当为sepo的时候,表示相减后,2个体有各自的独立边界。
keep1与keep2是询问相减后,保留哪个体?当第一个为keep时,保留nv1,都缺省的时候,操作结果最终只有一个体,比如:vsbv,1,2,sepo,,keep,表示执行1-2的操作,结果是保留体2,体1被删除,还有一个1-2的结果体,现在一共是2个体(即1-2与2),且都各自有自己的边界。
如vsbv,1,2,,keep,,则为1-2后,剩下体1和体1-2,且2个体在边界处公用。
同理,将v换成a 及l是对面和线进行减操作!mp,lab, mat, co, c1,…….c4 定义材料号及特性lab: 待定义的特性项目(ex,alpx,reft,prxy,nuxy,gxy,mu,dens) ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号(缺省为当前材料号)co: 材料特性值,或材料之特性,温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项,2次项,3次项,4次项的系数定义DP材料:首先要定义EX和泊松比:MP,EX,MA T,……MP,NUXY,MAT,……定义DP材料单元表(这里不考虑温度):TB,DP,MA T进入单元表并编辑添加单元表:TBDATA,1,CTBDATA,2,ψTBDATA,3,……如定义:EX=1E8,NUXY=0.3,C=27,ψ=45的命令如下:MP,EX,1,1E8MP,NUXY,1,0.3TB,DP,1TBDATA,1,27TBDATA,2,45这里要注意的是,在前处理的最初,要将角度单位转化到“度”,即命令:*afun,degVSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Type,是选择的方式,有选择(s),补选(a),不选(u),全选(all)、反选(inv)等,其余方式不常用Item, Comp 是选取的原则以及下面的子项如volu 就是根据实体编号选择,loc 就是根据坐标选取,它的comp就可以是实体的某方向坐标!其余还有材料类型、实常数等MIN, VMAX, VINC,这个就不必说了吧!,例:vsel,s,volu,,14vsel,a,volu,,17,23,2上面的命令选中了实体编号为14,17,19,21,23的五个实体VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体nv1:初始体号nv2:最终的体号ninc:体号之间的间隔kswp=0:只删除体kswp=1:删除体及组成关键点,线面如果nv1=all,则nv2,ninc不起作用其后面常常跟着一条显示命令VPLO,或aplo,nplo,这个湿没有参数的命令,输入后直接回车,就可以显示刚刚选择了的体、面或节点,很实用的哦!Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点(缺省)R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0” 使用正负号“1”仅用绝对值下面是单元生死第一个载荷步中命令输入示例:!第一个载荷步TIME,... !设定时间值(静力分析选项)NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿-拉夫森选项ESTIF,... !设定非缺省缩减因子(可选)ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点(不与活动单元相连的结点)D,ALL,ALL,0 !约束所有不活动的结点自由度(可选)NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SA VESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D, F,SF和BF命令得到更详细的解释。
Finish(退出四大模块,回到BEGIN层)/clear (清空内存,开始新的计算)1.定义参数、数组,并赋值.2./prep7(进入前处理)定义几何图形:关键点、线、面、体定义几个所关心的节点,以备后处理时调用节点号。
设材料线弹性、非线性特性设置单元类型及相应KEYOPT设置实常数设置网格划分,划分网格根据需要耦合某些节点自由度定义单元表3./solu加边界条件设置求解选项定义载荷步求解载荷步4./post1(通用后处理)5./post26 (时间历程后处理)6.PLOTCONTROL菜单命令7.参数化设计语言8.理论手册Finish(退出四大模块,回到BEGIN层)/clear (清空内存,开始新的计算)1.定义参数、数组,并赋值.dim, par, type, imax, jmax, kmax, var1, vae2, var3 定义数组par: 数组名type:array 数组,如同fortran,下标最小号为1,可以多达三维(缺省)char 字符串组(每个元素最多8个字符)tableimax,jmax, kmax 各维的最大下标号var1,var2,var3 各维变量名,缺省为row,column,plane(当type为table时) 2./prep7(进入前处理)2.1 设置单元类型及相应KEYOPTET, itype, ename, kop1……kop6, inopr 设定当前单元类型Itype:单元号Ename:单元名设置实常数Keyopt, itype, knum, valueitype: 已定义的单元类型号knum: 单元的关键字号value: 数值注意:如果,则必须使用keyopt命令,否则也可在ET命令中输入2.2 定义几个所关心的节点,以备后处理时调用节点号。
n,node,x,y,z,thxy, thyz, thzx 根据坐标定义节点号如果已有此节点,则原节点被重新定义,一般为最大节点号。
Ansys命令流大全1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9此命令用已知的一组关键点点(P1~P9)来定义面(Area),最少使用三个点才能围成面,同时产生转围绕些面的线。
点要依次序输入,输入的顺序会决定面的法线方向。
如果超过四个点,则这些点必须在同一个平面上。
Menu Paths:MainMenu>Preprocessor>Create>Arbitrary>Through KPs2、*ABBR,Abbr,String--定义一个缩略语.Abbr:用来表示字符串"String"的缩略语,长度不超过8个字符.String:将由"Abbr"表示的字符串,长度不超过60个字符.3、ABBRES,Lab,Fname,Ext-从一个编码文件中读出缩略语.Lab:指定读操作的标题,NEW:用这些读出的缩略语重新取代当前的缩略语(默认)CHANGE:将读出的缩略语添加到当前缩略语阵列,并替代现存同名的缩略语.Ext:如果"Fname"是空的,则缺省的扩展命是"ABBR".4、ABBSAV,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命令流解释大全A N S Y S命令流解释大全Document serial number【LGGKGB-LGG98YT-LGGT8CB-LGUT-一、定义材料号及特性mp,lab, mat, co, c1,…….c4lab: 待定义的特性项目(ex,alpx,reft,prxy,nuxy,gxy,mu,dens)ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号(缺省为当前材料号)c 材料特性值,或材料之特性,温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项,2次项,3次项,4次项的系数二、定义DP材料:首先要定义EX和泊松比:MP,EX,MAT,……MP,NUXY,MAT,……定义DP材料单元表(这里不考虑温度):TB,DP,MAT进入单元表并编辑添加单元表:TBDATA,1,CTBDATA,2,ψTBDATA,3,……如定义:EX=1E8,NUXY=,C=27,ψ=45的命令如下:MP,EX,1,1E8MP,NUXY,1,TB,DP,1TBDATA,1,27TBDATA,2,45这里要注意的是,在前处理的最初,要将角度单位转化到“度”,即命令:*afun,deg三、单元生死载荷步!第一个载荷步TIME,... !设定时间值(静力分析选项)NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿-拉夫森选项ESTIF,... !设定非缺省缩减因子(可选)ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点(不与活动单元相连的结点)D,ALL,ALL,0 !约束所有不活动的结点自由度(可选)NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SAVESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SF和BF命令得到更详细的解释。
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"的缩略语,长度不超过8个字符.String:将由"Abbr"表示的字符串,长度不超过60个字符.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为半径。
9、*AFUN,Lab在参数表达式中,为角度函数指定单位.Lab:指定将要使用的角度单位.有3个选项.RAD:在角度函数的输入与输出中使用弧度单位(默认)DEG:在角度函数的输入与输出中使用度单位.STAT:显示该命令当前的设置(即是度还是弧度).10、Agen, itime,na1,na2,ninc,dx,dy,dz,kinc,noelem,imove !面积复制命令。
1模型简介模型以壳单元shell63为主,横梁采用beam44单元,吊杆采用link10单元,横梁、吊杆与拱肋的连接采用节点耦合的方式。
不考虑纵坡和横披。
111213141516X Y ZROTCPACEL图1 三维图111213141516X Y Z 图2 正立面图2建模命令流finish/clear/config,nproc,2/prep7!双核运行!材料!钢材EX,1,2.06e11PRXY ,1,0.3DENS,1,7800MPDA TA,ALPX,1,,1.2e-5!吊杆钢材EX,2,2.06e11PRXY ,2,0.3DENS,2,7800MPDA TA,ALPX,2,,1.2e-5!桥面钢材包括铺装层重量EX,3,2.06e11PRXY,3,0.3DENS,3,7800*3MPDA TA,ALPX,3,,1.2e-5!单元ET,1,plane82ET,2,shell63!设置shell63关键字KEYOPT,2,1,0KEYOPT,2,2,0KEYOPT,2,3,0KEYOPT,2,5,2KEYOPT,2,6,0KEYOPT,2,7,0KEYOPT,2,8,0KEYOPT,2,9,0KEYOPT,2,11,0et,3,beam189!设置beam189关键字KEYOPT,3,1,1KEYOPT,3,2,0KEYOPT,3,4,2KEYOPT,3,6,3KEYOPT,3,7,2KEYOPT,3,8,3KEYOPT,3,9,3KEYOPT,3,10,0KEYOPT,3,11,0KEYOPT,3,12,0ET,4,solid45et,5,link10r,1,0.010r,2,0.016r,3,0.02r,4,0.04!建立各种截面!钢横梁!HL1/HL3a=1.0b=0.680th=0.02ts=0.016RECTNG,-a/2,a/2,-b/2,b/2, RECTNG,-a/2+ts,a/2-ts,-b/2+th,b/2-th,RECTNG,-a/2+ts,-a/2+ts+0.16,-0.01,0.01, RECTNG,a/2-ts,a/2-ts-0.16,-0.01,0.01, RECTNG,-0.01,0.01,-b/2+0.02,-b/2+0.02+0.16, RECTNG,-0.01,0.01,b/2-0.02-0.16,b/2-0.02, AGLUE,allwprot,0,90,0wpoff,0,0,-b/2+thASBW,allwpoff,0,0,b/2-th-0.01ASBW,allwpoff,0,0,0.02ASBW,allwpoff,0,0,b/2-th-0.01ASBW,allwprot,0,0,90wpoff,0,0,-a/2+tsASBW,allwpoff,0,0,a-2*tsASBW,allwpoff,0,0,-a/2+ts+0.01ASBW,allwpoff,0,0,-0.02ASBW,allESIZE,0.3,0,amesh,allSECWRITE,sec3SECTYPE,3,BEAM,MESH, SECOFFSET,CENT,,,SECREAD,'sec3','SECT',' ',MESHaclear,allADELE,all, , ,1WPCSYS,-1,0!HL2!吊杆SECTYPE, 4, BEAM, CTUBE, , 0 SECOFFSET, CENTSECDA TA,0.032,0.04,,0,0,0,0,0,0,0!建立跨中截面所有线k,1,0,-2.058,0+5.45k,2,0,-2,0+5.45k,,0,-1.6,0+5.45k,,0,-0.8,0+5.45k,,0,-0.4,0+5.45k,,0,-1.6,0.108+5.45k,,0,-1.2,0.108+5.45k,,0,-0.8,0.108+5.45k,,0,-0.4,0.108+5.45k,,0,-2.541661,0.200339+5.45 k,,0,-2.37,0.576+5.45k,,0,-1.6,0.576+5.45k,,0,-1.2,0.576+5.45k,,0,-0.8,0.576+5.45k,,0,-0.4,0.576+5.45k,,0,-2.742,0.684+5.45k,,0,-2.37,0.684+5.45k,,0,-2,0.684+5.45k,,0,-1.6,0.684+5.45k,,0,-1.2,0.684+5.45k,,0,-0.8,0.684+5.45k,,0,-0.4,0.684+5.45k,,0,0,0.684+5.45k,,0,0,0+5.45*do,i,1,5,1l,i,i+1*enddol,6,25*do,i,17,24,1l,i,i+1*enddo*do,i,3,6,1l,i,i+4*enddo*do,i,20,23,1l,i,i-7*enddol,18,12larc,1,17,11l,2,19!建立拖拉的路径线k,26,-40,0,0+5.45l,25,26!拖拉成箱梁*do,i,1,25,1ADRAG,i, , , , , ,26*enddoallselNUMMRG,KP, , , ,LOWwprot,0,0,90*do,i,1,15,1wpoff,0,0,-2.5ASBW,all*enddo!通过循环生成横隔板,选择封闭的为轮廓线来生成!生成矩形部分*do,i,0,16,1lsel,s,loc,x,-2.5*i-0.001,-2.5*i+0.001lsel,r,loc,y,-2.001,0.001lsel,u,loc,z,0.001+5.45,0.2+5.45lsel,u,loc,z,0.683+5.45,0.6+5.45al,all*enddo!生成扇形部分*do,i,0,16,1lsel,s,loc,x,-2.5*i-0.001,-2.5*i+0.001lsel,r,loc,y,-1.999,-2.749lsel,u,loc,z,0.683+5.45,0.6+5.45al,all*enddoallsel!多次分割WPCSYS,-1,0asel,u,loc,z,0.001+5.45,0.2+5.45asel,u,loc,z,0.683+5.45,0.6+5.45wprot,0,90,0*do,i,1,6,1KWPA VE,iasbw,all*enddoKWPA VE,18asbw,allallselWPA VE,0,0,0wpoff,0,0,0.2asbw,allwpoff,0,0,1.58asbw,allallselWPCSYS,-1,0KWPA VE,10asbw,allKWPA VE,16asbw,allWPCSYS,-1,0wpoff,0,0,0.2+5.45asbw,allwpoff,0,0,0.268asbw,all!删除需要空心的地方*do,j,0,3*do,i,1,3,1asel,s,loc,x,-2.5*(i+4*j)-0.001,-2.5*(i+4*j)+0.001 asel,r,loc,z,0.2+5.45,0.2+0.268+5.45asel,r,loc,y,-0.2,-0.2-1.580ADELE,all, , ,1*enddo*enddoallselaglue,allNUMMRG,KP, , , ,LOWNUMCMP,KPNUMCMP,LINENUMCMP,AREAWPCSYS,-1,0!建立壳单元拱肋!上拱肋跨中断面关键点k,10100,0,-3.25-0.5,-0.35+10.1k,10101,0,-3.25-0.5,0+10.1k,10102,0,-3.25-0.5,0.35+10.1k,10103,0,-3.25-0.34,0+10.1k,10104,0,-3.25-0.25,-0.35+10.1k,10105,0,-3.25-0.25,-0.19+10.1k,10106,0,-3.25-0.25,0.19+10.1k,10107,0,-3.25-0.25,0.35+10.1k,10108,0,-3.25-0.16,0+10.1k,10109,0,-3.25+0,-0.35+10.1k,10110,0,-3.25+0,0+10.1k,10111,0,-3.25+0,0.35+10.1k,10112,0,-3.25+0.16,0+10.1k,10113,0,-3.25+0.25,-0.35+10.1k,10114,0,-3.25+0.25,-0.19+10.1k,10115,0,-3.25+0.25,0.19+10.1 k,10116,0,-3.25+0.25,0.35+10.1 k,10117,0,-3.25+0.34,0+10.1 k,10118,0,-3.25+0.5,-0.35+10.1 k,10119,0,-3.25+0.5,0+10.1 k,10120,0,-3.25+0.5,0.35+10.1 !下拱肋跨中断面关键点k,10200,0,-3.25-0.5,-0.35+5 k,10201,0,-3.25-0.5,0+5k,10202,0,-3.25-0.5,0.35+5k,10203,0,-3.25-0.34,0+5k,10204,0,-3.25+0,-0.35+5k,10205,0,-3.25+0,-0.19+5k,10206,0,-3.25+0,0.19+5k,10207,0,-3.25+0,0.35+5k,10208,0,-3.25+0.34,0+5k,10209,0,-3.25+0.5,-0.35+5 k,10210,0,-3.25+0.5,0+5k,10211,0,-3.25+0.5,0.35+5!连成线NUMCMP,LINELSTR, 10118, 10113 LSTR, 10113, 10109 LSTR, 10109, 10104 LSTR, 10104, 10100 LSTR, 10100, 10101 LSTR, 10101, 10102 LSTR, 10102, 10107 LSTR, 10107, 10111 LSTR, 10111, 10116 LSTR, 10116, 10120 LSTR, 10120, 10119 LSTR, 10119, 10118 LSTR, 10111, 10110 LSTR, 10110, 10109 LSTR, 10116, 10115 LSTR, 10113, 10114 LSTR, 10119, 10117 LSTR, 10110, 10112 LSTR, 10107, 10106 LSTR, 10104, 10105 LSTR, 10101, 10103 LSTR, 10110, 10108LSTR, 10209, 10204LSTR, 10204, 10200LSTR, 10200, 10201LSTR, 10201, 10202LSTR, 10202, 10207LSTR, 10207, 10211LSTR, 10211, 10210LSTR, 10210, 10209LSTR, 10210, 10208LSTR, 10201, 10203LSTR, 10204, 10205LSTR, 10207, 10206!通过旋转生成面!建立上拱肋圆心关键点k,10300,0,0,-27.141k,10301,0,-3,-27.141!建立下拱肋圆心关键点k,10302,0,0,-67.525+5k,10303,0,-3,-67.525+5 WPCSYS,-1,0!上拱肋*do,i,2499,2520,1AROTA T, i,,,,,,10300,10301,43.21447, *enddo!下拱肋*do,i,2521,2532,1AROTA T, i,,,,,,10302,10303,22.18745, *enddo!拱脚的处理ASEL,S,LOC,Y,-3+0.251,-10 ALLSEL,BELOW,AREA NUMMRG,KP, , , ,LOW NUMCMP,KPNUMCMP,LINENUMCMP,AREANUMMRG,KP, , , SELE,LOW ASEL,R,LOC,Z,0,4.5wprot,0,90,0wpoff,0,0,3.00asbw,allwpoff,0,0,0.50asbw,allASEL,S,LOC,Y,-3+0.251,-10WPCSYS,-1,0wprot,0,0,90*do,i,1,5,1wpoff,0,0,-3asbw,all*enddoWPA VE,0,0,0*do,i,1,2,1wpoff,0,0,-10asbw,all*enddoasel,r,loc,x,-25,-20AOVLAP,1593,1625KWPLAN,-1, 1441, 1440, 1112ASBW,allasel,r,loc,x,-25,-23.5ADELE,all, , ,1ASEL,S,LOC,Y,-3+0.251,-10ALLSEL,BELOW,AREANUMMRG,KP, , , SELE,LOWWPCSYS,-1,0!拱脚加劲板LSTR, 1409, 1431al,2710,2703,2533LSTR, 1432, 1405al,3177,2536,2992!上拱肋横隔板*do,i,0,5,1lSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5.001+0.35,10.101+0.35lsel,r,loc,x,-3*i-0.001,-3*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-3*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-3*i)**2)** 0.5)+0.05lsel,r,loc,y,-2.749,-3.251al,alllSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5.001+0.35,10.101+0.35lsel,r,loc,x,-3*i-0.001,-3*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-3*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-3*i)**2)** 0.5)+0.05lsel,r,loc,y,-3.249,-4al,all*enddoi=1lSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5.001,10.101+0.35lsel,r,loc,x,-10*i-0.001,-10*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-10*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-10*i)**2) **0.5)+0.05lsel,r,loc,y,-2.749,-3.251al,alllSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5.001,10.101+0.35lsel,r,loc,x,-10*i-0.001,-10*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-10*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-10*i)**2) **0.5)+0.05lsel,r,loc,y,-3.249,-4al,alli=2lSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,3.8,10.101+0.35lsel,r,loc,x,-10*i-0.001,-10*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-10*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-10*i)**2) **0.5)+0.05lsel,r,loc,y,-2.749,-3.251al,alllSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,3.8,10.101+0.35lsel,r,loc,x,-10*i-0.001,-10*i+0.001lsel,u,loc,y,-3.001,-2.999lsel,u,loc,y,-3.501,-3.499lsel,u,loc,z,10.1-(37.241-(37.241**2-(-10*i)**2)**0.5)-0.05,10.1-(37.241-(37.241**2-(-10*i)**2) **0.5)+0.05lsel,r,loc,y,-3.249,-4al,all!下拱肋横隔板*do,i,0,5,1lSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5-(67.525-(67.525**2-(-3*i)**2)**0.5)-0.6,5-(67.525-(67.525**2-(-3*i)**2)**0.5)+0. 6lsel,r,loc,x,-3*i-0.001,-3*i+0.001lsel,u,loc,y,-3.251,-3.249lsel,u,loc,z,5-(67.525-(67.525**2-(-3*i)**2)**0.5)-0.05,5-(67.525-(67.525**2-(-3*i)**2)**0.5)+ 0.05al,all*enddo*do,i,10,20,10lSEL,S,LOC,Y,-3+0.251,-10lsel,r,loc,z,5-(67.525-(67.525**2-(-i)**2)**0.5)-0.6,5-(67.525-(67.525**2-(-i)**2)**0.5)+0.6 lsel,r,loc,x,-i-0.001,-i+0.001lsel,u,loc,y,-3.251,-3.249lsel,u,loc,z,5-(67.525-(67.525**2-(-i)**2)**0.5)-0.05,5-(67.525-(67.525**2-(-i)**2)**0.5)+0.05 al,all*enddoallsel!删除上拱肋多余部分asel,s,loc,y,-3.251,-10asel,r,loc,x,-12,-50KWPA VE, 1411asbw,allasel,r,loc,z,kz(997),20ADELE,all, , ,1asel,s,loc,y,-3.251,-10asel,r,loc,z,5.5,20ADELE,all, , ,1asel,s,loc,y,-2.749,-10NUMMRG,KP, , , SELE,LOW!建立新吊杆线LSTR, 1174, 1152LSTR, 1215, 1190LSTR, 1252, 1229LSTR, 1289, 1266LSTR, 1326, 1303LSTR, 1363, 1340FLST,5,6,4,ORDE,2FITEM,5,2501FITEM,5,-2506LSEL,R, , ,P51Xcm,diaogan,lineallsel!建立横梁线allselk,10001,0,0,kz(1158)k,10002,0,-2.73,kz(1158) k,10003,-10,0,kz(1387) k,10004,-10,-2.73,kz(1387) k,10005,-20,0,kz(1410) k,10006,-20,-2.73,kz(1410)l,10001,10002l,10003,10004l,10005,10006lsel,s,loc,z,-10,5.01lsel,r,loc,y,0,-2.73 WPCSYS,-1,0KWPA VE, 10001 wprot,0,90,0wpoff,0,0,1.6lsbw,allcm,hengliang,line!划分!划分横梁type,3mat,1real,1SECNUM,3esize,0.5lmesh,hengliang!划分吊杆allselTYPE,5mat,2real,1SECNUM,4esize,20lmesh,diaoganallsel!划分横隔板type,2mat,1*do,i,0,40,2.5asel,s,loc,x,-i+0.001,-i-0.001 asel,r,loc,y,0,-2.75 ALLSEL,BELOW,AREA LESIZE,all,0.25, , , , , , ,1 aclear,allamesh,all*enddo!划分纵隔板allselLESIZE,all, , ,-1, ,1asel,s,loc,y,-0.001,0.001 asel,a,loc,y,-2.001,-1.999 ALLSEL,BELOW,AREA LESIZE,all,0.4, , , , , , ,1 type,2mat,1real,3amesh,all!划分加劲肋asel,s,loc,y,-2.370,-2.370*do,i,0.4,1.6,0.4asel,a,loc,y,-i-0.001,-i+0.001 *enddoALLSEL,BELOW,AREA LESIZE,all, , ,-1, ,1 LESIZE,all,0.4, , , , , , ,1 type,2mat,1real,1amesh,allallsel!划分箱体外包面ASEL,S,MA T,,0asel,r,loc,y,0.001,-2.75+0.001 LESIZE,all, , ,-1, ,1 LESIZE,all,0.4, , , , , , ,1 real,2amesh,allallsel!拱肋划分allselLESIZE,all, , ,-1, ,1mat,1real,3esize,0.5asel,s,loc,y,-2.749,-10 amesh,allallsel!镜像生成全桥LSEL,S,TYPE,,3LSEL,a,TYPE,,5 LSYMM,Y,all, , , ,0,0 asel,s,loc,y,-0.001,-10 ARSYM,Y,all, , , ,0,0LSEL,S,TYPE,,3LSEL,a,TYPE,,5lsel,r,loc,x,-50,-0.001 LSYMM,x,all, , , ,0,0 asel,s,loc,x,-0.001,-50 ARSYM,x,all, , , ,0,0allselNUMMRG,NODE, , , ,LOW NUMMRG,KP, , , ,LOW NUMMRG,NODE,,, ,LOW !统一所有单元坐标allselcsys,0EMODIF,P51X,ESYS, ,!耦合横梁和箱梁!跨中nsel,s,loc,x,-0.001,0.001 nsel,r,loc,y,-1.601,-1.599n1=node(0,-1.60,0)nsel,u,,,n1n2=node(0,-1.60,0)nsel,a,,,n1cp,1,ux,n1,n2cp,next,uy,n1,n2cp,next,uz,n1,n2nsel,s,loc,x,-0.001,0.001 nsel,r,loc,y,1.601,1.599n1=node(0,1.60,0)nsel,u,,,n1n2=node(0,1.60,0)nsel,a,,,n1cp,next,ux,n1,n2cp,next,uz,n1,n2!H2横梁处nsel,s,loc,x,-10.001,-9.999 nsel,r,loc,y,-1.601,-1.599n1=node(0,-1.60,0) nsel,u,,,n1n2=node(0,-1.60,0) nsel,a,,,n1cp,next,uy,n1,n2cp,next,uz,n1,n2nsel,s,loc,x,-10.001,-9.999 nsel,r,loc,y,1.601,1.599n1=node(0,1.60,0)nsel,u,,,n1n2=node(0,1.60,0)nsel,a,,,n1cp,next,uz,n1,n2!H2'横梁处nsel,s,loc,x,10.001,9.999 nsel,r,loc,y,-1.601,-1.599n1=node(0,-1.60,0) nsel,u,,,n1n2=node(0,-1.60,0) nsel,a,,,n1cp,next,uy,n1,n2cp,next,uz,n1,n2nsel,s,loc,x,10.001,9.999 nsel,r,loc,y,1.601,1.599n1=node(0,1.60,0)nsel,u,,,n1n2=node(0,1.60,0)nsel,a,,,n1cp,next,uz,n1,n2!H1横梁处nsel,s,loc,x,-20.001,-19.999 nsel,r,loc,y,-1.601,-1.599n1=node(0,-1.60,0) nsel,u,,,n1n2=node(0,-1.60,0) nsel,a,,,n1cp,next,uy,n1,n2cp,next,uz,n1,n2nsel,s,loc,x,-20.001,-19.999 nsel,r,loc,y,1.601,1.599n1=node(0,1.60,0)nsel,u,,,n1n2=node(0,1.60,0)nsel,a,,,n1cp,next,uz,n1,n2!H1'横梁处nsel,s,loc,x,20.001,19.999 nsel,r,loc,y,-1.601,-1.599n1=node(0,-1.60,0) nsel,u,,,n1n2=node(0,-1.60,0) nsel,a,,,n1cp,next,uy,n1,n2cp,next,uz,n1,n2nsel,s,loc,x,20.001,19.999 nsel,r,loc,y,1.601,1.599n1=node(0,1.60,0)nsel,u,,,n1n2=node(0,1.60,0)nsel,a,,,n1cp,next,uz,n1,n2allsel!耦合横梁与拱肋!跨中横梁n1=node(0,-2.72,5)n2=node(-0.5,-2.75,5)n3=node(0.5,-2.75,5)n4=node(0,-2.75,5)n5=node(0,-2.75,5+0.35) cp,next,all,n1,n4n1=node(0,2.72,5)n2=node(-0.5,2.75,5)n3=node(0.5,2.75,5)n4=node(0,2.75,5)n5=node(0,2.75,5+0.35) cp,next,all,n1,n4!H2n1=node(-10,-2.72,4.255) n2=node(-10.5,-2.75,4.255) n3=node(-9.5,-2.75,4.255)n4=node(-10,-2.75,4.255)n5=node(-10,-2.75,4.255+0.35) cp,next,all,n1,n4n1=node(-10,2.72,4.255)n2=node(-10.5,2.75,4.255)n3=node(-9.5,2.75,4.255)n4=node(-10,2.75,4.255)n5=node(-10,2.75,4.255+0.35) cp,next,all,n1,n4n1=node(10,-2.72,4.255)n2=node(10.5,-2.75,4.255)n3=node(9.5,-2.75,4.255)n4=node(10,-2.75,4.255)n5=node(10,-2.75,4.255+0.35) cp,next,all,n1,n4n1=node(10,2.72,4.255)n2=node(10.5,2.75,4.255)n3=node(9.5,2.75,4.255)n4=node(10,2.75,4.255)n5=node(10,2.75,4.255+0.35) cp,next,all,n1,n4!H1n1=node(-20,-2.72,4.274)n2=node(-20.5,-2.75,4.274-0.35) n3=node(-20.5,-2.75,4.274+0.35) n4=node(-19.5,-2.75,4.274-0.35) n5=node(-19.5,-2.75,4.274+0.35) n6=node(-20,-2.75,4.274)n7=node(-20,-2.75,4.274+0.35) cp,next,all,n1,n6n1=node(-20,2.72,4.274)n2=node(-20.5,2.75,4.274-0.35) n3=node(-20.5,2.75,4.274+0.35) n4=node(-19.5,2.75,4.274-0.35) n5=node(-19.5,2.75,4.274+0.35) n6=node(-20,2.75,4.274)n7=node(-20,2.75,4.274+0.35) cp,next,all,n1,n6n1=node(20,-2.72,4.274)n2=node(20.5,-2.75,4.274-0.35) n3=node(20.5,-2.75,4.274+0.35) n4=node(19.5,-2.75,4.274-0.35) n5=node(19.5,-2.75,4.274+0.35) n6=node(20,-2.75,4.274)n7=node(20,-2.75,4.274+0.35)cp,next,all,n1,n6n1=node(20,2.72,4.274)n2=node(20.5,2.75,4.274-0.35)n3=node(20.5,2.75,4.274+0.35)n4=node(19.5,2.75,4.274-0.35)n5=node(19.5,2.75,4.274+0.35)n6=node(20,2.75,4.274)n7=node(20,2.75,4.274+0.35)cp,next,all,n1,n6! 施加约束nsel,s,loc,z,5,50D,node(-40,-1.60,0), , , , , ,UY,UZ, , , ,D,node(-40,1.60,0), , , , , ,Uz,, , , ,D,node(-30,-1.60,0), , , , , ,UY,UZ, , , ,D,node(-30,1.60,0), , , , , ,Uz,, , , ,D,node(30,-1.60,0), , , , , ,UY,UZ, , , ,D,node(30,1.60,0), , , , , ,Uz,, , , ,D,node(40,-1.60,0), , , , , ,UY,UZ, , , ,D,node(40,1.60,0), , , , , ,Uz,, , , ,nsel,s,loc,x,-25.366,-50nsel,r,loc,z,-10,0.256D,all, , , , , ,all,, , , ,nsel,s,loc,x,25.366,50nsel,r,loc,z,-10,0.256D,all, , , , , ,all,, , , ,allsel!修改桥面单元,使之重量包括铺装层重量asel,s,loc,z,6.133,6.135asel,r,loc,y,2.73,-2.73ALLSEL,BELOW,AREAEMODIF,all,MA T,3,allsel3影响线求解(静力法)!在给一排单元上施加-1000kN/m^2的均布荷载,每排单元在纵向上的长度为2.5/7m。