2D线性弹性模拟器 返回
材料力学·有限元法

2D线性弹性模拟器 — 平面应力与平面应变

对比「平面应力」与「平面应变」两种2维简化方法的工具。输入面内应力σx、σy、τxy,从胡克定律实时求出应变分量、面外分量、von Mises应力,观察薄板与厚物的结果差异。

参数设置
杨氏模量 E
GPa
材料的弹性模量。钢约200GPa
泊松比 ν
横向收缩程度。钢约0.3
应力 σx
MPa
x方向法向应力(正为拉伸)
应力 σy
MPa
y方向法向应力(正为拉伸)
剪切应力 τxy
MPa
面内剪切应力
分析模式
切换面外方向的处理方式
计算结果
法向应变 εx (×10⁻³)
法向应变 εy (×10⁻³)
剪切应变 γxy (×10⁻³)
面外分量
von Mises应力 (MPa)
剪切弹性模量 G (GPa)
受应力材料单元 — 变形动画

方形材料单元受σx、σy、τxy作用。虚线为原始形状,填充为夸大后的变形形状(颜色表示von Mises应力大小)。

莫尔圆 — 面内应力状态
应变分量(×10⁻³)
理论·主要公式

平面应力(σz=0、薄板)的胡克定律:

$$\varepsilon_x=\frac{\sigma_x-\nu\sigma_y}{E},\quad \varepsilon_y=\frac{\sigma_y-\nu\sigma_x}{E},\quad \gamma_{xy}=\frac{\tau_{xy}}{G}$$

σx、σy:面内法向应力,τxy:剪切应力,E:杨氏模量,G:剪切弹性模量。平面应力下产生面外应变εz=-ν(σx+σy)/E。

平面应变(εz=0、厚物)的面外应力:

$$\varepsilon_z=0,\quad \sigma_z=\nu(\sigma_x+\sigma_y)$$

面外应变被约束为零,但作为约束反力产生面外应力σz。剪切弹性模量在两种模式下通用:G=E/(2(1+ν))。

2D线性弹性介绍

🙋
材料力学里出现「平面应力」和「平面应变」,听起来很相似,搞不清区别。到底有什么不同?
🎓
两者都是「把本来3维的问题简化为2维来求解」的假设。关键区别在于「如何处理厚度方向(面外、z方向)」。平面应力是「把z方向应力设为零」,平面应变是「把z方向应变设为零」。简单说,薄物用平面应力,厚物用平面应变。
🙋
为什么要根据薄厚来区分呢?
🎓
薄板在面内受拉时,因为厚度方向很薄,z方向可以自由地缩小。所以z方向应力几乎为零=平面应力。反过来,大坝的截面那么大,z方向沿着深度延伸很远,周围的材料死死地按住中央的单元,z方向无法伸缩。所以z方向应变为零=平面应变,这就是区别。
🙋
明白了。那右上角的「分析模式」切换后,同样的σx、σy却得到不同结果,就是因为这个原因?
🎓
正是。平面应力中σz=0,但泊松效应让板在厚度方向自动缩小,所以εz=-ν(σx+σy)/E不为零。平面应变中反过来,εz=0被强制约束,反作用力产生面外应力σz=ν(σx+σy)。σz有没有,由此连von Mises应力都会变。切换模式,看「面外分量」卡片和von Mises数值的对比吧。
🙋
von Mises应力是用来判断屈服的指标吧。2D情况下公式会变吗?
🎓
公式形式会变。平面应力因为σz=0,简化为√(σx²-σxσy+σy²+3τxy²)这么清爽的式子。平面应变的σz不为零,得用3轴通用公式√(0.5((σx-σy)²+(σy-σz)²+(σz-σx)²)+3τxy²)。实务中常说「这零件是薄的,用平面应力2D模型算」,但假设错了就应力见鬼。所以一开始最好两个都算算,摸清底细。
🙋
所以在做FEM的2D分析之前,用这个工具先预演一遍,才能有把握。
🎓
完全同意。FEM的2D单元库里有「平面应力单元」和「平面应变单元」两种,选错了就是解错问题。提前用这个工具知道σx、σy、τxy时对应的应变和von Mises值,要是FEM结果跟预期差十倍,一下子就能反应过来「哦,单元类型选错了」。下面的莫尔圆也是检查面内主应力、最大剪切应力的经典工具。

常见问题

