outcomes

  1. 理解在求解一般形状微元的FE矩阵的复杂计算,这为在自然坐标系下建立有限元提供了条件。
  2. 在自然坐标系下构建1D曲线元和四边形元,并推导FE矩阵。
  3. 在上述构建过程中学习转换的雅可比行列式。
  4. 理解数值整合来计算FE积分,尤其是高斯勒让德正交。(Gauss Legendre Quadrature)

普通形状的有限元

2D四边形:
FEA 笔记5
3D六面体:
FEA 笔记5

FE公式-物理坐标VS自然坐标

形状函数依赖于元素的几何形状.
积分限不容易定义,因此FE积分的计算变得很繁琐。
积分限制取决于元素形状
FEA 笔记5
FEA 笔记5
FEA 笔记5
FEA 笔记5
FEA 笔记5

自然坐标下的FE公式

积分范围一般为[0,1]or[-1,1]

FEA 笔记5

自然坐标下的FE公式-1D曲线

假设问题域是二维空间中的一条曲线,有限元法可以沿着线性单元求解,但由于线性单元无法精确匹配曲率,在精度上有所妥协。在这种情况下,二次元(在几何中)可能更合适。
FEA 笔记5
通常,FE的积分由以下确定:
FEA 笔记5
现在,假设可以定义从物理坐标到局部坐标的映射函数
FEA 笔记5
首先,坐标s由曲线的位置向量(x, y)定义。
FEA 笔记5

回想一下,shape函数用于定义字段变量
FEA 笔记5
同样的形状函数集也可以用于几何图形。因此,对曲线元素使用二次形状函数,
FEA 笔记5
从几何微积分中可以看出,初等长度为
FEA 笔记5
其中
FEA 笔记5
基本长度变成
FEA 笔记5
雅可比行列式(Jacobian of transformation)导数项:
FEA 笔记5
FE积分:
FEA 笔记5

1D积分计算

对于积分为多项式函数的FE积分,最有效的数值积分方法是高斯-勒让德积分法。
FEA 笔记5

对于高斯-勒让德求积,仅用高斯点就可以求出多项式函数(2 - 1)阶的精确解。
最后,利用高斯-勒让德积分法求出了有限元积分。
FEA 笔记5

相关文章:

  • 2021-12-04
  • 2022-03-03
  • 2021-11-16
  • 2022-01-12
  • 2022-12-23
  • 2022-12-23
  • 2022-12-23
  • 2022-12-23
猜你喜欢
  • 2021-07-18
  • 2021-11-28
  • 2021-06-30
  • 2021-06-02
  • 2021-06-16
  • 2021-07-02
相关资源
相似解决方案