用户登录

你在这里

为什么在非线性有限元中使用速率方程?

Ramdas chennamsetti的照片

嗨!

我有一个非常基本的问题。

在非线性有限元公式中,我们使用速率方程(虚功),但在线性有限元中,我们不使用速率方程。为什么? ?

是不是因为非线性解是迭代解(时间可能是虚时间)?

我要求那些有想法的人给出一些解释。

提前感谢,

问候,

——Ramdas

评论

黄睿的照片

前段时间我也有同样的问题。我相信你是对的,速率方程是用来迭代求解非线性方程的,例如,用共轭梯度法。我现在面临另一个关于速率方程的问题。有几种不同速率的定义(Lie速率、Jaumann速率等)。我们用哪一个?为什么?

谢谢。

RH

arash_yavari的照片

亲爱的Ramdas和Rui:

我不是这方面的专家,但我有一些意见可能会有所帮助。

1)正如Biswajit已经提到的,连续介质力学中现有的客观时间导数都是李氏导数(相对于空间速度场)。事实上,不同的目标导数只不过是同一个李氏导数的不同表示(使用度规张量来提高和降低指标)。众所周知,长期以来,没有“更多”的客观压力率。

需要注意的是,目标时间导数的范围比李氏导数大。下面两篇论文对此进行了详细的解释。

i) F. Bampi和A. Morro,连续介质物理的客观性和客观时间导数,物理基础10(1980),905-920。
ii)Jean-Luc Thiffeault,动力学系统的协变时间导数,物理学报34(2001),5875-5885。

2)非耗散固体的种类比超弹性固体大得多。下面这篇最近的文章讨论了其中一些微妙的问题。

王志强,非耗散固体的响应研究,中国机械工程学报,2007,33(4),357-367。

问候,

有关压力率的更多信息,请参阅我在Wikiversity网站上发布的笔记

http://en.wikiversity.org/wiki/Continuum_mechanics/Objective_stress_rates

使用哪种汇率并不是特别重要。重要的是,你的材料参数应该使你使用的速率反映现实。也就是说,在确定参数时应牢记特定的速率,而不仅仅是与小应变弹性(例如)中使用的参数相同。你要记住,一阶低弹性不是与历史无关的。

jaaumann速率之所以受欢迎,是因为它易于实现。

荣天的画

我张贴了一个问题,只是发现这个线程
讨论。我再次提出我的问题。

http://万博manbetx平台m.limpotrade.com/node/4725

在上面的帖子中,我们展示了这一点特鲁斯戴尔率可以
通过假设F =,简化为格林-纳赫迪速率。R
并且可以
进一步简化为若曼速率,假设W =。(.) R (T),
在那里,=。意思是近似相等

在拉伸主导变形时,三种速率给出不同的速率
应力率。
这通常被解释为我们需要一个不同的
不同物镜速率下的切向模量。
然而,它是
当我们使用a时,很难理解为什么我们需要改变“材料”模量
不同的“数学”形式的客观比率,因为他们都应该是
“等效”。

因此,一个简单的解释可能是,Jaumann和Green-Naghdi比率
当拉伸变形占主导地位时是“不准确的”。
我们都知道
在剪切变形中,若曼速率随剪切应力的变化而变化
Truesdell给出了确切的答案(线性变化)(参见Ted Belytschko, Wing Kam Liu和Brian)
莫兰。连续体与结构的非线性有限元。换句话说,
这三种速率形式可能不相等,特别是在拉伸时
变形优势种

一个简单的问题是,为什么我们经常用乔曼速率,而不是
大变形有限元分析中的Truesdell率?从…的观点来看
执行时,无论哪个速率我们都不应该有太大的差别
选择。

请参阅http://万博manbetx平台m.limpotrade.com/files/stress%20rate.pdf
三个比率之间的关系。

荣添

Ramdas chennamsetti的照片

谢谢你! !我下载了你的笔记。

在这里,我想了解为什么会出现速率方程(而不是哪个速率方程)?它的起源是什么?

为什么要用单位时间表示虚功?为什么我们不能用没有速率的非线性有限元公式,用没有速率的虚功?

请解释一下。

提前感谢,

问候,

——Ramdas

Siva P V Nadimpalli的照片

你好Ramdas,

我之前也有同样的问题,但我找到了一些有用的信息。

我没有在(公式)编码等方面做过广泛的工作,因此我不知道从计算的角度来看上述观点的任何原因。然而,如果我们从实验的角度来看,这是有道理的。

