首页/文章/ 详情

走近弹性力学物理方程

精品
作者优秀平台推荐
详细信息
文章亮点
作者优秀
优秀教师/意见领袖/博士学历/特邀专家/独家讲师
平台推荐
主编推荐
2年前浏览2468


图片


导读:3年前,力学酒吧公 众 号是我最初学习力学理论和有限元分析的科普平台之一,现在我仍然是张伟伟老师的粉丝。张老师不仅才华横溢、文思泉涌,而且持之以恒的在力学理论和科普分享这条路越走越宽,受到了万千用户的好评,让我仰慕已久。

在2022仿真秀主办的“力学与有限元学习月”来临之际(见上图),我们获得了张老师的授权,将在仿真秀分享他的原创力学科普文章,希望能引发理工科学子和学习型工程师朋友们共鸣。并期待张老师在仿真秀平台给用户带来一场面对面的力学科普公开课。

图片

图片来网络

以下是正文:

我们知道,弹性力学基本方程包括平衡方程、几何方程和物理方程。其中,物理方程表征了材料在发生变形时应力与应变之间所满足的关系,因此,也被称为应力-应变关系。因其又反应材料的本质特性,也被称为本构关系,在弹性力学中简化为广义胡克定律。弹性力学旨在构建固体变形的一般规律,本文将从数学连续性的角度出发,考察在三维空间一般情形下的应力-应变关系,并讨论如何由一般情况导出广义胡克定律。


一、一般情形
在三维空间一般情况下,考虑剪应力互等定理,一点处的应力状态和应变状态可由6个独立的应力分量和6个独立的应变分量唯一确定,它们分别为
图片
认为应力分量为应变分量的一般函数,设
图片
这就将每一个应力分量写成了关于6个应变分量的多元函数。我们以σ1为例来考虑,只要函数f1为光滑连续函数,利用多元函数泰勒级数展开,可得
图片
图片
图片
因此,式(1)可近似写成一次函数
图片
为了方便描述,我们做如下标记:
图片
图片
图片
类似地,所有的应力分量都利用泰勒级数展开,并忽略二阶以上的高阶无穷小,都可以写成式(2)的形式。将应力分量和应变分量各自排成一列成为列阵,即
图片
因此,三维一般情况下,应力-应变关系可以表示为矩阵形式
图片
这里,[C]为物理方程的系数矩阵,共36个常数。可以看出,当系数矩阵中所有常数均不为0时,所有6个应变分量都会影响任何一个应力分量,这将意味着剪切应变不仅会引起切应力的改变,还会引起到正应力的改变;同样,线应变的改变不仅会引起正应力的改变,也会引起切应力的改变。反过来,应力对应变的影响也类似的耦合在一起。具有这样特性的材料被称为极端各向异性材料,这样的材料在任意方向上都将表现出不同的力学性能。
不过,即便是极端各向异性材料,系数矩阵的36个常数也并非完全独立的。我们知道,当压缩一个线性弹簧时,外力做功等于弹性势能的增加量,为
图片
图片
上式中假定了材料为线弹性材料。并且e表示应变能密度。引入张量记号,上式可写为
图片
在张量体系中,具有某一项中有两个相同的下标时,称之为哑标,表示进行求和计算,将式(5)展开后,即为式(4)。再回到式(3),针对于每一个应力分量,利用哑标可写为
图片
这里,有一对下标j,表示需要求和。即上式展开后,可写成
图片
当i从1变化到6,所有的应力分量均可表出。现在把式(6)代入到式(5)中,得

图片

二、具有一个弹性对称面的情形
所谓弹性对称面具有这样的性质:关于对称面对称的点处力学性能完全一致。设某材料具有一个弹性对称面,将该对称面设为xoy坐标平面,这样在材料上建立图1所示的两个对称坐标系,在两个坐标系中求解力学量,并不会因为坐标系的不同而不同。


图片


