Workbench DS非线性分析

合集下载

Workbench12.0力学非线性_Diag

Workbench12.0力学非线性_Diag

ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.
7-4
April 30, 2009 Inventory #002659
Workbench Mechanical - Diagnostics
... Solution Information
Training Manual
Recall from Solver Section 2 that the Contact Region name and the contact set ID are listed. From that example, we know that the above contact set 9 is part of “Teeth 3” region.
Training Manual
The Solver Output can provide detailed text output about the solution. It is useful to become familiar with how to read this file.
1) The beginning of the Solver Output simply shows the ANSYS license used (in this case, ANSYS Multiphysics) and the version number.
A. Solution Information
Training Manual
In Chapter 2, the Solution Information branch was introduced
– Recall that with the Solution Information branch, the detailed Solver Output from ANSYS can be reviewed, and convergence graphs, such as the Force Convergence behavior, can be plotted.

ANSYS第三讲_Standard 中的非线性分析

ANSYS第三讲_Standard 中的非线性分析
– 把新的残差R2与容差进行比较, 察看在u2 处是否得到收敛解 。
L3.16
ABAQUS/analysis_ ABAQUS/Standard中的非线性分析
求解平衡方程
– 该过程将一直重复,直到力的残差在允许的容差之内。每次迭代i需要: 1. 形成切线刚度Ki。 2. 求解系统方程组,得到位移修正ci+1 。 • 修正位移的估计值: ui+1 = ui + ci+1。 3. 基于ui+1 计算内力向量Ii+1。 4. 进行平衡收敛判断: • 是否Ri+1 在容差之内?
非线性输入文件
• 分析步和过程输入
*STEP, NLGEOM, INC=25
– NLGEOM: 包括所有由以下原因引起的几何非线性效应: • 大挠度、大旋转、大变形。 • 预载荷(初始应力)。
• 载荷刚度。
– 如果上面列出的项不重要,应用NLGEOM选项得到的结果同没有应用 NLGEOM选项得到的结果类似,但是求解的费用更高。 – INC=25: 在本例中允许的最大增量为25: • 如果在施加全部载荷之前达到了最大增量数量,程序将会中止。 • 保证程序的运行时间不会太长—用户可以重新启动分析。 • 默认值为100。
(Eq. 3.3)
L3.10
ABAQUS/analysis_ ABAQUS/Standard中的非线性分析
求解平衡方程
– 略去高阶项,方程可以写为
Ki ci1 P(ui ) I (ui ),
其中 K i
I (ui ) P(ui ) 为切线刚度。 u u
– 解的下一次近似为
*MODAL DYNAMIC (只允许固定时间增量)
*COUPLED THERMAL-ELECTRIC – 物理(真实)时间或名义时间的度量取决于过程和时间相关过程或率相关行为。

【ANSYS非线性分析】4-非线性分析方法

【ANSYS非线性分析】4-非线性分析方法

01112121222y y d N d d R d M d d R ελφ⎧⎫
⎧⎫⎡⎤⎧⎫=∆+⎨⎬⎨⎬⎨⎬⎢⎥
⎣⎦⎩⎭
⎩⎭⎩⎭ 改写为,
11112021222y y d N R d d d d M R d d εφλ-⎡⎤⎧⎫⎧⎫
⎧⎫=-⎨⎬⎨⎬⎨⎬⎢⎥
-⎩⎭⎩⎭⎩⎭
⎣⎦ 求解过程中,可控制d φy 的值,求出相应的0d ε及荷载增量比例因子d λ。

由于ij d 与截面应变平面有关,需要迭代才能使截面补平衡力12,R R 趋近于零。

图4-9 位移控制法 在结构分析中控制指定位移增量,则P —δ曲线的下降段不难求得。

将底端固定顶端自由的柱,在柱顶端施加水平荷载,将柱的加载点处换为支座,而分析时控制该支座位移并求出该支座的反力,图4—9表示了得到的全过程分析P-δ曲线。

对于一般结构,将刚度矩阵重新排列,使得选择的控制位移排到最后,将原矩阵分块表示成以下形式,
111211121
22222K K du P R K K du P R ⎡⎤⎧⎫⎧⎫⎧⎫=∆+⎨⎬⎨⎬⎨⎬⎢⎥⎣⎦⎩⎭⎩⎭⎩⎭
λ 改写方程为,
11
11121221
2222K P R K du du K P R K -⎡⎤⎧⎫⎧⎫⎧⎫
=-⎨⎬⎨⎬⎨⎬⎢⎥-∆⎩⎭⎣⎦⎩⎭⎩⎭
λ 需要指出的是,改写以后的系数矩阵是不对称的,也不是带状的,求解时需要较多的存储单元。

§4.5.4 修正完善后的弧长法 1.弧长法的基本原理
仍从结构增量平衡方程:{}{}{}11i i i i K w P g --=-∆∆λ∆。

04 ANSYS13.0 Workbench 结构非线性培训 一般过程解析

04 ANSYS13.0 Workbench 结构非线性培训 一般过程解析

... 建立非线性模型
• 对任何结构单元, DOF(自由度)求解Du 是对节点求解 • 应力和应变是在积分点计算. 由DOF推导而来.
– 例如, 可由位移确定应变 ,经:
Training Manual
Dε BDu
s, e u
– 这里 B 称为 应变-位移矩阵
• 右图所示的一 4节点四边形单元有 2x2个积分点, 红点为积分点. • 在后处理结果中, 积分点的应力/应变值经外插值或复制到节点位置。
Training Manual
– 仅适用于高阶单元. – 当一部件厚度方向只有一个单元时,强制使用完全积分有助于提高精确度.
2-6
Workbench Mechanical - General Nonlinear Procedures
... 建立非线性模型
• WB Mechanical 默认采用高阶单元(有中节点)来划分网格.
第二章
一般过程
Workbench – Mechanical 结构非线性
2-1
Workbench Mechanical - General Nonlinear Procedures
章节概述