我想大家都知道,在不同的加载速率下,大多数工程材料在非线性部分(即塑性变形范围)的材料响应都发生了显著的变化。材料的弹性性能不受加载速率的影响(粘弹性除外)。这可能是由于一些微观机制。

因此,在弹性区域中可能不需要速率方程来模拟材料的响应。

我希望这些信息对你有所帮助。

问候,

湿婆

你好湿婆,

你在帖子中提到的是关于材料的本构行为。

但是,我们应该注意的是,本构方程的速率形式与速率相关不是同一个概念。

对于非线性塑性本构方程,加载速率也不影响塑性行为。

但速率表仍被广泛使用。

这是因为使用速率形式可以从增量描述中得到完整的应力-应变路径,这更加方便。

假设材料经历循环加载或棘轮,在这种情况下,用速率形式可以很容易地描述循环应力-应变行为,特别是塑性加载-弹性卸载的转变。但是使用总形式可能是一个大问题,因为在相同的应变下,你可能有多个应力值。你必须付出很多努力来规划这种情况。

速率依赖(粘塑性/粘弹性)的事情完全是另一个话题。也许我们可以开始另一个话题来讨论它。

问候,

Xiaoteng王

典型的棘轮s-s曲线,

Siva P V Nadimpalli的照片

你好
Xiaoteng王

谢谢你的解释,我想我没有理解清楚这个问题。作为一个做实验的人,当提到速率时,我是从本构行为的角度来看的。你说得完全正确,现在我明白问题的主要内容了。

问候,

湿婆

谢谢你的网站我在这里找到了很多有用的信息

这里有两个问题:

1)虚功速率(或虚功原理)

速率形式本构方程

这两种形式对于非线性有限元分析都不是必需的,尽管速率形式被广泛使用(特别是在显式代码中)。

Simo和Pister的“Computational Inelasticity”采用了不使用虚功原理的速率形式的方法。那本书讨论了弹性的速率形式本构方程的问题,并最终提出了全形式超弹性是最好的方法。由于各种原因,费率表用于依赖于历史的行为(其中一些由另一位评论者引用)。

如果你喜欢Belytschko、Liu、Moran的虚功原理的速率形式“连续体和结构的非线性有限元”,他给出了大量的例子。

速率形式的本构方程通常比总形式更容易构造(并且可能更容易用增量线性形式表示)。

我希望本构关系方面的专家能给我们提供一些关于使用速率本构方程的更深层次原因的见解。

——Biswajit

黄睿的照片

我认为,在求解非线性平衡方程的过程中,即使对于全形式超弹性,也需要用速率方程来定义切模张量。这是我的问题,因为切模量的定义取决于你使用的应力率。例如,在ABAQUS UMAT中,假设了一个对称的切模张量,但我不清楚切模是如何定义的,以及如何使用它来求解非线性方程。

RH

ABAQUS调用用户定义的切模量(或雅可比矩阵)来表示初始刚度矩阵K。

然后ABAQUS采用迭代法求解非线性方程f=Ku。

对称切模张量具有更好的迭代稳定性。

为了表述雅可比矩阵,必须首先写出应力增量张量(ds)和总应变增量张量(de)之间的关系。这种关系不容易表述,但从本构方程中我们最终可以推导出来。

一阶导数d(ds)/d(de)是材料雅可比矩阵j,为了实现这个目标,你需要一个位张量分析。

幸运的是,我们不需要一个非常精确的J,我们只需要一个近似于真实值的J来帮助求解器解出方程,而不需要这么多的迭代。

一旦你定义了自己的本构模型,你就可以用j来表示。因为你已经定义了塑性应变增量、总应变增量和应力增量之间的关系。

1)切线模量为dsigma/depsilon。如果是封闭形式那么正切可以不用速率形式来计算。

2)你是对的,在推导一致的切线模量时经常使用速率方程。然而,该速率方程可以通过在物质框架中对超弹性关系进行时间导数,然后在需要时将其推进到空间框架中来生成。这相当于取李氏导数,但与弹性材料的行为一致(没有历史依赖)。你可以在http://en.wikiversity.org/wiki/Nonlinear_finite_elements/Rate_form_of_hyperelastic_laws

3)你可以得到一个精确的关系,一致的切模量(除非弹性模量依赖于压力)后一点代数创伤。

这与低弹性模型不同——无论你选择什么样的客观应力率。

