![简明弹性力学及有限元法](https://wfqqreader-1252317822.image.myqcloud.com/cover/246/40937246/b_40937246.jpg)
第四节 物理方程——应力与应变的关系
上面已经导出了平面问题的平衡方程和几何方程。但是,仅有这两组方程,还不能求解问题。这是因为式(2-1)、式(2-3)只含应力分量,而式(2-6)只含应变分量。我们必须将它们联系起来,表示应力分量与应变分量的关系式,称为物理方程。
在完全弹性的各向同性体内,应变分量与应力分量之间的关系已在材料力学中根据广义胡克(R·hooke)定律导出如下:
![](https://epubservercos.yuewen.com/781C32/21277075708309006/epubprivate/OEBPS/Images/17323_19_1.jpg?sign=1739566065-2zL2ujMEgFADhBkg6fPbJ9rORPy9qipF-0-1c38233735fc48d0535cc3bce59e8ccd)
式中 E——拉压弹性模量,又简称为弹性模量;
μ——侧向收缩系数,又称为泊松系数;
G——剪切弹性模量,又称为刚度模量。
这三个弹性常数之间有如下的关系:
![](https://epubservercos.yuewen.com/781C32/21277075708309006/epubprivate/OEBPS/Images/17323_19_2.jpg?sign=1739566065-eR0LqZWHVb5OphfxIYPnznm16aUnSZ1b-0-fb1f504f96387574a08d4a892bfce6a0)
因此,只要知道这三个弹性常数中的任何两个,我们就可以根据式(2-9)求出第三个。
因为假定物体是线性弹性的、均匀的、各向同性的,所以,这些弹性常数不随应力或应变的大小而变,不随位置坐标而变,也不随方向而变。
在平面问题中,式(2-9)所示的应力与应变关系可以得到简化,下面分两种情况来讨论。
一、平面应力问题
在平面应力问题中,σz=τzy=τzx=0。在式(2-8)的第一式及第二式中删去σz,并将式(2-9)代入式(2-8)中的第三式,得
![](https://epubservercos.yuewen.com/781C32/21277075708309006/epubprivate/OEBPS/Images/17323_19_3.jpg?sign=1739566065-5M2v7URHPKBUiYiwCHQPhytPsMNj9LDa-0-866bd8a7bbc99b159a48e932aff4c222)
这就是平面应力问题的物理方程。
此外,式(2-8)中的第三式成为
![](https://epubservercos.yuewen.com/781C32/21277075708309006/epubprivate/OEBPS/Images/17323_20_1.jpg?sign=1739566065-zZ0KCVMnBSVgEANJGqWy170RRc5Lu5Lq-0-690c83ba3081a193924c3f6616d1cafa)
由式(2-11)可见,只要求出了σx和σy,也就可以求得εz,因此,εz不是独立的量。在以后的讨论中,我们不把它看作基本未知量,但可以由式(2-11)的εz求得薄板厚度的改变。
最后,由式(2-8)的第五、第六式可得γzy=γzx=0。
式(2-10)是由应力求应变的物理方程。反之,可由式(2-10)解出应力,即可得出由应变求应力的物理方程,可表示为
![](https://epubservercos.yuewen.com/781C32/21277075708309006/epubprivate/OEBPS/Images/17323_20_2.jpg?sign=1739566065-2RISHSVZbwHEL4ZGoToBynXHqBn2MwDq-0-0acb2cf617fc5cdc81bc5851c426ccc8)
为了便于以后的应用,式(2-12)可写成矩阵形式:
![](https://epubservercos.yuewen.com/781C32/21277075708309006/epubprivate/OEBPS/Images/17323_20_3.jpg?sign=1739566065-9jxmSklUp8fnzD4ACRGFc20nkTyCtwpP-0-d810a34b71ce8c6fac917d8136886917)
其中
[D]矩阵的元素只与弹性常数E及μ有关,故称为平面应力问题的弹性矩阵,它是一个对称矩阵。
二、平面应变问题
在平面应变问题中,因为物本的所有各点都不沿z方向移动(即ω=0),所以z方向的线段都没有伸缩,因而得εz=0。于是由式(2-8)中的第三式,得
![](https://epubservercos.yuewen.com/781C32/21277075708309006/epubprivate/OEBPS/Images/17323_20_4.jpg?sign=1739566065-JkM3JWEAnubj5AHPGkbLeuHNCz2OITHn-0-64b26f635711a5a9c1e86473b37e576f)
此式说明只要求出了σx和σy,也就随之确定了σz,因此,σz也不作为基本未知量。所以,在分析平面应变问题时,只有σx、σy和τxy三个独立的应力分量。将式(2-15)代入式(2-8)中的第一式及第二式,并注意式(2-10)中的第三式仍然适用,得
![](https://epubservercos.yuewen.com/781C32/21277075708309006/epubprivate/OEBPS/Images/17323_20_5.jpg?sign=1739566065-74BsSX47YKsgJ2PEKsOVYfN6YuXnOa4M-0-f15bb93f8c748afc13fdea95bc38ea1b)
式(2-16)就是平面应变问题的物理方程。此外,在平面应变问题中,因τzy=0,τzx=0,所以也有γzy=0,γzx=0。
如将两类平面问题的物理方程式(2-11)和式(2-16)加以比较,则可看出两者虽不一样,但是,若把平面应力问题的物理方程式(2-11)中的E换成E1、μ换成μ1,即
![](https://epubservercos.yuewen.com/781C32/21277075708309006/epubprivate/OEBPS/Images/17323_21_1.jpg?sign=1739566065-mycyQpoalTn5easKJfAWtFmiDLXU6Zdh-0-16aaa0248f2ffd2130f3c7eac4f16d58)
就得到平面应变问题的物理方程式(2-16)。其中的第三式也并不例外,读者由E1和μ1的定义不难证明
![](https://epubservercos.yuewen.com/781C32/21277075708309006/epubprivate/OEBPS/Images/17323_21_2.jpg?sign=1739566065-5Hu6aZUATa7VBh90kh6I61dsem92xF88-0-11bbfbb46fb67a0b63363d74f998e788)
应用式(2-17)的关系,将式(2-16)改写成
![](https://epubservercos.yuewen.com/781C32/21277075708309006/epubprivate/OEBPS/Images/17323_21_3.jpg?sign=1739566065-iiMUTVgcOogRdm1npqHzx2wEZ9tYJpQi-0-d6ca1240d9e8a7d2e6f0619de09f0be1)
由式(2-19)看出,两类平面问题可以按照同样的方式进行分析,只要利用转换关系式(2-17),就可将平面应力问题的物理方程转换成平面应变问题的物理方程。为了便于以后的应用,平面应当变问题的物理方程也可写成式(2-13)的矩阵形式{σ}=[D]{ε},只是其中的弹性矩阵[D],应由式(2-17),将式(2-14)转换成如下形式
![](https://epubservercos.yuewen.com/781C32/21277075708309006/epubprivate/OEBPS/Images/17323_21_4.jpg?sign=1739566065-PfgrwRZhAW566ZuCA3Ni7h3zQvBh9CAT-0-e9f7dd38492cdec6105bebb61e04c8a8)
两类平面问题的物理方程虽用统一的式(2-13),但应用时需注意:平面应力问题的弹性矩阵使用式(2-14);平面应变问题的弹性矩阵使用式(2-20)。