ANSYS在玻璃钢-混凝土复合材料结构中的计算分析

殷波

（扬州大学水利与建筑工程学院土木工程系，扬州225009）

摘要：混凝土结构由于受荷载变化、材料、施工质量等因素影响，会造成结构的强度、刚度不足，玻璃钢-混凝土复合材料结构则改善其性能。本文通过ansys有限元软件计算，分析说明了玻璃钢-混凝土复合材料结构将有力的提高结构的强度和刚度。

关键词：ansys，混凝土、玻璃钢-混凝土、有限元、复合材料

ANSYS’s calculation in glass fibre reinforced plastics- concrete

composite material structure

YIN BO

(Dept.of Civil Engin,Hydr and Civil Engin Coll, Yangzhou University, Yangzhou, 225009, China)

Abstract: As the variety of load、material、constructional quality and so on ,concrete structure may be insufficient in intensity and rigidity. Glass fibre reinforced plastics- concrete composite material structure can improve its capacity. With the calculation of ansys, this paper indicate that glass fibre reinforced plastics-concrete composite material structure will raise the intensity and rigidity.

Key words: ansys ,concrete, glass fibre reinforced plastics- concrete, finite element; composite material

1．前言

钢筋混凝土结构合理地利用钢筋和混凝土两种材料的力学性能，因而具有整体性、耐久性等优点。但同时也存在着自重大等缺点，在受弯构件中较显著。自重太大，对于设计大跨度结构、高层建筑以及抗震结构都很不利。针对普通钢筋混凝土梁以及经由玻璃钢加固的钢筋混凝土复合梁，通过ansys有限元软件计算分析说明了玻璃钢-混凝土复合材料结构将提高其强度和刚度。

2．用ANSYS进行有限元计算

有限元法是目前工程技术领域中实用性最强、应用最为广泛的数值计算方法。它的基本思想是将问题的求解域划分为一系列单元，单元之间靠节点连接。单元内部点的待求物理量可由单元节点物理量通过选定的函数关系插值求得。由于单元形状简单，易于由平衡关系或能量关系建立节点量之间的方程式，然后由单元方程再形成总体代数方程组，加入边界条件后即可对方程组求解。以ANSYS为代表的有限元分析软件具有以下优点:减少设计成本、缩短设计和分析的循环周期、增加产品和工程的可靠性、采用优化设计、降低材料的消耗和成本、在产品制造或工程施工前预先发现潜在的问题、可以进行模拟实验分析、进行机械事故分析,查找事故原因。ANSYS软件是集结构、热、流体、电磁、声学于一体的大型通用有限元分析软件,可广泛地应用于土木工程、

交通、水利、铁道、石油化工、航空航天、机械制造、国防军工、电子、造船、生物医学、地矿、日用家电等一般工业及科学研究。

本文采用钢筋混凝土梁为单筋矩形梁，截面为 26kN 26kN

100ｍｍ×150ｍｍ，底部配有2φ12钢筋。混凝土抗压 0.353kN/m

强度ｆｃ=30ＭＰａ。梁受荷载：0.353ｋＮ/ｍ作用和

三分点处各受两个集中荷载26ｋＮ作用如图1所示。

图1 梁受荷载图

2.1 钢筋混凝土梁的计算

钢筋混凝土结构比较特殊，是由两种材料所组成，可分别建立其模型。也有较方便的方法，本文即是采用的solid65钢筋混凝土单元。先建立实体模型，设置solid65单元属性和材料特性。然后对其进行网格划分，再施加荷载和位移约束，即可对该钢筋混凝土梁进行有限元计算，如图2所示。

图2 钢筋混凝土梁计算Sx图

2.2 玻璃钢加固的钢筋混凝土梁的计算

在上述的钢筋混凝土梁上、下各加一层玻璃纤维-环氧复合材料，厚度均为2mm。建立实体模型，利用solid65对钢筋混凝土进行网格划分。利用叠层复合材料实体单元、solid46单元对玻璃纤维进行网格划分。设置复合材料特性，施加荷载和位移约束，得玻璃钢加固的钢筋混凝土梁计算结果如图3所示。

