一简述ANSYS软件的发展史。 1970年，Doctor John Swanson博士洞察到计算机模拟工程应该商品化，于是创立了ANSYS公司，总部位于美国宾夕法尼亚州的匹兹堡。30年来，ANSYS 公司致力于设计分析软件的开发，不断吸取新的计算方法和技术，领导着世界有限元技术的发展，并为全球工业广泛接受，其50000多用户遍及世界。 ANSYS软件的第一个版本仅提供了热分析及线性结构分析功能，像当时的大多数程序一样，它只是一个批处理程序，且只能在大型计算机上运行。 20世纪70年代初。ANSYS软件中融入了新的技术以及用户的要求，从而使程序发生了很大的变化，非线性、子结构以及更多的单元类型被加入到子程序。70年代末交互方式的加入是该软件最为显著的变化，它大大的简化了模型生成和结果评价。在进行分析之前，可用交互式图形来验证模型的几何形状、材料及边界条件；在分析完成以后，计算结果的图形显示，立即可用于分析检验。 今天软件的功能更加强大，使用更加便利。ANSYS提供的虚拟样机设计法，使用户减少了昂贵费时的物理样机，在一个连续的、相互协作的工程设计中，分析用于整个产品的开发过程。ANSYS分析模拟工具易于使用、支持多种工作平台、并在异种异构平台上数据百分百兼容、提供了多种耦合的分析功能。 ANSYS公司对软件的质量非常重视，新版的必须通过7000道标准考题。业界典范的质保体系，自动化规范化的质量测试使ANSYS公司于1995年5月在设计分析软件中第一个通过了ISO9001的质量体系认证。 ANSYS公司于1996年2月在北京开设了第一个驻华办事机构，短短几年的时间里发展到北京、上海、成都等多个办事处。ANSYS软件与中国压力容器标准化技术委员会合作，在1996年开发了符合中国JB4732-95国家标准的中国压力容器版。作为ANSYS集团用户的铁路机车车辆总公司，在其机车提速的研制中，ANSYS软件已经开始发挥作用。 二节点﹑单元﹑单元类型的基本概念。 节点：几何模型通过划分网格，转化为有限元模型，节点构成了网格的分布和形状，是构成有限元模型的基本元素。 单元：有限元模型的组成元素，主要有点、线、面、体。 单元类型：根据实体模型划分网格时所要确定的单元的形状，是单元属性的一部分，单元类型决定了单元的自由度，包括线单元（梁、杆、弹簧单元）、壳单元（用于薄板或曲面模型）、二维实体单元、三维实体单元、线性单元、二次单元和P–单元。 三用ANSYS软件进行分析的一般过程。 1建立有限元模型 （1）指定工作文件名和工作标题。 该项工作并不是必须要求做的，但是做对多个工程问题进行分析时推荐使用工作文件名和工作标题。
“有限元分析及应用”课程有限元分析软件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
有限元 学院：机电学院 专业： 姓名： 学号：
一、问题描述 如图所示的平面，板厚为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 在弹出的对
1 E c ； / E c lE s _2卜+僅 12 (5-30) 通过大型有限元软件ANSYS 对简支梁进行模拟计算 下面以钢筋混凝土简支梁的 ANSYS ①程序数值模拟的应用实例，对ANSYS ⑧程序的应用方法及 模拟效果进行验证，梁的尺寸、配筋及荷载如图5-9所示。钢筋采用H 级钢，混凝土强度等级为 C30。 2.1单元类型 i )混凝土单元：采用 ANSYS ①程序单元库中 SOLID65单元。 (ii )纵向钢筋：PIPE20 (iii )横向箍筋:PIPE20 2.2材料性质 i )、混凝土材料 [16~ 19] 混凝土立方体抗压强 度 f cu ( N / mm 2) 弹性模量E c 2 (N/mm ) 泊松 比 V 单轴抗压强度f c ' 2 (N/mm ) 单轴抗拉强度f r (N/mm ) 裂缝间剪力 传递 系数P t 张开 闭合 30 24000 0.20 25.0 3.1125 0.35 0.75 ?单轴受压应力-应变曲线(二-；曲线) 在ANSYS @程序分析中，需要给出混凝土单轴受压下的应力应变曲线。在本算例中，混凝土单 轴受压下的应力应变采用 Sargin 和Saenz 模型[17,18]: ①22①22 E20 ①22 RCBEAM-03 图5-9 2①82①82①8 2①82①82①8 ① 8@75@75@75 2①22①22①22 150 150 150 150 RCBEAM-01 150150150 RCBEAM-02 (b )、梁断面图 梁尺寸、配筋及荷载示意图 f ①24 ①24 ①22 150 150 ■4- ------------- P P 125 1200 600 (a )、梁的几何尺寸及荷载示意图 600
/一个周边简支的圆盘，其中心受到一个冲杆的周期作用(假定冲杆是刚性的)，需要进行圆盘在冲杆的周期作用下的塑性分析。本实例的模型简图如图19.1所示，材料特性如下所示，塑性时的应力－应变关系如表19.1，载荷历史如表19.2所示。 弹性模量：EX=70000，泊松比：NUXY=0.325 /PREP7 /TITLE,Circular Plate Loaded by a Circular Punch - Kinematic Hardening !* 下面定义建模分析时需要的参数 EXX=70000 RPL=65 RPU=5 H=6.5 STS1=55 STN1=STS1/EXX STS2=112 STN2=0.00575 STS3=172 STN3= 0.02925 STS4=241 STN4= 0.1 NEX=15 NET=2 NEX1=nint(0.8*NET) NEX2=NEX-NEX1 !* ET,1,42,,,1 ！定义单元PLANE42，设置为轴对称 !* MP,EX,1,EXX ！定义材料属性 MP,NUXY,1,0.325 !* TB,KINH,1,1,4, ！定义多线性随动强化准则 TBPT,,STN1,STS1 TBPT,,STN2,STS2 TBPT,,STN3,STS3 TBPT,,STN4,STS4 !* 创建节点 N,1,RPL,,,,,, N,2,0,,,,,, N,3,,H/2,,,,, !* 创建关键点 K,1,,-(H/2),, K,2,RPU,-(H/2),, K,3,RPL,-(H/2),, KGEN,2,ALL, , , ,H, ,3,0 ！复制并平移关键点
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程序）。
第四章子结构 什么是子结构？ 子结构就是将一组单元用矩阵凝聚为一个单元的过程。这个单一的矩阵单元称为超单元。在ANSYS分析中，超单元可以象其他单元类型一样使用。唯一的区别就是必须先进行结构生成分析以生成超单元。子结构可以在ANSYS/Mutiphysics，ANSYS/Mechanical和ANSYS/Structural中使用。 使用子结构主要是为了节省机时，并且允许在比较有限的计算机设备资源的基础上求解超大规模的问题。原因之一如a)非线性分析和带有大量重复几何结构的分析。在非线性分析中，可以将模型线性部分作成子结构，这样这部分的单元矩阵就不用在非线性迭代过程中重复计算。在有重复几何结构的模型中（如有四条腿的桌子），可以对于重复的部分生成超单元，然后将它拷贝到不同的位置，这样做可以节省大量的机时。 子结构还用于模型有大转动的情况下。对于这些模型，ANSYS假定每个结构都是围绕其质心转动的。在三维情况下，子结构有三个转动自由度和三个平动自由度。在大转动模型中，用户在使用部分之前无须对子结构施加约束，因为每个子结构都是作为一个单元进行处理，是允许刚体位移的。 另外一个原因b)一个问题就波前大小和需用磁盘空间来说相对于一个计算 1
机系统太庞大了。这样，用户可以通过子结构将问题分块进行分析，每一块对于计算机系统来说都是可以计算的。 如何使用子结构 子结构分析有以下三个步骤： ●生成部分 ●使用部分 ●扩展部分 生成部分就是将普通的有限元单元凝聚为一个超单元。凝聚是通过定义一组主自由度来实现的。主自由度用于定义超单元与模型中其他单元的边界，提取模型的动力学特性。图4-1是一个板状构件用接触单元分析的示意。由于接触单元需要迭代计算，将板状构件形成子结构将显著地节省机时。本例中，主自由度是板与接触单元相连的自由度。 图4-1 子结构使用示例 2
输气管道受力分析（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拾取菜
（ii ）纵向钢筋：PIPE20 （iii ）横向箍筋：PIPE20 2.2 材料性质 （i ）、混凝土材料 表5-4 混凝土材料的输入参数一览表[16~19] ·单轴受压应力-应变曲线（εσ-曲线） 在ANSYS ○R 程序分析中，需要给出混凝土单轴受压下的应力应变曲线。在本算例中，混凝土单轴受压下的应力应变采用Sargin 和Saenz 模型[17,18]： 2 21??? ? ??+???? ??-+= c c s c c E E E εεεεε σ （5-30）
式中取4' 4')108.0028.1(c c c f f -=ε；
断面图配筋图断面图配筋图断面图配筋图RCBEAM-01 RCBEAM-02 RCBEAM-03 图5-12 各梁FEM模型断面图 （a）单元网格图（b）钢筋单元划分图 图5-13 算例（一）的FEM模型图 2.4 模型求解 在ANSYS○R程序中，对于非线性分析，求解步的设置很关键，对计算是否收敛关系很大，对于混凝土非线性有限元分析，在计算时间容许的情况下，较多的求解子步（Substeps）或较小的荷载步和一个非常大的最大子步数更容易导致收敛[2]。在本算例中，设置了100个子步。最终本算例收敛成功，在CPU为P41.6G、内存为256MB的微机上计算，耗时约为8小时。 2.5 计算结果及分析 2.5.1 荷载—位移曲线 图5-14为ANSYS○R程序所得到的各梁的荷载-跨中挠度曲线，从图中可以看出： （i）、梁RCBEAM-01：曲线形状能基本反映钢筋混凝土适筋梁剪切破坏的受力特点，而且荷载-跨中挠度曲线与钢筋混凝土梁的弯剪破坏形态非常类似，即当跨中弯矩最大截面的纵筋屈服后，由于裂缝的开展，压区混凝土的面积逐渐减小，在荷载几乎不增加的情况下，压区混凝土所受的正应力和剪应力还在不断增加，当应力达到混凝土强度极限时，剪切破坏发生，荷载突然降低。
平面问题斜支座的处理 如图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 中使用函数加载的一个简单例子 本文将通过一个具体实例说明在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 所示：
龙源期刊网 https://www.doczj.com/doc/629352708.html, 格构梁的ANSYS有限元模拟分析实例运用作者：张少剑刘真 来源：《城市建设理论研究》2013年第10期 摘要:本文通过一工程实例运用ansys模拟计算。针对格构梁的研究，合理地简化模型，取出1.5米宽的土体、梁和面层单元，两边加对称约束，从而达到模拟空间结构梁的目的。本文还模拟了基坑的开挖过程的时空效应，共分七步，土体在自重应力作用下的沉降为第一步，梁与面层的激活、力的施加和土层杀死共分六步。梁的最大受力状态并不发生在最后一步完成后，而是在第六工况。 关键词:格构梁有限元分析模拟分析 中图分类号：K826.16 文献标识码：A 文章编号： 1 土体、梁、锚索和混凝土面层共同作用 基坑支护的受力机理是土体的土压力作用在格构梁和混凝土面层上，混凝土面层的力传递到格构梁上，格构梁再把它受到的力传递到和它相连的锚索上，锚索则和被支护土体嵌固为一体，格构梁和混凝土面层除起到承受土压力外，格构梁还起到平均弯矩和变形的作用，喷射混凝土面层则有保护土体表面，防止土体表面非格构梁作用部位坍塌的作用。 2模型简化及技术处理 根据基坑开挖深度，根据实际的土体性质建立土体模型。格构梁的作用是承受弯矩的，可以选用Beam4梁单元，考虑到钢筋混凝土格构梁中有钢筋的作用，其弹性模量、泊松比等设置有所调整。在建模时，如果混凝土面层的长宽与厚度的比都大于5，所以在有限元分析中采用板壳单元可以全面地反映其变形特征和应力分布规律。混凝土面层用Shell63单元模拟，其参数的取值和梁单元相同。 由于格构梁的受力性状，锚索的模拟对格构梁的受力影响较小，本模型忽略考虑锚索的模拟。预应力锚索的作用简化为作用在纵横梁交点处的集中力。 对于格构梁和土体、混凝面层之间的接触，模型采用节点耦合，以实现共同变形和受力。 3.1ANSYS有限元模拟计算 3.1.1模型的参数 1．土体的参数见下表：
ANSYS 结构稳态(静力)分析之经典实例-命令流格式.txt两人之间的感情就像织毛衣，建立 的时候一针一线，小心而漫长，拆除的时候只要轻轻一拉。。。。/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生 一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT
/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
过盈配合在机械产品的装配中使用的相当普遍。比如轴与轴承、轴与轴瓦、汽车的制动盘等，都是通过一定的过盈量来使两个装配部件紧密地连接起来。 下面讨论如何在ANSYS 中正确地模拟过盈配合。 过盈配合在有限元分析中是一种典型的非线性接触行为。在有限元分析中设定了接触，从本质上来讲就是对相互接触的两个部件施加了某种约束，不同的接触算法对于接触约束的处理方法有所不同。接触约束的理论算法的选择，在ANSYS 中是通过设置contact 单元的KEOPT(2) 选项来实现的。 在ANSYS 中目前主要有5 种接触约束算法： KEYOPT(2)=0 Augmented Lagrangian - 加强的拉格朗日算法，这是ANSYS 的缺省选择； KEYOPT(2)=1 Penalty function - 罚函算法； KEYOPT(2)=2 Multipoint constraint (MPC) - 多点约束算法； KEYOPT(2)=3 Lagrange multiplier on contact normal and penalty on tangent - 接触法向采用拉格朗日乘子，接触切向采用罚函数的综合算法。 KEYOPT(2)=4 Pure Lagrange multiplier on contact normal and tangent - 法向和切向均采用拉格朗日乘子算法。 各种不同的约束算法各有其优缺点，各有各自最适用的场合，具体情况需要具体对待。大部分情况下，默认选择KEYOPT(2)=0 就够用了。 过盈配合所致的接触分析的难点在于如何确定初始接触状态。初始
接触状态设置得不对，会导致错误的计算结果或者不准确的计算结果，下面举两个例子来说明。ANSYS仿真计算代做：模态分析，瞬态动力学，谐响应分析和谱分析、械结构的疲劳、损伤，CFD流体；结构的强度评估和优化；企鹅：690294845 例1．两个圆柱体在几何上是刚好接触，划分网格后有限元模型有间隙。如图1 所示。 这两个圆柱体，在几何上是刚好相切的，即处于几何上刚好接触的初始状态。划分网格后，由于在圆周上用小段直线代替了弧线，两个圆柱体之间产生了一定的间隙，两个圆柱体的有限元模型的初始状态不再是接触的。此时，如果接触参数设置不当，就会因为初始约束不足，圆柱体出现刚体位移，得到错误的结果。
ANSYS中的超单元 从8.0版开始，ANSYS中增加了超单元功能，本文通过一些实际例子，探讨了ANSYS 中超单元的具体使用。 1 使用超单元进行静力分析 根据ANSYS帮助文件，使用超单元的过程可以划分为三个阶段(称为Pass)： (1) 生成超单元模型(Generation Pass) (2) 使用超单元数据(Use Pass) (3) 扩展模型(Expansion Pass) 以下摘自htbbzzg邹老师博客，请勿乱传! 下面以一个例子加以说明： 一块板，尺寸为20×40×2，材料为钢，一端固支，另一端承受法向载荷。 首先生成原始模型se_all.db，即按照整个结构进行分析，以便后面与超单元结果进行比较： 首先生成两个矩形，尺寸各为20×2。然后定义单元类型shell63； 定义实常数1为： 2 (板厚度)。 材料性能：弹性模量E=201000；波松比μ=0.3；密度ρ=7.8e-9； 单位为mm-s-N-MPa。 采用边长1划分单元；一端设置位移约束all，另一端所有(21个)节点各承受Z向力5。计算模型如下图：
静力分析的计算结果如下：
为了后面比较的方便，分别给出两个area上的结果：
超单元部分，按照上述步骤操作如下： (1) 生成超单元 选择后半段作为超单元，前半段作为非超单元(主单元)。 按照ANSYS使用超单元的要求，超单元与非超单元部分的界面节点必须一致(重合)，且最好分别的节点编号也相同，否则需要分别对各节点对建立耦合方程，操作比较麻烦。 实际上，利用ANSYS中提供的mesh200单元，对超单元和非超单元的界面实体，按照同一顺序，先于所有其它实体划分单元，很容易满足界面节点编号相同的要求。对于多级超单元的情况，则还要结合其它操作(如偏移节点号等)以满足这一要求。 对于本例，采用另一办法，即先建立整个模型，然后再划分超单元和非超单元。即：将上述模型分别存为se_1.db (超单元部分)和se_main.db (非超单元部分)两个文件，然后分别处理。 对于se_1.db模型，按照超单元方式进行处理。由于模型及边界条件已建立，只需删除前半段上的划分，结果就是超单元所需的模型。 然后直接进入创建超单元矩阵的操作，首先说明一下创建超单元矩阵的一般步骤： A进入求解模块： 命令：/Solu GUI：Main menu -> Solution B设置分析类型为“子结构或部件模态综合“
焊接模拟ansys实例 !下面的命令流进行的是一个简单的二维焊接分析, 利用ANSYS单元生死和热-结构耦合分析功能进 !行焊接过程仿真, 计算焊接过程中的温度分布和应力分布以及冷却后的焊缝残余应力。 finish /clear /filnam,1-2D element birth and death /title,Weld Analysis by "Element Birth and Death" /prep7 /unit,si !采用国际单位制 !****************************************************** et,1,13,4 !13号二维耦合单元, 同时具有温度和位移自由度 et,2,13,4 !1号材料是钢 !2号材料是铝 !3号材料是铜 !铝是本次分析中的焊料, 它将钢结构部分和铜结构部分焊接起来 !下面是在几个温度点下, 各材料的弹性模量 mptemp,1,20,500,1000,1500,2000 mpdata,ex,1,1,1.93e11,1.50e11,0.70e11,0.10e11,0.01e11 mpdata,ex,2,1,1.02e11,0.50e11,0.08e11,0.001e11,0.0001e11 mpdata,ex,3,1,1.17e11,0.90e11,0.30e11,0.05e11,0.005e11 !假设各材料都是双线性随动硬化弹塑性本构关系 !下面是各材料在各温度点下的屈服应力和屈服后的弹性模量 tb,bkin,1,5 tbtemp,20,1 tbdata,1,1200e6,0.193e11 tbtemp,500,2 tbdata,1, 933e6,0.150e11 tbtemp,1000,3 tbdata,1, 435e6,0.070e11 tbtemp,1500,4 tbdata,1, 70e6,0.010e11 tbtemp,2000,5 tbdata,1, 7e6,0.001e11 ! tb,bkin,2,5 tbtemp,20,1 tbdata,1,800e6,0.102e11 tbtemp,500,2 tbdata,1,400e6,0.050e11 tbtemp,1000,3 tbdata,1, 70e6,0.008e11