ANSYS中混凝土的本构关系

ANSYS中混凝土的本构关系
ANSYS中混凝土的本构关系

一、关于模型

钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种,即分离式、分布式和组合式模型。考虑钢筋和混凝土之间的粘结和滑移,则采用引入粘结单元的分离式模型;假定混凝土和钢筋粘结很好,不考虑二者之间的滑移,则三种模型都可以;分离式和分布式模型适用于二维和三维结构分析,后者对杆系结构分析比较适用。裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型,后者目前尚处研究之中,主要应用的是前两种。离散裂缝模型和分布裂缝模型各有特点,可根据不同的分析目的选择使用。随着计算速度和网格自动划分的快速实现,离散裂缝模型又有被推广使用的趋势。

就ANSYS而言,她可以考虑分离式模型(solid65+link8,认为混凝土和钢筋粘结很好,如要考虑粘结和滑移,则可引入弹簧单元进行模拟,比较困难!),也可采用分布式模型(带筋的solid65)。而其裂缝的处理方式则为分布裂缝模型。

二、关于本构关系

混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类,其中研究最多的是非线性弹性和弹塑性本构关系,其中不乏实用者。混凝土破坏准则从单参数到五参数模型达数十个模型,或借用古典强度理论或基于试验结果等,各个破坏准则的表达方式和繁简程度各异,适用范围和计算精度差别也比较大,给使用带来了一定的困难。就ANSYS而言,其问题比较复杂些。

1 ANSYS混凝土的破坏准则与屈服准则是如何定义的

采用tb,concr,matnum则定义了W-W破坏准则(failure criterion),而非屈服准则(yield criterion)。W-W破坏准则是用于检查混凝土开裂和压碎用的,而混凝土的塑性可以另外考虑(当然是在开裂和压碎之前)。理论上破坏准则(failure criterion)和屈服准则(yield criterion)是不同的,例如在高静水压力下会发生相当的塑性变形,表现为屈服,但没有破坏。而工程上又常将二者等同,其原因是工程结构不容许有很大的塑性变形,且混凝土等材料的屈服点不够明确,但破坏点非常明确。

定义tb,concr matnum后仅仅是定义了混凝土的破坏准则和缺省的本构关系,即W—W 破坏准则、混凝土开裂和压碎前均为线性的应力应变关系,而开裂和压碎后采用其给出的本构关系。但屈服准则尚可另外定义(随材料的应力应变关系,如tb,MKIN,则定义的屈服准则是Von Mises,流动法则、硬化法则也就确定了)。

2 定义tb,concr后可否定义其它的应力应变关系

当然是可以的,并且只有在定义tb,concr后,有些问题才好解决。例如可以定义tb,miso,输入混凝土的应力应变关系曲线(多折线实现),这样也就将屈服准则、流动法则、硬化法则等确定了。

这里可能存在一点疑问,即ANSYS中的应力应变关系是拉压相等的,而混凝土材料显然不是这样的。是的,因为混凝土受拉段非常短,认为拉压相同影响很小,且由于定义的tb,concr中确定了开裂强度,所以尽管定义的是一条大曲线,但应用于受拉部分的很小。

三、具体的系数及公式

1 定义tb,concr时候的两个系数如何确定

一般的参考书中,其值建议先取为~(江见鲸),原话是“在没有更仔细的数据时,不妨先取~进行计算”,足见此~值的可用程度。根据我的经验和理由,建议此值取大些,即开裂的剪力传递系数取,(定要>)闭合的剪力传递系数取。支持此说法的还有现行铁路桥规的抗剪计算理论,以及原公路桥规的容许应力法的抗计剪计算。

2 定义混凝土的应力应变曲线

单向应力应变曲线很多,常用的可参考国标混凝土结构规范,其中给出的应力应变曲线是二次曲线+直线的下降段,其参数的设置按规范确定即可。当然如有实测的应力应变曲线更好了。

四、关于收敛的问题

ANSYS混凝土计算收敛(数值)是比较困难的,主要影响因素是网格密度、子步数、收敛准则等,这里讨论如下。

