用ABAQUS玩转液压油缸模拟：基于CEL算法的加载模型

xEurCjvwu

234人浏览 · 2026-03-25 10:00:00

xEurCjvwu · 2026-03-25 10:00:00 发布

ABAQUS有限元模型：基于CEL算法的液压油缸加载模型。使用ABAQUS有限元软件，基于CEL算法，模拟了液压油缸在荷载作用下，结构的受力和内部液体压强变化，其中油缸采用拉格朗日体，内部液体使用欧拉体，分析了加压过程的分析状态，效果较好。包括模型文件和教程。

最近在搞液压油缸的有限元模拟，想看看用ABAQUS能不能把这个过程模拟得逼真一些。于是，我决定尝试一下基于连续欧拉朗格朗日（CEL）算法的模型。说实话，刚开始接触 CEL 算法的时候，感觉挺复杂的，但一步步搞下来，发现其实也还好，效果还不错。

一、为什么选择CEL算法？

在液压油缸的模拟中，液体和固体的相互作用是一个关键点。传统的拉格朗日法可能会因为网格畸变导致计算困难，而欧拉法在处理大变形时又可能会出现精度问题。于是，CEL算法来了，它把固体用拉格朗日网格表示，液体用欧拉网格表示，这样两全其美。液体的运动由欧拉网格计算，而固体的变形由拉格朗日网格处理，两者通过接触算法相互作用。

听起来很高大上，对吧？不过实际操作起来，配置还是得仔细点。

二、模型搭讪：从结构到液体

首先是模型的建立。液压油缸的主体结构是一个金属圆筒，活塞在内部运动。按照CEL算法的要求，我需要把金属圆筒和活塞作为拉格朗日体，而内部的液压油则作为欧拉体。

具体来说，模型大致这样：

油缸主体：一个圆柱体，金属材料，用四面体单元（C3D4）划分网格。
活塞：一个圆盘，同样用四面体单元划分，固定在油缸的一端。
液压油：填充在油缸内部，用欧拉体（FLUENT流体单元）划分网格。

这样配置的好处是，金属部分的网格保持不变形，而液体则可以通过欧拉网格流动。两者之间的接触面需要定义接触对，这样液体才能带动活塞运动。

三、代码部分：配置接触和边界条件

说到配置，代码部分当然少不了。这里我写了一个简单的ABAQUS脚本，用于定义模型的基本设置。

# 模型建立
model = mdb.Model(name='HydraulicCylinderCEL')
part_cylinder = model.Part(name='Cylinder', dimensionality=THREE_D)
# 创建油缸主体几何
# ...创建圆柱体和活塞几何...

# 定义材料
model.Material(name='Steel')
steel = model.materials['Steel']
steel.Elastic(table=((200e9, 0.3),))
# 定义液压油材料
model.Material(name='Oil')
oil = model.materials['Oil']
oil.Fluid(table=((800, 1.0, 0.0),))  # 流体密度、粘度、膨胀系数

# 网格划分
part_cylinder.setElementType(elements=C3D4)
part_cylinder.generateMesh()

# 定义接触属性
contact_prop = model.ContactProperty('ContactProp')
contact_prop.ContactControls(tangential_stiffness=1e5)
contact_pair = model.ContactPair(name='FluidStructure',
                                  surface_1=part_cylinder.faces.findAt((0,0,0)),
                                  surface_2=part_cylinder.faces.findAt((0,0,1)))

# 边界加载
# 对活塞施加位移载荷
displacement = 0.1  # 单位：米
step = model.steps['InitialStep']
step.DisplacementBC(amplitude=displacement,
                     region=part_cylinder.sets['PistonSet'],
                     dof=3)  # 沿z轴位移

从这段代码可以看出，我定义了两种材料：金属和液压油。其中，金属部分用四面体单元，而液压油用流体单元。接触属性的定义也很重要，它决定了液体和固体之间的相互作用。

四、模拟过程：液体压强和结构应力

当活塞被施加位移载荷时，内部的液压油就会被压缩，导致压强变化。通过模拟，可以实时观察到液体压强的变化以及结构的应力分布。

有意思的是，在模拟中，液体的流动和固体的变形是同时进行的。通过ABAQUS的后处理，可以看到活塞运动过程中，液体在油缸内的流动情况，以及油缸壁的应力分布。

五、结果可视化：压力场和位移云图

模拟完成后，导出的结果让人眼前一亮。压力场的分布清晰可见，活塞运动引起的液体压强变化也一目了然。同时，结构的位移云图显示了油缸壁的变形程度。

# 后处理
odb = session.openOdb('HydraulicCylinderCEL.odb')
# 显示压力场
pressure = odb Field: pressure
session.viewports['Viewport: 1'].odbDisplay.setPrimaryVariable(variable=pressure)
# 显示活塞位移
displacement = odb Field: displacement
session.viewports['Viewport: 2'].odbDisplay.setPrimaryVariable(variable=displacement)

从结果来看，活塞运动过程中，液体的压强变化非常剧烈，尤其是在活塞快速移动的时候。同时，油缸壁的应力主要集中在活塞接触区域，这和实际应用中的情况非常吻合。