浅谈“塑性力学”教学中的Lode应力参数拓展

作者: 王宏伟 付佳卉

浅谈“塑性力学”教学中的Lode应力参数拓展0

[摘 要] 从课堂教学的角度分析了当前塑性力学教材中Lode应力参数的局限性,给出了其计算公式在π平面内的拓展形式。结果显示,拓展后的Lode应力参数有6种不同的计算公式,传统Lode应力参数计算公式只是主应力次序是σ1≥σ2≥σ3时的一种情况。当主应力次序不再为σ1≥σ2≥σ3时,Lode应力参数的计算公式随之发生改变。同时,由于Lode应力参数是建立材料屈服条件的基础,所以采用研究结果可以在不确定主应力大小次序的情况下绘制屈服曲线。所探讨的结果拓展了Lode应力参数的计算公式形式,完善了π平面内屈服曲线的绘制途径。

[关键词] 塑性力学;Lode应力参数;π平面;主应力

[基金项目] 2021年度中国矿业大学(北京)本科教育教学改革与研究项目(J210614);2020年度中国矿业大学(北京)本科教育教学改革与研究项目(J200706);2018年度中国矿业大学(北京)“越崎青年学者”(2018QN13)

[作者简介] 王宏伟(1981—),男,山西朔州人,博士,中国矿业大学(北京)力学与建筑工程学院教授,博士生导师,主要从事工程力学、岩土工程和采矿工程等方面的理论及教学研究;付佳卉(1997—),女,辽宁大连人,中国矿业大学(北京)力学与建筑工程学院2019级工程力学专业硕士研究生,研究方向为岩石强度理论。

[中图分类号] TP028.8 [文献标识码] A [文章编号] 1674-9324(2022)12-0173-04 [收稿日期] 2021-07-26

在“塑性力学”理论教学中,主应力空间中π平面内的应力分解会引出一组重要的参数,即Lode参数,包括Lode应力参数μσ、Lode应力角θσ和偏应力矢量的模rσ[1,2]。这组参数是建立材料屈服条件的基础,借助这组参数可以直观地绘制屈服曲线和屈服面[3]。

当前塑性力学教材中,Lode参数推导过程是在π平面内建立笛卡尔正交坐标系,并将坐标系的y轴和主应力σ2轴的投影轴σ2'轴的正向重合,而且主应力顺序是σ1≥σ2≥σ3,由此才可得到教材中的Lode应力参数μσ和Lode应力角θσ[4-7]。但是学生在学习过程中总会提出:为什么π平面内的正交坐标系y轴要和σ2'轴正向重合;如果y轴和其他投影轴重合将得到什么不同的结果和物理意义;尽管σ1≥σ2≥σ3是常采用的主应力次序,但是在一般情况下,往往无法事先判断材料内各点的三个主应力大小次序[8],那么如果主应力不按此排序,Lode参数有无变化。

为了解开这些疑问,笔者尝试将π平面内的正交坐标系的y轴分别和主应力轴在π平面内的投影轴的正向和负向分别重合,考察将得到何种结果。首先简单复盘Lode参数的推导过程,并且分析结果的局限性。

一、Lode参数的推导过程

(一)Lode应力角和Lode应力参数

(二)现有结果的局限性

当前的塑性力学教材中,如果主应力大小次序按照σ1≥σ2≥σ3排序,则可以由公式(1)和(2)确定Lode应力角和Lode应力参数的取值范围,即-30°≤θσ≤30°,-1≤μσ≤1,也就是图2(a)所示的π平面内60°范围内的偏应力状态变化。然而,当主应力大小次序不按照上述排序时,公式(1)和(2)就无法确定Lode应力角和Lode应力参数的取值范围,π平面内除上述60°的其他范围内的偏应力状态也无法回答。这导致学生不能充分理解Lode应力参数在π平面上的物理意义。所以,需要进一步推导Lode参数的其他计算公式,解释其物理意义。笔者接下来将π平面内的y轴分别沿着σ1、σ2和σ3轴的正向和负向,分析得到结果。

二、Lode参数的拓展

(一)Lode应力参数的其他形式

在教学过程中,为了让学生能充分理解Lode应力参数,笔者将π平面内的正交坐标系的y轴分别与重合于σ1',σ2'和σ3'轴的正向和负向[8]。

与前述推导相比,变化的只是矢量OP'在坐标系x轴和y轴上的分量,因此其他情况下的Lode参数推导过程略去。表1给出了当y轴沿着π平面内各投影轴时的Lode应力角和Lode应力参数的结果。结果显示,Lode应力参数有6个不同的公式,而公式(2)只是y轴沿着σ2'轴的正向时的一种情况。这6种不同的情况恰好可以将Lode应力角[-30°,30°]的范围旋转360°后覆盖了整个π平面,如图2(b)所示。

(二)主应力和主偏应力大小次序的几种情况

图3是塑性力学教材中讲解到Tresca屈服曲线时给出的主应力大小次序分布情况[1],但是教材中并未给出明确解释,下面将根据Lode应力参数的拓展公式来解释这一问题。