1 网格密度:网格密度适当能够收敛。不是网格越密越好,当然太稀也不行,这仅仅是就收敛而言的,不考虑计算费用问题。但是究竟多少合适,没有找到规律,只能靠自己针对情况慢慢试算。

2 子步数:NSUBST的设置很重要,设置太大或太小都不能达到正常收敛。这点可以从收敛过程图看出,如果F范数曲线在[F]曲线上面走形的很长,可考虑增大nsubst。或者根据经验慢慢调正试算。

3 收敛精度:实际上收敛精度的调正并不能彻底解决收敛的问题,但可以放宽收敛条件以加速吧。一般不超过5%(缺省是%),且使用力收敛条件即可。

4 混凝土压碎的设置:不考虑压碎时,计算相对容易收敛;而考虑压碎则比较难收敛,即便是没有达到压碎应力时。如果是正常使用情况下的计算,建议关掉压碎选项;如果是极限计算,建议使用concr+MISO且关闭压碎检查;如果必设压碎检查,则要通过大量的试算(设置不同的网格密度、NSUBST)以达到目的,但也很困难。

5 其他选项:如线性搜索、预测等项也可以打开,以加速收敛,但不能根本解决问题。

6 计算结果:仅设置concr,不管是否设置压碎,其一般P-F曲线接近二折线;采用concr+miso则P-F曲线与二折线有差别,其曲线形状明显是曲线的。

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)需要参数: 单轴抗拉强度,单轴,双轴抗压强度,围压压力,在围压作用下双轴,单轴抗压强度

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分析

摘要 本文介绍ANSYS 模拟钢筋混凝土梁的过程,讨论了有限元模型的建立以及在 ANSYS 中的实现,给出了用分离式配筋方法对混凝土梁的分析的一般过程。并给出了详细的命令流过程。并在此基础上对混凝土梁进行了分析,讨论了在力的作用下混凝土梁的塑形变形和裂缝的发展过程。 关键词 Ansys 混凝土梁 分离式配筋 The analysis of mechanics of a reinforced concrete based on ANSYS Abstract This paper introduces ANSYS simulation of the reinforced concrete beam process, discusses the establishment of the finite element model and the realization, and gives the ANSYS reinforcement method with separate the analysis of concrete beams of the general process. And gives the detailed command flow process. Based on the analysis of concrete beams, and discussed the concrete beam under the action of forces of the body deformation and fracture process. Keywords Ansys concrete beams reinforced separated 1 引言 由于钢筋混凝上材料性质复杂,使其表现出明显的非线性行为[1]。长期以来采用线弹性理论的设计方法来研究钢筋混凝上结构的应力或内力,显然不太合理,尽管有此理论是基于人量试验数据上的经验公式,还是不能准确反映混凝上的力学性能,特别是受力复杂的重要结构,必须采用三维钢筋混凝上非线性有限元方法才能很好地掌握其力学性能。利用ANSYS 对钢筋混凝上结构弹塑性的仿真分析,可以对结构自开始受荷载直到破坏的全过程进行分析,获得不同阶段的受力性能。本文将以混凝土梁的弹塑性分析为例,介绍在Ansys 中分析材料非线性问题的具体实现方法。 2 问题介绍 如图所示的钢筋混凝土梁[2],横截面尺寸为200400b h mm mm ?=?,梁的跨度为3.0L m =,支座宽度为250mm 采用C20混凝土,梁内受拉纵筋3φ20,架立筋采用2φ12, 箍筋采用φ6@150,钢筋保护层厚度为25mm 。如图一。 图一 对于梁中所采用的所有钢筋,弹性模量为5 2.110MPa ?,抗拉强度设计值210MPa , 密度33 7.810/kg m ?,泊松比为0.3。

有限元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中混凝土的本构关系

