用户登录

你在这里

为什么非线性有限元中的速率方程?

Ramdas chennamsetti的照片

嗨!

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

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

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

我请有想法的同学解释一下。

提前谢谢你,

问候,

——Ramdas

评论

黄睿的照片

我也曾想过同样的问题。我相信你是对的,迭代求解非线性方程需要速率方程,例如,用共轭梯度法。我现在要面对另一个关于速率方程的问题。有几种不同的速率的定义(谎言率,Jaumann率等)。我们用哪一个?为什么?

谢谢。

RH

arash_yavari的照片

亲爱的Ramdas和Rui:

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

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

人们应该注意到,目标时间导数的类别比李导数更大。下面的两篇论文对此进行了详细的解释。

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

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

K.R. Rajagopal和A.R. Srinivasa,关于非耗散固体的响应,皇家学会学报A 463(2007), 357-367。

问候,

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

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

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

乔曼比率很受欢迎,因为它很容易实施。

容天的照片

我发了个问题,发现了这条
讨论。我再发一遍我的问题。

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

在上面的文章中,我们展示了这一点Truesdell率可以通过
假设F .=,简化为Green-Naghdi速率。R
可以是
假设W .=,进一步简化为Jaumann速率。(.) R (T),
在那里,=。均值近似相等

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

一个简单的解释就是Jaumann和Green-Naghdi速率
当拉伸变形占主导时是“不准确的”。
正如我们所知
在剪切变形中,Jaumann速率给出了一个随剪切应力变化的“sin”,而
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)虚功率(或虚权原理)

2)速率形成本构方程

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

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

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

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

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

——Biswajit

黄睿的照片

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

RH

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

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

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

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

那么一阶导数d(ds)/d(de)就是物质雅可比矩阵j。为了达到这个目标,你需要一个位张量分析。

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

一旦你定义了你自己的本构模型,你总是可以制定j,因为塑性应变增量,总应变增量和应力增量之间的关系已经由你自己定义。

1)切线模量为d / desilon。如果你有封闭形式的,那么切线可以不用速率形式来计算。

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

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

对于低弹性模型来说,情况并非如此——无论你选择什么样的客观应力率。

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

Govindjee教授和Acharya教授都是这些问题的专家。让我们希望他们能抽出时间来发表意见。我发现Simo 1988年在cmae上的论文特别有用(尽管非常密集)。

——Biswajit

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

1)使用汇率形式没有什么基础。一个人可以
我们将使用非费率形式(对于正在考虑的情况)。然而,有时人们发现从比率的角度考虑更容易。

2)在连续介质力学中,并非所有的客观速率都是如此
Lie导数的简单变化——尽管很多可能是
表示为李导数的合理简单变分
(在旁观者的眼中)。真正重要的是:你能做得更多吗
简单地表达你的材料的行为在一个速率或
另一个?


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

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

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

Sanjay Govindjee教授
加州大学伯克利分校

孙怀志的照片

亲爱的Govindjee教授:

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

谢谢,

WaiChing

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

高中学历|家庭学校文凭

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


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

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

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

对于不依赖于时间速率的材料,由于其始终是路径相关的,因此必须使用“增量形式”。在这种情况下,在迭代的每一步,也应该使用虚拟功原理的“增量线性形式”。

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

希望这能有所帮助。

致以最亲切的问候
X.Guo

黄睿的照片

非常感谢Xu、Xiaoteng和Biswajit在高弹性理论中的切线模量方面的帮助。Biswajit在维基学院的文章对详细的张量代数特别有帮助。

从高弹性的任意应变能密度函数出发,我能够得到有限元的弱形式的全形式非线性方程。为了用迭代法解这个非线性方程,我必须以增量形式线性化它。这里我遇到了切线模量的问题,因为基尔霍夫应力的增量由两部分组成。一个可以写成c:d其中c是第四个弹性张量Biswajit写的.此外,增量还取决于当前状态下的真应力和旋转张量,旋转张量一般是不对称的。通过使用李导数,第二部分消失了。但是我相信这两部分都是解非线性方程所需要的。在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看做一个自变量。

2)速率本构关系通常是这样的

dS/dt = A(E) dE/dt

正如前面的评论所讨论的那样,这种关系可以有各种用途。

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

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

——Biswajit

Ramdas chennamsetti的照片

谢谢你! !

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

谢谢你!

与问候,

——Ramdas

Ramdas chennamsetti的照片

mr . Govindji教授,

非常感谢你的解释。感谢Biswajitda博士和其他所有参与这篇文章的人。

与问候,

——Ramdas

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

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

索志刚的照片

我写了我关于框架冷漠的讲座笔记然后发到网上。有几个人发表了意见。评论和注释都可能有助于讨论。

注意一个建设性的介绍在[1]中,基于固体力学公式作为一个简单的拉格朗日系统的客观时间导数。这使我们能够区分变形率(原则上是李氏导数)和应力率(实际上是沿曲线表示变形过程的协变导数)。此外,还强调了Daleckii-Krein公式在理解广义应变理论中的作用,并特别注意了对数时间导数。

[1]菲亚拉,Z.:客观时间导数修正。ZAMP 71,文章编号:4 (2020)

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

最近的评论

更多的评论

辛迪加

订阅辛迪加