Training Manual
这章介绍一般工具和程序,不是对特殊来源非线性的详细介绍, 但介绍了 达到收敛的有用措施和后处理结果:
A. B. C. D.
建立非线性模型 分析设置 非线性结果后处理 作业
2-2
Workbench Mechanical - General Nonlinear Procedures
A. 建立非线性模型
什么是建立非线性模型与线性模型的不同?
Training Manual
• 某些情况,它们没有不同!

Solidworks的线性和非线性分析方法与技巧

Solidworks的线性和非线性分析方法与技巧

Solidworks的线性和非线性分析方法与技巧Solidworks是一款广泛使用的三维计算机辅助设计(CAD)软件,它提供了丰富的分析工具,使工程师能够进行线性和非线性分析。

本文将介绍Solidworks中常用的线性和非线性分析方法与技巧。

一、线性分析方法与技巧1. 静态分析:静态分析用于研究物体在外力作用下的静止行为。

在Solidworks 中进行静态分析时,需要定义边界条件、材料特性和加载条件。

通过分析结果,可以获得物体的变形、内应力和应变等信息。

2. 模态分析:模态分析用于研究结构的固有频率和振型。

它对于预测结构的共振问题和自由振动问题非常有用。

在Solidworks中进行模态分析时,可以选择求解结构的前几个固有频率和相应的振型。

3. 热传导分析:热传导分析用于研究热量在物体内部的传导过程。

在Solidworks中进行热传导分析时,需要定义边界条件、材料的热传导性质和温度加载条件。

通过分析结果,可以获得物体的温度分布和热传导通量等信息。

4. 疲劳分析:疲劳分析用于研究物体在交变载荷下的寿命和破坏位置。

在Solidworks中进行疲劳分析时,需要定义材料的疲劳特性、加载条件和疲劳强化因素。

通过分析结果,可以获得物体的寿命预测和疲劳破坏位置等信息。

在进行线性分析时,以下是一些Solidworks中常用的技巧:- 合理使用边界条件:在定义边界条件时,需要根据实际情况选择合适的固定支撑、约束和加载类型。

合理的边界条件能够提高分析结果的准确性。

- 网格划分:在进行线性分析之前,需要对物体进行网格划分。

合理的网格划分能够提高分析的精度和计算效率。

较小的网格将更准确地捕捉结构中的应变和应力变化。

- 结果后处理:Solidworks提供了丰富的结果后处理工具,如变形云图、应力云图和应变云图等。

通过仔细观察分析结果,可以发现潜在的问题并进行进一步的优化设计。

二、非线性分析方法与技巧1. 大变形分析:大变形分析用于研究物体在承载过程中的形状变化。

AnsysWorkbench工程应用之——结构非线性(上):几何非线性(1)

AnsysWorkbench工程应用之——结构非线性(上):几何非线性(1)

AnsysWorkbench工程应用之——结构非线性(上):几何非线性(1)在上一篇文章中,我们已经对非线性有了初步的了解,本文将详细介绍几何非线性,以及非线性的计算原理。

本文较长,将分为上下两篇。

1 大变形选项在线性计算中,结构的刚度矩阵是不变的,而且不应用非线性收敛准则。

大变形Large Deflection只是一种算法,它考虑了结构变形后的刚度矩阵重建。

理论上来说,开启大变形后计算精度更高,但是将消耗更多的计算资源和时间。

所以在小变形、小转动等问题中,无需开启大变形选项。

假设构件转角为α,小转动时,构成结构刚度矩阵的三角函数因子cosα≈1,所以无需开启大变形选项。

而当大变形、大转动时,cosα与1相差较大,此时重建刚度矩阵是必要的。

如果出现了大变形、大转动,程序计算后会跳出警告信息,提示用户开启大变形选项。

一般工程经验上,对于普通精度要求的问题,变形超过5%或转角超过10°时,建议打开大变形开关。

2 非线性控制2.1收敛求解的原理在非线性求解中,真实的载荷-位移曲线是未知的,不能直接使用一组线性方程得到结果,而需要使用一系列的线性方程迭代,逼近非线性解。

有限元计算中使用的迭代过程为Newton-Raphson方法,简称牛顿法,达到收敛的迭代称为平衡迭代。

牛顿法的原理如下:在牛顿法中,第一次迭代,施加总载荷Fa。

得到位移结果x1。

根据位移,算出内力F1 。

如果Fa≠F1,系统不平衡。

因此,根据当前的条件,计算新的刚度矩阵(虚线的斜率)。

Fa-F1的值称为不平衡力或残余力。

残余力达到足够小时,求解收敛。

重复以上过程,直到Fa=Fi。

在这个例子之中,四次迭代之后,系统达到平衡,求解收敛。

实际计算中,残差Fa-Fi不可能正好等于0,所以规定只要残差小于一个规定的微小量,就认为计算已经收敛了,这个微小量就是力收敛准则,此处以[R]表示,即Fa-Fi<[R],则达到收敛。

每一次迭代中,当残差小于[R]表现为收敛,大于[R]表现为发散。

ANSYSWorkbench基础教程与工程分析详解第七章结构非线性分析

ANSYSWorkbench基础教程与工程分析详解第七章结构非线性分析

前面的内容属于线性问题,其符合虎克定律(Hooke),满足公式:F=kx。

其中,k表示刚度矩阵常量,力与位移呈线性关系。

实际工程中多数结构的力与位移是呈非线性关系的,出现非线性行为,即载荷能够引起结构刚度的显著改变。

引起结构刚度变化的原因有:应变超出弹性极限,即产生塑性变形;大挠度,如钓鱼竿受力变形的过程;接触,物体之间的接触变形。

本章所要学习的内容包括:¾了解结构非线性基础¾熟悉ANSYS Workbench软件大变形分析的步骤¾了解结构非线性分析的应用场合¾理解非线性分析的计算结果¾了解非线性分析与其他分析的不同之处7.1 结构非线性分析基础7.1.1 引起非线性的原因结构在承受大变形时,几何形状发生变化会导致结构的非线性变化,如悬臂杆一端受力使杆发生弯曲,力臂明显减少,从而使得杆端的刚度不断增大,这是大挠度引起的非线性响应。