图3 玻璃钢加固的钢筋混凝土梁计算Sx图

2.3 两种梁计算结果比较

对上述两种梁，通过ansys有限元软件计算，得出梁的最大挠度、横截面上的正应力及von-mises应力数值如表1所示。

表1 两种梁计算结果比较

梁类型最大挠度（m）SMN (Pa) SMX (Pa) Von-mises(Pa)

钢筋混凝土梁.650E-03 -.140E+08 .135E+08 .161E+08

玻璃钢加固钢筋混凝土梁 .457E-04 -.239E+08 .501E+07 .456E+08 3．结论

1．钢筋混凝土构件及玻璃钢加固的钢筋混凝土构件均为多种不同性质材料所组成，理论计算较为复杂，而应用ansys有限元软件计算则很行之有效。

2．通过ansys计算表明：玻璃钢加固的钢筋混凝土梁挠度大为减小，这对减轻钢筋混凝土构件自重具有指导意义。

3．玻璃钢加固的钢筋混凝土梁的横截面上的正应力及Von-mises应力数值都较普通钢筋混凝土梁的数值低，这同时说明提高了梁的强度，增加了构件抵抗破坏的能力。

参考文献

[1]陈精一蔡国忠.电脑辅助工程分析 ANSYS使用指南[M].北京:中国铁道出版社，2001

[2] 扬庆生.复合材料细观结构力学与设计[M].北京: 中国铁道出版社，2000

[3] 李顺林等. 复合材料工作手册[M].北京:航空工业出版社，1988

ANSYS软件介绍与实例讲解

一简述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）指定工作文件名和工作标题。该项工作并不是必须要求做的，但是做对多个工程问题进行分析时推荐使用工作文件名和工作标题。

有限元ansys静力分析的一个小例子

有限元学院：机电学院专业：姓名：学号：

一、问题描述如图所示的平面，板厚为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中简支梁的模拟计算

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

ANSYS新手入门学习心得

(1) 如果你模拟结构体中裂缝扩展过程的模拟，在Ansys中可以用全解耦损伤分析方法来近似模拟裂缝扩展，我曾用Ansys软件中提供的可以定义10,000个材料参数和单元ekill/alive 功能完成了层状路面体中表面裂缝和反射裂缝在变温作用下的扩展过程的模拟。我模拟的过程相对来说比较简单，模拟过程中我们首先要知道裂缝的可能扩展方向，这样在裂缝可能扩展的带内进行网格加密处理，加密到什么程度依据计算的问题来确定。 (2) 如果采用断裂力学理论计算含裂缝结构体的应力强度因子，建模时只需在裂尖通过命令kscon生成奇异单元即可。Ansys模块中存在的断裂力学模块可以计算I、II、III型应力强度因子(线弹性断裂力学)和J积分(弹塑性断裂力学)，在Ansys中verification里面有一个计算I型应力强度因子的例子vm143，参见该例子就可以了。 (3) 如果通过断裂力学模拟裂缝的扩展过程，需要采用动态网格划分，这方面我没有做，通过Ansys的宏命令流应该可以实现。技术参考可参阅文献：杨庆生、杨卫.断裂过程的有限元模拟.计算力学学报，1997,14(4). (4) 我现在做动荷载作用下路面结构体中应力强度因子的分布规律，我是通过位移插值得到不同时间点处的应力强度因子。如果想这样做，可参阅理论参考中关于应力强度因子计算说明。 1. 讨论两种Ansys求极限荷载的方法（1）力加载可以通过对应的方法（比如说特征值屈曲）估计结构的极限荷载的大致范围，然后给结构施加一个稍大的荷载，打开自动荷载步二分法进行非线性静力分析，最后计算会因不收敛终止，则倒数第二个子步对应的就是结构的极限荷载；另外，也可以选择弧长法，采用足够的子步（弧长法可以一直分析到极限承载力之后的过程）同样可以从绘制的荷载位移曲线或计算结果中找出结构的极限荷载。（2）位移加载给结构施加一个比较大的位移，打开自动荷载步二分法进行非线性分析，保证足够的子步数，这样也可以分析到极限荷载以后，通过绘制荷载位移曲线或查看相应结果文件也可知道结构的极限荷载。希望众高手讨论一下（1）弧长法求极限荷载的收敛性问题，如何画到荷载位移曲线的下降段？（2）位移法求极限荷载的具体步骤？ 2. 需要注意的问题 1. 由于SOLID 65单元本身是基于弥散裂缝模型和最大拉应力开裂判据，因此在很多情况下会因为应力集中而使混凝土提前破坏，从而和试验结果不相吻合，因此，在实际应用过程中应该对单元分划进行有效控制，根据作者经验，当最小单元尺寸大于5cm 时，就可以有效避免应力集中带来的问题； 2. 支座是另一个需要注意的问题。在有限元分析中，很多时候约束是直接加在混凝土节点上，这样很可能在支座位置产生很大的应力集中，从而使支座附近的混凝土突然破坏，造成求解失败。因此，在实际应用过程中，应该适当加大支座附近单元的尺寸或者在支座上加一些弹性垫块，避免支座的应力集中；

