结构力学数值方法：有限元法(FEM)：接触问题的有限元分析.pdf

结构力学数值方法：有限元法(FEM)：接触问题的有限元分

析

1绪论

1.1有限元法的历史和发展

有限元法（FiniteElementMethod,FEM）起源于20世纪40年代末，最初

由工程师们在解决结构力学问题时提出。1943年，R.Courant在解决弹性力学

问题时首次使用了有限元的概念。然而，直到1956年，当O.C.Zienkiewicz和

Y.K.Cheung在《工程计算》杂志上发表了一篇关于有限元法的文章后，这一方

法才开始在工程界广泛传播。自那时起，FEM迅速发展，成为解决复杂工程问

题的强有力工具，其应用领域从最初的结构力学扩展到流体力学、热力学、电

磁学等多个领域。

1.2接触问题在工程中的重要性

接触问题是指两个或多个物体在接触面上的相互作用，这种作用可能包括

摩擦、磨损、粘附等现象。在工程设计和分析中，接触问题的准确模拟对于预

测结构的性能、优化设计、减少故障和提高安全性至关重要。例如，在机械设

计中，齿轮、轴承、连接件等部件的接触分析可以帮助工程师理解并控制接触

应力，避免过早的磨损和失效。在土木工程中，地基与结构之间的接触分析对

于评估结构的稳定性和承载能力同样重要。

1.3接触问题的分类

接触问题根据接触面的性质和接触条件的不同，可以分为以下几类：

1.点接触：两个物体在单点接触，常见于尖锐物体与平面的接触。

2.线接触：接触发生在一条线上，如齿轮的齿面接触。

3.面接触：接触发生在两个物体的面上，如压板与工件的接触。

4.滑动接触：接触面之间存在相对滑动，涉及摩擦力的计算。

5.粘附接触：接触面之间存在粘附力，常见于微纳米尺度的接触问

题。

6.分离接触：接触面之间可以分离，涉及间隙和接触状态的判断。

每种类型的接触问题都有其特定的分析方法和挑战，FEM通过建立接触面

的离散模型，可以有效地模拟这些复杂现象。

1

1.3.1示例：使用Python和FEniCS进行面接触问题的有限元分析

假设我们有两个弹性体，一个为圆柱形压头，另一个为平面工件，我们想

要分析压头在工件上施加压力时的接触应力分布。以下是一个简化版的Python

代码示例，使用FEniCS库进行有限元分析。

fromfenicsimport*

importmatplotlib.pyplotasplt

#创建网格

mesh=RectangleMesh(Point(0,0),Point(1,1),10,10)

#定义函数空间

V=VectorFunctionSpace(mesh,Lagrange,2)

#定义边界条件

defboundary(x,on_boundary):

returnon_boundary

bc=DirichletBC(V,Constant((0,0)),boundary)

#定义变分问题

u=TrialFunction(V)

v=TestFunction(V)

f=Constant((0,-1))#压力

g=Constant((0,0))#边界力

#弹性体的材料属性

E=1e3#弹性模量

nu=0.3#泊松比

mu=E/(2*(1+nu))

lmbda=E*nu/((1+nu)*(1-2*nu))

#应力应变关系

defsigma(u):

returnlmbda*tr(eps(u))*Identity(2)+2.0*mu*eps(u)

#应变位移关系

defeps(u):

returnsym(nabla_grad(u))

#定义变分形式

a=inner(sigma(u),eps(v))*dx

L=inner(f,v)*dx+inner(g,v)*ds

2

#求解变分问题

u=Function(V)

solve(a==L,u,bc)

#可视化结果

plot(u)

plt.show()

1.3.2代码解释

1.创建网格：我们使用Rectangl