图1 关于xoy平面对称两个对称坐标系
仍考虑应变能密度e,将式(7)展开后,有
图片
假设材料上任意一点A,设其在坐标系(a)中的坐标为(x,y,z), 则在坐标系(b)中坐标为(x,y,-z)。假设A点在坐标系(a)中z方向的位移为w,则在坐标系(b)中z方向的位移为-w。考察几何方程,
图片
图片
图片
系数矩阵可记为
图片
这样,对于物理方程系数矩阵,其独立的21个常数又减少了8个,共有13个独立常数。
三、具有三个弹性对称面的情形
如果某一材料有三个弹性对称平面,如图2所示,设空间任意点A在坐标系(a)中的坐标为(x,y,z),在坐标系(b)中的坐标为(x,-y,z),在坐标系(c)中的坐标为(-x,y,z),在坐标系(d)中的坐标为(x,y,-z),因此,对坐标求导后,在两对称的坐标系内,互为相反数。位移分量u,v,w在各自的对称坐标系下也会反号。


图片

图2  关于三个坐标面对称的坐标系
(a)与(b)关于xoz对称;(a)与(c)关于zoy对称;(c)与(d)关于xoy对称
图片
四、横观各向同性情形
图片

图片

图3 微元体模型

图片

图4  纯剪应力状态
图片
图片
图片
图片
五、各向同性情形
各向同性材料表示材料任意方向的性质都相同,因此在x,y,z方向施加载荷,弹性系数相同,即
图片
图片
图片
图片
根据拉梅给出的物理方程
图片
图片
六、最后一个问题

在应用物理方程时,我们经常听说物理方程的其它名称,如本构关系,应力-应变关系,广义胡克定律。实际上,这个名称并不完全相同,或者它们各有侧重。

称其为物理方程时,是相对于平衡方程、几何方程而言的。平衡方程表达了微元体受力方面要满足的平衡条件,几何方程表达了微元体发生变形时需要满足的几何条件。这两类方程的建立不涉及任何有关材料的性质,物理方程则补充了材料的性质,把特定材料的力和变形联系起来,反应了材料物理变形特性。

称其为本构关系时,有一种专门把材料的力与变形性能单独拿出来进行研究的意味。如材料的本构关系研究中,就不涉及材料的平衡方程、几何方程,仅仅是研究材料的力学性能。更具体一些,本构关系研究主要在于判定某种材料的本构关系含有多少个独立常数,并证明它们的完备性,以及通过实验确定出这些常数的值。

应力-应变关系是一个不严谨的称呼。例如,一种新材料在某一方向上进行了拉伸实验,得到了加载条件下的应力-应变曲线,也被称为材料的应力-应变关系。但是,这条曲线无论如何也不能称之为本构关系。为了获得本构关系,需要从不同的加载方向、以及不同的加载形式(拉伸、压缩、扭转、剪切等)对材料做全面的力学性能实验,最终获得完备的、不可约的本构常数(本构关系的系数矩阵)。这时才能称之为本构关系。

可见,应力-应变关系,本构关系,物理方程,三个概念之间有一种递进关系。在实验中通过加载测试,得到的应力-应变曲线都可以成为应力-应变关系,但还不能称之为本构关系;当对材料进行各向加载、获得全面的力学性能后,通过分析、加工可获得本构关系,算是较“全面地”掌握了材料的力学性能;本构关系可以作为物理方程,与平衡方程、几何方程一起,组成封闭的方程组,对弹性力学问题进行求解。

至于通常所说的广义胡克定律,则是对各向同性材料本构关系在线弹性条件下的一个简化。在力学中,只要提到胡克定律,就指力与变形成正比例关系。“广义”突出了材料在复杂应力状态下,以区别于单向应力状态。广义胡克定律从形式上可称为应力-应变关系,它表征了各向同性材料的本构关系,可以与平衡方程、几何方程一起,作为物理方程求解弹性力学问题。

参考文献
张少实, 庄茁. 《复合材料与粘弹性力学》第2版. 机械工业出版社. 2011.6

徐芝纶. 《弹性力学》第5版. 高等教育出版社. 2016.3

(完)

:张伟伟 仿真秀专栏作者

声明:原文首发力学酒吧公 众号,本文已授权转载。部分图片和内容源自网络,如有不当请联系我们,欢迎分享,禁止私自转载,转载请联系我们。
科普理论通用瞬态动力学静力学结构基础其他软件
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2022-03-08
最近编辑:2年前
仿真圈
技术圈粉 知识付费 学习强国
获赞 10083粉丝 21544文章 3539课程 219
点赞
收藏
作者推荐
未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