ANSYS塑性变形模拟例子

/一个周边简支的圆盘，其中心受到一个冲杆的周期作用(假定冲杆是刚性的)，需要进行圆盘在冲杆的周期作用下的塑性分析。本实例的模型简图如图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学习心得

一学习ANSYS需要认识到的几点相对于其他应用型软件而言，ANSYS作为大型权威性的有限元分析软件，对提高解决问题的能力是一个全面的锻炼过程，是一门相当难学的软件，因而，要学好ANSYS，对学习者就提出了很高的要求，一方面，需要学习者有比较扎实的力学理论基础，对ANSYS分析结果能有个比较准确的预测和判断，可以说，理论水平的高低在很大程度上决定了ANSYS使用水平；另一方面，需要学习者不断摸索出软件的使用经验不断总结以提高解决问题的效率。在学习ANSYS的方法上，为了让初学者有一个比较好的把握，特提出以下五点建议：（1）将ANSYS的学习紧密与工程力学专业结合起来毫无疑问，刚开始接触ANSYS时，如果对有限元，单元，节点，形函数等《有限元单元法及程序设计》中的基本概念没有清楚的了解话，那么学ANSYS很长一段时间都会感觉还没入门，只是在僵硬的模仿，即使已经了解了，在学ANSYS之前，也非常有必要先反复看几遍书，加深对有限元单元法及其基本概念的理解。作为工程力学专业的学生，虽然力学理论知识学了很多，但对许多基本概念的理解许多人基本上是只停留于一个符号的认识上，理论认识不够，更没有太多的感性认识，比如一开始学ANSYS时可能很多人都不知道钢材应输入一个多大的弹性模量是合适的。而在进行有限元数值计算时，需要对相关参数的数值有很清楚的了解，比如材料常数，直接关系到结果的正确性，一定要准确。实际上在学ANSYS时，以前学的很多基本概念和力学理论知识都忘得差不多了，因而遇到有一

定理论难度的问题可能很难下手，特别是对结果的分析，需要用到《材料力学》，《弹性力学》和《塑性力学》里面的知识进行理论上的判断，所以在这种情况下，复习一下《材料力学》，《弹性力学》和《塑性力学》是非常有必要的，加深对基本概念的理解，实际上，适当的复习并不要花很多时间，效果却很明显，不仅能勾起遥远的回忆，加深理解，又能使遇到的问题得到顺利的解决。在涉及到复杂的非线性问题时（比如接触问题），一方面，不同的问题对应着不同的数值计算方法，求解器的选择直接关系到程序的计算代价和问题是否能顺利解决；另一方面，需要对非线性的求解过程有比较清楚的了解，知道程序的求解是如何实现的。只有这样，才能在程序的求解过程中，对计算的情况做出正确的判断。因此，要能对具体的问题选择什么计算方法做出正确判断以及对计算过程进行适当控制，对《计算方法》里面的知识必须要相当熟悉，将其理解运用到ANSYS的计算过程中来，彼此相互加强理解。要知道ANSYS是基于有限元单元法与现代数值计算方法的发展而逐步发展起来的。因此，在解决非线性问题时，千万别忘了复习一下《计算方法》。此外，对《计算固体力学》也要有所了解（一门非常难学的课），ANSYS对非线性问题处理的理论基础就是基于《计算固体力学》里面所讲到的复杂理论。作为学工程力学的学生，提高建模能力是非常急需加强的一个方面。在做偏向于理论的分析时，可能对建模能力要求不是很高，但对于实际的工程问题，有限元模型的建立可以说是一个最重要的问题，而后