4)对于复杂的材料模型,如具有非关联流动规律的土壤塑性模型,切线模量往往是不对称的。我希望ABAQUS允许不对称切模,但我不太了解这个软件。

戈文德吉教授和阿查里亚教授都是其中一些问题的专家。让我们希望他们能抽出时间参与进来,给出他们的建议。我发现Simo 1988年关于came的论文特别有用(尽管非常密集)。

——Biswajit

不管怎样,既然你问了。以下是我的两点心得:

1)使用速率表没有什么根本的问题。一个人可以平等地
我们将使用非比率表格(对于正在考虑的情况)。然而,有时人们发现用利率来考虑更容易。

2)在连续介质力学中并非所有客观速率都是如此
李氏导数的简单变体——尽管可以有很多
表示为李氏导数的合理的简单变化
(情人眼里出西施)。真正重要的是:你能做得更多吗
用单速率或单速率轻松表达材料的行为
另一个?


力学和计算力学中速率形式的历史
“特别”更多的是为了方便。书面比率
形成有限变形的本构关系看起来更像
小应变公式。这对人们来说更容易
想想模型。不幸的是,当计算出现时,这也给人们一种错误的印象,认为他们知道如何在有限变形下实现它们
仿真代码。注意,我说的是思考!他们实际上并没有这么做
从纯粹主义者的观点来看是正确的。即使在今天,许多商业法规也是如此
做得不正确(有弹性或无弹性);他们省略了条款
这导致了违反热力学第二定律的本构评价。
西莫和皮斯特对此有一篇很好的论文。这里也参考了一份关于这个问题某些方面的小文件:

S. Govindjee,《准确性》
纺丝过程中Jaumann应力速率方程的积分及稳定性
Engng尸体。”比较,14,14-30(1997)。

4)速率形式仅有助于表示所需的切矩阵
在有限元计算中,当速率形式被精确地积分进去时
本构评估子程序。否则全局Newton-Raphson迭代需要的是算法切矩阵;也就是切矩阵
该算法采用相关联的积分率形式。

桑杰·戈文杰教授
加州大学伯克利分校

孙外庆的照片

亲爱的Govindjee教授:

我曾听到有人评论说,即使在小变形情况下,非关联塑性模型也可能违反热力学第二定律。你介意就这个话题给我们一些参考吗?

谢谢,

WaiChing

我想我们都知道,当受到不同的加载速率时,
大多数工程材料显示家庭文凭材料的重大变化
非线性部分(即塑性变形范围)的响应。
材料的弹性性能不受拉伸速率的影响
加载(粘弹性除外)。这可能是由于一些微观的
机制。

高中毕业|家庭学校文凭

如果你喜欢速率形式的虚功原理“非线性”
连续体和结构的有限元”,作者:Belytschko, Liu, Moran
给出了大量的例子。论文写作帮助


似乎我们正在谈论的“速率”有两个含义:“速率”相对于实时和“速率”相对于一些参数(例如,对负载增量)。

如果材料性能是速率相关的,或者系统中有一些能量耗散取决于物体的速度(例如与速度相关的摩擦),则应使用与时间相关的“速率形式”虚功原理(虚功率原理)。在这种情况下,测试场通常被称为“虚拟速度”,基本的未知变量是结构中的TRUE速度场。

对于超弹性材料的分析,由于这是一个保守系统,平衡方程可以用“总形式”表示。但这种“总形式”方程在状态变量(如位移场等)方面是非线性的。因此,对其进行线性化和迭代求解是必要的。在迭代的每一步,都应该使用虚功原理的“增量线性形式”。在这种情况下,测试场通常被命名为“虚拟位移增量”,基本的未知变量是与当前负载阶跃相关联的结构中的TRUE(线性化意义上的)增量位移场。在此线性化过程中也得到了切线模量。

对于与时间速率无关的材料的塑性分析,由于它总是与路径相关,因此必须使用“增量形式”。在这种情况下,在迭代的每一步,也应该使用虚拟工作原理的“增量线性形式”。

从我的观点来看,对于与雷尔时间无关的问题,术语“增量线性形式”比“速率形式”更合适。当然,您也可以认为“速率”是相对于负载增量的速率。

希望这对你有所帮助。

致以最亲切的问候
X.Guo

黄睿的照片

我感谢Xu, Xiaoteng和Biswajit帮助我解决了超弹性中的切线模量问题。比斯瓦吉特在维基大学的文章对详细的张量代数特别有帮助。

