弹性力学仿真软件：FEMAP：FEMAP网格生成技术.pdfVIP

弹性力学仿真软件：FEMAP：FEMAP网格生成技术

1弹性力学基础

1.1弹性力学基本概念

在工程和物理学中，弹性力学是研究物体在外力作用下变形和应力分布的

学科。它基于材料的弹性性质，即当外力去除后，物体能够恢复其原始形状。

弹性力学的核心概念包括：

弹性体：能够发生弹性变形的物体。

应力：单位面积上的内力，通常用张量表示，分为正应力和剪应

力。

应变：物体变形的程度，也是用张量表示，分为线应变和剪应变。

弹性模量：描述材料弹性性质的物理量，包括杨氏模量、剪切模

量和泊松比。

1.2应力与应变分析

1.2.1应力分析

应力是弹性力学中的关键概念，它描述了物体内部的力分布。在三维空间

中，应力可以用一个3x3的对称张量表示，包括三个正应力分量和三个剪应力

分量。

1.2.2应变分析

应变是物体变形的度量，它与应力通过胡克定律（Hooke’sLaw）相关联。

在弹性范围内，应变与应力成线性关系，由弹性模量决定。

1.2.3胡克定律示例

假设有一个立方体，其杨氏模量为，泊松比为。当立方体受到均匀的拉

伸力时，我们可以使用胡克定律来计算其应变。

#胡克定律计算应变示例

#定义材料属性

杨氏模量，单位：帕斯卡（）

E=200e9#Pa

nu=0.3#泊松比

#定义应力

单位：帕斯卡（）

sigma_x=100e6#Pa

1

#计算应变

epsilon_x=sigma_x/E#线应变

epsilon_y=epsilon_z=-nu*epsilon_x#剪应变

#输出结果

print(f线应变：{epsilon_x:.6f})

print(f剪应变（y方向）：{epsilon_y:.6f})

print(f剪应变（z方向）：{epsilon_z:.6f})

1.3有限元法原理

有限元法（FiniteElementMethod,FEM）是一种数值方法，用于求解复杂

的弹性力学问题。它将连续体离散成有限数量的单元，每个单元用简单的数学

模型表示，然后通过组合这些单元来近似整个结构的行为。

1.3.1有限元法步骤

1.结构离散化：将结构划分为多个小的单元。

2.单元分析：为每个单元建立力学模型，通常使用矩阵表示。

3.整体分析：将所有单元的力学模型组合成一个整体的系统方程。

4.求解：使用数值方法求解系统方程，得到结构的应力和应变分布。

1.3.2有限元法示例

以下是一个使用Python和SciPy库求解简单梁的有限元法示例。假设梁的

长度为1米，截面为矩形，宽度和高度均为0.1米，材料的杨氏模量为200GPa，

泊松比为0.3。

importnumpyasnp

fromscipy.sparseimportdiags

fromscipy.sparse.linalgimportspsolve

#定义梁的属性

length=1.0#梁的长度，单位：米

width=0.1#梁的宽度，单位：米

height=0.1#梁的高度，单位：米

杨氏模量，单位：帕斯卡（）

E=200e9#Pa

nu=0.3#泊松比

#定义网格

n_elements=10#单元数量

n_nodes=n_elements+1#节点数量

dx=length/n_elements#单元长度

#定义节点坐标

2

x=np.linspace(0,length,n_nodes)

#定义单元刚度矩阵

k=(E*width*height)/dx#单元刚度

K=diags([k,-k],[0