ANSYS二次开发概述
标准ANSYS程序是一个功能强大、通用性好的有限元分析程序，同时它还具有良好的开放性，用户可以根据自身的需要在标准ANSYS版本上进行功能扩充和系统集成，生成具有行业分析特点和符合用户需要的用户版本的ANSYS程序。开发功能包括四个组成部分：
⑴．参数化程序设计语言（APDL）
⑵．用户界面设计语言（UIDL）
⑶．用户程序特性（UPFs）
⑷．ANSYS数据接口
APDL所能实现的功能通俗的说来应该是次于UPF而强与UIDL，但实际上是由于三者具体侧重点不同造成的：UIDL主要控制GUI界面的各类二次开发方法，涉及的分析部分就要少一些，APDL可以称其为和分析部分频繁打交道的一组小型工具，功能强大，但不和UIDL一样能够非常具体的针对某一两方面的二次开发处理，通常情况下融合在分析的角角落落中。UPF是三者之间的最强者，能完成最复杂的二次开发工作，比如说构建新单元，复杂数据库交互，外围命令定制等，但UPF在很多情况下也借助了APDL命令来完全实现其功能。同样也能在UIDL中嵌入APDL命令，来构建比较复杂的GUI二次开发工作。
UIDL、APDL和UPF三者各有所长，密不可分。结合使用三者，就能够实现任何强大的分析功能。
5.2 Ansys的开发功能组成部分Ansys的开发功能由三个部分组成：参数化程序设计语言（APDL）、用户界面设计语言（UIDL）、用户程序特性（UPFs）
5.2.1 参数化程序设计语言（APDL）
参数化程序设计语言(APDL-ANSYS Parametric Design Language)实质上由类似于FORTRAN77的程序设计语言部分和1000多条ANSYS命令组成。其中，程序设计语言部分与其它编程语言一样，具有参数、数组表达式、函数、流程控制（循环与分支）、重复执行命令、缩写、宏以及用户程序等。标准的ANSYS 程序运行是由1000多条命令驱动的，这些命令可以写进程序设计语言编写的程序，命令的参数可以赋确定值，也可以通过表达式的结果或参数的方式进行赋值。从ANSYS命令的功能上讲，它们分别对应ANSYS分析过程中的定义几何模型、划分单元网格、材料定义、添加载荷和边界条件、控制和执行求解和后处理计算结果等指令。很多情况下，APDL主要用在优化设计或者自适应网格划分中。
用户可以利用程序设计语言将ANSYS命令组织起来，编写出参数化的用户程序，从而实现有限元分析的全过程，即建立参数化的CAD模型、参数化的网格划分与控制、参数化的材料定义、参数化的载荷和边界条件定义、参数化的分析控制和求解以及参数化的后处理。
宏是具有某种特殊功能的命令组合，实质上是参数化的用户小程序，可以当作ANSYS的命令处理，可以有输入参数或没有输入参数。
缩写是某条命令或宏的替代名称，它与被替代命令或宏存在一一对应的关系，在ANSYS中二者是完全等同的，但缩写更符合用户习惯，更易于记忆，减少敲击键盘的次数。ANSYS工具条就是一个很好的缩写例子。
5.2.2 用户界面设计语言（UIDL）
标准ANSYS交互图形界面可以驱动ANSYS命令，提供命令的各类输入参数接口和控制开关，用户在图形驱动的级别上进行有限元分析，整个过程变得直观轻松。用户图形界面设计语言（UIDL-User Interface Design Language）就是编写或改造ANSYS图形界面的专用设计语言，GUI方面几乎全部的二次开发功能都将由它完成,主要完成以下三种图形界面的设计：
⑴．主菜单系统及菜单项
⑵．对话框和拾取对话框
⑶．帮助系统
通过用户界面设计语言（UIDL），用户可以在扩充ANSYS功能的同时建立起对应的图形驱动界面，如在主菜单的某位置增加菜单项，设计对应的对话框、拾取对话框，实现参数的输入和其它程序运行的控制，同时提供相应的联机帮助，使操作者能方便地获取系统帮助。主要功能有：
[1] 组织强大的菜单系统。
在Ansys中也能轻松做出可以和VC，VB之类主流GUI开发工具媲美的菜单响应效果。
[2] 构建功能繁复的对话框。
利用强大的UIDL工具，也能轻松架构起进行大型工程分析的实用对话框向导。
[3] 建立自己的联机帮助
Ansys中的联机帮助非常实用，构建自己完善的帮助系统，UIDL是这方面不可或缺的理想开发工具
一个完整的UIDL控制文件结构大致如图5-1：
控制文件头
结构块结构
图5-1 UIDL控制文件结构图
任何一个UIDL控制文件开头都是一个控制文件头，其后接一个或多个结构块结构。控制文件头第一行必须有:F filename，filename是UIDL控制文件名。控制文件头第二行必须有:D description，description是对本文件的一些说明。注意到description中有时能带%E%扩展，但只有当拥有类似SCCS的系统（含一
源码控制系统），ANSYS才能有效的进行%E%扩展，否则只有手动把这些说明替代%E%写入description中。
控制文件第三行必须有:I 0, 0, 0，各个0必须出现在第9、18、27行。只要在这些位置填入0，ANSYS在调用该文件后会自动在这些位置填入GUI界面的位置信息。
结构块结构是一个UIDL控制文件的核心，它涵盖了菜单信息，命令信息，以及帮助文件信息，按照其不同的类型可划分为菜单结构块，命令结构块和帮助结构块。一般来说函数结构块还都伴随着构建一个对话框结构。
结构块结构基本框架如图5-2。
头部分
数据控制部分
尾部分
图5-2 UIDL结构块结构图
ANSYS在调用GUI界面时会自动调用menulistxx.ans文件（由版本不同而改变），该文件中描述了UIDL前处理器到哪里去寻找UIDL控制文件。ANSYS 在其docu/目录中有一个基本的menulistxx.ans文件和对应的基本UIDL控制文件。默认情况下，ANSYS就调用这一文件。
5.2.3用户程序特性（UPFs）
用户程序特性（UPFs）向用户提供丰富的FORTRAN77用户程序开发子程序和函数，用户利用它们从开发程序源代码的级别上扩充ANSYS的功能。使用这些子程序和函数，编写用户功能的源代码程序，在与ANSYS版本要求匹配的FORTRAN或C++编译器上重新编译和连接，生成用户版本的ANSYS程序。另外，还提供了外部命令功能，允许用户创建ANSYS可以利用的共享库。用户可以开发下列方面的功能程序：
⑴．开发用户子程序实现从ANSYS数据库中提取数据或将数据写入ANSYS 数据库。该种子程序可以编译连接到ANSYS中，此时ANSYS提供了10个数据库操作命令；如果作为外部命令处理，可以在ANSYS的任何模块中运行；
⑵．利用ANSYS提供的子程序定义各种类型的载荷，其中包括BF或BFE 载荷、压力载荷、对流载荷、热通量和电荷密度等；
⑶．利用ANSYS提供的子程序定义各种材料特性，包括塑性、蠕变、膨胀、粘塑性、超弹、层单元失效准则等；
⑷．利用ANSYS提供的子程序定义新单元和调整节点方向矩阵；
⑸．利用ANSYS提供的子程序修改或控制ANSYS单元库中的单元；
(6).利用UEROP创建用户优化程序；ANSYS程序作为子程序在用户程序中调用。
一个典型的UPFs包括下列步骤：
1.在FORTRAN77中编制用户程序。在ANSYS中所有的用户程序源代码都是公开的。大部分完成至少一个简单的功能，因此在编制程序前应列出一份完整的可用程序表。
2.编译并将用户程序连接到ANSYS程序中。在ANSYS Installation and Configuration Guide中有这方面的详细描述。
3.用户可能要验证自己做的改动是否影响其他ANSYS标准功能的使用。可以通过做几个ANSYS Verification Manual中的例题来验证。在ANSYS中也有输入这些问题的方法。
4.用用户认为可以满足要求的过程检验用户子程序。请牢记本步由用户完成，用户要对子程序负责。
ANSYS程序将在用户使用时自动激活一些子程序，如用户单元等。例如，要激活一个用户单元，仅仅在模型建立时指定它为其中的一个单元类型（ET命令或MainMenu>Preprocessor>Element Type>Add/Edit/Delete），将单元类型属性指针设置好，然后用该单元进行实体划分[AMESH,VMESH等]或直接生成单元[ET等命令]即可。
对于其他一些UPFs，用户需要做一些别的操作来激活它们。即，需要输入USRCAL命令或选择下列菜单：Main Menu>Preprocessor>Loads>-Load Step Opts-Other>User Routines或Main Menu>Solution>-Load Step Opts-Other>User Routines。如果没有做这一步，将执行标准的ANSYS功能。例如，如果想施加对流载荷，即便连接了用户对流子程序，缺省时仍是使用标准的ANSYS功能。因此用户必须使用USRCAL命令（或相应路径）激活正确的用户子程序才能使用用户过程。请参阅USRCAL命令描述得到命令可以影响的用户过程列表。使用NSVR命令定义在用户过程中要存储的另外的变量数目。（NSVR命令无相应的路径）。
另一个有用的命令是/UCMD，可以使用户在用户过程之外生成自己的命令。（/UCMD命令无GUI路径。）假定用户连接抛物线形分布压力的子程序。如果子程序名为USERnn（nn=1到10），可以用自己的命令调用程序：
/UCMD,PARAB,1
PARAB将成为合法的ANSYS命令，用来调用用户子程序USER01。可以用命令调用10个子程序。将/UCMD命令嵌入启动文件（START.ANS）中，用户就可以在所有ANSYS过程中使用用户子程序。
第24卷第5期辽宁工学院学报V o l.24 N o.5 2004年10月JOU RNAL O F L I AON I N G I N ST ITU T E O F T ECHNOLO GY O ct.2004① 基于AN SYS的二次开发技术的实现方法 吴 鹏1,曾 红1,韩 迈2 (1.辽宁工学院,辽宁锦州 121001;2.鞍山广播电视大学,辽宁鞍山 114000) 摘 要:基于大型通用有限元分析软件AN SYS8.0环境,对AN SYS二次开发技术进行了探讨,并对AN SYS 三种开发工具进行了详细的介绍。论述了采用二次开发方法设计产品的必要性和重要性,证实了以AN SYS为平台开发专业模块的可行性,提高了工作效率,缩短了产品的开发研制周期。 关键词:AN SYS;二次开发;A PDL;U I DL;U PF s 中图分类号:T P391.72 文献标识码:B 文章编号:100521090(2004)0520025205 Realization of Secondary D evelop m en t of Technology Based on ANS Y S W U Peng1,ZEN G Hong1,HAN M ai2 (1.L iaoning Institute of T echno logy,J inzhou121001,Ch ina;2.A nshan R adi o&TV U niversity,A nshan114000,Ch ina) Key words:AN SYS;Secondary developm en t;A PDL;U I DL;U PF s Abstract:T he m ethod of secondary developm en t of techno logy on the basis of large-scale fin ite elem en t analysis softw are—AN SYS is described and app roached,w h ich details th ree k inds of de2 velop ing too ls of AN SYS.It dem on strates the necessity and i m po rtance of the m ethod of sec2 ondary developm en t of techno logy.T he feasib ility of develop ing p rofessi onal m odu le on the AN2 SYS p latfo r m is verified,w o rk ing efficiency i m p roved,and the developm en t cycle of the p roducts sho rtened. 从20世纪70年代以来,随着计算技术的飞速发展,结构分析有了很大的突破,国外相继出现了许多大型通用有限元分析程序,如AN SYS, ABAQU S,M A RC和M SC NA STRAN等,这些程序具有良好的界面、方便的前后处理和强大的计算分析功能以及开放的二次开发系统。 AN SYS软件是融热、电、磁、流体、结构、声学于一体的大型通用有限元分析软件。具有强大的求解器和前、后处理功能,为解决复杂、庞大的工程项目提供了一个强有力的工具。然而,正是由于AN2 SYS的通用性特点,使其对不同行业的专业性模块的分析不具有针对性,复杂的英文界面和繁琐的分析步骤都给从事有限元分析的技术人员造成了很大的障碍。另外,虽然AN SYS有较强大的前、后处理功能,但使用者必须具有较高的相关力学知识和丰富的分析经验,在几何建模简化和力学建模等前处理方面需要花费很多时间和精力。因此,基于这些不便因素,在熟练应用AN SYS软件的基础上,结合具体各行业的实践经验,利用AN SYS内部提供的二次开发工具,用户可在AN SYS系统中开发出具有中文界面的、特定功能的专用模块,可以有效地提高设计的效率和质量,充分体现了专业化、用户化、便 ①收稿日期:2004206228 基金项目:辽宁省教育厅科研资助项目(20032086)作者简介:吴鹏(19792),男,辽宁盘锦人,硕士生。
万方数据
第３期张庆峰等：基于ＡＮＳＹＳ二次开发的管系结构应力分析系统—－７９—．大，计算结果可靠。但它要求使用者具有一定的有限元知识背简单。 景，并同时具有较强专业知识水平、较强的结构分析能力和扎实 的英语基础。鉴于上述特点，使其对压力管系的有限元分析不 具有针对性。复杂的英文界面和繁琐的分析步骤都给从事压力 管系有限元分析的技术人员造成了很大的障碍。因此，基于这 些不便因素，为适用不同层次的用户使用，利用ＡＮＳＹＳ内部提 供的二次开发工具。把ＡＮＳＹＳ作为结构分析工具，建立了特别 适用于结构应力分析的中文界面环境、菜单和工具杆的管系结 构分析系统模块。此模块以向导的方式来进行每一步骤，各步 骤附有帮助文件，充分体现了专业化、用户化、便捷化的特点。 如图１所示。 图３管系图 图１绘制管系图 ４应用实例 利用在役压力管道系统的应力分析模块对某厂核反应器再循环装置管线进行应力分析，如图２所示。 图２核反应器再循环装置回路管线图 ４．１核反应器再循环装置回路管线概况 下面是一个应用该软件对在役核压力回路管线进行应力结构分析的简例。如图３所示，假定核反应器再循环装置的回路管线中发现了二处裂纹。这些裂纹可能是由于在生产或制造过程没有操作经验或某种晶间应力腐蚀所引起的。这两个裂纹，①和②，存在于旁路与核反应再循环装置回路管线主管路相连的焊接部位，它们可认为是复合缺陷。旁路管线的内径是２８２ｍｍ，主管路的内径是４５０ｒａｍ，厚度是３１．７６ｍｍ。这些管路和弯管是ＳＡ３３３ＧＲ６型材料，弹性模量是１８８ＧＰａ。 ４．２管系的结构分析 借助ＡＮＳＹＳ的二次开发功能，在开发“含缺陷压力管系风险分析系统”时。在结构应力分析模块中，选择了国际著名的ＡＮＳＹＳ有限元分析软件作为结构应力分析工具，并为适用不同层次的用户的需要，针对ＡＮＳＹＳ的管路系统模块的特征，对ＡＮＳＹＳ进行了二次开发，建立了专用程序的同时建立起对应的图形驱动界面，使得前处理建模、计算和后处理操作等变得十分 图４管系应力分布云图 ５结论 通过开发以ＡＮＳＹＳ为平台的管系应力分析系统，证实了运用ＡＮＳＹＳ内部提供的ＡＰＤＬ语言和ＵＩＤＬ语言进行开发专业模块的可行性，并且达到了界面简洁、易操作的预期功能。 利用建立在ＡＮＳＹＳ二次开发基础上强大的管道结构应力分析模块，可以在制定管道的检修计划时，方便地确定出管道高度应集中部位，有针对性地选择焊缝并进行射线探伤，使管线的安全状况分析更加准确。有针对性地选择焊缝并进行射线探伤，使得管道的安全状况分析更为准确。同时，也可以利用该系统为分析工具，制定出旨在降低失效风险的管道结构改进措施，优化管道结构。以较低的成本提高管道的完整性水平。 因此，该系统的推广应用，对提高企业的压力管道管理水平，保障安全生产和技术进步具有重要意义。 参考文献 １ＡＮＳＹＳＡＰＤＬＰｒｏｇｒａｍｍｅｒ’ｓＧｕｉｄｅＲｅｌｅａｓｅ５．５．ＡＮＳＹＳ。Ｉｎｅ． ２ＴｈｅＵＩＤＬＰｒｏｇｒａｍｍｅｒ’ｓＧｕｉｄｅＲｅｌｅａｓｅ５．５．ＡＮＳＹＳ．Ｉｎｅ． ３谢禹钧，蔺永诚，等．含缺陷压力关系失效风险分析系统（Ｉ）【Ｊ】．石油化工设备，２００２，３１（４）：４—６． ４谢禹钧，蔺永诚等．含缺陷压力关系失效风险分析系统（ｎ）【Ｊ】．石油化工设备，２００２，３１（５）：４～６． ５程进，江见鲸等．基于ＡＮＳＹＳ的程序界面设计及应用。四川建筑科学研究。２００２，２８． ６沈士明，在役压力管道安全评定研究的现状与发展。中国机械工程。 １９９７．８． ７ＡＰＤＬ参数化有限元分析技术及其应用实例，中国水利水电出版社， ２００３． 万方数据
ANSYS的二次开发技术 ANSYS 的二次开发技术 ANSYS 提供的二次开发工具有三个：参数化设计语言(ANSYS Parametric Design Language，APDL)，用户界面设计语言(User Interface Design Language，UIDL)以及用户可编程特性(User Programmable Features，UPFs)。其中，前两种可归类为标准使用特性，后一种为非标准使用特性。 ANSYS 参数化设计语言(APDL)APDL 扩展了传统有限元分析范围之外的能力，提供了建立标准化零件库、序列化分析、设计修改、设计优化以及更高级 的数据分析处理能力，包括灵敏度研究等。ANSYS 用户可编程特性(UPFs)利 用UPFs，用户可以开发下列方面的功能程序：(1) 开发用户子程序实现从ANSYS 数据库中提取数据或将数据写入ANSYS 数据库。该种子程序可以编译 连接到ANSYS 中，此时ANSYS 提供了10 个数据库操作命令；如果作为外部 命令处理，可以在ANSYS 的任何模块中运行；(2) 利用ANSYS 提供的子程序 定义各种类型的载荷，其中包括BF 或BFE 载荷、压力载荷、对流载荷、热通 量和电荷密度等；(3) 利用ANSYS 提供的子程序定义各种材料特性，包括塑性、蠕变、膨胀、粘塑性、超弹、层单元失效准则等；(4) 利用ANSYS 提供 的子程序定义新单元和调整节点方向矩阵，ANSYS 最多可以有6 个独立的新 单元USER100-USER105；( 5) 利用ANSYS 提供的子程序修改或控制ANSYS 单元库中的单元；(6) 利用UEROP 创建用户优化程序，可以用自己的算法和 中断准则替换ANSYS 优化过程。(7) ANSYS 程序作为子程序在用户程序中调用，如用户自定义的优化算法。ANSYS 软件本身是通过FORTRAN 和C 语言 开发的。使用UPFs 进行二次开发，在安装ANSYS 的基础上，还需要Compaq Visual FORTRAN 和MS Visual C++的支持。
有限元 学院：机电学院 专业： 姓名： 学号：
一、问题描述 如图所示的平面，板厚为0.01m，左端固定，右端作用50kg的均布载荷，对其进行静力分析。弹性模量为210GPa，泊松比为0.25. 二、分析步骤 1．启动ansys，进入ansys界面。 2．定义工作文件名 进入ANSYS/Multiphsics的的程序界面后，单击Utility Menu菜单下File中Change Jobname的按钮，会弹出Change Jobname对话框，输入gangban为工作文件名，点击ok。 3．定义分析标题 选择菜单File-Change Title在弹出的对话框中，输入Plane Model作为分析标题，单击ok。 4．重新显示 选择菜单Plot-Replot单击该按钮后，所命令的分析标题工作文件名出现在ANSYS 中。 5．选择分析类型 在弹出的对话框中，选择分析类型，由于此例属于结构分析，选择菜单Main Menu：Preferences，故选择Structural这一项，单击ok。 6．定义单元类型 选择菜单Main Menu-Preprocessor-Element Type-Add/Edit/Delete单击弹出对话框中的Add按钮，弹出单元库对话框，在材料的单元库中选Plane82单元。即在左侧的窗口中选取Solid单元，在右侧选择8节点的82单元。然后单击ok。 7．选择分析类型 定义完单元类型后，Element Type对话框中的Option按钮被激活，单击后弹出一个对话框，在Elenment behavior中选择Plane strs w/ thk，在Extra Element output 中，选择Nodal stress，单击close，关闭单元类型对话框。 8．定义实常数 选择菜单Main Menu-Preprocessor-Real Constants Add/Edit/Delete执行该命令后，在弹出Real Constants对话框中单击Add按钮，确认单元无误后，单击ok，弹出Real Constants Set Number 1，for Plane 82对话框，在thickness后面输入板的厚度0.01单击ok，单击close。 9．定义力学参数 选择菜单Main Menu-Preprocessor-Material Props-Material Model 在弹出的对
ansys二次开发教程+实例 第3章ANSYS基于VC++6.0的二次开发与相互作用分析在ANSYS中的实现 3.1 概述 ANSYS是一套功能十分强大的有限元分析软件，能实现多场及多场耦合分析；是实现前后处理、求解及多场分析统一数据库的 一体化大型FEA软件；支持异种、异构平台的网络浮动，在异种、异构平台上用户界面统一、数据文件全部兼容，强大的并行计算功能 支持分布式并行及共享内存式并行。该软件具有如下特点： (1) 完备的前处理功能 ANSYS不仅提供了强大的实体建模及网格划分工具，可以方便地构造数学模型，而且还专门设有用户所熟悉的一些大型通用有 限元软件的数据接口（如MSC／NSSTRAN，ALGOR，ABAQUS等），并允许从这些程序中读取有限元模型数据，甚至材料特性和边 界条件，完成ANSYS中的初步建模工作。此外，ANSYS还具有近200种单元类型，这些丰富的单元特性能使用户方便而准确地构建出 反映实际结构的仿真计算模型。 (2) 强大的求解器 ANSYS提供了对各种物理场量的分析，是目前唯一能融结构、热、电磁、流体、声学等为一体的有限元软件。除了常规的线性、 非线性结构静力、动力分析外，还可以解决高度非线性结构的动力分析、结构非线性及非线性屈曲分析。提供的多种求解器分别适用于 不同的问题及不同的硬件配置。 (3) 方便的后处理器 ANSYS的后处理分为通用后处理模块（POST1）和时间历程后处理模块（POST26）两部分。后处理结果可能包括位移、温度、应力、应变、速度以及热流等，输出形式可以有图形显示和数据列表两种。 (4) 多种实用的二次开发工具 ANSYS除了具有较为完善的分析功能外，同时还为用户进行二次开发提供了多种实用工具。如宏（Marco）、参数设计语言（APDL）、用户界面设计语言（UIDL）及用户编程特性（UPFs），其中APDL（ANSYS Parametric Design Language）是一种非常类似于Fortran77的参数化设计解释性语言，其核心内容为宏、参数、循环命令和条件语句，可以通过建立参数化模型来自动完成一些通用性强的任务；UIDL（User Interf ace Design Language）是ANSYS为用户提供专门进行程序界面设计的语言，允许用户改变ANSYS的图形用户界面(GUI)中的一些组项，提供了一种允许用户灵活使用、按个人喜好来组织设计ANSYS图形用户界面的强有力工具；UPFs(User Programmable Features)提供了一套Fortran77函数和例程以扩展或修改程序的功能，该项技术充分显示了ANSYS的开放体系，用户 不仅可以采用它将ANSYS程序剪裁成符合自己所需的任何组织形式（如可以定义一种新的材料，一个新的单元或者给出一种新的屈服 准则），而且还可以编写自己的优化算法，通过将整个ANSYS作为一个子程序调用的方式实现。 鉴于上述特点，近几年来，ANSYS软件在国内外工程建设和科学研究中得到了广泛的应用。但这些应用大多局限于直接运用ANSYS软件进行实际工程分析，对利用ANSYS提供的二次开发工具进行有限元软件设计却很少涉及。本文首次利用ANSYS软件的二次开发功能，以VC++6.0为工具，运用APDL语言，对ANSYS进行二次开发，编制框筒结构－桩筏基础－土相互作用体系与地震反应分析程序。 3.2 程序设计目标 针对某一实际工程问题，ANSYS所提供的APDL语言可对ANSYS软件进行封装。APDL语言即ANSYS软件提供的参数化设计 语言，它的全称是ANSYS Parametric Design Language。使用APD L语言可以更加有效地进行分析计算，可以轻松地进行自动化工作（循环、分支、宏等结构），而且，它是一种高效的参数化建模手段。使用APDL语言进行封装的系统可以只要求操作人员输入前处理 参数，然后自动运行ANSYS进行求解。但完全用APDL编写的宏还存在弱点。比如用APDL语言较难控制程序的进程，虽然它提供了 循环语句和条件判断语句，但总的来说还是难以用来编写结构清晰的程序。它虽然提供了参数的界面输入，但功能还不是太强，交互性 不够流畅。针对这种情况，本文用VC++6.0开发框筒结构－桩筏基础－土相互作用有限元分析程序（简称LW S程序）。
１．ＡＮＳＹＳ分析的原理和步骤 ＡＮＳＹＳ的热分析［１］包括稳态和瞬态两种，如果系统的温度场与时间无关，则称该系统处于稳定的热状态，简称稳态；如果系统的温度场随时间发生变化，则称系统处于瞬态。显然，大体积混凝土的浇筑过程属于瞬态分析，也属于非线性分析。 我们不仅要进行混凝土温度场的模拟还要进行应力场的模拟，所以要用到ＡＮＳＹＳ中耦合分析，ＡＮＳＹＳ提供了两种分析耦合场的方法：直接耦合与间接耦合。 直接耦合法的耦合单元包含所有必须的自由度，仅仅通过一次求解就能得出耦合场分析结果；间接耦合法是以特定的顺序求解单个物理场的模型，通过把第一次场分析的结果作为第二次场分析的载荷来实现两种场的耦合。如我们用到的热－应力耦合分析就是将热分析得到的节点温度作为载荷施加在后序的应力分析中来实现耦合的。基本步骤如下： 第一步：进行热分析，可选择ＳＯＬＩＤ７０单元； 第二步：重新进入前处理器，转换单元类型；将热单元转换为相应的结构单元，原来的ＳＯＬＩＤ７０单元将自动转换为ＳＯＬＩＤ４５单元，其对应的命令是ＥＴＣＨＧ，ＴＴＳ。 第三步：设置结构分析中的材料属性； 第四步：读入热分析结果并将其作为载荷；可采用命令ＬＤＲＥＡＤ读入热分析的节点温度，或点击ＭａｉｎＭｅｎｕ＞Ｓｏｌｕｔｉｏｎ＞ＬｏａｄＡｐｐｌｙ＞Ｔｅｍｐｅｒａｔｕｒｅ＞ＦｒｏｍＴｈｅｒｍａｌＡｎａｌｙｓｉｓ。注意，结果文件的扩展名为＊．ｒｔｈ。 第五步：指定参考温度；在参考温度处，热应力值为零。 第六步：求解及后处理。 ２．温度场的求解 ２．１三种基本传热方式 （１）热传导，遵循傅里叶定律（导热基本定律）：ｑ″＝－λｄＴ ｄｘ ，式中ｑ″为热流密度（Ｗ／ｍ２），λ为导热系数（Ｗ／ｍ?℃），“－”表示热量流向温度降低的方向。 （２）热对流，用牛顿冷却方程来描述：ｑ″＝β（ＴＳ－ＴＢ），式中β为对流换热系数，ＴＳ为固体表面的温度，ＴＢ为周围流体的温度。 （３）热辐射，指物体发射电磁能，并被其它物体吸收转变为热的热量交换过程。 ２．２边界条件 （１）第一类边界条件是指混凝土表面温度Ｔ是时间τ的已知函数，即 Ｔ（ｘ，ｙ，ｚ，τ）＝Ｔｂ（τ） （２）第二类边界条件是指混凝土表面的热流量是时间的已知函数，即 －λ$Ｔ $ｎ ＝Ｔ′（τ） 式中λ—— —导热系数，Ｗ／ｍ?℃或ｋＪ／ｍ?ｈ?℃，Ｗ／ｍ?℃＝３．６ｋＪ／ｍ?ｈ?℃； ｎ—— —表面外法线方向，若表面是绝热的，有：$Ｔ $ｎ ＝０。 （３）第三类边界条件假定经过混凝土表面的热流量与混凝土表面温度Ｔ和气温Ｔａ之差成正比，即 －λ$Ｔ $ｎ ＝β（Ｔ－Ｔａ） 式中β—— —表面放热系数，也称对流系数，Ｗ／ｍ２?℃。其数值与风速ｖａ（ｍ／ｓ）有密切的关系，固体表面在空气中的放热系数可用以下两式计算，单位是ｋＪ／ｍ２?ｈ?℃。 粗糙表面：β＝２３．９＋１４．５０ｖａ（１）光滑表面：β＝２１．８＋１３．５３ｖａ（２）当有模板和保温层时，可按下式计算：β＝１ ∑ δ ｉ λ ｉ ＋１ β ｑ （３）式中δｉ—— —各种保温材料的厚度（ｍ）； λｉ—— —各种保温材料的导热系数（Ｗ／ｍ?Ｋ），可按表１取值［２］； βｑ—— —空气的传热系数，可取２３（Ｗ／ｍ２?Ｋ）。 表１各种保温材料的导热系数λ值（Ｗ／ｍ?Ｋ） （４）当两种条件不同的固体接触时，如接触良好，则在接触面上温度和热流量都是连续的，即Ｔ１＝Ｔ２，λ１（ !Ｔ １ !ｎ ）＝λ２（!Ｔ２ !ｎ ）。 混凝土与空气接触（包括有养护条件）的边界可按照第三类边界条件处理： ＮＳＥＬ，，，！选择与空气接触的表面节点 ＳＦ，ＡＬＬ，ＣＯＮＶ，β，Ｔａｉｒ，！加载表面散热系数和环境温度 混凝土与地基或基岩的边界可以按照第四类边界条件处理，通过定义两种材料的导热系数和初始温度即可。 ２．３热学参数取值基本参数较容易获得，也可参考下表： 表２材料的基本热学参数 ２．３．１水化热的施加在ＡＮＳＹＳ中，混凝土的水化热是通过生热率ＨＧＥＮ来施加的。顾名思义，生热率就是单位时间内混凝土的生热量，即所产生的热量对时间的导数，用表达式表示为： ｈｇｅｎ＝ｄＱ ｄｔ （４）式中：Ｑ—— —混凝土中产生的热量； ｈｇｅｎ—— —混凝土生热率。 混凝土的水化热放热过程与混凝土的绝热温升过程具有一致性，若取指数经验式： ＡＮＳＹＳ模拟大体积混凝土浇筑过程的参数分析 赵英菊王社良康宁娟 （西安建筑科技大学土木工程学院陕西西安７１００５５） 摘要：建筑工程中的大体积混凝土结构越来越多，利用有限元程序ＡＮＳＹＳ进行施工过程的模拟仿真可以形象地给出温度场和应力场的分布情况，同时能考虑各参数随时间的变化。时变参数的选取及其在程序中的实现是仿真分析中的重点和难点，特总结归纳，并给出解决的方法供参考。 关键词：ＡＮＳＹＳ；混凝土；浇筑；时变参数 材料名称λ材料名称λ 木模０．２３黏土砖０．４３ 钢模５８油毡０．０５ 草袋０．１４沥青矿棉０．０９～０．１２ 木屑０．１７沥青玻璃棉毡０．０５ 矿渣０．４７泡沫塑料制品０．０３～０．０５ 黏土１．３８～１．４７泡沫混凝土０．１０ 干砂０．３３水０．５８ 湿砂１．３１空气０．０３ 名称数值单位名称数值单位 混凝土的密度２４００ｋｇ／ｍ３混凝土的导热系数２．７１０Ｗ／ｍ?℃ 土壤的密度１７５０ｋｇ／ｍ３土壤的导热系数０．５８６Ｗ／ｍ?℃ 混凝土的比热０．９６３ｋＪ／ｋｇ?℃混凝土的线膨胀系数１０×１０－６℃ 土壤的比热１．００５ｋＪ／ｋｇ?℃混凝土的导温系数０．００４２ｍ２／ｈ９６
基于VB语言的EXCEL、CST以及HFSS等的二次开发 代码1：vb创建xls表，并写入内容 Set ExcelApp = CreateObject("Excel.Application") '创建EXCEL对象Set ExcelBook = ExcelApp.Workbooks.Add Set ExcelSheet = ExcelBook.Worksheets(1) '添加工作页ExcelSheet.Activate '激活工作页 ExcelApp.DisplayAlerts = False https://www.doczj.com/doc/d02794745.html,="sheet1" ExcelSheet.Range("A1").Value = 100 '设置A1的值为100 ExcelBook.SaveAs "d:\test.xls" '保存工作表 msgbox "d:\test.xls创建成功!" ExcelBook.close set excelApp=nothing set ExcelBook=nothing set ExcelSheet=nothing 将以上代码copy到记事本存为"writexls.vbs"文件，可运行测试 代码2：读execel文件 Set ExcelApp = CreateObject("Excel.Application") '创建EXCEL对象Set ExcelBook = ExcelApp.Workbooks.open("d:\test.xls") Set ExcelSheet = ExcelBook.Worksheets(1) msgbox ExcelSheet.Range("A1").Value 将以上代码copy到记事本存为"readxls.vbs"文件，可运行测试 代码3：上述代码联合调试 Dim ExcelApp,ExcelBook,ExcelSheet Set ExcelApp = CreateObject("Excel.Application") '创建EXCEL对象Set ExcelBook = ExcelApp.Workbooks.Add Set ExcelSheet = ExcelBook.Worksheets(1) '添加工作页ExcelSheet.Activate '激活工作页 ExcelApp.DisplayAlerts = False https://www.doczj.com/doc/d02794745.html,="sheet1" ExcelSheet.Range("A1").Value = 100 '设置A1的值为100 ExcelBook.SaveAs "D:\Study\VBS\Book1.xls" '保存工作表 msgbox "d:\Book1.xls创建成功!" ExcelBook.close set excelApp=nothing set ExcelBook=nothing set ExcelSheet=nothing 'ExcelApp.WorkBooks.Close 'ExcelApp.Quit
“有限元分析及应用”课程有限元分析软件ANSYS6.xed 上机指南 清华大学机械工程系 2002年9月
说明 本《有限元分析软件ANSYS6.1ed:上机指南》由清华大学机械工程系石伟老师组织编写，由助教博士生孔劲执笔, 于2002年9月完成，基本操作指南中的所有算例都在相应的软件系统中进行了实际调试和通过。 本上机指南的版权归清华大学机械工程系所有，未经同意，任何单位和个人不得翻印。
目录 Project1 简支梁的变形分析 (1) Project2 坝体的有限元建模与受力分析 (3) Project3 受内压作用的球体的应力与变形分析 (5) Project4 受热载荷作用的厚壁圆筒的有限元建模与温度场求解 (7) Project5 超静定桁架的有限元求解 (9) Project6 超静定梁的有限元求解 (11) Project7 平板的有限元建模与变形分析 (13)
Project1 梁的有限元建模与变形分析 计算分析模型如图1-1 所示, 习题文件名: beam。 NOTE:要求选择不同形状的截面分别进行计算。 梁承受均布载荷：1.0e5 Pa 图1-1梁的计算分析模型 梁截面分别采用以下三种截面（单位：m）： 矩形截面：圆截面：工字形截面： B=0.1, H=0.15 R=0.1 w1=0.1，w2=0.1，w3=0.2, t1=0.0114，t2=0.0114，t3=0.007 1.1进入ANSYS 程序→ANSYSED 6.1 →Interactive →change the working directory into yours →input Initial jobname: beam→Run 1.2设置计算类型 ANSYS Main Menu: Preferences →select Structural →OK 1.3选择单元类型 ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete… →Add… →select Beam 2 node 188 →OK (back to Element Types window)→Close (the Element Type window) 1.4定义材料参数 ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural→Linear→Elastic→Isotropic→input EX:2.1e11, PRXY:0.3→OK 1.5定义截面 ANSYS Main Menu: Preprocessor →Sections →Beam →Common Sectns→分别定义矩形截面、圆截面和工字形截面：矩形截面:ID=1,B=0.1，H=0.15 →Apply →圆截面：ID=2,R=0.1 →Apply →工字形截面：ID=3,w1=0.1，w2=0.1，w3=0.2，t1=0.0114，t2=0.0114，t3=0.007→OK
ansys二次开发 1概述 ANSYS是一套功能十分强大的有限元分析软件，能实现多场及多场耦合分析；是实现前后处理、求解及多场分析统一数据库的一体化大型FEA软件；支持异种、异构平台的网络浮动，在异种、异构平台上用户界面统一、数据文件全部兼容，强大的并行计算功能支持分布式并行及共享内存式并行。该软件具有如下特点：(1)完备的前处理功能 ANSYS不仅提供了强大的实体建模及网格划分工具，可以方便地构造数学模型，而且还专门设有用户所熟悉的一些大型通用有限元软件的数据接口（如MSC／NSSTRAN，ALGOR，ABAQUS等），并允许从这些程序中读取有限元模型数据，甚至材料特性和边界条件，完成ANSYS中的初步建模工作。此外，ANSYS还具有近200种单元类型，这些丰富的单元特性能使用户方便而准确地构建出反映实际结构的仿真计算模型。 (2)强大的求解器 ANSYS提供了对各种物理场量的分析，是目前唯一能融结构、热、电磁、流体、声学等为一体的有限元软件。除了常规的线性、非线性结构静力、动力分析外，还可以解决高度非线性结构的动力分析、结构非线性及非线性屈曲分析。提供的多种求解器分别适用于不同的问题及不同的硬件配置。 (3)方便的后处理器 ANSYS的后处理分为通用后处理模块（POST1）和时间历程后处理模块（POST26）两部分。后处理结果可能包括位移、温度、应力、应变、速度以及热流等，输出形式可以有图形显示和数据列表两种。 (4)多种实用的二次开发工具 ANSYS除了具有较为完善的分析功能外，同时还为用户进行二次开发提供了多种实用工具。如宏（Marco）、参数设计语言（APDL）、用户界面设计语言（UIDL）及用户编程特性（UPFs），其中APDL（ANSYS Parametric Design Language）是一种非常类似于Fortran77的参数化设计解释性语言，其核心内容为宏、参数、循环命令和条件语句，可以通过建立参数化模型来自动完成一些通用性强的任务；UIDL（User Interface Design Language）是ANSYS为用户提供专门进行程序界面设计的语言，允许用户改变ANSYS的图形用户界面(GUI)中的一些组项，提供了一种允许用户灵活使用、按个人喜好来组织设计ANSYS图形用户界面的强有力工具；UPFs(User Programmable Features)提供了一套Fortran77函数和例程以扩展或修改程序的功能，该项技术充分显示了ANSYS的开放体系，用户不仅可以采用它将ANSYS程序剪裁成符合自己所需的任何组织形式（如可以定义一种新的材料，一个新的单元或者给出一种新的屈服准则），而且还可以编写自己的优化算法，通过将整个ANSYS作为一个子程序调用的方式实现。 鉴于上述特点，近几年来，ANSYS软件在国内外工程建设和科学研究中得到了广泛的应用。但这些应用大多局限于直接运用ANSYS软件进行实际工程分析，对利用ANSYS提供的二次开发工具进行有限元软件设计却很少涉及。本文首次利用ANSYS软件的二次开发功能，以VC++6.0为工具，运用APDL语言，对ANSYS进行二次开发，编制框筒结构－桩筏基础－土相互作用体系与地震反应分析程序。2程序设计目标 针对某一实际工程问题，ANSYS所提供的APDL语言可对ANSYS软件进行封装。APDL语言即ANSYS软件提供的参数化设计语言，它的全称是ANSYS Parametric
输气管道受力分析（ANSYS建模） 任务和要求： 按照输气管道的尺寸及载荷情况,要求在ANSYS中建模，完成整个静力学分析过程。求出管壁的静力场分布。要求完成问题分析、求解步骤、程序代码、结果描述和总结五部分。所给的参数如下： 材料参数：弹性模量E=200Gpa; 泊松比0.26；外径R?=0.6m；内径R?=0.4m；壁厚t=0.2m。输气管体内表面的最大冲击载荷P为1Mpa。 四.问题求解 （一）．问题分析 由于管道沿长度方向的尺寸远大于管道的直径，在计算过程中忽略管道的端面效应，认为在其长度方向无应变产生，即可将该问题简化为平面应变问题，选取管道横截面建立几何模型进行求解。 （二）．求解步骤 定义工作文件名 选择Utility Menu→File→Chang Jobname 出现Change Jobname对话框，在[/FILNAM] Enter new jobname 输入栏中输入工作名LEILIN10074723，并将New log and eror file 设置为YES，单击[OK]按钮关闭对话框 定义单元类型 1)选择Main Meun→Preprocessor→Element Type→Add/Edit/Delte命令，出现Element Type 对话框，单击[Add]按钮，出现Library of Element types对话框。 2）在Library of Element types复选框选择Strctural、Solid、 Quad 8node 82,在Element type reference number输入栏中出入1，单击[OK]按钮关闭该对话框。 3. 定义材料性能参数 1）单击Main Meun→Preprocessor→Material Props→Material models出现Define Material Behavion 对话框。选择依次选择Structural、Linear、Elastic、Isotropic选项，出现Linear Isotropic Material Properties For Material Number 1对话框。 2）在EX输入2e11，在Prxy输入栏中输入0.26，单击OK按钮关闭该对话框。 3）在Define Material Model Behavion 对话框中选择Material→Exit命令关闭该对话框。 4.生成几何模型、划分网格 1）选择Main Meun→Preprocessor→Modeling→Create→Areas→Circle→Partail→Annulus出现Part Annulus Circ Area对话框，在WP X文本框中输入0，在WP Y文本框中输入0，在Rad1文本框中输入0.4，在Theate-1文本框中输入0，在Rad2文本框中输入0.6，在Theate-2文本框中输入90，单击OK按钮关闭该对话框。 2）选择Utility Menu→Plotctrls→Style→Colors→Reverse Video，设置显示颜色。 3）选择Utility Menu→Plot→Areas，显示所有面。 4) 选择Main Menu→Preprocessor→Modeling→Reflect→Areas，出现Reflect Areas拾取菜
ANSYS 中使用函数加载的一个简单例子 本文将通过一个具体实例说明在ANSYS 中如何使用函数加载，后续将通过该实例在分析过程中遇到的一个问题提出自己的一点看法。 实例的具体说明： 一个1/4 圆柱，内半径30 mm，外半径42 mm，长度100mm，如图1 所示： 所用材料为双线性弹塑性材料，其机械性能为： 弹性模量 E = 201000 Mpa；泊松比μ=0.3
屈服应力σ= 200 Mpa；切线模量Et = 2010使用单元类型solid185 (8 节点六面体单元)。 取整体单元边长4 mm，然后可以直接对该几何模型划分MAP 网格，划分网格结果如图2： 约束条件为： 轴向两个截面为对称边界条件；一个端面约束轴向位移Uz。 载荷条件为： 在外表面施加变化的压力载荷，载荷函数为： P (y) = 8e7 + 7E7 * (Y/42)
即： X = 0 ，Y = 42 (最高点) 时，P = 15E7； X = 42，Y = 0 （最低点）时，P = 8E7。 我们采用函数方式来施加这一压力载荷，首先定义函数： 在Solution 模块中，点击菜单路径： Solution > Define Loads > Apply > Functions > Define/Edit 将会弹出一个函数编辑器，可以在其中定义所需的函数。 在函数编辑器中，函数类型选择为Single equation，即单值函数；计算函数值时使用的插值坐标系( (x,y,z) interpreted in CSYS) 选择0，即总体直角坐标系，如图3 所示：
然后，在函数编辑器中间位置的“Result = “ 小窗口中输入要定义的函数表达式，如果表达式中有x, y, z, time 等变量(供定义函数时使用的“自变量”)，可以用{X},{Y},{Z},{TIME} 等的形式输入；或者点击下面一个小窗口右边的小箭头，会出现一个下拉列表，列出可以选择的变量，然后从该列表中选择某个自变量，则该自变量会按照上述格式写入函数中，如图5 所示：
平面问题斜支座的处理 如图5-7所示，为一个带斜支座的平面应力结构，其中位置2及3处为固定约束，位置4处为一个45o的斜支座，试用一个4节点矩形单元分析该结构的位移场。 (a)平面结构(b)有限元分析模型 图5-7 带斜支座的平面结构 基于ANSYS平台，分别采用约束方程以及局部坐标系的斜支座约束这两种方式来进行处理。 (7) 模型加约束 左边施加X,Y方向的位移约束 ANSYS Main Menu: Solution →Define Loads →Apply →-Structural→Displacement On Nodes →选取2，3号节点→OK →Lab2: All DOF(施加X,Y方向的位移约束) →OK 以下提供两种方法处理斜支座问题，使用时选择一种方法。 ?采用约束方程来处理斜支座 ANSYS Main Menu：Preprocessor →Coupling/ Ceqn →Constraint Eqn ：Const :0, NODE1:4, Lab1: UX,C1:1,NODE2:4,Lab2:UY,C2:1→OK 或者?采用斜支座的局部坐标来施加位移约束 ANSYS Utility Menu：WorkPlane →Local Coordinate System →Create local system →At specified LOC + →单击图形中的任意一点→OK →XC、YC、ZC分别设定为2，0，0,THXY：45 →OK ANSYS Main Menu：Preprocessor →modeling →Move / Modify →Rotate Node CS →To active CS → 选择4号节点 ANSYS Main Menu：Solution →Define Loads →Apply →Structural →Displacement On Nodes →选取4号节点→OK →选择Lab2:UY(施加Y方向的位移约束) →OK 命令流； !---方法1 begin----以下的一条命令为采用约束方程的方式对斜支座进行处理 CE,1,0,4,UX,1,4,UY,-1 !建立约束方程(No.1): 0=node4_UX*1+node_UY*(-1) !---方法1 end --- !--- 方法2 begin --以下三条命令为定义局部坐标系，进行旋转，施加位移约束 !local,11,0,2,0,0,45 !在4号节点建立局部坐标系 !nrotat, 4 !将4号节点坐标系旋转为与局部坐标系相同 !D,4,UY !在局部坐标下添加位移约束 !--- 方法2 end
将ANSYS作为子程序调用 对于优化或参数化设计，可以在VC或FORTRAN中将ANSYS作为子程序调用。具体调用方法如下： 1.在VC中调用ANSYS ::WinExec("d:/ANSYS57/BIN/INTEL/ANSYS57 -b -p ansys_product_feature -i input_file -o output_file",SW_SHOWNORMAL); 2.在FORTRAN中调用ANSYS LOGICAL(4) result RESULT=SYSTEMQQ('d:\ANSYS57\BIN\INTEL\ANSYS57 -b -p ansys_product_feature -i input_file -o output_file') 3.说明 1和2中，input_file为用APDL语言编写的ANSYS输入文件。 ansys_product_feature为你的ANSYS产品特征代码。 需要注意的是，在VC中调用ANSYS时，需要加一条判断语句，以确定ANSYS 已经执行完毕。 在ANSYS中当然也可以以VC或FORTRAN作为子程序调用。可以参看有关ANSYS二次开发方面的资料。 这个方法应该是与系统无关的。 在FORTRAN中不需要判断，FORTRAN会等ANSYS执行完毕才继续执行下一条语句。 在VC中，我没有找到与FORTRAN类似的函数，只好加一条循环判断语句。 如果谁能找着这样的函数，请告诉我，谢谢！ 判断方法很简单，只需判断错误文件file.err是否可写就可以了。 因为当ANSYS在运行时，file.err是不可写的，只有当它运行完毕，此文件才可写。 好贴，请教如何在VB中调用ansys？ 这是我早期的帖子，请参考： https://www.doczj.com/doc/d02794745.html,/dispbb ... ID=2923&ID=2923 在第二页中已经对VC调用ANSYS的方法进行了更新，如下所示： VC调用ANSYS的示例程序。 //Test.cpp #include "stdio.h" #include"process.h" void main() { int result; printf("Solving..."); result=system("d:/ANSYS57/BIN/INTEL/ANSYS57 -b -p ansysul -i test.txt -o test.out"); //不用::WinExec，就用不着等待语句，可以实现用FORTRAN调用一样的效果。 printf("Solution finished..."); }
/PREP7 /TITLE,Steady-state thermal analysis of pipe junction /UNITS,BIN ! 英制单位；Use U. S. Customary system of units (inches) ! /SHOW, ! Specify graphics driver for interactive run ET,1,90 ! Define 20-node, 3-D thermal solid element MP,DENS,1,.285 ! Density = .285 lbf/in^3 MPTEMP,,70,200,300,400,500 ! Create temperature table MPDATA,KXX,1,,8.35/12,8.90/12,9.35/12,9.80/12,10.23/12 ! 指定与温度相对应的数据材料属性；导热系数；Define conductivity values MPDATA,C,1,,.113,.117,.119,.122,.125 ! Define specific heat values（比热） MPDATA,HF,2,,426/144,405/144,352/144,275/144,221/144 ! Define film coefficient;除144是单位问题，上面的除12也是单元问题 ! Define parameters for model generation RI1=1.3 ! Inside radius of cylindrical tank RO1=1.5 ! Outside radius Z1=2 ! Length RI2=.4 ! Inside radius of pipe RO2=.5 ! Outside pipe radius Z2=2 ! Pipe length CYLIND,RI1,RO1,,Z1,,90 ! 90 degree cylindrical volume for tank WPROTA,0,-90 ! 旋转当前工作的平面；从Y到Z旋转-90度；；Rotate working plane to pipe axis CYLIND,RI2,RO2,,Z2,-90 ! 角度选择在了第四象限；90 degree cylindrical volume for pipe WPSTYL,DEFA ! 重新安排工作平面的设置；另外WPSTYL,STAT to list the status of the working plane；；Return working plane to default setting BOPT,NUMB,OFF ! 关掉布尔操作的数字警告信息；Turn off Boolean numbering warning VOVLAP,1,2 ! 交迭体；Overlap the two cylinders /PNUM,VOLU,1 ! 体编号打开；Turn volume numbers on /VIEW,,-3,-1,1