力学与实践, 2021, 43(1): 139-143 DOI: 10.6052/1000-0879-20-181

教育研究

弹性力学应力状态理论求解器开发及应用研究1)

何峰,2), 任天娇, 杨松, 赵娜

辽宁工程技术大学力学与工程学院,辽宁阜新 123000

DEVELOPMENT AND APPLICATION OF STRESS STATE THEORY SOLVER FOR ELASTICITY 1)

HE Feng,2), REN Tianjiao, YANG Song, ZHAO Na

Department of Mechanics and Engineering, Liaoning Technical University, Fuxin 123000, Liaoning, China

通讯作者: 2)何峰,副教授,从事工程力学及岩土工程等方面的教学和研究工作。E-mail:12041827@qq.com

责任编辑: 胡漫

收稿日期: 2020-05-6   修回日期: 2020-07-22   网络出版日期: 2021-02-08

基金资助: 1)国家自然科学基金.  51774167
校金课项目.  2019030

Received: 2020-05-6   Revised: 2020-07-22   Online: 2021-02-08

作者简介 About authors

摘要

论文基于吴家龙教授编著的《弹性力学》教材内容和课程"三难"特点,以应力状态理论部分常见题型求解任意微分斜面应力矢量、正应力、切应力、主应力和应力主方向及最大切应力为例,应用Matlab-GUI模块编制应力状态理论求解器,通过界面物理参数输入和结果显示过程,再现弹性力学求解思维,验证学生数理求解结论;把繁琐抽象的力学公式和微分斜面结果进行图形可视化,使学生的抽象思维和力学思维得到兼容;通过弹性力学应力状态理论求解器开发及应用,可为学生力学学习及力学同行教学提供借鉴。

关键词: 弹性力学 ; 应力状态求解器 ; 微分斜面 ; 可视化

Abstract

Based on the textbook Elasticity by Wu Jialong, and the "trilemma" characteristics of the course, taking the common questions in the stress state theory as examples, the stress vector, the normal stress, the shear stress, the principal stress, the principal direction of the applied force and the maximum shear stress on any differential inclined plane are included. The Matlab-GUI module is used to compile the stress state theory solver, By processing the inputting physical parameters and displaying results, the thinking process of the elasticity is reproduced and the mathematical solution is verified, the complex abstract mechanical formulas and the differential slope results are clarified in graphic visualization, so that the students' abstract thinking and mechanical thinking are linked. The stress state theory solver of elasticity can provide a reference for students' mechanics learning and for mechanics teaching.

Keywords: elasticity ; stress state solver ; differential slope ; visualization

PDF (6640KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

何峰, 任天娇, 杨松, 赵娜. 弹性力学应力状态理论求解器开发及应用研究1). 力学与实践[J], 2021, 43(1): 139-143 DOI:10.6052/1000-0879-20-181

HE Feng, REN Tianjiao, YANG Song, ZHAO Na. DEVELOPMENT AND APPLICATION OF STRESS STATE THEORY SOLVER FOR ELASTICITY 1). MECHANICS IN ENGINEERING[J], 2021, 43(1): 139-143 DOI:10.6052/1000-0879-20-181

弹性力学课程是普通高校力学类、土木类等理工科本科专业核心课,是从力学基本理论过渡到工程实际应用的桥梁;具有教师难教、学生难学、工程思维和力学思维难兼容的"三难"课程特点。基于课程重要性和特点,众多学者进行了弹性力学教学改革,如张伟伟等[1]提出了弹性力学三段式教学方法, 即围绕每一知识点, 按照工程背景、数理基础、力学原理进行划分, 在教学实施中遵循先工程, 后数理,再力学的讲解顺序。潘东辉等[2]将MATLAB/PDE工具箱引入弹性力学教学,实现弹性力学的可视化教学。刘杨等[3]提出了"互联网$+$项目指派"式教学模式在《弹性力学及有限单元法》的教学中的实践方式和考核办法。陈小亮等[4]利用Maple软件探索了弹性力学应力函数逆解法的计算机求解规范流程。卢小雨等[5]利用Maple来求解弹性力学中的一些具体问题。邢静忠[6]将MAPLE引入力学教学,引导和培养学生利用数学工具的习惯和能力,强化算法设计和程序的通用性和灵活性,为处理复杂问题提供帮助。丁洲祥等[7]详细介绍了计算机代数系统MAPLE在土力学与基础工程研究型教学中的具体应用。

