首页/文章/ 详情

连续介质力学的基本概念

11月前浏览1169

连续介质力学理论是一种在宏观上模拟固体和流体的理论,它忽略了不均匀性,如分子、颗粒或晶体结构。因此,在这个理论中,固体和流体的行为可以是空间变量的拥有属性平滑函数。连续介质力学的主题包括以下基本课题: (1)不考虑原因的运动和变形的研究(运动学) ,(2)内力的研究(动力学) ,(3)保守方程或平衡原理,表明有某些重要的物理性质,如质量,动量和能量必须保守,和(4)本构模型,提供运动学和动力学变量之间的关系。它是通过一个本构模型,在连续介质力学中,我们可以区分固体和流体,橡胶和岩石等等。

本节目录如下:

1.1 运动与变形

1.2 应变度量

1.3 应力度量

1.4 客观应力率

1.5 守恒方程

1.6 本构模型

1.1 运动与变形

在连续介质力学中,物体B被认为是由无限组材料点形成的,这些材料点被赋予了一定的力学性质。物体初始未变形构型中材料点的位置矢量表示为X,相对于某个坐标基。X被命名为材料或拉格朗日坐标。同一材料点在变形配置中的位置由称为空间坐标或欧拉坐标的x表示。

固体的运动(变形)用函数φ(X,t)来描述。空间坐标和材料坐标之间的关系可以建立如下:

 

物体的位移场、速度场和加速度场是描述物体运动的主要运动学场。材料点X的位移,用u(X,t)表示,是其当前位置φ(X,t)和初始位置φ(X,0)之间的差。所以

 

一个物质点 X 的速度,用 v (X,t)表示,定义为该物质点位置的变化率,即

 

这就是拉格朗日速度场。该速度场存在欧拉形式,但由于MPM采用拉格朗日描述,因此此处不再讨论

材料点X的加速度是其速度的变化率,即速度的材料时间导数,

 

变形梯度张量F是有限变形连续力学中的一个关键量,因为所有的变形量都是从它导出的。它是一个线性映射算子,将参考配置中的每个无穷小线性元素dX映射为当前配置中的无穷小线性元dX。它被定义为

 

其次,介绍了物质时间导数的概念。要理解这一重要概念,请考虑以下情况。假设我们在物体上定义了一个特定的场(标量、矢量或张量),我们想知道在给定材料点X的变化率。这被称为φ的材料时间导数。这个概念有两个定义,分别对应于物质描述和空间描述

1. 拉格朗日描述。在拉格朗日描述中,自变量是材料坐标X和时间t。因此,我们所要做的就是取给定场的偏导数。对于材料场φ(X,t),其材料时间导数为

 

其中前两个方程指示材料时间导数的标准表示法。由于MPM采用拉格朗日描述,因此物质时间导数非常简单

2. 欧拉式的描述。所考虑的字段为φ(x,t)。这种情况要复杂得多,因为不仅时间变化,而且考虑的粒子的空间位置x也变化。我们必须计算材料描述的关于时间的偏导数,保持x固定

 

利用链式规则,我们得到了欧拉标量场的物质时间导数的重要公式

 

项dφ/dt被称为空间-时间导数,项φ,jvj是对流项,也称为输运项。

1.2 应变度量

在非线性连续介质力学中,有许多度量应变和应变速率的方法。对于任何刚体运动,特别是刚体转动,应变度量必须消失。在此,我们回顾一些常用于非线性连续介质力学的应变度量方法。它们是右柯西-格林变形张量 C,格林应变张量 E,和变形速率张量 D

右柯西-格林变形张量写为

 

其中上标T表示转置运算符。格林应变张量由下式给出

 

这个应变张量度量 dX 和 dX 长度的平方的差

速度的空间梯度或速度梯度张量 L 定义为速度的空间梯度

 

速度梯度 L 允许变形梯度 F 的材料时间导数写为

 

在第二个等式中,我们使用了拉格朗日场的物质时间导数随物质梯度而交换的事实。注意到,这一事实并不普遍适用于欧拉领域。

将速度梯度张量 L 分解为对称和反对称两部分

 

其是二阶张量的标准分解。变形率张量D被定义为L的对称部分,而自旋张量被定义为斜对称部分。使用这些定义,可以编写

 

其中 D 描述了拉伸和剪切的速率。

1.3 应力度量

由于存在不同的应变度量,因此也存在与之工作共轭的不同的应力度量。最常用的应力张量是(1)柯西应力,(2)基尔霍夫应力。这些应力张量之间的关系如表1所示。柯西应力是与变形体积的变形率D的真应力和功共轭。基尔霍夫应力也称为加权柯西应力系与初始体积的变形张量率的功共轭。第一次PK应力是非对称的,与变形梯度的速率共轭功。

第二皮奥拉-基尔霍夫应力是一个完全材料对称的应力张量,与格林应变速率张量是功共轭的。请注意,一些作者,例如Belytschko等人(2000),更喜欢使用标称应力,其转置为第一PK应力。

表1 不同应力测量之间的关系

 

1.4 客观应力率