七个ansys经典入门实例

“有限元分析及应用”课程有限元分析软件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二次开发及实例

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心得

1. 讨论两种Ansys求极限荷载的方法（1）力加载可以通过对应的方法（比如说特征值屈曲）估计结构的极限荷载的大致范围，然后给结构施加一个稍大的荷载，打开自动荷载步二分法进行非线性静力分析，最后计算会因不收敛终止，则倒数第二个子步对应的就是结构的极限荷载；另外，也可以选择弧长法，采用足够的子步（弧长法可以一直分析到极限承载力之后的过程）同样可以从绘制的荷载位移曲线或计算结果中找出结构的极限荷载。（2）位移加载给结构施加一个比较大的位移，打开自动荷载步二分法进行非线性分析，保证足够的子步数，这样也可以分析到极限荷载以后，通过绘制荷载位移曲线或查看相应结果文件也可知道结构的极限荷载。希望众高手讨论一下（1）弧长法求极限荷载的收敛性问题，如何画到荷载位移曲线的下降段？（2）位移法求极限荷载的具体步骤？ 2. 需要注意的问题 1. 由于SOLID 65单元本身是基于弥散裂缝模型和最大拉应力开裂判据，因此在很多情况下会因为应力集中而使混凝土提前破坏，从而和试验结果不相吻合，因此，在实际应用过程中应该对单元分划进行有效控制，根据作者经验，当最小单元尺寸大于5cm 时，就可以有效避免应力集中带来的问题； 2. 支座是另一个需要注意的问题。在有限元分析中，很多时候约束是直接加在混凝土节点上，这样很可能在支座位置产生很大的应力集中，从而使支座附近的混凝土突然破坏，造成求解失败。因此，在实际应用过程中，应该适当加大支座附近单元的尺寸或者在支座上加一些弹性垫块，避免支座的应力集中； 3. 六面体的SOLID 65 单元一般比四面体的单元计算要稳定且收敛性好，因此，只要条件允许，应该尽量使用六面体单元； 4. 正确选择收敛标准，一般位移控制加载最好用位移的无穷范数控制收敛，而用力控制加载时可以用残余力的二范数控制收敛。在裂缝刚刚出现和接近破坏的阶段，可以适当放松收敛标准，保证计算的连续性； 3. 关于下降段的问题 1）在实际混凝土中都有下降段，但是在计算的时候要特别小心下降段的问题。 2）下降段很容易导致计算不收敛，有时为了计算的收敛要避免设置下降段，采用rush模型。 3）利用最大压应变准则来判断混凝土是否破坏。 4. Solid65单元中的破坏准则 1）采用Willam&Warnke五参数破坏准则 2）需要参数：单轴抗拉强度，单轴，双轴抗压强度，围压压力，在围压作用下双轴，单轴抗压强度 5. 近来我对混凝土单元进行了一点思考，有一些想法，贴在下面，共同探讨： 1）分析混凝土结构，选择合理的材料特性是建立模型的关键，所以有必要弄清混凝土的材料特性。混凝土是脆性材料,并具有不同的拉伸和压缩特性。典型混凝土的抗拉强度只有抗压强度的8％-15%。在ANSYS中，对于混凝土单元，材料特性ANSYS要求输入以下数据（为了清楚起见，我将几个系数均译为了中文）：弹性模量、泊松比、张开与闭合滑移面的剪切强度缩减系数、抗拉与抗压强度、极限双轴抗压强度、周围静水应力状态、静水应力状态下单轴与双轴压缩的