一、关于模型 钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种,即分离式、分布式和组合式模型。考虑钢筋和混凝土之间的粘结和滑移,则采用引入粘结单元的分离式模型;假定混凝土和钢筋粘结很好,不考虑二者之间的滑移,则三种模型都可以;分离式和分布式模型适用于二维和三维结构分析,后者对杆系结构分析比较适用。裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型,后者目前尚处研究之中,主要应用的是前两种。离散裂缝模型和分布裂缝模型各有特点,可根据不同的分析目的选择使用。随着计算速度和网格自动划分的快速实现,离散裂缝模型又有被推广使用的趋势。 就ANSYS而言,她可以考虑分离式模型(solid65+link8,认为混凝土和钢筋粘结很好,如要考虑粘结和滑移,则可引入弹簧单元进行模拟,比较困难!),也可采用分布式模型(带筋的solid65)。而其裂缝的处理方式则为分布裂缝模型。 二、关于本构关系 混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类,其中研究最多的是非线性弹性和弹塑性本构关系,其中不乏实用者。混凝土破坏准则从单参数到五参数模型达数十个模型,或借用古典强度理论或基于试验结果等,各个破坏准则的表达方式和繁简程度各异,适用范围和计算精度差别也比较大,给使用带来了一定的困难。 就ANSYS而言,其问题比较复杂些。 1 ANSYS混凝土的破坏准则与屈服准则是如何定义的? 采用tb,concr,matnum则定义了W-W破坏准则(failure criterion),而非屈服准则(yield criterion)。W-W破坏准则是用于检查混凝土开裂和压碎用的,而混凝土的塑性可以另外考虑(当然是在开裂和压碎之前)。理论上破坏准则(failure criterion)和屈服准则(yield criterion)是不同的,例如在高静水压力下会发生相当的塑性变形,表现为屈服,但没有破坏。而工程上又常将二者等同,其原因是工程结构不容许有很大的塑性变形,且混凝土等材料的屈服点不够明确,但破坏点非常明确。 定义tb,concr matnum后仅仅是定义了混凝土的破坏准则和缺省的本构关系,即W—W破坏准则、混凝土开裂和压碎前均为线性的应力应变关系,而开裂和压碎后采用其给出的本构关系。但屈服准则尚可另外定义(随材料的应力应变关系,如tb,MKIN,则定义的屈服准则是Von Mises,流动法则、硬化法则也就确定了)。 2 定义tb,concr后可否定义其它的应力应变关系 当然是可以的,并且只有在定义tb,concr后,有些问题才好解决。例如可以定义tb,miso,输入混凝土的应力应变关系曲线(多折线实现),这样也就将屈服准则、流动法则、硬化法则等确定了。 这里可能存在一点疑问,即ANSYS中的应力应变关系是拉压相等的,而混凝土材料显然不是这样的。是的,因为混凝土受拉段非常短,认为拉压相同影响很小,且由于定义的tb,concr 中确定了开裂强度,所以尽管定义的是一条大曲线,但应用于受拉部分的很小。 三、具体的系数及公式 1 定义tb,concr时候的两个系数如何确定? 一般的参考书中,其值建议先取为0.3~0.5(江见鲸),原话是“在没有更仔细的数据时,不妨先取0.3~0.5进行计算”,足见此0.3~0.5值的可用程度。根据我的经验和理由,建议此值取大些,即开裂的剪力传递系数取0.5,(定要>0.2)闭合的剪力传递系数取1.0。支持此说法的还有 现行铁路桥规的抗剪计算理论,以及原公路桥规的容许应力法的抗计剪计算。

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 钢筋混凝土建模