本构方程常用速率表示,即应力速率与变形速率之间的关系。在大转动情况下,简单地使用应力张量的材料导数,例如 Cauchy 应力 Dσ/Dt 的速率,是错误的,因为它不能在叠加刚体运动下正确地转化为张量。我们讨论了三种常用的客观应力率: Jaumman 率、 Truesdell 率和 Green-Naghdi 率。一个本构模型可以根据这些目标应力率中的任何一个来制定,并且从一个应力率转换到另一个应力率需要重新制定本构模型。用方程式可以表示为:

 

1.5 守恒方程

连续介质力学中一组重要的方程式是守恒方程式或平衡方程式。对于热力学系统,守恒定律包括

1. 质量守恒

2. 线性动量守恒

3. 角动量守恒定律

4. 能量守恒

质量守恒。质量守恒定律用质量守恒方程来描述,或者通常称为连续性方程

 

如果密度不变,即材料是不可压缩的,因此密度的材料时间导数消失,连续性方程变为vi,i=0,这是众所周知的不可压缩性条件。对于拉格朗日描述,质量守恒的一个更简单的代数方程由下式给出

 

线性动量守恒。这个定律要求线性动量在时间上的变化等于作用在物体上的所有外力(体积和表面力)的总和。它用所谓的动量方程来描述。

 

这个定律要求柯西应力是一个对称张量角动量守恒定律。

能量守恒。这条定律指出,物体内总能量(包括内能和动能)的变化率等于外力所做功的速率加上热通量 q 和能源所提供的功的速率。

 

其中 e 是特定的内能; ρs 表示单位体积的热源。

1.6 本构模型

前面给出的所有方程都与材料无关:它们对固体和流体都有效。为了对材料行为进行建模,需要一个本构方程或本构关系——与运动量(如应变)相关的运动量(例如应力)之间的关系。正是通过本构方程,人们才能区分流体与固体、混凝土与橡胶等。第一个本构方程(本构定律)由Robert Hooke提出,被称为Hooke定律。它涉及线性弹性材料的情况。从那时起,已经开发了大量的构造模型来表征各种各样的天然和工程材料。值得注意的是,数值模拟只能与所使用的材料模型一样精确。

来源:STEM与计算机方法
非线性UM理论材料数字孪生人工智能
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-11-01
最近编辑:11月前
江野
博士 等春风得意,等时间嘉许。
获赞 58粉丝 96文章 153课程 0
点赞
收藏
作者推荐

模拟固体开裂的数值方法:断裂相场法(fracture phase-field method)