ANSYS分析指南精华：子结构

第四章子结构什么是子结构？子结构就是将一组单元用矩阵凝聚为一个单元的过程。这个单一的矩阵单元称为超单元。在ANSYS分析中，超单元可以象其他单元类型一样使用。唯一的区别就是必须先进行结构生成分析以生成超单元。子结构可以在ANSYS/Mutiphysics，ANSYS/Mechanical和ANSYS/Structural中使用。使用子结构主要是为了节省机时，并且允许在比较有限的计算机设备资源的基础上求解超大规模的问题。原因之一如a)非线性分析和带有大量重复几何结构的分析。在非线性分析中，可以将模型线性部分作成子结构，这样这部分的单元矩阵就不用在非线性迭代过程中重复计算。在有重复几何结构的模型中（如有四条腿的桌子），可以对于重复的部分生成超单元，然后将它拷贝到不同的位置，这样做可以节省大量的机时。子结构还用于模型有大转动的情况下。对于这些模型，ANSYS假定每个结构都是围绕其质心转动的。在三维情况下，子结构有三个转动自由度和三个平动自由度。在大转动模型中，用户在使用部分之前无须对子结构施加约束，因为每个子结构都是作为一个单元进行处理，是允许刚体位移的。另外一个原因b)一个问题就波前大小和需用磁盘空间来说相对于一个计算 1

机系统太庞大了。这样，用户可以通过子结构将问题分块进行分析，每一块对于计算机系统来说都是可以计算的。如何使用子结构子结构分析有以下三个步骤： ●生成部分 ●使用部分 ●扩展部分生成部分就是将普通的有限元单元凝聚为一个超单元。凝聚是通过定义一组主自由度来实现的。主自由度用于定义超单元与模型中其他单元的边界，提取模型的动力学特性。图4-1是一个板状构件用接触单元分析的示意。由于接触单元需要迭代计算，将板状构件形成子结构将显著地节省机时。本例中，主自由度是板与接触单元相连的自由度。图4-1 子结构使用示例 2

ansys工程实例(4经典例子)解析

输气管道受力分析（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应用实例：钢筋混凝土简支梁数值模拟

（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：曲线形状能基本反映钢筋混凝土适筋梁剪切破坏的受力特点，而且荷载-跨中挠度曲线与钢筋混凝土梁的弯剪破坏形态非常类似，即当跨中弯矩最大截面的纵筋屈服后，由于裂缝的开展，压区混凝土的面积逐渐减小，在荷载几乎不增加的情况下，压区混凝土所受的正应力和剪应力还在不断增加，当应力达到混凝土强度极限时，剪切破坏发生，荷载突然降低。

ANSYS分析报告

《大型结构分析软件的应用及开发》学习报告学院：建筑工程学院专业班级：工程力学141 姓名：付贤凯指导老师：姚激学号：201411012111

1．模型介绍如下图所示的一桁架结构，受一集中力大小为800N的作用，杆件的弹性模量为200GPa，泊松比为0.3。杆件的截面为正方形达长为1m，横截面面积为1m2。现求它的变形图与轴力图。图1 桁架模型与受力简图（单位：mm） 2．建模与划分网格利用大型有限元软件ANSYS，采用Link，2Dspar 1的单元进行模拟，通过网格的划分得到如图2所示的有限元模型。图2 有限元模型

结合有限元模型中的约束条件为左侧在X与Y方向铰支固定，荷载条件为最右侧处施加向下的集中力P=800N。施加约束与荷载后的几何模型如图4所示。图3 施加荷载与约束的几何模型 3．位移与轴力图因在Y方向受力，所以主要做Y方向的位移图，又因为杆件在轴线方向有变形，故在X 方向仍有一定的位移。则图5为变形前后的板件形状。图6为模型沿Y方向的位移图，图7为模型沿X方向的位移图，图8为模型的总位移图。图4 桁架变形前后形状图

图5 Y方向位移图图6 X方向位移图

