用户登录

你在这里

雅可比矩阵

kapornith的照片

亲爱的专家,

我对雅可比矩阵在有限元分析中的解释有点困惑。雅可比矩阵在这里表示什么?这只是两组坐标系之间的变换矩阵,还是网格元素的变形(拉伸/剪切)的度量?

有什么特别的参考文献吗?

问候

Kapil沙玛

免费的标签:

通常,在有限元环境中,雅可比行列式用于有限元的等参公式。如果描述其几何形状的插值函数与描述位移分布的插值函数相同,则称为等参有限元。在等参公式中,一个元素的坐标系被转换成另一个称为自然坐标系的坐标系。在这种坐标系的变换中,需要控制元素的畸变,确定逆变换的存在性(从自然坐标系到初始坐标系),确定元素中可以计算导数的区域等。因此,为了实现这些,将雅可比矩阵定义为nxn矩阵(下面用MATLAB语言编写):

J = [diff (x1,ξ1),diff (x1,ξ2),……diff (x1,ξn);

差异(x2,ξ1),差异(x2,ξ2),……差异(x2,ξn);

diff (xn,ξ1),diff (xn,ξ2),……diff (xn,ξn))

它的行列式叫做雅可比行列式。你可以参考标准的有限元教科书了解更多细节。如果你只关注雅可比矩阵,而不考虑FEM,你可以认为雅可比矩阵和行列式与n个变量的n个函数相关联,通常用于坐标系的变换。

最好的问候,

乔治Papazafeiropoulos

__________________________

希腊空军中尉

土木工程师,硕士,博士研究生

kapornith的照片

嗨,乔治,

谢谢你的回复。这确实是一个很好的解释。我发现这里面的数学比我想象的要多。

非常感谢笑

Kapil

SivaSrinivasKolukula的照片

有限元中刚度矩阵、质量矩阵等的计算都涉及到积分。根据形状函数中使用的插值,这种集成可能简单也可能复杂。为了便于积分,我们使用等参数公式。(当然,使用等参数公式还有其他优点)。在等参公式中,将物理坐标转换为自然坐标。这个变换是由雅可比矩阵完成的。自然坐标介于+1和-1之间,与物理坐标长度无关。这使得我们的集成变得简单。因为我们有+1和-1的极限。数值积分技术可以很容易地使用。 SO in this formulation we are reducing the length of dimensions of the coordinates for simplyfing our integration. After integration is done, the length of dimesnions should be replaced, it is done by the Jacobian. Multiplying the integration obtained with natural coordinates with det(Jacobian) to get the integration in physical coordinates.

因此雅可比矩阵可以看作是用自然坐标的长度来描述物理坐标长度的尺度因子。

一维雅可比矩阵给出长度,二维雅可比矩阵给出面积,三维雅可比矩阵给出体积。

例:看笛卡尔坐标和极坐标的关系

dx。Dy = r.dr.dθ

这里r正好扮演雅可比矩阵的角色。这里J = r。

关于

Sreenu

kapornith的照片

谢谢Sreenu,

很好的解释:)

问候

Kapil

Jayadeep U. B。的图片

亲爱的Kapil,

你可能会对《力学》中雅可比矩阵的详细讨论感兴趣。万博manbetx平台请参阅下面的链接:

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

问候,

Jayadeep

附:致iMechanica的管理员:我们能不能想个办法避免同样万博manbetx平台的博客在iMechanica上重复出现?我建议你在发布博客之前提供一个类似主题的链接(我知道现在有一个“类似链接”选项;在发布新博客的过程中将其作为一个步骤,并使用一些基于相关性的排序机制,可能会完成这项工作)。我觉得最好避免重复,除非真的需要……谢谢。

索志刚的照片

亲爱的Jayadeep:好建议。当我们更新软件时,我们将研究选项。完全独立于软件,它有助于用户将新帖子链接到iMechanica上的相关帖子,以及其他在线和离线资源。万博manbetx平台谢谢你提供这样的链接。

kapornith的照片

亲爱的Jayadeep,

谢谢你的回复。这个查询更多的是一个有限元软件/有限元理论的雅可比项的比较。我最近得到了一个回复,我认为可能值得与大家分享。

在有限元法中,一个单元的雅可比矩阵将自然坐标空间中的量与实空间中的量联系起来。与理想形状单元相比,变形越大,物理量从自然空间到现实空间的转换就越差。实际上,这意味着元件形状的畸变会在从自然空间到实空间的数学转换中引入误差。

在有限元软件中,雅可比矩阵(也称为雅可比比)是对给定元素与理想形状元素偏差的度量。雅可比矩阵的取值范围是-1.0到1.0,其中1.0表示形状完美的元素。元素的理想形状取决于元素类型。检查是通过将参数坐标中的理想元素映射到全局坐标中定义的实际元素来执行的。例如,理想四元的角在参数坐标中的坐标为(-1,-1)、(1,-1)、(1,1)和(-1,1)。雅可比矩阵的行列式涉及到将其拟合到全局坐标空间所需的参数空间的局部拉伸。所附的图片说明了一个真实的QUAD4元素,经常在网格中发现,以及理想的QUAD4元素。

元素的雅可比矩阵是一个方阵,一维元素的维数为1x1,二维元素的维数为2x2,三维元素的维数为3x3。这个矩阵里的项在大多数情况下是依赖于参数坐标r st的函数。因此,一个元素的雅可比矩阵的行列式将取决于有限元包将在计算中使用哪些点。许多软件(例如HyperMesh)在每个元素的积分点(也称为高斯点)处计算雅可比矩阵的行列式。例如,一个有4个积分点的元素将有4个行列式值。有限元词汇中的雅可比矩阵定义为雅可比矩阵行列式的最小值与最大值之间的比值。因此,雅可比矩阵总是在0到1之间。

不同的求解器代码使用不同的集成点模式,并且相同的求解器可能对相同元素配置的不同表述使用不同的模式。有的求解器采用单元节点坐标计算雅可比矩阵,有的采用积分点坐标计算雅可比矩阵,有的采用积分点坐标计算雅可比矩阵,有的采用积分点坐标计算雅可比矩阵。也可以选择FE软件应该使用哪种类型的点(集成点或节点/角)。

这种止回阀包括一级和二级固体元件和板式元件;然而,一阶tria和tetra元素的雅可比矩阵值总是1.0,因为雅可比矩阵中的项是常量。


注释:有有限元软件将雅可比矩阵计算为雅可比矩阵行列式的最大值和最小值之间的比率。这是ANSYS和COSMOS/SolidWork的情况。在这些软件帮助中,有最大允许雅可比值的建议。这些建议是在软件公司的FEA专家进行详尽的相关测试后创建的。

有谁能告诉我如何在ABAQUS的UEL程序中找到二维4点四边形单元的雅可比矩阵和弹性矩阵吗?

提前感谢。

订阅评论查看

更多的评论

辛迪加

认购辛迪加