ANSYS 钢筋混凝土建模 一、简介 钢筋混凝土有限元建模的方法与结果评价(前后处理),是对钢筋混凝土结构进行数值模拟的重要步骤,能否把握模型的可行性、合理性,如何从计算结果中寻找规律,是有限元理论应用于实际工程的关键一环。Blackeage以自己做过的一组钢筋混凝土暗支撑剪力墙的数值模拟为例,从若干方面提出一些经验与建议。希望大家一起讨论、批评指正(wang.jian@https://www.360docs.net/doc/3312452591.html,)。 程序:ANSYS 单元:SOLID65、BEAM188 建模方式:分离 暗支撑剪力墙结构由北京工业大学曹万林所提出,简言之就是一种在普通钢筋配筋情况下,加配斜向钢筋的剪力墙结构。 二、单元选择 以前经常采用的钢筋混凝土建模方法是通过SOLID65模拟混凝土,通过SOLID65的实常数指定钢筋配筋率,后来发现这种整体式的模型并不理想,而且将钢筋周围的SOLID65单元选择出来,再换算一个等效的配筋率,工作量也并不小。最关键的是采用整体式模型之后,得不出什么有意义的结论,弄一个荷载-位移曲线出来又和实验值差距比较大。只有计算的开裂荷载与实验还算是比较接近,但这个手算也算得出来的东西费劲去装模作样的建个模型又有什么意义? 所以,这次我尝试采用分离式的模型,钢筋与混凝土单元分别建模,采用节点共享的方式。建模时发现,只要充分、灵活地运用APDL的技巧,处理好钢筋与混凝土单元节点的位置,效率还是很高的。 暗支撑剪力墙数值模型 看过很多的资料,分离式模型是用LINK8与SOLID65的组合方式,这样做到是非常直观,因为LINK8是spar类型的单元,每个节点有3个自由度,这与SOLID65单元单节点自由度数量是一致的。但是问题也就由此产生,当周围的混凝土开裂或是压碎时,SOLID65将不能对LINK8的节点提供足够地约束(如

七个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中混凝土的计算问题

ANSYS中混凝土的计算问题【精华】 最近做了点计算分析,结合各论坛关于这方面的讨论,就一些问题探讨如下,不当之处敬请指正。 一、关于模型 钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种,即分离式、分布式和组合式模型。考虑钢筋和混凝土之间的粘结和滑移,则采用引入粘结单元的分离式模型;假定混凝土和钢筋粘结很好,不考虑二者之间的滑移,则三种模型都可以;分离式和分布式模型适用于二维和三维结构分析,后者对杆系结构分析比较适用。裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型,后者目前尚处研究之中,主要应用的是前两种。离散裂缝模型和分布裂缝模型各有特点,可根据不同的分析目的选择使用。随着计算速度和网格自动划分的快速实现,离散裂缝模型又有被推广使用的趋势。 就ANSYS而言,她可以考虑分离式模型(solid65+link8,认为混凝土和钢筋粘结很好,如要考虑粘结和滑移,则可引入弹簧单元进行模拟,比较困难!),也可采用分布式模型(带筋的solid65)。而其裂缝的处理方式则为分布裂缝模型。 二、关于本构关系 混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类,其中研究最多的是非线性弹性和弹塑性本构关系,其中不乏实用者。混凝土破坏准则从单参数到五参数模型达数十个模型,或借用古典强度理论或基于试验结果等,各个破坏准则的表达方式和繁简程度各异,适用范围和计算精度差别也比较大,给使用带来了一定的困难。 就ANSYS而言,其问题比较复杂些。 1 ANSYS混凝土的破坏准则与屈服准则是如何定义的? 采用tb,concr,matnum则定义了W-W破坏准则(failure criterion),而非屈服准则(yield criterion)。W-W破坏准则是用于检查混凝土开裂和压碎用的,而混凝土的塑性可以另外考虑(当然是在开裂和压碎之前)。理论上破坏准则(failure criterion)和屈服准则(yield criterion)是不同的,例如在高静水压力下会发生相当的塑性变形,表现为屈服,但没有破坏。而工程上又常将二者等同,其原因是工程结构不容许有很大的塑性变形,且混凝土等材料的屈服点不够明确,但破坏点非常明确。 定义tb,concr matnum后仅仅是定义了混凝土的破坏准则和缺省的本构关系,即W—W破坏准则、混凝土开裂和压碎前均为线性的应力应变关系,而开裂和压碎后采用其给出的本构关系。但屈服准则尚可另外定义(随材料的应力应变关系,如tb,MKIN,则定义的屈服准则是V on Mises,流动法则、硬化法则也就确定了)。

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 轴力图

相关文档
最新文档