由于Lode应力参数的取值范围为-1≤μσ≤1,为了满足这个条件,不同的Lode应力参数计算公式的主应力的大小次序是不同的,如表1所示。例如在图2(a)中的简单应力状态下,从单向拉伸状态A点,到纯剪切状态B点,再到单向压缩状态C点所覆盖的范围,即为对应的Lode应力参数所应该满足的主应力取值范围,即一般规定的σ1≥σ2≥σ3的分布情况。另外,由于主偏应力的方向和主应力方向重合,因此主偏应力大小次序S1≥S2≥S3也在这个区域内。为了清晰的给出主应力和主偏应力的分布情况,图3绘制了6种主应力和主偏应力的大小次序对应的分布区域。

由图3可知,当前塑性力学教材中主应力大小次序σ1≥σ2≥σ3只是图中y轴沿着σ2'轴的正向时的一种情况。当Lode应力参数为其他公式时,屈服曲线位于其他范围内,主应力和主偏应力的大小次序应该是图3中的各种情况。这里需要说明的是,虽然在π平面内主要应该讨论的是偏应力状态,但是图3也给出了主应力的分布情况。这是因为Lode应力参数的引出主要为了建立屈服条件或者绘制屈服曲线,而屈服条件的建立还是应该体现主应力的大小和方向。这也是塑性力学教材中在讲解到π平面内Tresca和Mises屈服曲线时,不同范围内主应力大小次序不一致的原因,但是教材中并未给出解释。

三、主应力次序变换时的简单应力状态和屈服曲线

Lode参数是建立材料屈服条件的基础,一般地,屈服曲线的绘制方法是通过简单应力状态确定若干在π平面内的特殊点,然后通过屈服曲线的对称性进行绘制。所谓简单应力状态,是指施加的外力使材料只发生一种基本变形,一般为单向拉伸、单向压缩和纯剪切状态。

因为图2(a)中的简单应力状态对应的主应力次序是σ1≥σ2≥σ3,因此现有教材是以这种情况来计算π平面内的特殊点的。那么其他主应力次序下如何确定屈服曲线,本文的Lode参数拓展可以解决这一问题。

由此可见,当主应力次序不同时,Lode应力参数的计算公式改变,则简单应力状态所对应的主应力也不同。当主应力大小次序按照表1中6种不同情况排序时,均可得到简单应力状态在π平面内的位置,这也可以帮助学生更加清晰直观地理解整个π平面的物理意义,如图2(b)所示。虽然改变主应力次序得到的Lode应力参数的形式不同,但是各个参数所体现的物理意义是相通的。因此,如果不确定主应力大小次序的情况下,依然可以按照图2(b)来确定π平面内简单应力状态对应的点,进而通过对称性绘制屈服曲线,关于这一点现有的塑性力学教材是没有涉及的。

结语

“塑性力学”是固体力学的重要分支,是力学学科中理论性较强的一门专业核心课程,该理论重点阐述了材料在外力作用下的屈服条件和本构关系。而Lode参数是建立材料屈服条件和本构关系的基础,笔者通过在π平面内将主应力大小次序拓展后发现,Lode应力参数计算公式存在6种不同的情况,采用这6种Lode参数结果可以在不确定主应力大小次序的情况下绘制屈服曲线。本文所探讨的结果拓展了Lode应力参数的计算公式形式,完善了π平面内屈服曲线的绘制途径,在教学中能够展现出理论的全面性,同时也锻炼了学生的思维拓展能力。

参考文献

[1]徐远杰,楚锡华.中主应力对Mohr-Coulomb材料屈服性质的影响[J].力学与实践,2009,31(4):76-77.

[2]刘家信.关于塑性力学中的Lode参数[J].东北重型机械学院学报,1987(1):70-78.

[3]赵长财,于忠海,袁荣娟,等.关于应变Lode参数的研究[J].塑性工程学报,1994(4):19-24.

[4]夏志皋.塑性力学[M].上海:同济大学出版社,1991:7.

[5]尚福林.塑性力学基础[M].3版.西安:西安交通大学出版社,2018:4.

[6]徐秉业.简明弹塑性力学[M].北京:高等教育出版社,2011:1.

[7]宋卫东.塑性力学[M].北京:科学出版社,2017:5.

[8]王宏伟,祝捷,刘德军.塑性力学[M].徐州:中国矿业大学出版社,2021:10.

Discussion on the Expansion of Lode Stress Parameters in the Teaching of Plasticity Course

WANG Hong-wei, FU Jia-hui

(School of Mechanics and Civil Engineering, China University of Mining and Technology-Beijing, Beijing 100083, China)

Abstract: The limitation of Lode stress parameters in the current plasticity textbook is analyzed from the perspective of classroom teaching. The expanded form of Lode stress parameter calculation formula in π plane is deduced in this paper. It is suggested that there are six formulas for the calculation of Lode stress parameter. The formula in the current plasticity textbook is one of the six cases when the principal stress sequence is . The form of Lode stress parameter formula could change according to the principal stress sequence. Meanwhile, since the Lode stress parameter is the basis for establishing the material yield criteria in the plasticity, the results of this paper can be used to draw the yield curve without determining the sequence of principal stress. Therefore, the results discussed in this paper can develop the formulas form of Lode stress parameter and improve the way of drawing the yield curve in π plane.

Key words: Plasticity; Lode stress parameters; π plane; principal stress

经典小说推荐

杂志订阅