文一: 裂缝变分一致相场各向异性损伤模型摘要:在脆性和准脆性固体中断裂的相场建模中,表征材料的非对称拉伸/压缩行为至关重要。现有的相场模型通常采用自由能密度的直观划分而不正确地捕捉裂纹边界条件,或者采用在失去变分一致性的情况下的特设混合公式。为了解决这个问题,本工作在脆性断裂和准脆性破坏的统一相场理论框架内提出了一个可变一致的相场各向异性损伤模型Wu[1,2]。与断裂的变分方法一致,采用了能量范数Wu和Cervera[3]中有效应力的正/负投影,最大限度地减少了驱动裂纹演化的储存能量的拉伸部分。自然出现了一个圆形的兰金标准来真实地描述脆性和准脆性固体的局部失效,而不需要特别的假设。应变局部化后恢复了混合模式粘结区模型,相关参数由约束拉伸下软化棒的解析解确定。有代表性的数值例子表明,该模型可以独立于网格离散化和长度尺度参数来捕捉固体中任意裂纹的扩展。值得注意的是,伪应力锁定并没有表现出来,这是出了名的伴随着经典的各向异性损伤模型。 图:具有尖锐裂纹/界面和几何正则化的开裂固体。在这里,眼睛和鼻子被视为嵌入熊面部(基质相)的不同相(未损坏的弹性夹杂物)。嘴部和不同相之间的界面都在几何上正则化为具有有限厚度的窄局部化带。 图:剪切下的单边缺口板:位移u*=0.02 mm时的损伤轮廓文二: 断裂相场模拟摘要:断裂是工程材料最常见的失效模式之一。因此,防止开裂引起的失效是结构设计中的一个主要问题。裂缝的计算建模是一种不可或缺的工具,不仅可以预测裂缝结构的失效,还可以深入了解混凝土、岩石、陶瓷、金属和生物软组织等许多材料的裂缝过程。本章从以下方面对所谓的相场断裂/损伤模型(PFM)的文献进行了广泛的概述,特别是脆性和准脆性材料的准静态和动态断裂。PFMs是多种断裂方法的正则化版本,该方法推广了Griffith的脆性断裂理论。它们可以通过非常简单的实施来处理拓扑复杂的裂缝,如二维和三维的起始、交叉和分支裂缝。我们的目标之一是证明PFM越来越受欢迎。为此,讨论了理论和计算两个方面,并提出了用PFMs成功和失败解决的广泛的基准问题(准静态和动态脆性/内聚断裂)。还记录了有待进一步调查的未解决问题。 图:具有两个平行缺口的岩石样品:荷载-位移曲线以及裂纹模式与“45-2a-2a”几何形状实验的比较 图:动态裂纹分支: 不同时刻的裂纹模式对于牵引力 σ01=42.5 MPa文三: 固体多物理裂缝的相场内聚区建模及其在Comsol Multiphysics中的开源实现摘要尽管纯力学问题中裂缝的相场模型很受欢迎,但它们在多物理问题中裂缝建模中的应用却很少被报道。这一方面可能归因于多物理现象所涉及的理论复杂性,另一方面也可能归因于这些耦合模型在自制平台中的繁琐实现。在这项工作中,相场内聚区模型(PF-CZM)被用作原型模型,以解决各种多物理问题中的断裂,例如热机械、化学机械、化学热机械、机电等。相关的理论和数值方面被归类为模块化结构,并详细介绍了在软件平台Comsol Multiphysics中的开源实现。为了验证PF-CZM在多物理问题中的断裂及其数值实现,考虑了一些具有代表性的基准实例。不仅将定性的裂纹模式,而且将定量的整体响应与可用的实验测试数据进行了比较。研究发现,在所有考虑的多物理问题中,裂缝的典型特征都得到了很好的捕捉。此外,与纯力学对应物一样,预测的裂纹模式和整体响应对相场长度尺度不敏感,这使得PF-CZM有望在其他更复杂的多物理问题中对裂纹进行建模。 图:加热陶瓷板的淬火试验:对不同相场长度尺度参数的裂纹模式进行数值预测。 图:约束混凝土环试验:不同长度尺度参数下不同时刻裂纹相场的演变。文四: 相场断裂的块体结构自适应网格细化和强形式弹性方法及其在分层、裂纹分支和裂纹偏转中的应用摘要断裂是大多数复合材料工程结构中普遍存在的现象,也是导致灾难性失效的主要机制。在过去的几十年里,出现了许多建模和预测裂纹失效的方法。断裂的相场方法使用替代损伤场来模拟裂纹扩展,消除了对明确裂纹网格划分的艰巨需求。在这项工作中,提出了一种新的数值框架来实现非均质材料中的混合相场断裂。所提出的方法基于“无回流”方法,用于在块结构自适应网格细化(BSAMR)网格上以强形式求解线性弹性方程。BSAMR的使用实现了高效和可扩展的重新分级,便于使用时间子循环进行显式时间积分,并允许以最小的计算成本在裂纹边界进行超高细化。该方法适用于各种简单的非均质结构:层压板、波状界面和圆形夹杂物。在每种情况下,都会进行无量纲参数研究,以确定行为区域,从而改变问题的几何形状和相对断裂能量释放率。在层压板和波状界面的情况下,分层和断裂区域对应于简单的分析预测。对于圆形夹杂物,夹杂物的模量比以及分层能量释放率和问题几何形状都是变化的。在这种情况下,观察到了各种各样的行为,包括偏转、分裂、分层和纯断裂。 图:BSAMR多级细化方案示意图 图:内含物:选择裂纹行为的可视化。在图像网格中,每列三个图像对应于图9上的一个点。顶部图像对应于0.5的模量比(软夹杂物),中间的模量比为1.0(相同材料),底部的模量比是2.0(硬夹杂物)。方框轮廓颜色对应于分类(夹杂物断裂、分支、分层、基质断裂)。文五: 超弹性材料大变形断裂的变节点自适应相场模型摘要在这项工作中,提出了一种自适应相场方法来模拟超弹性材料在大变形下的断裂。可变节点元素有助于自适应网格细化,可变节点元素可以灵活地用作所采用的四叉树网格中的过渡元素。为了控制自适应过程,我们提出了一种相场和能量网格细化的组合准则,其中能量部分利用了本工作中使用的AT1相场模型的应变能阈值。考虑了可压缩和不可压缩的新胡克模型,后者由平面应力条件强制执行,以简化实现。与可用的数值和实验参考数据以及固定的局部预细化网格相比,研究了几个具有代表性的例子来验证所提出的自适应相位场方法的准确性和效率。仿真结果表明,当相场启动时,网格细化准则的能量部分可以有效地防止延迟损伤演化。最后,模拟了含有夹杂物的超弹性复合材料的断裂过程,以证明所提出的方法在大变形下再现复杂失效现象的能力。 图:在具有过渡元素的四叉树网格方案下,由相场和能量准则相结合控制的自适应网格细化:(a)在网格细化之前;(b) 网格细化后;(c) 分层四叉树结构。 图:张力作用下的竹筒形预裂试样。不同加载阶段的相位场轮廓:(a)-(e)在未变形配置中,具有自适应细化网格图;(f) –(j)在没有网格图的相应变形配置中。请注意,以下变形配置中不使用放大倍数。 图:含有刚性夹杂物的超弹性复合材料。不同加载阶段的相位场轮廓:(a)-(e)在未变形配置中,具有自适应细化网格图;(f) –(j)在没有网格图的相应变形配置中。(k) –(o)实验结果来源:STEM与计算机方法

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