本文基于吴家龙教授编著的《弹性力学》教材内容和课程特点,以应力状态理论部分的求解斜面应力分量及主应力和最大切应力为例,应用Matlab-GUI模块,通过编制程序把繁琐抽象的力学公式结果进行图形可视化,从而提升学生学习弹性力学兴趣。

1 应力状态理论[8]

弹性力学的研究对象为三维弹性体,由于满足连续性和均匀性假设,因此从微分单元体分析入手,确定一点不同截面的应力变化规律称为应力状态分析。研究思路是首先确定应力状态描述方法,包括应力矢量定义及其分解为正应力、切应力和应力坐标轴分量,应力张量表示一点应力状态;其次从静力学观点出发,建立平衡微分方程和应力边界条件;任意截面的应力分量确定;应力分量转轴公式;一点的特殊应力确定:主应力和主平面、最大切应力。该部分教学中学生和老师常见问题有:如何直观体现已知方向余弦?如何确定应力主方向、最大切应力微分面方位或已知微分面方程求方向余弦,同时可视化微分斜面?

2 求解器运行原理

(1)微分面三元一次方程系数与其方向余弦关系

$F(x,y,z)=Ax+By+Cz+D $
$\left. {\begin{array}{l} l=\dfrac{A}{\sqrt {A^{2}+B^{2}+C^{2}} } \\ m=\dfrac{B}{\sqrt {A^{2}+B^{2}+C^{2}} } \\ n=\dfrac{C}{\sqrt {A^{2}+B^{2}+C^{2}} } \\ \end{array}} \right\} $
$l^{2}+m^{2}+n^{2}=1 $

式(1)$\sim$式(3)中,$F(x,y,z)$为三元一次函数; $A$, $B$,$C$,$D$为方程系数;$l$,$m$,$n$为方向余弦。

(2)斜面上应力分量(图1)

$\left. {\begin{array}{l} f_{vx} =l\sigma_{x} +m\tau_{xy} +n\tau_{zx} \\[1mm] f_{vy} =l\tau_{xy} +m\sigma_{y} +n\tau_{yz} \\[1mm] f_{vz} =l\tau_{zx} +m\tau_{yz} +n\sigma_{z} \\[1mm] \end{array}} \right\} $
$f_{v} =\sqrt {f_{vx}^{2}+f_{vy}^{2}+f_{vz}^{2}} $
$\sigma_{v} =lf_{vx} +mf_{vy} +nf_{vz} $
$\tau_{v}^{2}=f_{v}^{2}-\sigma_{v}^{2} $

式(4)$\sim$式(7)中,$f_{vx}$, $f_{vy}$, $f_{vz} $为沿坐标轴投影应力分量;$l,m,n$为方向余弦;$\sigma_{x}$,$\sigma_{y}$,$\sigma_{z}$,$\tau_{xy}$,$\tau_{yz}$,$\tau_{zx}$为正应力和切应力分量;$f_{v} $为斜面应力矢量;$\sigma_{v} $为斜面正应力;$\tau_{v} $为斜面切应力。

图1

图1   微分斜面上应力分量


(3)主应力、应力主方向和最大切应力及方位

$\sigma^{3}-I_{1} \sigma^{2}+I_{2} \sigma -I_{3} =0 $
$\left. {\begin{array}{l} I_{1} =\sigma_{x} +\sigma_{y} +\sigma_{z} \\ I_{2} =\sigma_{x} \sigma_{y} +\sigma_{x} \sigma_{z} +\sigma_{y} \sigma_{z} -\\\qquad \tau_{xy}^{2}-\tau_{yz}^{2}-\tau_{zx}^{2} \\ I_{3} =\sigma_{x} \sigma_{y} \sigma_{z} -\sigma_{x} \tau_{yz} ^{2}-\sigma_{y} \tau_{zx}^{2}-\\\qquad \sigma_{z} \tau_{xy}^{2}+2\tau _{xy} \tau_{yz} \tau_{zx} \\ \end{array}} \right\} $

式(8)$\sim$式(9)中,$I_{1}$, $I_{2}$, $I_{3} $为应力张量不变量;$\sigma_{x}$, $\sigma_{y}$, $\sigma_{z}$,$\tau_{xy}$, $\tau _{yz}$, $\tau_{zx} $为正应力和切应力分量;$\sigma $为主应力。