此外,钓鱼竿也是常见的几何非线性,如图7-1所示。

几何非线性主要有大应变、大挠度、应力刚化引起的非线性响应。

非线性应力-应变关系是典型的材料非线性。

影响材料应力-应变关系的因素有加载历史、环境问题、加载的时间总量等。

材料非线性如图7-2所示。

图7-1 钓鱼竿大变形图7-2 材料非线性接触是一种很普遍的非线性行为,是状态变化非线性类型中一个特殊且很重要的部分。

当两个接触物体相互接触或者分离时会发生刚度的突然变化,此时也会出现非线性。

在非线性静力分析中,刚度矩阵[K ]依赖于位移矩阵[x ]:[k(x)](x)={F}. 式中,力与位移的关系是非线性的,同样可参考图7-2。

Contact (接触类型) Iterations (迭代次数) Normal Behavior (法向分离) Tangential Behavior (切向滑移) Bonded (绑定) 1 Closed (无间隙) Closed (不能滑移) No Separation (不分离) 1 Closed (无间隙) Open (允许滑移) Frictionless (光滑) Multiple (多次) Open(允许有间隙) Open (允许滑移) Rough (粗糙) Multiple (多次) Open(允许有间隙) Closed (不能滑移) Frictional (摩擦)Multiple (多次)Open(允许有间隙)Open (允许滑移)其中,Bonded 和No Separate 两种接触是最基础的线性行为,故仅需要迭代一次,所以计算速度非常快。

AnsysWorkbench非线性分析牛顿辛普森法过程 ppt课件

AnsysWorkbench非线性分析牛顿辛普森法过程 ppt课件

2. 材料非线性
非线性应力应变关系是非线性结构行为的普遍原因.
应力
应力
应变

应变
橡胶
3. 状态变化非线性 导致刚度突然变化的状态改变是非线性行为的另一个
普遍原因.
在此例中, 随着载荷的 增加, 接触状态从“开” 变为 “闭合”, 从而引 起刚度变化.
载荷步、子步、平衡迭代 a. 在静态分析中允许设置多个分析步,并一步一步
的求解 b.对于静态分析,终止时间被用作确定载荷步和载
荷子步的追踪器
终止时间
终止时间
载荷步与子步的关系
子步与平衡迭代的关系
载荷的施加方式
设置多步载荷
Fa :节点载荷 外载荷和约束力 Fnr :节点力 节点受其附近的单元施加的力
12 34
构件处于平衡状态,因
F2
此每个节点也平衡,即
F1 Fa
F4 节点力应等于节点载荷 :
导入装配体,软件自动检测接触面并生成接触对,也可 以手动添加接触对
接触单元就像是在接触区域面上的“皮肤”,提 供了部件之间的相互关系。
在“contact”分支点击某个接触对,
构成这个接触对的部件就会变成半透明
与该接触对无关的部件变成全透明 接触面用红色和蓝色表示 Contact-接触面 Target-目标面
计算膜片弹簧大端加载——位移特性曲线的 A-L公式为 :
F 16(1 E 2 1 )tlL (n R rl)2[h (R L r l1)h (2(R L rl)1)t2]
将膜片弹簧的尺寸代入 A-L公式为得大端载荷-位移曲线为 :
F 12.6 99 1 9 329 1 28 912 138
接触的类型:
Bonded(绑定): 这是AWE中关于接触的默认设置。如果 接触区域被设置为绑定,不允许面或线间有相对滑动或分离 。

ansys workbench ws05_非线性诊断

ansys workbench ws05_非线性诊断

…作业5 - 非线性诊断
• 后处理结果...
ANSYS Workbench – Simulation ANSYS Workbench – Simulation ANSYS Workbench – Simulation ANSYS Workbench – Simulation
Training Manual
Training Manual
ANSYS Workbench – Simulation ANSYS Workbench – Simulation ANSYS Workbench – Simulation ANSYS Workbench – Simulation
• •
ANSYS Workbench – Simulation ANSYS Workbench – Simulation ANSYS Workbench – Simulation ANSYS Workbench – Simulation
… 作业 5- 非线性诊断

Training Manual
为了解决这个问题,在与接触区域相关的曲面上插入网格控制为 0.5mm(总共有6个面). 局部的网格细化在这个关键部位增加了更多 的接触高斯点从而降低了收敛迭代的次数。同时也提高了结果质 量。
ANSYS Workbench – Simulation ANSYS Workbench – Simulation ANSYS Workbench – Simulation ANSYS Workbench – Simulation
Training Manual
作业5 非线性诊断
常见问题解答
作业 5- 非线性诊断
• 目的
Training Manual
ANSYS Workbench – Simulation ANSYS Workbench – Simulation ANSYS Workbench – Simulation ANSYS Workbench – Simulation

AnsysWorkbench工程应用之——结构非线性(中):材料非线性(1)弹塑篇

AnsysWorkbench工程应用之——结构非线性(中):材料非线性(1)弹塑篇

AnsysWorkbench工程应用之——结构非线性(中):材料非线性(1)弹塑篇本文结合材料知识与工程应用,从理论到实践,从书本到实操、从动脑到动手,保姆式手把手介绍非线性材料本构使用方法!这也可能是您在网上能找到的关于Ansys Workbench非线性材料蕞啰嗦(xiangxi)的一篇基础性文章。

材料的应力应变关系一般用材料本构来表示,本构模型又称材料的力学本构方程,或材料的应力-应变模型,是描述材料的力学特性(应力-应变-强度-时间关系)的数学表达式。

Ansys Workbench提供了丰富的非线性材料本构,用户也可基于试验数据定义自己的非线性材料。

材料的响应与载荷或变形施加的速率无关的材料称为率不相关材料,如弹塑性、超弹性(橡胶等)、混凝土等材料,大多数金属在低温(≤30%左右的熔点)和低应变速率时,为率无关材料,通常所说的塑性也就是率无关塑性。

