弹塑性弯曲直梁的回弹变分原理及应用.docxVIP

弹塑性弯曲直梁的回弹变分原理及应用

陈英杰，郭清超，周梦飞，张华

(燕山大学河北省土木工程绿色建造与智能运维重点实验室,秦皇岛066004；燕山大学建筑工程与力学学院,秦皇岛066004)

关于物体变形问题的研究由来已久，而应用变分原理研究物体的变形经过了漫长的时间，经研究发现，理想弹塑性材料的耦联系统塑性成形过程可以视为实际系统的变形，对应的反耦联系统的变形可以作为实际系统的回弹变形。因此，可以应用反耦联系统求解出弹塑性梁对应的回弹挠曲线方程。当前工程中涉及梁的变形及回弹问题还有很多是根据弹性理论或经验公式解决的，存在一定的误差，也在一定程度上造成材料浪费，为了限制弹塑性梁的变形，需要建立能够指导实际生产应用的理论。弹塑性梁的振动问题[1-2]和弯曲问题[3-4]一直被广大学者所关注，由于材料本构关系的原因，荷载作用下的梁会产生弯曲变形，从而截面会产生弹性区和塑性区，如何划分弹塑性区和弹塑性梁的变形是研究热点之一。受力构件卸载后，结构产生了反向变形，由于材料的回弹值很小，因此需要研究出有效、精确的方法解决回弹问题。中国学者对回弹问

题开展了大量研究。吴斌等[5]基于U形板热冲压构件的残余应力与回弹影响之

间的关系进行了研究；罗之军[6]建立了电子设备罩壳弹塑性材料有限元模型，对结构的回弹现象进行了研究；钱伟长研究发现，由于材料本构关系的原因，荷载作用下的梁会发生弯曲变形，从而会产生弹性区和塑性区，而弹性区和塑性区的分界面并不能事先给定，需要找出划分弹塑性区的分界线；前人采用边界待定的变分法[8-10]解决了塑性力学中理想弹塑性材料变形体弹性区和塑性区的划分问题，以及理想弹塑性材料的梁在不同条件下的变形及回弹问题，通过求解可以得到该类问题的精确解。对弹塑性直梁弯曲变形和回弹变形的研究具有重要的理论意义和工程应用价值。为了方便叙述，弯曲直梁变形均为有限变形。

弯曲直梁的反耦联系统和反耦联方程

首先引入弯曲直梁的反耦联系统和反耦联方程的概念，考虑材料、形状和尺寸相同，位移边界条件相同，但体力及边界力相反的两直梁系统。其中之一为成形直梁，该直梁在加载时发生弹塑性变形。另一直梁成为回弹直梁，该直梁在加载时只发生弹性变形。成形直梁和回弹直梁构成弯曲直梁-反耦联系统。成形直梁如图1所示。

成形直梁的平衡方程、静力边界条件和位移边界条件分别表示为

(1)

(2)

式中：M为成形直梁弯矩；Mx=0、Mx=l分别为成形直梁的梁左、右端弯矩；

Wx=0、Wx=l分别为成形直梁的梁左、右端挠度。

(3)

回弹直梁和相应的力学方程分别表示为

x为成形直梁从原点o到直梁上任意一点的距离；为梁左、右端弯矩；表示梁左、右端挠度；q为均布外荷载；l为跨度图1成形直梁Fig1Straightbeamforming

(4)

(5)

(6)

式中：M′为回弹直梁弯矩；M′x=0、M′x=l分别为回弹直梁的梁左、右端弯矩；W′x=0、W′x=l分别为回弹直梁的梁左、右端挠度。

式(1)～式(3)分别与式(4)～式(6)合并，则得方程为

(7)

(8)

(9)

式中：式(7)为成形体与回弹体构成的反耦联系统的反耦联平衡方程；式(8)为反耦联静力边界条件；式(9)为反耦联位移边界条件。

弯曲直梁的回弹变分原理

弯曲直梁的回弹最小势能原理

采用加权余量法对式(7)、式(8)进行计算，则有：

(10)

对式进行相应的变分运算。

(11)

则得

(12)

式中：Π′sp为回弹直梁的弹性势能；A(K′)为回弹势能密度；式中：δ为变分极值；w′为回弹直梁的容许挠度；Me为梁水平方向弯矩κe为弹塑性区与弹性区交界处的转角；κ1为弹塑性区的转角；ξ为梁跨内弹塑性区长度；E为杨氏模量；J为惯性矩；A(K′)为回弹势能密度；K′为回弹直梁弯曲曲率。

对式(13)进行变分运算。

(13)

则得

(14)

根据变分法基本预备定理，则得欧拉方程为

(15)

(16)

(17)

弯曲直梁的回弹最小余能原理

应用加权余量法，可得回弹总余能Π′sc为

(18)

因

(19)

即

(20)

式中：B(M′)为回弹余能密度；K为成形直梁弯曲曲率；w为成形直梁的容许挠度。

(21)

由于内力是容许的，故应满足平衡方程和静力边界条件：

(22)

δM′x=0=δM′x=l=0

(23)

于是式(21)变为

(24)

据变分法基本预备定理，则得欧拉方程为

(25)

w′x=0-wx=0=0

(26)

w′x=l-wx=l=0

(27)

式中：wx=0、wx=l分别为成形直梁的梁左、右端容许挠度；w′x=0、w′x=l分别为回弹直梁的梁左、右端容许挠度。

弯曲直梁的广义回弹变分原理

弯曲直梁的广义回弹势能原理

应用拉格朗