$\begin{eqnarray} \left. {\begin{array}{l} l_{1} \left( {\sigma_{x} -\sigma_{1} } \right)+m_{1} \tau_{xy} +n_{1} \tau_{zx} =0 \\ l_{1} \tau_{xy} +m_{1} \left( {\sigma_{y} -\sigma_{1} } \right)+n_{1} \tau_{yz} =0 \\ l_{1} \tau_{zx} +m_{1} \tau_{yz} +n_{1} \left( {\sigma_{z} -\sigma _{1} } \right)=0 \\ l_{1}^{2}+m_{1}^{2}+n_{1}^{2}=1 \\ \end{array}} \right\} \end{eqnarray} $

式(10)中,$l_{1}$, $m_{1}$, $n_{1} $为第一主应力方向余弦;$\sigma_{x}$, $\sigma_{y}$, $\sigma_{z}$, $\tau_{xy}$,$\tau_{yz}$, $\tau_{zx}$为正应力和切应力分量;$\sigma_{1} $为第一主应力。

$\begin{eqnarray} \left. {\begin{array}{l} l_{2} \left( {\sigma_{x} -\sigma_{2} } \right)+m_{2} \tau_{xy} +n_{2} \tau_{zx} =0 \\ l_{2} \tau_{xy} +m_{2} \left( {\sigma_{y} -\sigma_{2} } \right)+n_{2} \tau_{yz} =0 \\ l_{2} \tau_{zx} +m_{2} \tau_{yz} +n_{2} \left( {\sigma_{z} -\sigma _{2} } \right)=0 \\ l_{2}^{2}+m_{2}^{2}+n_{2}^{2}=1 \\ \end{array}} \right\} \end{eqnarray} $

式(11)中,$l_{2}$, $m_{2}$, $n_{2} $为第二主应力方向余弦;$\sigma_{x}$, $\sigma_{y}$, $\sigma_{z}$, $\tau_{xy}$, $\tau_{yz}$, $\tau_{zx}$为正应力和切应力分量;$\sigma_{2} $为第二主应力。

$\begin{eqnarray} \left. {\begin{array}{l} l_{3} \left( {\sigma_{x} -\sigma_{3} } \right)+m_{3} \tau_{xy} +n_{3} \tau_{zx} =0 \\ l_{3} \tau_{xy} +m_{3} \left( {\sigma_{y} -\sigma_{3} } \right)+n_{3} \tau_{yz} =0 \\ l_{3} \tau_{zx} +m_{3} \tau_{yz} +n_{3} \left( {\sigma_{z} -\sigma_{3} } \right)=0 \\ l_{3}^{2}+m_{3}^{2}+n_{3}^{2}=1 \\ \end{array}} \right\} \end{eqnarray} $

式(12)中, $l_{3}$,$m_{3}$,$n_{3} $为第二主应力方向余弦;$\sigma_{x}$, $\sigma_{y}$, $\sigma_{z}$,$\tau_{xy}$,$\tau_{yz}$,$ \tau_{zx}$为正应力和切应力分量;$\sigma_{3} $为第三主应力。

3 Matlab-GUI模块应力状态理论求解器及应用

编制基于Matlab-GUI模块应力状态理论求解器,程序流程如图2所示。通过该求解器(图3)可显示应力状态图和微分斜面,通常遇到题型有两种:一种题型为给出方向余弦可求出三元一次方程系数,进而显示微分面方程和斜面图,另一种题型若直接给出三元一次方程及系数也可求出方向余弦和显示微分斜面,同时可求出三个主应力及主方向和最大切应力及方向余弦,可通过主应力微分面按钮,如图4图5所示,分别显示主平面和最大切应力平面,验证主应力所处面切应力是为零,应力矢量和正应力相等,取得最大切应力微面,正应力不一定为零等结论。

图2

图2   应力状态理论求解器流程图


图3

图3   应力状态理论求解器


图4

图4   第一主应力微分斜面


图5

图5   最大切应力微分斜面


已知物体内一点的应力张量为

$\begin{eqnarray*} \left[ {{\begin{array}{*{20}c} {10} & 0 & {-10} \\ 0 & {-10} & {10} \\ {-10} & {10} & {20} \\ \end{array} }} \right]\mbox{MPa} \end{eqnarray*}$

求通过该点平面为$x+\sqrt 6 y+3z=1$的微面上沿$x$,$y$与$z$轴的应力分量、总应力、正应力、剪应力和主应力、应力主方向和最大切应力?

4 结论

