ansys非线性分析例子
ANSYS-6-非线性分析应用
第六章 钢筋混凝土结构非线性分析应用§6.1截面非线性分析例 1: 钢筋混凝土单筋矩形截面，混凝土和钢筋的应力-应变关系选自CEB 模型规范（1990），见下图6-1-1，图 6.1-1 截面和材料应力-应变关系极限弯矩 M u : 用弧长法对截面进行全过程分析，对给定的弯矩M y , 计算相应的截面应变平面({}[]T z y ϑϑεε0=).计算不平衡弯矩及相应的应变平面增量，直至满足收敛条件。
再增加弯矩∆M y , 计算相应的应变平面增量，等等，图6-1-2为截面弯矩-曲率关系曲线。
图 6.1-2 弯矩-曲率关系曲线 例2: 采用不同应力-应变关系（EC2规范, CEB 规范），钢筋混凝土矩形截面的几何尺寸和配筋同例1，非线性分析结果见图6-1-4。
力-应变关系随应变而逐渐的降低，截面刚度降低的也比较缓慢。
图 6.1-4 CEB 规范与EC2 规范建议的应力-应变关系截面分析结果比较例 3: 异形截面非线性分析. 此例Georg Knittel [32]计算过，Knittel 选择的材料应力应变关系取自德国规范DIN 1045(见图 6.1-5). 截面形状和尺寸见图6.1-6. Knittel 分析的截面极限承载力为，{}{}N M M y z T T=--005026000075... 相应的应变矢量为，{}{}{}TT z y 009343.0006976.0004359.00--==ϑϑεε. 用弧长法分析时取的参照荷载值为，{}{}N M M yz T T =--00050026000075... 截面极限荷载为，{}{}N M M y z T T =--004991490263211600076718...(a) DIN 1045建议的混凝土应力-应变关系 (b) DIN 1045建议的钢筋应力-应变关系图 6.1-5 DIN 1045规范建议的应力-应变关系图 6.1-6 钢筋混凝土柱截面图 6.1-7 极限状态时混凝土压应力分布图 6.1-8 弯矩-曲率(M y- y) 关系曲线§6.2 受弯和偏压构件非线性分析6.2.1 简化计算利用虚功原理计算荷载挠度曲线：设两点集中加载简支梁，弯矩图、曲率分布图如下，图6-2-1 梁内力与变形取支撑条件相同的简支梁为虚梁，拟求跨中挠度，在虚梁跨中施加单位荷载（求转角加单位力矩）。
ansys高级非线性分析-第九章几何不稳定性
失稳准则、米泽斯失稳准则和霍夫失稳准则等。
这些判据可以帮助我们确定结构的临界载荷和失稳模式，从而
03
采取相应的措施来提高结构的稳定性。
几何不稳定性的影响因素
材料性质
材料的弹性模量、泊松比、屈 服强度等都会影响结构的稳定
性。
结构形状和尺寸
结构的形状、尺寸、支承条件 等都会影响其稳定性。
外部载荷
外部载荷的大小、方向和分布 也会影响结构的稳定性。
案例二：高层建筑的几何不稳定性分析
总结词
高层建筑的几何不稳定性分析是确保高层建筑结构安全的重要环节。
详细描述
利用ANSYS的高级非线性分析功能，可以对高层建筑在不同风载、地震等载荷作 用下的结构响应进行模拟，评估其稳定性和安全性，为设计提供依据。
案例三：重型机械的几何不稳定性分析
总结词
重型机械的几何不稳定性分析是确保 重型机械在各种工况下安全运行的关 键。
02
几何不稳定性分析在复杂边界条件、多物理场耦合等方面的研究尚不够深入， 需要进一步拓展研究范围，完善分析方法。
03
随着计算机技术和数值计算方法的不断发展，几何不稳定性分析的计算效率和 精度将得到进一步提高，为工程实际提供更加准确和可靠的理论支持。
THANKS
感谢观看
现象。
在非线性分析中，需要考虑 结构在变形过程中形状和尺 寸的变化，以及由此引起的
力和位移的重新分布。
几何非线性行为通常出现在大 变形、应力刚化、旋转软化和
塑性流动等情况下。
几何不稳定性判据
01
几何不稳定性是指结构在某些条件下失去稳定性，发生屈曲或 失稳的现象。
02
判据是用来判断结构是否稳定的准则，常用的判据包括：欧拉
ANSYS材料非线性分析
【分享】ANSYS7.0超弹材料的定义－新的曲线拟合功能--摘自ansys用户专区几何非线性几何非线性不受敛主要原因1.网格质量，特别是warpage2.约束方程，少用刚性连接3.收敛准则，可适当加大容差4.荷载步设置，可适当加大步数最近碰到一个对我来说很意外的问题：如果确实如此希望大家以后小心大家知道定义接触后会自动生成一组实常数，前几天我碰到一个问题，需定义超过10组实常数，接触对很多，好像有20多处，按照常规步骤划分完所有网格，当时因为有一个实常数参数没确定，便预留了最后一组（第10组）实常数里面的参数为空，接下来就定义了所有的接触对，由于所有接触对里的设置一样，ANSYS在我保存db完重新打开后便把我所有的接触对综合成一个了！接下来我就把第十组实常数里面的参数补上了，但在求解时却提示我该实常数同时被两种单元（包括CNTACT单元）同时占用，出现错误！！检查了半天才发现自动生成的接触对实常数把第10组实常数也占用了！我实在没找到什么好的解决办法，只得把接触对删除了重新定义，那可是上百多个面的选取过程，痛苦不堪简直！ANSYS里接触对面的选取时还不能针对Component操作！ANSYS7.0超弹材料的定义－新的曲线拟合功能ANSYS7.0中的超弹材料模拟能力得到了很大的加强，在ANSYS6.1的超弹材料模型的基础上又增加了Gent, Yeoh, Blatz-Ko, and Ogden (Foam)四种超弹性材料模型，使得其超弹模拟能力得到了进一步扩展。
ANSYS7.0中对超弹能力最吸引人的增强还不在于此，而是在于其曲线拟合能力的大幅度扩展，不再像ANSYS6.1以前的版本一样曲线拟合仅仅局限于Mooney-Rivlin模型，而是将其扩展到所有的超弹模型，这样，用户可以利用实验得到的应力应变数据直接让程序自己拟合出任意一种超弹材料模型的参数，大大方便了用户的使用。
以下就ANSYS7.0的超弹拟合功能做一简单介绍。
Ansys Workbench非线性分析 牛顿辛普森法过程
然而, 相当多的结构在力和位移之间没有线性关系 因为此类结构的 F-u 图不是直线, 这样的结构称为非线性结构 . – 刚度不再是一个常数K; 它成为施加载荷的函数, KT (切 线刚度). 普通的例子是韧性金属的拉伸试验:
载荷位移曲线
Fa Fnr
R
下一次迭代用刚度矩阵
u0 △ u u1
u0: 所设初始位移值
u
位移
KT :切向刚度
u1: 下次迭代位移值
R F a F nr
收敛： Newton-Raphson法需要一个收敛的度量以决定 何时结束迭代。给定节点Fa，节点力Fnr ，在一个体 中，节点载荷必须与节点力平衡
高级接触选项包括： 自动探测尺寸 Auto detection dimension and slider 非对称接触Asymmetric contact 接触结果工具Contact results tool 接触算法contact formulations Pinball 控制
膜片弹簧接触设置
带支撑环的膜片弹簧边界条件及求解设置
3 2 F1 269.991 2981 1 9238 1
膜片弹簧大端载荷-位移曲线
3 2 F1 269.991 2981 1 9238 1
注：膜片弹簧载荷挠度曲线为非线性，是由于大变形引起 的几何非线性，但材料始终发生的是弹性变形，没有发生 塑性变形
回忆材料力学中低碳钢力学性能试验
ANSYS讲义非线性分析
t1
t2
时间 t
XJTU
自动时间步(续)
• 自动时间步算法是 非线性求解控制 中包含的多种算法的一种。
（在以后的非线性求解控制中有进一步的讨论。） • 基于前一步的求解历史与问题的本质，自动时间步算法或者增加
或者减小子步的时间步大小。
XJTU
5) 输出文件的信息
在非线性求解过程中，输出窗口显示许多关于收敛的信息。输出 窗口包括：
子步
时间 ”相关联。
“时间
两个载荷步的求解 ”
XJTU
在非线性求解中的 “ 时间 ”
• 每个载荷步与子步都与 “ 时间 ”相关联。 子步 也叫时间步。
• 在率相关分析（蠕变，粘塑性）与瞬态分析中，“ 时间 ”代表真实 的时间。
• 对于率无关的静态分析，“ 时间 ” 表示加载次序。在静态分析中， “ 时间 ” 可设置为任何适当的值。
最终结果偏离平衡。
u 位移
XJTU
1) Newton-Raphson 法
ANSYS 使用Newton-Raphson平衡迭代法 克服了增量
求解的问题。 在每个载荷增量步结束时，平衡迭代驱 使解回到平衡状态。
载荷
F
4 3 2
1
u 位移
一个载荷增量中全 Newton-Raphson 迭代 求解。（四个迭代步如 图所示）
XJTU
非线性分析的应用(续)
宽翼悬臂梁的侧边扭转失 稳
一个由于几何非线性造 成的结构稳定性问题
XJTU
非线性分析的应用(续)
橡胶底密封 一个包含几何非线 性（大应变与大变 形），材料非线性 （橡胶），及状态 非线性（接触的例 子。
XJTU
非线性分析的应用(续)
材料非线性接触设置实例_ANSYS Workbench有限元分析实例详解（静力学）_[共11页]
![材料非线性接触设置实例_ANSYS Workbench有限元分析实例详解（静力学）_[共11页]](https://img.taocdn.com/s3/m/9f18735f84868762cbaed53c.png)
第5章 非线性静力学分析– 424 –同理，右键点击Connections 插入Connection Group4，隐含core 、solid 、shell 三个零件，在Geometry 处选择剩下的13个零件，调整Tolerance Slider 为0，然后自动生成接触，如图5-3-83所示。
图5-3-83 接触设置4由于产生的接触对较多，不可避免会有接触对重复现象，所以生成完所有接触对之后，右键点击Connections →Search Connections for Duplicate Pairs ，软件将自动检查重复接触对，然后手动删除，多次检查，直到出现：no connections with duplicate pairs have been found 提示。
接触对全部定义完以后，可以再统一修改接触类型为No Separation 或Frictionless ；或者在定义接触之前，就修改Tool →Option →Mechanical →Connections →Type 为No Separation 或Frictionless ，这样定义的接触对默认为不分离或无摩擦。
4．小结对于复杂零件的接触设置，通过定义多个Connection Group ，定义不同的接触公差，可以有效地提高软件自定义接触的准确性。
如果是更加复杂的整件，通过External Model 模块可以装配有限元模型（支持主流有限元软件的网格文件，且不受版本限制）。
该模块可以保留网格文件中的命名选择、网格控制，如果是装配部件，还会保留接触对设置，可以对Solid 单元、Shell 单元的高阶、低阶单元模型进行装配，而且ANSYS 后续版本都在强化该模块装配后的智能操作。
5.3.8 材料非线性接触设置实例接触分析过程中，往往伴随着材料非线性特征，这两种非线性结合在一起，极易不收敛。
初学者在学习过程中，由于参照例子一步一步操作，知其然不知其所以然，造成面临实际不收敛问题时，往往不知所措。
ansys非线性接触分析中的接触行为
ansys非线性接触分析中接触行为接触是状态改变非线性，经典ANSYS版本中共提供了7种接触行为，每一种都有其特点及相应的应用范围，在选用的时候应该谨慎。
（1）标准接触行为（standard）该接触行为包括了法向接触闭合和分开行为，在该接触模式中既考虑粘着摩擦同时也考虑了滑动摩擦。
如图上，AB与BC本来是分开的，中间通过B点连接，当在A点施加力F，AB慢慢贴近BC，最终靠在一起。
但F撤销后，AB在恢复力的作用下慢慢回复到初始分开状态。
标准接触行为包括了分开状态→闭合状态→分开状态。
当AB与BC靠在一起时，既存在正压力，同时还有沿BC圆弧切线方向的摩擦力。
（2）粗糙接触行为（rough）该接触行为包括了法向接触闭合和分开行为，但滑动行为在此是不会发生的。
原因是所有参与接触的表面都被假定为非常粗糙，以致于可以认为摩擦力无穷大而不能够产生相对滑动。
在这种接触行为中，接触的两个物体或部件之间，除了存在正压力外，还有切向摩擦力，但是接触部分之间不可以产生相对滑动。
（3）绑定接触行为（bonded）是指一旦接触关系建立，那么目标面及接触面就被假定为粘结在一起（不可以分开）。
（4）绑定接触行为（始终）（bonded（always））任何初始时在许可接触容差范围内探测到的接触点或者是那些即将进入接触的点在后续的分析中将被绑定在一起。
这种接触行为的典型应用，如在组装分析中将两种不同网络的组件“加”在一起。
线性静态分析也可以用该种接触行为来解决，虽然由于有接触单元的存在，分析中将会提示为非线性分析，但往往只要一步迭代就完成了。
（5）绑定接触行为（初始接触）（bonded（initial））绑定仅发生在初始状态下就接触的面上，初始状态下没有接触的部分将继续保持分开。
典型的例子是通过焊接连接在一起的两个物体，焊接部分始终保持连接，没有焊接的部分保持分离状态。
（6）不分开型（no separation）一旦接触关系建立，目标面及接触面便被约束在一起了，但还是允许接触面之间有滑动。
ANSYS非线性接触问题分析汇总
ANSYS 分线性接触问题分析汇总接触非线性是一门复杂的学科，ANSYS 关于计算非线性接触的设置选项多只又多，很多人摸不到头脑，本文就基于ANSYS 模拟过的几个接触实例，研究了相关设置选项对接触结果的影响。
实例1：橡胶密封圈配合接触研究—非线性求解设置对结果的影响密封圈配合模型简图见图1，左右两端为刚体，中间圆部分为橡胶密封圈，将刚体2沿刚体1方面移动，从而实现橡胶圈密封作用，采用plane182单元，设置轴对称行为，建立橡胶密封圈与刚体接触模型，见图2。
图1 密封圈配合模型简图 图2 密封圈配合有限元模型图接触对采用默认设置，摩擦系数取0.10，研究非线性求解器设置对收敛方面的影响，大变形静态（Large Displacement Static ）效应打开，自动时间步长（Automatic time stepping ）打开，子步数（Number of substeps ）设置为50，线性搜索（Line search ）打开。
1 收敛准则对结果的影响此实例收敛准则默认采用力收敛结合力矩收敛准则（基于L2范数），收敛容差（Tolerance ）默认为0.001，工程上认为0.05的收敛容差足够满足要求。
表 1 收敛容差对计算结果的影响收敛容差 最大应力/ MPa报错与否？ 0.001 4.12364报错 0.05 4.12785 报错 0.14.12996报错查看报错信息，见图3，表示单元过于扭曲，建议提高子步数或降低时间步长，需要提高网格质量，也要考虑材料属性，接触对及约束方程的合理性，若在第一步迭代就如此，需要预先执行单元形状检查。
图3 报错信息刚体1刚体2密封圈橡胶密封圈配合Von Mises应力云图见图4。
图4 橡胶密封圈配合Von Mises应力2 子步数对结果的影响此实例子步数设置为50、100、200、500，收敛容差（Tolerance）默认为0.001，研究子步数对收敛的影响。
ANSYS材料非线性分析
弹簧单元模拟土壤
杆为平面应变分析
土壤单元选项及实常数
材料性能：d-p模型（土壤）
铁桩材料
布尔运算：叠合。
并压缩实体编号（方便划网）
选择所有
all
相加分离的线
delte 划网注意划网属性与面编号对应
定义映射的边再划分面
转换为柱坐标系
创建节点并复制
通过节点再创建单元，注意选设置好单元属性
选择r=180（圆边界节点）(from all)及外节点180*(2*1.732)/3(also select)进行全约束
转换为直角坐标
对桩进行压力加载500（Mpa）先选择节点（通过坐标值）
创建静态分析
并设置非线性收敛标准
求解保存
结果（应力、位移、剪力）。
ANSYS-1-非线性分析概述
第一章钢筋混凝土结构非线性分析概述1.1 钢筋混凝土结构的特性1.钢筋混凝土结构由两种材料组成，两者的抗拉强度差异较大，在正常使用阶段，结构或构件就处在非线性工作阶段，用弹性分析方法分析的结构内力和变形无法反映结构的真实受力状况；2.混凝土的拉、压应力-应变关系具有较强的非线性特征；3.钢筋与混凝土间的黏结关系非常复杂，特别是在反复荷载作用下，钢筋与混凝土间会产生相对滑移，用弹性理论分析的结果不能反映实际情况；4.混凝土的变形与时间有关：徐变、收缩；5.应力-应变关系莸软化段：混凝土达到强度峰值后有应力下降段；6.产生裂缝以后成为各向异形体。
混凝土结构在荷载作用下的受力-变形过程十分复杂，是一个变化的非线性过程。
1.2 混凝土结构分析的目的和主要内容《混凝土结构设计规范》（GB50010-2002）中新增的主要内容：(1)混凝土的本构关系和多轴强度：给出了单轴受压、受拉非线性应力-应变（本构）关系，混凝土二轴强度包络图、三轴抗压强度图和三轴应力状态破坏准则；（2）结构分析：规范概括了用于混凝土结构分析的5类方法，列入了结构非线性分析方法。
一、结构分析的基本目的：计算在各类荷载作用下的结构效应——内力、位移、应力、应变根据设计的结构方案确定合理的计算简图，选择不利荷载组合，计算结构内力，以便进行截面配筋计算和采取构造措施。
二、结构分析的主要内容：（1）确定结构计算简图：考虑以下因素：(a)能代表实际结构的体形和尺寸；(b)边界条件和连接方式能反映结构的实际受力状态，并有可靠的构造措施；(c)材料性能符合结构的实际情况；(d)荷载的大小、位置及组合应与结构的实际受力吻合；(e)应考虑施工偏差、初始应力及变形位移状况对计算简图进行适当修正；(f)根据结构受力特点，可对计算简图作适当简化，但应有理论或试验依据，或有可靠的工程经验；(g)结构分析结果应满足工程设计的精度要求。
（2）结构作用效应分析：根据结构施工和使用阶段的多种工况，分别进行结构分析，确定最不利荷载效应组合。
基于ANSYS的索—混凝土组合梁的非线性分析
ｍ ｎｓ ｅｔ建模 ， 不用划分 网格。６施加边界条件。需要注意的是 ： ） 在划
， … 】 … ， ’ ｏ … １ Ｉ， ； ， … ， … ， … ， ， ，
Ｂａ ｅ ｎ ＡＮＳ ｉ ｉｅ ｅ ｅ ｅ ｎａ ｙ ｉ ｆ ｔ ｅ ｌ ｉｎ ｏ ｕ ｓ ａ ｐ ｐｅ ｉ ｅ ｐ ｌｉ ｅ ｄ ｓｄｏ ＹＳ ｆｎ ｔ ｌ ｍ ｎｔａ ｌ ｓｓ ｏ ｈｅ ｒ ｃ ａ ｌ ｆｓ ｂ ｅ ｉ ｌｎ ｕ ｌｎｇ ｈ ａ ＇
基 于 ＡＮ Ｙ Ｓ Ｓ的索一 混凝 土 组 合 梁 的 非线 性 分析
范思冬
摘
秦乃兵
张
雷
要： 主要针 对索一 混凝土 组合 梁这 一新 型结 构， 用大 型有 限元分析 软件 Ａ ＹＳ进行 了非 线性分析 ， 利 ＮＳ 研究 了索一
混 凝 土 组 合 梁 在 一定 荷 载 作 用 下 的 受力 模 型 以及 其 部 件 的拉 压 情况 ， 通 过 对 比 证 实 了索一 混凝 土 组 合 梁 的优 越 性 。 并 关键词 ： 索一 混凝 土 组合 梁 ， 线 性 ， 限元 分 析 非 有
合梁结构 的静力特性 的非线 性 阶段进行 分析 。由于混凝 土承受
表 ２ 混 凝 土 的 基本 参数
张开裂缝剪切传递系数 闭合裂缝剪切传递系数 单轴极限抗拉强度／ ａ ＭＰ 单轴极限抗压强度／ ａ ＭＰ 等压双轴抗压强度 周围静水压力 静水压力下 双轴极限抗压强度 ０ ５ ．５ ０ ９ ．５ １４ ．３ １ ． ４３ 一１不考虑混凝土压碎的情况 ） （ ０ ０
泊松 比
０３ ．
构、 组合结构得 到较迅 速的发 展。随着 商业 和住 宅 的要 求提 高 ，
对建筑开间要求越来越 大。传统 的做 法是加 大梁高 或采用 预应 力梁 ， 也有采 用钢结 构的 。加大梁 高对空 间利用有影 响 ， 用预 采 应力 梁施工麻 烦 ， 用钢结 构造价较 高。因此 ， 采 采用传 统的钢筋 混凝土梁结构能否做成跨度大于 １ ５ｍ或更大 ， 而且梁高 又不大 ， 或根据建筑功能的要求 ， 在传 统梁 高范 围内布置 管线 （ 包括 空调
ANSYS,ADINA在钢筋混凝土结构非线性分析中的应用与算例分析
表1跨中截面的应变和应力计算 结果对比
结论与展望
结论与展望
通过以上分析和算例对比，我们可以得出以下结论：ANSYS和ADINA在钢筋混 凝土结构非线性分析中均具有较高的准确性和可靠性，二者的计算结果非常接近。 这些软件提供了丰富的建模功能、强大的求解器和灵活的材料本构模型选择，使 得它们能够广泛应用于各种复杂的钢筋混凝土结构非线性分析问题中。
ANSYS和ADINA在钢筋混凝土结构非线性分析中的应用步骤和常见问题
钢筋混凝土结构非线性分析基本原理和方法
ANSYS和ADINA都是广泛用于工程模拟的有限元软件，它们在钢筋混凝土结构 非线性分析中具有相似的应用步骤。具体而言，可以分为以下几个步骤：
钢筋混凝土结构非线性分析基本原理和方法
1、建立模型：根据实际结构，建立相应的计算模型，包括几何形状、材料属 性、边界条件等。
谢谢观看
ANSYS,ADINA在钢筋混凝土结构 非线性分析中的应用与算例分
析
目录
01 引言
03 表1跨中截面的应变 和应力计算结果对比
02 算例分析 04 结论与展望 Nhomakorabea 内容摘要
ANSYS和ADINA在钢筋混凝土结构非线性分析中的应用与算例分析
引言
引言
钢筋混凝土结构非线性分析在工程实践中具有重要意义，它能够帮助工程师 们更准确地预测结构的性能和行为。在众多非线性分析软件中，ANSYS和ADINA是 最常用的两个工具。本次演示将介绍这两个软件在钢筋混凝土结构非线性分析中 的应用原理、方法和算例，并对其优势和不足进行分析，以期为工程师们提供有 益的参考。
ANSYS线性与非线性分析比较
例：某漏斗状金属板长度为0.085m，两端宽度为0.022m，漏斗半径为0.008m，漏斗处宽度为0.008m，材料为Von Mises多线性各同性，弹性模型E=9.71E10Pa，泊松比V=0.3。
材料应力应变数据如表所示。
分析金属板在轴向0.25mm位移下的应力。
表材料应力应变数据应变（mm/mm）应力（MPa）0.000402 3.90E+070.002118 1.96E+080.003137 2.60E+080.004109 3.03E+080.004631 3.18E+080.005293 3.34E+080.006131 3.48E+080.006926 3.58E+080.007748 3.66E+080.008534 3.72E+080.009548 3.78E+080.010442 3.82E+080.011548 3.86E+080.012236 3.88E+080.014749 3.95E+080.016139 3.97E+080.018424 4.02E+080.020001 4.04E+08线性分析FINI/CLE/FILNAME,EX1-1-1/TITLE,EXAMPLE 1-1-1/UNITS,SI !使用国际单位制*AFUN,DEG !指定角度单位制为度!参数设定*SET,LL,0.085 !长*SET,BL,0.022 !宽*SET,RF,0.008 !漏斗半径*SET,AF,0.008 !漏斗处宽*SET,EXX,9.71E10 !弹性模量*SET,vv,0.3 !弹性模量!进入前处理器/PREP7ET,1,PLANE182 !定义单元类型为182号单元MP,EX,1,EXX !定义弹性模量MP,NUXY,1,vv !定义泊松比!建立四分之一几何模型RECTNG,0,RF+AF/2,0,BL/2 !建立矩形CYL4,,RF+AF/2,RF !建立圆面ASBA,1,2 !面相减得到漏斗RECTNG,RF+AF/2,LL/2,0,BL/2 !建立矩形平板AGLUE,1,3 !面粘接!划分网格MSHKEY,1 !打开映射网格开关MSHAPE,0 !使用四边形或者六面体网格/PNUM,LINE,1 !显示线的编号/PSYMB,LDIR,1 !显示线的方向LPLOT !绘制线LESIZE,9,,,6 !设置9号线划分份数为6 LESIZE,8,,,6LESIZE,11,,,12LESIZE,6,,,6LESIZE,7,,,6,3AMAP,2,8,1,7,9 !通过8,1,7,9关键点映射2号面LESIZE,3,,,12LESIZE,5,,,12LESIZE,4,,,6AMESH,1 !划分1号面!通过对称生成整体模型NSYM,X,175,ALL !通过X对称生成节点，节点增量175 ESYM,,175,ALL !生成单元EPLOTNSYM,Y,350,ALLESYM,,350,ALLEPLOTNUMMRG,ALL !合并所有元素NUMCMP,ALL !压缩编号FINI!进入求解器/SOLUNSEL,S,LOC,X,-LL/2 !选择左边界节点D,ALL,UX,0 !施加位移约束NSEL,R,LOC,Y,0 !选择左边界中轴线上节点D,ALL,ALL !全约束NSEL,S,LOC,X,LL/2D,ALL,UX,0.00025 !施加位移载荷ALLSELSOLVEFINI!进入后处理器/POST1 !进入后处理器非线性分析FINI/CLE/FILNAME,EX1-1-1/TITLE,EXAMPLE 1-1-1/UNITS,SI !使用国际单位制*AFUN,DEG !指定角度单位制为度!参数设定*SET,LL,0.085 !长*SET,BL,0.022 !宽*SET,RF,0.008 !漏斗半径*SET,AF,0.008 !漏斗处宽*SET,EXX,9.71E10 !弹性模量*SET,vv,0.3 !弹性模量!进入前处理器/PREP7ET,1,PLANE182 !定义单元类型为182号单元MP,EX,1,EXX !定义弹性模量MP,NUXY,1,vv !定义泊松比KEYOPT,1,3,0 !定义平面应力问题!定义应力应变数据点TB,MISO,1,1,25 !定义本构关系MISO TBTEMP,0 !定义温度TBPT,,0.000402,3.90E7 !在应力应变曲线上定义点TBPT,,0.002118,1.96E8TBPT,,0.003137,2.60E8TBPT,,0.004109,3.03E8TBPT,,0.004631,3.18E8TBPT,,0.005293,3.34E8TBPT,,0.006131,3.48E8TBPT,,0.006926,3.58E8TBPT,,0.007748,3.66E8TBPT,,0.008534,3.72E8TBPT,,0.009548,3.78E8TBPT,,0.010442,3.82E8TBPT,,0.011548,3.86E8TBPT,,0.012236,3.88E8TBPT,,0.014749,3.95E8TBPT,,0.016139,3.97E8TBPT,,0.018424,4.02E8TBPT,,0.020001,4.04E8TBPLOT,MISO,1 !图形显示应力应变曲线!建立四分之一几何模型RECTNG,0,RF+AF/2,0,BL/2 !建立矩形CYL4,,RF+AF/2,RF !建立圆面ASBA,1,2 !面相减得到漏斗RECTNG,RF+AF/2,LL/2,0,BL/2 !建立矩形平板AGLUE,1,3 !面粘接!划分网格MSHKEY,1 !打开映射网格开关MSHAPE,0 !使用四边形或者六面体网格/PNUM,LINE,1 !显示线的编号/PSYMB,LDIR,1 !显示线的方向LPLOT !绘制线LESIZE,9,,,6 !设置9号线划分份数为6LESIZE,8,,,6LESIZE,11,,,12LESIZE,6,,,6LESIZE,7,,,6,3AMAP,2,8,1,7,9 !通过8,1,7,9关键点映射2号面LESIZE,3,,,12LESIZE,5,,,12LESIZE,4,,,6AMESH,1 !划分1号面!通过对称生成整体模型NSYM,X,175,ALL !通过X对称生成节点，节点增量175 ESYM,,175,ALL !生成单元EPLOTNSYM,Y,350,ALLESYM,,350,ALLEPLOTNUMMRG,ALL !合并所有元素NUMCMP,ALL !压缩编号FINI!进入求解器/SOLUNSEL,S,LOC,X,-LL/2 !选择左边界节点D,ALL,UX,0 !施加位移约束NSEL,R,LOC,Y,0 !选择左边界中轴线上节点D,ALL,ALL !全约束OUTRES,ALL,ALL !设定输出所有的结果TIME,1 !指定求解时间NSUBST,10 !指定载荷子步NLGEOM,ON !打开大变形选项NSEL,S,LOC,X,LL/2D,ALL,UX,0.00025 !施加位移载荷ALLSELSOLVEFINI!进入后处理器/POST1 !进入后处理器。
Ansys第21例非线性屈曲分析实例
Ansys第21例⾮线性屈曲分析实例第21例⾮线性屈曲分析实例—悬臂梁本例通过计算悬臂梁的临界载荷，介绍了利⽤ANSYS进⾏⾮线性屈曲分析的⽅法、步骤和过程。
21.1⾮线性屈曲分析过程1．建⽴模型⾮线性屈曲分析的建模过程与其他分析相似，包括选择单元类型、定义单元实常数、定义材料特性、定义横截⾯、建⽴⼏何模型和划分⽹格等。
2．求解(1)进⼊求解器。
(2)指定分析类型。
⾮线性屈曲分析属于⾮线性静⼒学分析。
(3)定义分析选项。
激活⼤变形效应。
(4)施加初始⼏何缺陷或初始扰动。
可以先进⾏线性屈曲分析，将分析所得到的屈曲模态形状乘以⼀个较⼩的系数后作为初始扰动施加到结构上，本例即采⽤该⽅法。
(5)施加载荷。
所施加的载荷应⽐预测值⾼10%⼀21%。
(6)定义载荷步选项。
(7)设置弧长法。
(8)求解。
3.查看结果在POST26时间历程后处理器中，建⽴载荷和位移关系曲线，从⽽确定结构的临界载荷。
21.2问题描述及解析解图21-1 (a)所⽰为⼀悬臂梁，图21-1 (b)为梁的横横截⾯形状，分析其在集中⼒P作⽤下的临界载荷。
已知截⾯各尺⼨分别为H=50mm、h=43mm、B=35mm、b=32mm，梁的长度L=1m。
钢的弹性模量E=2xl011N/m2，泊松⽐p=0.3。
图21-1⼯⼦悬臂梁21.3分析步骤21.3.1改变任务名拾取菜单Utility Menu→File→Change Jobname，弹出如图21-2所⽰的对话框，在“[/FJLNAM]”⽂本框中输⼊EXAMPLE21，单击“OK”按钮。
21.3.2选择单元类型拾取菜单Main Menu→Preprocessor→Element Type→Add/Edit/Delete，弹出如图21-3所⽰的对话框，单击“Add.”按钮，弹出如图21-4所⽰的对话框，在左侧列表中选“Structural Beam”，在右侧列表中选“3 node 189”，单击“OK”按钮，返回到如图21-3所⽰的对话框，单击“Close”按钮。
《有限元教程》20例ANSYS经典实例
《有限元教程》20例ANSYS经典实例有限元方法在工程领域中有着广泛的应用，能够对各种结构进行高效精确的分析和设计。
其中，ANSYS作为一种强大的有限元分析软件，被广泛应用于各个工程领域。
下面将介绍《有限元教程》中的20个ANSYS经典实例。
1.悬臂梁的静力分析：通过加载和边界条件，研究悬臂梁的变形和应力分布。
2.弯曲梁的非线性分析：通过加载和边界条件，研究受弯曲梁的非线性变形和破坏。
3.柱体的压缩分析：研究柱体在压缩载荷作用下的变形和应力分布。
4.钢筋混凝土梁的受弯分析：通过添加混凝土和钢筋材料属性，研究梁的受弯变形和应力分布。
5.圆盘的热传导分析：根据热传导方程，研究圆盘内部的温度分布。
6.输电线杆的静力分析：研究输电线杆在风载荷和重力作用下的变形和应力分布。
7.轮胎的动力学分析：通过加载和边界条件，研究轮胎在不同路面条件下的变形和应力分布。
8.支架的模态分析：通过模态分析，研究支架的固有频率和振型。
9.汽车车身的碰撞分析：通过加载和边界条件，研究汽车车身在碰撞中的变形和应力分布。
10.飞机翼的气动分析：根据飞机翼的气动特性，研究翼面上的气压分布和升力。
11.汽车车身的优化设计：通过参数化建模和优化算法，寻找最佳的车身结构设计。
12.轮毂的疲劳分析：根据材料疲劳寿命曲线，研究轮毂在不同载荷下的寿命。
13.薄膜材料的热应力分析：根据热应力理论，研究薄膜材料在不同温度下的应变和应力。
14.壳体结构的模态分析：通过模态分析，研究壳体结构的固有频率和振型。
15.地基基础的承载力分析：通过加载和边界条件，研究地基基础的变形和应力分布。
16.水坝的稳定性分析：根据水力和结构力学，研究水坝的稳定性和安全性。
17.风机叶片的动态分析：通过加载和边界条件，研究风机叶片在不同风速下的变形和应力分布。
18.圆筒容器的蠕变分析：根据蠕变理论，研究圆筒容器在持续加载下的变形和应力。
19.桥梁结构的振动分析：通过模态分析，研究桥梁结构的固有频率和振型。
非线性2【ANSYS非线性分析】
1几何非线性分析随着位移增长，一个有限单元已移动的坐标可以以多种方式改变结构的刚度。
一般来说这类问题总是是非线性的，需要进行迭代获得一个有效的解。
大应变效应一个结构的总刚度依赖于它的组成部件（单元）的方向和单刚。
当一个单元的结点经历位移后，那个单元对总体结构刚度的贡献可以以两种方式改变变。
首先，如果这个单元的形状改变，它的单元刚度将改变。
（看图2─1(a)）。
其次，如果这个单元的取向改变，它的局部刚度转化到全局部件的变换也将改变。
（看图2─1（b)）。
小的变形和小的应变分析假定位移小到 足够使所得到的刚度改变无足轻重。
这种刚度不变假定意味着使用基于最初几何形状的结构刚度的一次迭代足以计算出小变形分析中的位移。
（什么时候使用“小”变形和应变依赖于特定分析中要求的精度等级。
相反，大应变分析说明由单元的形状和取向改变导致的刚度改变。
因为刚度受位移影响，且反之亦然，所以在大应变分析中需要迭代求解来得到正确的位移。
通过发出NLGEOM ，ON （GUI 路径Main Menu>Solution>Analysis Options)，来激活 大应变效应。
这效应改变单元的形状和取向，且还随单元转动表面载荷。
（集中载荷和惯性载荷保持它们最初的方向。
）在大多数实体单元（包括所有的大应变和超弹性单元），以及部分的壳单元中大应变特性是可用的。
在ANSYS/Linear Plus 程序中大应变效应是不可用的。
图1─11 大应变和大转动大应变处理对一个单元经历的总旋度或应变没有理论限制。
（某些ANSYS单元类型将受2到总应变的实际限制──参看下面。
）然而，应限制应变增量以保持精度。
因此，总载荷应当被分成几个较小的步，这可以〔NSUBST ，DELTIM ，AUTOTS 〕，通过GUI 路径 Main Menu>Solution>Time/Prequent)。
无论何时当系统是非保守系统，来自动实现如在模型中有塑性或摩擦，或者有多个大位移解存在，如具有突然转换现象，使用小的载荷增量具有双重重要性。
2-非线性本构关系【ANSYS非线性分析】
第二章材料本构关系§2.1本构关系的概念本构关系：应力与应变关系或内力与变形关系结构的力学分析，必须满足三类基本方程:（1）力学平衡方程:结构的整体或局部、静力荷载或动力荷载作用下的分析、精确分析或近似分析都必须满足；（2）变形协调方程:根据结构的变形特点、边界条件和计算精度等，可精确地或近似地满足；（3）本构关系:是连接平衡方程和变形协调方程的纽带，具体表达形式有：材料的应力-应变关系，截面的弯矩-曲率关系，轴力-变形（伸长、缩短）关系，扭矩-转角关系，等等。
所有结构（不同材料、不同结构形式和体系）的力学平衡方程和变形协调方程原则上相同、数学形式相近，但本构关系差别很大。
有弹性、弹塑性、与时间相关的粘弹性、粘塑性，与温度相关的热弹性、热塑性，考虑材料损伤的本构关系，考虑环境对材料耐久性影响的本构关系，等等。
正确、合理的本构关系是可靠的分析结果的必要条件。
混凝土结构非线性分析的复杂性在于：钢筋混凝土---复杂的本构关系：有限元法---结构非线性分析的工具：非线性全过程分析---解决目前结构分析与结构设计理论矛盾的途径：§2.2 一般材料本构关系分类1．线弹性(a) 线性本构关系； (b) 非线性弹性本构关系图2-1 线弹性与非线性弹性本构关系比较在加载、卸载中，应力与应变呈线性关系：}]{[}{εσD = （图2-1a ） 适用于混凝土开裂前的应力-应变关系。
2． 非线性弹性在加载、卸载中，应力与应变呈非线性弹性关系。
即应力与应变有一一对应关系，卸载沿加载路径返回，没有残余变形（图2-1b ）。
}{)]([}{εεσD = 或 }{)]([}{εσσD =适用于单调加载情况结构力学性能的模拟分析。
3． 弹塑性图2 – 2 弹塑性本构关系(a)典型弹塑性；(b)理想弹塑性；(c)线性强化；(d)刚塑性典型的钢筋拉伸应力、应变曲线 （图2-2（a ））包含弹性阶段（OA ）、流动阶段（AB ）及硬化阶段（BC ）。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
/BATCH /COM,ANSYS RELEASE 12.1 UP20091102 /input,menust,tmp,'',,,,,,,,,,,,,,,,1 /GRA,POWER /GST,ON /PLO,INFO,3 /GRO,CURL,ON /CPLANE,1 /REPLOT,RESIZE WPSTYLE,,,,,,,,0 *SET,B,0.005 *SET,H,0.005 *SET,H , *SET,H,0.005 *SET,L,1 *SET,A,B*H *SET,I,B*H**3/12 /PREP7 !* ET,1,BEAM23 !* !* R,1,A,I,H, !* !*
4.9 选择Main Menu>Solution>Load Step Opts>Write LS File,在对话框中输入 3，单击OK,如下图所 示
4.10 选择Main Menu>Solution>Solve >From LS Files，弹 出如下对话框，并进 行如下设置单击ok
5.1 选择Main Menu>General Postproc>Read Results>First Set,读取第一个载荷步的求解结果 5.2 Main Menu>General Postproc>plot Results>Contour Plot>Nodal Solu,弹出如下对 话框，并进行如下页的设置
5.11单击ok,此时ANSYS窗口将显示载荷为14N的杆件变形形状， 结果如下的图所示
5.12 选择Main Menu>General Postproc>List Results,弹出对话框进行 如下左的设置，点OK，显示载荷为14N的杆件节点位移，如下右图. 5.13 点击右上角的“X”关闭按钮，弹出对话框，选择save everything， 然后点击ok。则退出了ansys,并保存了该分析的所有文件。
ห้องสมุดไป่ตู้
2.4 选择Main Menu>Preprocessor> Material Props>Material Models, 弹出对话框如左上图，在 Material Prop Available 栏中依次展开 Structural>linear>isotr opic,弹出如左下所示对 话框，EX输入8e10， prxy输入0.3单击OK关闭 该对话框，选择Material >Exit关闭左上对话框。
4.1 选择Main Menu>Solution>Analysis Type>New Analysis,在弹出 的对话框中，选择分析类型 Static,单击OK，如下左图所示
4.2 选择Main Menu>solution>Analys is Type>Sol’n Control 命令，弹出对话框，参照 下右图进行设置
4.6 选择Main Menu>Solution>Define Loads>Apply>Structural>Force/Moment>On Nodes,弹出如左拾取菜单,在其输入 10，单击ok，弹出对话框，对其进行如下右上图的设置，单击apply,在弹出如左 拾取菜单，输入10，单击ok对其进行如下右下图的设置，单击OK。
/GO D,P51X, , , , , ,ALL, , , , , FLST,2,1,1,ORDE,1 FITEM,2,10 !* /GO F,P51X,FY,-12 FLST,2,1,1,ORDE,1 FITEM,2,10 !* /GO F,P51X,FX,0.05 LSWRITE,1, FLST,2,1,1,ORDE,1 FITEM,2,10 !* /GO F,P51X,FY,-13 FLST,2,1,1,ORDE,1 FITEM,2,10 !* /GO F,P51X,FX,0.05 LSWRITE,2, FLST,2,1,1,ORDE,1 FITEM,2,10 !*
4.7 选择Main Menu>Solution>Load Step Opts>Write LS File,在对话框中输入2，单击OK,如下图所示
4.8 选择Main Menu>Solution>Define Loads>Apply>Structural>Force/Moment>On Nodes,弹出如左拾取菜单,在其输入 10，单击ok，弹出对话框，对其进行如下右上图的设置，单击apply,在弹出如左 拾取菜单，输入10，单击OK对其进行如下右下图的设置，并单击OK。
1.1 单击 开始>所有 程序 >ansys12.1>mench anical APDL Product Launcher窗口如左图。 在working Directory 栏中设置工作目录， 在job name中输入如 左图所示的工作名。 单击Run,进入ansys。
1.2 选择 Main Menu>file>Change Title,弹 出如下图所示的对话框在其输入工作标题为 slender rod。
2.2 选择Main Menu>Preprocessor>Element Type>add/Edit/Delete弹出对话框，单击add， 选择StructualBeam中的2D plastic23单元，如图，单击ok关闭对话框，单击close关闭Element Type对话框
2.3 选择Main Menu>preprocessor > real constants>add/edit /Delete，在弹出如左 上图的对话框单击OK, 并弹出如左下图所示的 对话框并对其进行如左 所示的参数设置，单击 OK。
2.1 选择Main Menu>Parameters>Scalar Parameters,弹出对话框，在selection中输入 B=0.005,单击accept;输入H=0.005,单击accept; 输入L=1,单击accept;输入A=B*H,单击accept;输 入I=B*H**3/12,单击accept,如下图
PRNSOL,U,COMP SET,FIRST SET,FIRST !* /EFACET,1 PLNSOL, U,SUM, 1,1.0 FINISH ! /EXIT,ALL
/GO F,P51X,FY,-14 FLST,2,1,1,ORDE,1 FITEM,2,10 !* /GO F,P51X,FX,0.05 LSWRITE,3, LSSOLVE,1,3,1, FINISH /POST1 SET,FIRST !* /EFACET,1 PLNSOL, U,SUM, 1,1.0 !* PRNSOL,U,COMP SET,NEXT !* /EFACET,1 PLNSOL, U,SUM, 1,1.0 !* PRNSOL,U,COMP SET,LAST !* /EFACET,1 PLNSOL, U,SUM, 1,1.0 !*
5.9 选择Main Menu>General Postproc>Read Results>Last Set,读取第 3个载荷步(14N)的求解结果
5.10 Main Menu>General Postproc>plot Results>Contour Plot>Nodal Solu,弹出如下对话框，并进行如下的设置
5.3 单击ok,此时ANSYS窗口将显示载荷为12N的杆件变形 形状，结果如下页的图所示
（载荷1）位移云图结果显示图
5.8 选择Main Menu>General Postproc>List Results,弹出对 话框进行如下左的设置，点OK，显示载荷为13N的杆件节点位移， 如下右图
18:49:38
06/16/2014
MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,8E10 MPDATA,PRXY,1,,0.3 N,1,,,,,,, N,10,0,L,0,,,, !* FILL,1,10,8, , ,1,1,1, !* FLST,2,2,1 FITEM,2,1 FITEM,2,2 E,P51X FLST,4,1,2,ORDE,1 FITEM,4,1 EGEN,9,1,P51X, , , , , , , , , , , FINISH /SOL !* ANTYPE,0 ANTYPE,0 NLGEOM,1 NSUBST,1,0,0 AUTOTS,1 FLST,2,1,1,ORDE,1 FITEM,2,1 !*
3.1选择Main Menu>Preprocessor> modeling>Create>Nod es>In Active CS,弹出对 话框，在NODE node number中输入1，单击 apply; 在NODE node number中输入10，输入 如下参数，单击ok.
3.2选择 Uility Menu>PlotCtrls>Numberin g,弹出对话框，使NODE node number复选框，使其off变为 on,如下所示,单击ok.