两者都是将3维弹性问题简化为2维的假设。平面应力假设面外方向应力为零(σz=0),适用于厚度方向自由变形的薄板。面外应变εz=-ν(σx+σy)/E不为零。平面应变假设面外方向应变为零(εz=0),适用于厚度方向变形受约束的大坝或长尺挤出材料。此时产生面外应力σz=ν(σx+σy)。相同的面内应力在两种模式下计算的应变与von Mises应力会有差异。
平面应力中:εx=(σx-νσy)/E,εy=(σy-νσx)/E,γxy=τxy/G。其中G是剪切弹性模量,G=E/(2(1+ν))。平面应变中考虑面外应力的贡献:εx=((1-ν²)σx-ν(1+ν)σy)/E,εy=((1-ν²)σy-ν(1+ν)σx)/E。剪切应变γxy=τxy/G在两种模式下相同。本工具根据选择的模式切换公式,同时显示三个应变分量与面外分量。
von Mises相当应力是多轴应力状态的单一标量指标,用于屈服判定。平面应力(σz=0)中:σ_vm=√(σx²-σxσy+σy²+3τxy²)。平面应变中σz=ν(σx+σy)不为零,使用通用公式:σ_vm=√(0.5((σx-σy)²+(σy-σz)²+(σz-σx)²)+3τxy²)。相同的面内应力在两种模式下的von Mises应力因σz的存在而改变,屈服倾向也会不同。
根据物体厚度与面内尺寸的比例判断。厚度远小于面内尺寸且面外方向自由变形时选择平面应力,如齿轮齿、薄板金部件、受面内荷载的面板等。厚度大且面外变形受周围约束时选择平面应变,如大坝截面、隧道、长挤出材料中部截面等。对于中等厚度可计算两种模式并比较,采用安全侧结果。

实际应用

薄板·板金部件设计(平面应力):受面内荷载的薄金属板、支架、齿轮齿轮廓、连杆机构的板件等都用平面应力求解。板厚相对于面内尺寸很小,面外方向应力可以忽略。穿孔板的应力集中、切口周边应力分布等是机械设计中最常见的领域,FEM也以平面应力单元为标准。

土建·地基结构物(平面应变):混凝土坝截面、隧道周边地基、挡土墙、长堤等厚度方向长期保持相同截面的结构物用平面应变分析。截面前后延伸很远,周围材料紧紧约束中央单元,导致面外变形受限。这时产生的面外应力σz如果忽视,应力评估会严重偏低。

长尺机械部件(平面应变):长管中央截面、挤出型材、旋转轴的中部、厚壁圆筒截面等也都是平面应变的典型。远离两端的中央区域,截面沿轴向一致变形,使得εz≈0成立。厚壁压力容器的应力分析也常从中央截面的平面应变假设出发。

CAE建模的模式选择检查:FEM建2D模型时最先要定的是「平面应力单元还是平面应变单元」。本工具能对单点的应力状态瞬间比较两种模式下的应变与von Mises,直观验证建模合理性。详细计算前用手算初步把握,若FEM结果相差十倍就立即怀疑单元类型选错,可作为可靠性检验工具。

常见误解与注意事项

最大误解是「平面应力则面外应变也为零,平面应变则面外应力也为零」的错误印象。实际恰恰相反:平面应力中σz=0但面外应变εz=-ν(σx+σy)/E不为零(板在厚度方向会缩小);平面应变中εz=0被强制约束,但作为反作用力产生面外应力σz=ν(σx+σy)。混淆「应力」与「应变」到底哪个为零,会导致完全曲解计算结果。本工具的「面外分量」卡片会根据选定模式切换显示εz还是σz,帮助理解。

其次「面内应力相同则von Mises应力也相同」的误判。平面应力和平面应变中σx、σy、τxy完全一致,但von Mises应力却不同。因为平面应变中面外应力σz=ν(σx+σy)非零,3轴通用公式算出的值与平面应力的值会偏离。特别当σx、σy同号且数值大时,σz也大,差异就明显。屈服判定必须用对应模式的von Mises应力。

最后要记住线性弹性本身的应用限界。本工具基于胡克定律,仅适用线性弹性范围。当应力超过材料的屈服应力(钢约200~400MPa)时,材料进入塑性变形,应力与应变不再成比例。von Mises应力超过屈服值的区域里,这里计算的应变会比实际值偏小。应变超过百分之几时的大变形、座屈、接触等问题也超出线性弹性范围。线性弹性仅在「小变形、屈服前」的前提下成立,时刻不能忘记这一点。

使用指南

  1. 输入材料的杨氏模量E(GPa)和泊松比ν(0~0.5)。例如软钢:E=200GPa、ν=0.3
  2. 设置X方向应力σx、Y方向应力σy,在平面应力(σz=0)或平面应变(εz=0)条件下选择
  3. 模拟器自动计算法向应变εx、εy、剪切应变γxy、von Mises应力,展示平面应力与平面应变的结果差

具体计算示例

厚度3mm的铝合金薄板(E=70GPa、ν=0.33)受到σx=150MPa、σy=50MPa的二轴应力。在平面应力条件下εx≈2.14×10⁻³、εy≈-0.36×10⁻³;同样荷载下厚壁铸铁(E=100GPa、ν=0.28)在平面应变条件受约束时,面外应力σz产生,εx≈1.43×10⁻³下降。两种条件的von Mises应力分别约155MPa与160MPa,在应力集中和疲劳评估中会导致不同判定结果

实务注意事项