材料的响应与载荷或变形施加的速率相关的材料称为率相关材料,如蠕变、黏弹性材料等。

材料的应力应变曲线也称为材料的响应曲线,是通过材料试验得到的,主要材料试验有单轴试验、等双轴试验、平面剪切试验、体积试验、松散试验等。

最常见的为单轴试验,可以测试拉伸也可以测试压缩,下图为某些塑形材料单轴拉伸试验的工程应力应变曲线。

1 率无关塑形1.1 基本理论1.1.1 比例极限与屈服极限结构的塑性响应基于单轴实验结果获得。

通过单轴应力-应变实验,可以得到材料的比例极限、屈服极限(或弹性极限)、应变强化。

对于塑形材料,当应力小于比例极限时,材料呈现线性; 当应力小于屈服点时,材料呈弹性,载荷卸除后,所有应变可以完全恢复;当应力大于屈服点时,材料呈塑性,载荷卸除后,应变不能完全恢复。

由于比例极限和屈服点非常接近,有限元软件假设两者值相等。

1.1.2 应力应变的工程曲线与真实曲线您一定很好奇,为什么材料力学课本中的塑形材料σ-ε曲线有下降段,而有限元分析软件中设置塑形材料不定义下降段,这不是因为忽然误差,而是因为材料力学课本上使用的是工程σ-ε曲线,也称名义应力-应变曲线,而有限元计算中往往使用切线斜率直线代替真实σ-ε曲线,他们的关系如下图。

day5 Workbench DS模态分析及屈曲分析

day5 Workbench DS模态分析及屈曲分析