从超弹性的任意应变能密度函数开始,我可以得到有限元弱形式的全形式非线性方程。为了通过迭代求解这个非线性方程,我必须以增量形式将其线性化。这里我遇到了切线模量的问题,因为基尔霍夫应力的增量由两部分组成。一个可以写成c:d其中c是第四个弹性张量比斯瓦吉特写了。此外,增量还取决于当前状态下的真实应力以及旋转张量,通常是非对称的。第二部分用李氏导数就消失了。然而,我相信这两个部分都需要解决非线性方程。在ABAQUS UMAT中,似乎只需要第一部分。我猜第二部分是由非线性求解器自动处理的,因为它并不真正依赖于材料。但我不确定。我将感谢任何ABAQUS专家的评论。

RH

你好RH,

在UMAT中,默认的公式是低弹性的,而不是超弹性的。如果激活大变形选项,则采用准弹性公式,求解器自动计算并传递旋转张量R。这意味着ABAQUS在默认的准弹性结构中计算并传递状态变量,并自动将应力速率旋转到Jaumann速率。你之前提到的第二部分是由ABAQUS计算的,但可能与其他公式不兼容....我n your case you write your constitutive model based on hyperelastic hypothesis, thus you have to be careful to define your own stress tensor and store it to the state variable matrix and you'd better not use the stress increment tensor given by ABAQUS.

问候

Xiaoteng王

没有速率方程就不能表示精确的路径。

虚功W=int(x(y)dy),其中x是载荷,通常是y的隐函数,dy是虚位移。

上面的方程是一个积分过程。

要得到相同的W有上千种可能的路径。

速率方程隐式地描述了确切的路径。

在线性有限元公式中,x与y的关系是简单的线性关系,因此我们可以直接写出虚功。

同样在非线性有限元中,如果你能明确地得到积分函数,你可以不用速率形式来表达它。

Ramdas chennamsetti的照片

嗨!

谢谢大家的讨论。我穿过了这条线。

如果本构定律依赖于速率,我理解了速率方程的使用。我真的不明白为什么速率无关本构定律需要速率方程。

比如说,我在研究几何非线性或接触分析,我需要用速率方程吗?请简单地解释一下。

谢谢你!

问候,

——Ramdas

1)即使本构律与速率相关,也不需要速率方程。

例如S = A E + B dE/dt

并不是一个真正的“速率方程”,尽管dE/dt出现在右边。我们可以把dE/dt看作一个自变量。

比率本构关系通常是这样的形式

dS/dt = A(E) dE/dt

如前面讨论过的,这种关系可以有多种用途。

3)你不“需要”速率方程,除非速率出现在公式的某些部分。准静态问题尤其如此。再一次看到和芮的讨论。

4)最简单的方法是生成你自己的接触公式,看看需要什么本构模型。Belytschko的书使用了一种虚拟权力的方法。这不是必须的。你可以用位移和功的变化来代替速度的变化。

——Biswajit

Ramdas chennamsetti的照片

谢谢你! !

正确的。我在看Belytschko的书。这种方法就是虚拟权力。你的意思是说,空间速度梯度的变化必须用变应变来代替吗?那么我们就需要相应的修改本构法。我说的对吗?

谢谢你!

与问候,

——Ramdas

Ramdas chennamsetti的照片

Govindji教授,

非常感谢你的解释。感谢Biswajitda博士和所有其他参与这个主题的人。

与问候,

——Ramdas

你不“需要”速率方程,除非速率出现在公式的某些部分准静态问题尤其如此再一次看到和芮的讨论。

区分本构方程的速率形式与速率相关的不是同一个概念。
当本构方程为非线性塑性时载荷也不影响塑性行为。但是速率形式至今仍被广泛使用。这是因为使用速率形式可以从增量描述中得到完整的应力-应变路径

索志刚的照片

我写下了关于框架无关性的讲义然后发布到网上。有几个人发表了评论。评论和注释都可能对这个讨论线索有所贡献。

注意一个建设性的介绍基于固体力学作为简单拉格朗日系统在b[1]中的表述的目标时间导数。这可以区分变形率和应力率,前者原则上是李导数,后者实际上是代表变形过程的曲线上的协变导数。此外,强调了Daleckii-Krein公式在理解广义应变理论中的作用,并特别注意了对数时间导数。

[10]费亚华,张志强。目标时间导数修正。ZAMP 71,文章编号:4 (2020)

订阅评论“为什么在非线性有限元中使用速率方程?”

最近的评论

更多的评论

辛迪加

认购辛迪加