图7总位移图分析所有的位移图可以看出从以看出左端变形最小，为零，右端变形最大。从总位移图可以看出最大的位移在左下点处，大小为0.164×10?5m。从X方向位移图可以看出，左下点处在X方向位移最大为0.36×10?6。从Y方向位移图可以看出最大位移在左下点处为0.164×10?5。都符合实际情况，图9为模型的轴力图。图8 轴力图

ANSYS 中使用函数加载的一个简单例子

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 所示：

格构梁的ANSYS有限元模拟分析实例运用

龙源期刊网 https://www.360docs.net/doc/8316255890.html, 格构梁的ANSYS有限元模拟分析实例运用作者：张少剑刘真来源：《城市建设理论研究》2013年第10期摘要:本文通过一工程实例运用ansys模拟计算。针对格构梁的研究，合理地简化模型，取出1.5米宽的土体、梁和面层单元，两边加对称约束，从而达到模拟空间结构梁的目的。本文还模拟了基坑的开挖过程的时空效应，共分七步，土体在自重应力作用下的沉降为第一步，梁与面层的激活、力的施加和土层杀死共分六步。梁的最大受力状态并不发生在最后一步完成后，而是在第六工况。关键词:格构梁有限元分析模拟分析中图分类号：K826.16 文献标识码：A 文章编号： 1 土体、梁、锚索和混凝土面层共同作用基坑支护的受力机理是土体的土压力作用在格构梁和混凝土面层上，混凝土面层的力传递到格构梁上，格构梁再把它受到的力传递到和它相连的锚索上，锚索则和被支护土体嵌固为一体，格构梁和混凝土面层除起到承受土压力外，格构梁还起到平均弯矩和变形的作用，喷射混凝土面层则有保护土体表面，防止土体表面非格构梁作用部位坍塌的作用。 2模型简化及技术处理根据基坑开挖深度，根据实际的土体性质建立土体模型。格构梁的作用是承受弯矩的，可以选用Beam4梁单元，考虑到钢筋混凝土格构梁中有钢筋的作用，其弹性模量、泊松比等设置有所调整。在建模时，如果混凝土面层的长宽与厚度的比都大于5，所以在有限元分析中采用板壳单元可以全面地反映其变形特征和应力分布规律。混凝土面层用Shell63单元模拟，其参数的取值和梁单元相同。由于格构梁的受力性状，锚索的模拟对格构梁的受力影响较小，本模型忽略考虑锚索的模拟。预应力锚索的作用简化为作用在纵横梁交点处的集中力。对于格构梁和土体、混凝面层之间的接触，模型采用节点耦合，以实现共同变形和受力。 3.1ANSYS有限元模拟计算 3.1.1模型的参数 1．土体的参数见下表：

ANSYS分析报告分析

有限元与CAE分析报告专业：班级：学号：姓名：指导教师： 2016年 1 月 2 日

简支梁的静力分析一、问题提出长3m的工字型梁两端铰接中间1.5m位置处受到6KN的载荷作用，材料弹性模量E=200e9，泊松比0.28，密度7850kg/㎡二、建立模型 1.定义单元类型依次单击Main Menu→Preprocessor→Elementtype→Add/Edit/Delete，出现对话框如图，单击“Add”，出现一个“Library of Element Type”对话框，在“Library of Element Type”左面的列表栏中选择“Structural Beam”，在右面的列表栏中选择3 node 189,单击“OK”。

2设置材料属性依次单击Main Menu→Preprocessor→MaterialProps>Material Modes，出现“Define Material ModelBehavior”对话框，在“Material Model Available”下面的对话框中，双击打开“Structural→Linear→Elastic→Isotropic”,出现对话框，输入弹性模量EX=2E+011,PRXY=0.28，单击“OK”。依次单击Main Menu→Preprocessor→MaterialProps>Material Modes，出现“Define Material ModelBehavior”对话框，在“Material Model Available”下面的对话框中，双击打开“Structural→Density”弹出对话框，输入DENS为7850 3.创建几何模型 1)设定梁的截面尺寸

ANSYS-结构稳态(静力)分析之经典实例-命令流格式

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