ANSYS BASIC TRAINING
模态分析基础
• 对于一个模态分析, 固有圆周频率 ωi 和 振型 φi 都能从矩阵方程式 对于一个模态分析, 里得到: 里得到:
([K ] − ω [M ]){φ } = 0
2 i i
在某些假设条件下的结果与分析相关: 在某些假设条件下的结果与分析相关: – [K] 和 [M] 是常量: 是常量: • 假设为线弹性材料特性 • 使用小挠度理论, 不包含非线性特性 使用小挠度理论, • [C] 不存在, 因此不包含阻尼 不存在, • {F} 不存在, 因此假设结构没有激励 不存在, • 根据物理方程 , 结构可能不受约束(rigid-body modes 结构可能不受约束(rigid (rigid或者部分/ present) ,或者部分/完全的被约束住 • 记住这些在Design Simulation中进行模态分析的假设是非常重要的. 记住这些在Design Simulation中进行模态分析的假设是非常重要的 中进行模态分析的假设是非常重要的.
ANSYS BASIC TRAINING
预应力模态分析
• 在某些情况下,进行模态分析时,需要考虑预应力效果. 在某些情况下,进行模态分析时,需要考虑预应力效果. – 在一个静态载荷(static) 的作用下, 结构的应力状态可能影响到 在一个静态载荷(static) 的作用下, 它的固有频率. 这一点是非常重要的, 它的固有频率. 这一点是非常重要的,尤其是对于那些在某一个 或两个尺度上很薄的结构. 或两个尺度上很薄的结构. – 现在,考虑一个吉它弦被调节的情况 – 当轴向载荷增加(拉紧)的 现在, 当轴向载荷增加(拉紧) 时候, 横向频率也随之相应的增加. 这是一个应力硬化的例子. 时候, 横向频率也随之相应的增加. 这是一个应力硬化的例子.

ANSYSWorkbench模型几何缺陷的引入与非线性屈曲分析

ANSYSWorkbench模型几何缺陷的引入与非线性屈曲分析

ANSYSWorkbench模型几何缺陷的引入与非线性屈曲分析在一些特定的非线性分析场合,如:非线性屈曲分析,需要为模型引入一个几何缺陷。

在较新的Workbench版本中已经可以实现这一功能,其实现过程分为两个阶段:第一阶段是进行特征值分析,得到变形模式;第二个阶段是引入小扰动作为结构的几何缺陷。

具体操作方法为:新建一个特征值分析的流程(系统A和B),再新建一个运行后续非线性屈曲分析的静力分析系统C,如下图所示。

共享材料Data数据B2:Engineering Data到C2:Engineering Data,传递变形结果数据B6:Solution到C4:Model,注意到此时C4:Model自动变为C3:Model,如下图所示。

选择单元格B6:Solution,通过Workbench菜单View>Properties,在Properties of Schematic B6:Solution中设置Scale Factor和引入缺陷的屈曲模态号Mode,如下图所示。

下面对一个典型计算实例进行简要的介绍。

悬臂的方形钢板,如下图所示。

钢板尺寸为1000mm×1000mm×10mm,底部(A端)为固定约束,顶部悬臂端(B端)作用有面内的荷载。

如果悬臂钢板所采用钢材弹性模量2.0×105MPa,屈服强度为300MPa,屈服后切线模量为2.0×103MPa。

计算此悬臂钢板的极限受压承载能力。

下图为材料非线性应力-应变曲线。

在Mechanical组件中进行特征值分析,得到一阶屈曲荷载约为4.35kN,引入基于第一阶屈曲变形最大侧向变形为10mm的几何缺陷,进行非线性屈曲分析。

经过计算后得到板的X方向(加载方向)反力与板自由端中间节点面外位移之间关系如下图所示。

由上述的载荷-变形曲线可知,悬臂板的极限承载力为曲线顶点,其数值约为41415N,低于特征值屈曲分析的临界承载力。

Workbench DS非线性分析

Workbench DS非线性分析

[( ) ( ) ( )] ( ) σo =
1 2
σx −σ y 2 + σ y −σz 2 + σz −σx
2
+
6
τ2 xy

2 yz
+
τ
2 xz
ANSYS BASIC TRAINING
Mises 屈服准则
• 在主应力空间,von Mises 屈服面是圆柱面.
• 材料在屈服面内呈弹性弹性行为。注意,多轴应力状态可以存在于圆柱 面内的任何一点。在圆柱面边上会产生屈服。没有应力状态会存在于圆 柱面以外。硬化规则会描述圆柱面随响应的变化。
– 在绑定的接触中,纯粹的罚函数法可以想象为在接 触面间施加了十分大的刚度系数来阻止相对滑动。 这个结果是在接触面间的相对滑动可以忽略的情况 下得到的。
– MPC 方程当中对接触面间的相对运动定义了约束 方程,因此没有相互的滑动。这个方程经常作为罚 函数法的最好的替代。
ANSYS BASIC TRAINING
In this case, the gap between the two parts is bigger than the pinball region, so no automatic gap closure will be performed.
ANSYS BASIC TRAINING
载荷方向
Fb2
Fb1
Fa
Fa1
ANSYS BASIC TRAINING
xa
xb
非线性分析设置
• 非线性静力分析的过程与线性静力分析很相似,因此本节中不再详述各 个步骤。黄色斜体字的步骤包含了非线性分析中的特有选项。 – 导入几何模型 – 指定材料属性(如果需要:金属塑性) – 定义接触选项 (如果需要) – 定义网格控制 (可选择) – 施加载荷或约束 – 选择需要查看的结果 – 设置非线性选项 – 求解模型 – 查看结果

14_Workbench_DS疲劳分析

14_Workbench_DS疲劳分析
• 材料的延展性, 材料的加工工艺 • 几何形状信息,包括表面光滑度、残余应力以及存在的应力集中 • 载荷环境, 包括平均应力、温度和化学环境
• 例如,压缩平均应力比零平均应力的疲劳寿命长,相反,拉伸平均应力 比零平均应力的疲劳寿命短. • 对压缩和拉伸平均应力,平均应力将分别提高和降低S-N曲线.
• 在材料特性的工作列表中,可以定义下列类型和输入的S-N曲线 – 插入的图表可以是线性的(“Linear”) 、半对数的(“SemiLog”即linear for stress, log for cycles)或双对数曲线( “Log-Log”)
– 记得曾提到的,S-N曲线取决于平均应力。如果S-N曲线在不同的 平均应力下都可适用的,那么也可以输入多重S-N曲线 • 每个S-N曲线可以在不同平均应力下直接输入 • 每个S-N曲线也可以在不同应力比下输入
max min
应力-寿命曲线
• 载荷与疲劳失效的关系,采用的是应力-寿命曲线或S-N曲线来表示: – 若某一部件在承受循环载荷, 经过一定的循环次数后,该部件裂纹或 破坏将会发展,而且有可能导致失效 – 如果同个部件作用在更高的载荷下,导致失效的载荷循环次数将减少 – 应力-寿命曲线或S-N曲线,展示出应力幅与失效循环次数的关系
应力定义
• 考虑在最大最小应力值min 和 max作用下的比例载荷、恒定振幅的情 况 – 应力范围 D 定义为 (max- min) – 平均应力 m 定义为 (max+min)/2 – 应力幅或交变应力 a是D/2 – 应力比 R 是 min/ max – 当施加的是大小相等且方向相反的载荷时,发生的是对称循环载荷. 这就是m = 0 ,R = -1的情况. – 当施加载荷后又撤除该载荷,将发生脉动循环载荷. 这就是 m = max/2 , R = 0的情况.

02 ANSYS13.0 Workbench 结构非线性培训 作业 小变形与大变形

02 ANSYS13.0 Workbench 结构非线性培训 作业 小变形与大变形

Workshop Supplement
• 改变 Solution Output 到力收敛 Force Convergence以便观察 Newton-Raphson 迭代的历史.
– 注意运行中发生的二分. 二分是从最后成功收敛子步开始以更小的增量 加载 (在指定范围使用更多载荷步 ).
WS2A-14
Workbench Mechanical – Structural Nonlinearities
WS2A-8
Workbench Mechanical – Structural Nonlinearities
… 作业 2A: 大变形
• 从初ement
WS2A-9
Workbench Mechanical – Structural Nonlinearities
已经建立了一端固定另一端压力载荷的弹簧模型.
– 高亮显示固支和压力载荷来确认模型适当约束和加载并准备求解.
WS2A-6
Workbench Mechanical – Structural Nonlinearities
… 作业2A: 大变形
• 注意求解设置
– Auto Time Stepping = Program Controlled – Large Deflection = Off
Workbench-Mechanical 结构非线性
WS2A-1
Workbench Mechanical – Structural Nonlinearities
作业 2A: 大变形
• 目标
Workshop Supplement
对比针对同一模型相同载荷和边界条件所进行的小变形、大变形分析结果.
• 模型描述
• • 高亮显示 Model 并改变名字为 “Large Deflection… 高亮显示 分析设置

工字梁Workbench非线性屈曲分析

工字梁Workbench非线性屈曲分析

工字梁Workbench非线性屈曲分析本文由Workbench小学授权转载1.失稳介绍根据失稳的性质,结构稳定性问题分为三类:①线性屈曲问题:当结构承受载荷达到某个值后,结构的平衡状态可能存在的同时,出现了第二个平衡状态。

压杆稳定性也属于这类问题。

②极限失稳点问题:当结构承受达到极限值后,结构的平衡状态可能存在的同时,不会出现新的平衡点③跳跃失稳问题:当结构承受达到极限值后,结构的平衡状态出现新的平衡状态。

本次仿真的主题是一个长度为1m的工字钢非线性屈曲分析。

2.建立模块关联3.建立模型注意:材料默认结构钢,也可以在材料库更改为存在双线性的结构钢。

4.划分网格注意:通过Mesh中的Edge sizing一Number of division选项将整个梁划分为100份。

5.边界条件添加边界条件说明:选择为节点过滤器,上端施加向下200N的力与远端约束,下端设置简支约束。

6.远端约束细节7.求解,后处理8.屈曲后处理查看其中在Eigenvalue Buckling的后处理中可以看到上图中的工字钢的屈曲载荷因子。

屈曲状态一般发生在低阶模态中,高阶不考虑该问题。

得到工字钢的屈曲载荷极限:屈曲载荷因子(5.29)x施加的力(200N)=1058N,退出并保存文件。

9.复制静力学模块注意:选择Duplicate选项进行复制,改变边界条件,将Force设置为1200N,并插入command设置初始缺陷,其他不变,设置1200N的原因是需要得到该工字钢是否在1058N处发生极限跳跃失稳。

注意:这里的初始缺陷施加为0.0001,具体值根据实际工程应用确定,command中的地址为源文件的存放地址。

10.分析设置注意:为了更清晰地捕捉到失稳点,将子步设置尽量大,并打开大变形与稳定性开关,如下图所示。

11.求解,后处理设置为总形变,Z方向形变与下端点的支座反力总变形图:最后通过chart功能建立变形与支座反力的图表,如下图所示:结论:从上图结果可以看出,工字型钢在1020N的逐渐开始失稳。

【ANSYSWorkbench仿真】非线性静力学分析(三):材料非线性分析

【ANSYSWorkbench仿真】非线性静力学分析(三):材料非线性分析

【ANSYSWorkbench仿真】非线性静力学分析(三):材料非线性分析contents•o材料非线性o▪超弹性本构模型▪▪Engineering Data 设置▪▪points▪材料本构参数▪蠕变▪相关单元技术▪单元中节点▪▪Kept(保留)OR Dropped(取消)材料非线性在应力水平低于比例极限时,应力应变关系为线性关系,超过这一极限后,应力应变关系为非线性关系,表现为非弹性或塑性(应变不可恢复状态)某种材料的应力应变曲线,主要材料试验有:单轴试验、等双轴试验、平面剪切试验、体积试验、松弛试验等。

超弹性本构模型超弹性(hyperelastic) 是指材料存在一个弹性势能函数,该函数是应变张量的标量函数,其对应变分量的导数是对应的应力分量,在卸载时应变可自动恢复的现象。

应力和应变不再是线性对应的关系,而是以弹性能函数的形式一一对应。

一般应有三种试验数据:单轴拉伸、双轴拉伸及平面拉伸Engineering Data 设置points•建立至少两倍于需要计算材料本构参数数目的数据点,同时需要考虑滞后等效应;•如果仅有单轴拉伸数据,则不能生成承受大剪切的模型;•试验数据应该包括全部关注的应变范围,如果只有0~100%应变的测试数据,则不能生成承受150%应变的模型,外推数据会产生极大的误差。

材料本构参数•Mooney-Rivlin 2 Parameter:拉应变为100%,压应变为30%,对于大应变,越高阶选项,精度越高。

•Neo-Hookean:应变能最简单的形式,可用应变范围为20%~30%。

•Arruda-Boyce:应变范围为300%。

Gent:应变范围为300%。

Blatz-Ko:描述可压缩泡沫橡胶材料的最简单形式。

•Polynomial 1st Order:与Mooney-Rivlin本构类似,等效于Mooney-Rivlin 2 Parameter•Yeoh 1st Order:一种缩减多项式超弹性本构,一阶等效于Neo-Hookean。

8_Workbench_DS热分析

8_Workbench_DS热分析

ANSYS BASIC TRAINING
稳态传热基础
• 对于一个DS中的稳态热分析, 温度 {T} 是由如下的矩阵求解:
K T T QT
这就导致了如下的一些假设: – 在稳态热分析中不考虑任何瞬态效应; – [K] 可以是常量或是温度的函数; • 每种材料属性中都可输入温度相关的热传导率; – {Q} 也可是常量或是温度的函数; • 在对流边界条件中可以输入温度相关的对流传热膜系数
ANSYS BASIC TRAINING
热分析使用的单元
• 在DS中, 可使用如下一些单元: – 实体通常用10节点四面体单元或 20节点六面体单元划分 • SOLID87和 SOLID90 – 面通常用 4节点 四边形壳单元划分 • SHELL57 使用实常数 • (SHELL131或SHELL132 目前已不使用) – 线通常用 2节点 线单元划分 • LINK33 使用实常数 • LINK33 还需要DM中定义的等效截面积 – 对热-应力分析,不使用耦合场单元。热-应力分析顺序进行,因 此可使用上面的热单元, 然后把温度场读进相应的结构单元中。
Inertial Loads
– 这是一些没有任何类比的载荷类型 • 对于像旋转速度、加速度之类的惯性载荷,没有热当量 • 对流边界条件的类比是一个结构分析中的“基础刚度”支撑, 类似于接地弹簧
ANSYS BASIC TRAINING
结果和后处理
• 可得到各种结果用于后处理: – 温度 – 热通量 – “反作用” 热流率
CT T K T T Qt, T
这里, t 是时间,{T} 是温度矩阵, [C] 是比热矩阵(热容), [K] 是热传导矩阵,{Q} 是热流率载荷向量。 • 在稳态热分析中, 所有时间相关的项都不考虑, 但非线性现象还可 能存在。

4_Workbench_DS简介

4_Workbench_DS简介
– 计算刚体和柔体组成的结构在时间载荷下的响应
• 稳态磁场分析
– 计算在稳态下空间的磁场分布
ANSYS BASIC TRAINING
与ANSYS求解器结合
• 用户还可以使用常规的ANSYS功能 ,但不是从DS中直接进入
– 用户可以在DS 环境下直接用前后处 理器(高级前后处理除外)。信息全 部保存在“branch”中,可以用到每 一个分析中
– 网格划分和加载也很容易地转换到常 规的ANSYS截面,对于习惯ANSYS界面 的用户是很有用的。
ANSYS BASIC TRAINING
Preliminary Decisions
Preprocessing
Solution
Postprocessing
ANSYS BASIC TRAINING
Design Simulation功能
• 线性应力分析:
– 基于承受静载荷下材料强度理论,计算应力、误差估计、安全系数等.
• 模态分析与预应力模态分析:
ANSYS BASIC TRAINING
基本过程
• 每一个分析都包括四个主要步骤: – 初步决定 • 是什么样的分析:静力、模态等等? • 什么样的模型: 零件或是组件? • 选什么样的单元: 平面的或是实体的? – 前处理 • 导入模型 • 给零件定义材料属性 • 划分网格 • 施加载荷和约束 • 要求结果 – 求解 – 后处理 • 观看结果 • 检验结果的正确性
• 形状优化(拓扑优化):
– 通过使用拓扑优化技术,对受载荷的零件体积优化给出预测
• 非线性结构分析:
– 计算静力载荷作用下的非线性系统变形和应力,包括材料非线性、 几何非线性和接触非线性.
• 柔体动力学分析

AnsysWorkbench工程应用之——结构非线性(序):概述

AnsysWorkbench工程应用之——结构非线性(序):概述

AnsysWorkbench工程应用之——结构非线性(序):概述各位道友,在时隔半年后,我又回来了,虽然由于水平有限,我依然无法给各位答疑解惑,但是我可以和大家一起学习学习一下基本理论知识。

从本文开始,图惜和大家开始共同学习结构非线性。

本文作为非线性知识的序篇,我们主要通过以下问题来学习结构非线性基本概念:问题1:什么是非线性?问题2:非线性的类型?1 什么是非线性前面介绍的许多内容都是结构线性问题,即满足胡克定律其中刚度矩阵[K]是一个常量,通俗地说,如果力F增大一倍,位移u也将增大一倍。

然而,实际工程中很多结构的力和位移的关系不呈线性关系,称之为非线性结构。

结构刚度不再是常数,而是随着载荷的变化而发生变化。

KT(切向刚度)代表了经过载荷位移曲线上的某一点,该曲线切线的斜率。

非线性问题分为三类,三类往往交叉出现:1.几何非线性,如大应变、大挠度大转动、应力钢化、旋转软化等。

2.材料非线性,如弹塑性、超弹性、蠕变、黏弹性等。

3.状态非线性,也称为边界条件非线性,如接触、生死单元等。

#2 非线性分析基础2.1几何非线性物体受载荷后,内部会发生变形,当应变远小于1%时,可忽略物体前后变形的形状和位置,简化为线性分析。

但是当结构承受大变形时,变形的几何形状可能会引起结构非线性响应。

一般几何非线性有大应变,大挠度,应力钢化等,它们的关系如下在Ansys Workbench中,如果要使用几何非线性功能,只需要打开分析设置中的大变形选项:Large Deflection=On,程序将考虑大应变,大挠度,应力钢化、旋转软化等效应。

(1)大应变,结构刚度由网格单元刚度和方向决定,单元的形状发生变化,从而最终引起结构的非线性响应。

所有的几何非线性现象几乎最终都会导致网格单元的大应变。

有限应变也属于这一类型,例如金属冷作成型过程中的有限塑性变形。

值得注意的是,大应变不一定导致大应力,如橡胶、海绵的大变形大变形不一定导致大应力。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
• 非线性特性有三个主要来源: – 几何非线性: 如果某个结构出现了大变形, 其变化 的几何外形会导致非线性行为。 – 材料非线性: 非线性的应力-应变关系, 如右图所示金 属的塑性, 也是非线性的另一个来源。 – 接触非线性: 接触效应是一种 “状态改变” 非线性, 当两接触体间互相接触或分离时会发生刚度 的突然变化,此时也会出现非线性。
Fb2 Fb1 Fa Fa1
ANSYS BASIC TRAINING
xa
xb
非线性分析设置
非线性静力分析的过程与线性静力分析很相似,因此本节中不再详述各 个步骤。黄色斜体字的步骤包含了非线性分析中的特有选项。 – 导入几何模型 – 指定材料属性(如果需要:金属塑性) – 定义接触选项 (如果需要) – 定义网格控制 (可选择) – – – – – 施加载荷或约束 选择需要查看的结果 设置非线性选项 求解模型 查看结果
Fa 3 F1 1 2 4
x1
x
ANSYS BASIC TRAINING
非线性求解
• 了解载荷管理 – 载荷步在通用加载中是变化的. • Simulation 通常在一个载荷步中求解所有非线性模型。但是,在有螺栓 预紧力载荷时,DS采用两个在载荷步,首先施加螺栓预紧力载荷, 然后 施加其它所有载荷. 这些载荷步可以认为是Fa和 Fb. – 子步以增量形式施加载荷 • 由于复杂的响应, 有可能需要按增量形式施加载荷。例如, Fa1 大约为Fa 的一半. 当 Fa1 收敛后,再施加全部的 Fa 载荷。在这个例子中 Fa 有 2个 子步而Fb有 3 个子步。 – 平衡迭代是修正求解以得到收敛子步 • 在右边的例子中,白色虚线间的迭代就是平衡迭代. Fb
ANSYS BASIC TRAINING
装配体 – 实体接触
对于高级用户,接触的另外一些选项可以进行修 改
– 方程式可以从 “Pure Penalty” 修改到 “Augmented Lagrange” 或者 “MPC”. • “MPC” 仅仅适用于绑定的接触 • “Augmented Lagrange” 应用于规则的 ANSYS模型中 – 在绑定的接触中,纯粹的罚函数法可以想象为在接 触面间施加了十分大的刚度系数来阻止相对滑动。 这个结果是在接触面间的相对滑动可以忽略的情况 下得到的。 – MPC 方程当中对接触面间的相对运动定义了约束 方程,因此没有相互的滑动。这个方程经常作为罚 函数法的最好的替代。

– 绑定的和不分离的接触是最基础的线性行为,仅仅需要一次迭代 – 无摩擦以及粗糙接触是非线性行为,需要多次迭代。但是, 需要注意的是仍然利用了小变形理论的假设。 • 当需要利用这些选项时,可以在相应的菜单下设 定 “Actual Geometry (and Specified Offset)” 或 “Adjusted to Touch”,其中允许用户调整ANSYS 模型的间隙到 ‘刚刚接触’ 的位置。
ANSYS BASIC TRAINING
装配体 – 实体接触
高级选项 (续):
– pinball region可以自己定义和显示出来 • pinball region定义了近距离开放式接触的位置 。而超出pinball region 范围之外的为远距离 开放式接触。 • 最初,pinball region 作为十分有效的接触探 测器使用,但是它也用于其它方面,例如绑定 接触等。 • 对于绑定或者不分离的接触,假如间隙或者渗 透小于pinball region,则间隙/渗透自动被删 除。
σ1 σ
• •
塑性
σy
弹性
ε σ3
σ2
主应力空间
单轴应力-应变
ANSYS BASIC TRAINING
硬化规则
硬化规则: • • 硬化规则描述屈服面的变化(大小,中心,形状)而导致塑性变形 硬化规则决定什么时候材料会随加载或卸载再次屈服 – 这和弹性、完全理想塑性材料形成鲜明对比,因为这种材料不会出 现硬化现象。就是说,屈服面保持固定。
σo =
1 2 2 2 (σ x − σ y )2 + (σ y − σ z )2 + (σ z − σ x )2 + 6 τ xy + τ yz + τ xz 2
[
(
)]
ANSYS BASIC TRAINING
Mises 屈服准则
在主应力空间,von Mises 屈服面是圆柱面. 材料在屈服面内呈弹性弹性行为。注意,多轴应力状态可以存在于圆柱 面内的任何一点。在圆柱面边上会产生屈服。没有应力状态会存在于圆 柱面以外。硬化规则会描述圆柱面随响应的变化。

ANSYS BASIC TRAINING
装配体 – 实体接触
四种接触类型可供选择:
Contact Type Bonded No Separation Frictionless Rough Iterations 1 1 Multiple Multiple Normal Behavior (Separation) Tangential Behavior (Sliding) Closed Closed Closed Open Open Open Open Closed
塑性 加载厚度的屈服面 弹性 初始屈服面
ANSYS BASIC TRAINING
各向同性硬化
各向同性硬化: • 各向同性硬化表明在塑性变化时屈服表面均匀膨胀. 术语 ‘各向同性’ 这 里是指屈服面在所有方向上均匀膨胀,这有别于“各向同性”屈服准则( 从材料的角度).
σ1
变形后的屈服 面
σ
σ' σy 2σ'
ANSYS BASIC TRAINING
非线性分析背景
在非线性静态分析中, 刚度 [K] 依赖于位移 {x},不再是常量:

[K (x )]{x} = {F }
– 从而力与位移的曲线将是非线性的, 如右图所示,因此当力加倍时, 位移和应力不一定会加倍。 – 非线性分析是迭代求解,因为载荷(F)和位移,响应(x)间的关系 之前并不知道。 – 不会考虑与时间相关的效应.

ANSYS BASIC TRAINING
装配体 – 实体接触
பைடு நூலகம்• 以下是ANSYS中接触算法的总结:
KEYOPT(2) Normal Direction 0 Augmented Lagrange 1 Penalty 2 MPC 3 Lagrange Multiplier 4 Lagrange Multiplier Tangential Direction Penalty Penalty MPC or None Penalty Lagrange Multiplier Normal Stiffness Yes Yes Tangential Stiffness Yes Yes Yes Contact Behavior Any Any Bonded, No Separation Any Any

ANSYS BASIC TRAINING
金属塑性
什么是塑性? • 当韧性材料承受的应力超过其弹性极限时就会屈服产生永久变形. – 塑性是只材料响应超过屈服极限. – 塑性对金属成形非常重要. – 塑性作为结构在服务中的能量吸收机制通样重要. • 在小塑性变形时就会破坏的材料是脆性材料. • 韧性响应在多数情况下比脆性响应安全.
[K ]{x} = {F }
F K x
– 由于假定 [K] 是常量, 因此本质上只允许线性行为。 如右图所示, 如果力加倍, 在线性分析中,位移(和应 力)也会加倍。 但在现实世界中很多情况下, 这种小位移理论是无效的. 此时就需要非线性分析.
ANSYS BASIC TRAINING
非线性分析背景
初始屈服面
ε σ3
σ2
ANSYS BASIC TRAINING
各向同性硬化
绘制应力-应变曲线可以方便的了解在加载和反向加载过程中的变化:

σ σ’ σy ε 2σ’
注意在随后的压缩屈服中,最 大应力等于拉伸阶段的最大应 力. 各向同性硬化通常用于大应变 和比例加载的仿真分析。它通 常不用于循环载荷的分析。
ANSYS BASIC TRAINING
应力-应变的曲线表示
曲线形状 • 2种不同类型的曲线表示:
σ
σ ε
ε
双线性
多线性
ANSYS BASIC TRAINING
仿真中的塑性总结
在仿真中金属塑性可以作为模型的一部分。但应记住以下几点: – 金属塑性处理弹性和非弹性(永久)变形. 当应力高于屈服强度时 就会发生非弹性或塑性变形。 当卸载后总会有可恢复的应变(弹性应 变). – 基于标量数据的应力-应变曲线通常来自单轴试验。而系统有可能 承受多轴应力状态, 因此在仿真中利用Mises 屈服准则通过标量试验 数据来近似多轴应力状态。在这中情况下,因该提供真实的应力应 变数据. – 当屈服发生后,由于应变硬化屈服点有可能会提高。这就改变了屈 服面,同时它在仿真中的发展取决于各向同性硬化假设 – 应力-应变曲线可以表现为双线性或多线性曲线

In this case, the gap between the two parts is bigger than the pinball region, so no automatic gap closure will be performed.
ANSYS BASIC TRAINING
载荷方向
DS非线性分析
ANSYS BASIC TRAINING
非线性分析
• 本章将会介绍下列内容: – 非线性理论的背景介绍 – 非线性分析设置 – 金属塑性 – 求解非线性模型 – 观察结果
ANSYS BASIC TRAINING
线性分析的背景
• 线性分析进行线性静态结构分析时的一些假设和限制 – 求解的矩阵方程是虎克定律:
相关文档
最新文档