论文结合弹性力学课程特点,力求突破"三难" 瓶颈;以弹性力学应力状态理论常见题型为}例,Matlab--GUI模块编制应力状态理论求解器,通过界面物理参数输入和结果显示过程,体现弹性力学中数理求解思路,简化了繁杂的计算过程;同时可形象直观再现和验证学生弹性力学数理计算结果,使得数学力学公式可视化,学生的抽象思维和具象思维达到兼容。最终达到提升学生学习弹性力学兴趣的目的。

参考文献

张伟伟, 田锦邦.

弹性力学的三段式教学方法

力学与实践, 2017,39(2):191-195

[本文引用: 1]

Zhang Weiwei, Tian Jinbang.

A three stages teaching method of elastic mechanics course

Mechanics in Engineering, 2017,39(2):191-195 (in Chinese)

[本文引用: 1]

潘东辉, 马崇武.

MATLAB/PDE在弹性力学可视化教学中的应用

力学与实践, 2014,36(4):500-504

DOI      URL     [本文引用: 1]

MATLAB 的PDE(partial differential equations)工具箱为平面问题求解结果提供了彩色可视化图形,便于弹性力学的学与教. 工具箱采用椭圆型偏微分方程求解平面问题的结果实质上是弹性力学位移法的数值解,该文推导并整理出PDE 的边界方程参数,用于描述弹性力学中常见的边界条件. 结合弯矩和均布载荷作用下的简支梁、受内外压的圆环板和受拉的带孔矩形板3 个例子,说明PDE 求解平面问题的方法与可视化的结果.

Pan Donghui, Ma Chongwu.

Application of MATLAB/PDE in elastic mechanics visualization teaching

Mechanics in Engineering, 2014,36(4):500-504 (in Chinese)

[本文引用: 1]

刘杨, 孟庆宇.

"互联网+项目指派"式教学模式在《弹性力学及有限单元法》教学中的应用

机械设计, 2018(S2):137-139

[本文引用: 1]

Liu Yang, Meng Qingyu.

The application of "internet + project assignment "teaching mode in the teaching of Elastic Mechanics and Finite Element Method

Journal of Machine Design, 2018(S2):137-139 (in Chinese)

[本文引用: 1]

陈小亮, 郑恒伟, 丁剑平.

弹性力学应力场的数学美感探讨

科学技术创新, 2017(22):103-104

[本文引用: 1]

Chen Xiaoliang, Zheng Hengwei, Ding Jianping.

Aesthetic stress field of elastic mechanics by maple

Scientific and Technological Innovation, 2017(22):103-104 (in Chinese)

[本文引用: 1]

卢小雨, 董春亮, 经来旺.

Maple在弹性力学教学中的应用

大众科技, 2013,15(9):104-105

[本文引用: 1]

Lu Xiaoyu, Dong Chunliang, Jing Laiwang.

Maple in elasticity teaching

Popular Science & Technology, 2013,15(9):104-105 (in Chinese)

[本文引用: 1]

邢静忠.

代数系统Maple在力学教学中的应用探讨

力学与实践, 2010,32(4):96-101

URL     [本文引用: 1]

计算机代数系统在力学教学中应用非常普遍. 摘录Maple在理论力学、材料力学、弹性力学和有限元中的应用, 以简便地完成模型求解和结果处理. 并扩展到备课、讲课和课后作业演练等环节, 提供最大限度地计算支持.将其引入力学教学, 引导和培养学生利用数学工具的习惯和能力, 强化算法设计和程序的通用性和灵活性, 为处理复杂问题提供帮助.

Xing Jingzhong.

Application of algebraic system Maple in mechanics teaching

Mechanics in Engineering, 2010,32(4):96-101 (in Chinese)

[本文引用: 1]

丁洲祥, 李涛, 白冰 .

MAPLE在土力学与基础工程研究型教学中的应用

力学与实践, 2013,35(6):87-89

DOI      URL     [本文引用: 1]

平面体系几何组成分析是结构力学的重要基础. 本文从广义三角形规则出发,提出了由“数字”引导进行平面体系几何组成分析的新思路,同时结合现有的分析方法归纳出更加完善的分析步骤,解决了初学者选择刚片、链杆以及铰时无从下手的难题.

Ding Zhouxiang, Li Tao, Bai Bing, et al.

Application of MAPLE in soil mechanics and basic engineering research-based teaching

Mechanics in Engineering, 2013,35(6):87-89 (in Chinese)

[本文引用: 1]

吴家龙. 弹性力学, 第3版. 北京: 高等教育出版社, 2016

[本文引用: 1]

Wu Jialong. Elasticity, 3rd edition. Beijing: Higher Education Press, 2016 (in Chinese)

[本文引用: 1]

/