结构力学第六章平面应力问题的有限单元法

1、结构力学第六章平面应力问题的有限单元法第1页，共32页，2022年，5月20日，19点0分，星期二 有限元方法已被公认为应力分析的有效工具，受到工程界科学研究者的广泛重视。现在的有限元软件种类繁多，较为突出的有：MSC/NASTRAN、ANSYS、MARC、Algor、ABAQUS、等。 本章主要通过板的平面应力问题介绍有限元基本方法。引言：第2页，共32页，2022年，5月20日，19点0分，星期二引言：第3页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 1.平面应力问题t/2t/2Oxyzy图6-1第4页，共32页，2022年，5月20日，19点0

2、分，星期二6.1平面应力问题及其基本方程式 因为板面上(z= t/2 )不受力，所以有应力分量 由于板很薄，故可以认为在整个薄板上所有各点上都有： 根据剪应力互等定律：第5页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 这样，板中任一点的九个应力分量就只剩下三个应力分量，即第6页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 因而这种问题称为平面应力问题。同时，由于板很薄，所以这三个应力分量，以及分析问题时须考虑的三个应变分量x 、y 、xy及和两个位移分量u，v，都可以认为沿厚度不变化。这就是说，它们只是坐标x

3、和y的函数，不随坐标z的变化而变化。 在平面应力问题中，可用如下三个向量分别表板中任一点的应力、应变和位移；(6-1)(6-2)(6-3)第7页，共32页，2022年，5月20日，19点0分，星期二 在船体结构中，很多问题可以简化为平面应力问题处理。例如甲板开口、舷侧、横梁开孔和肘板的强度问题等等。6.1平面应力问题及其基本方程式 2.基本方程式 基本方程式包括平衡微分方程式、几何方程式和物理方程式，此外还有边界条件方程式。下面依次到处平面应力问题的这些方程式 第8页，共32页，2022年，5月20日，19点0分，星期二xyxyCXYxy6.1平面应力问题及其基本方程式 (1)平衡微分方程式

4、根据平衡条件导出的各应力分量之间的微分关系就是平衡微分方程式图6-2第9页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 首先以通过中心C，并平行于z轴的直线为矩轴，列力矩平衡方程Mc=0。第10页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 其次，以x轴为投影轴，列出力投影的平衡方程Fx=0：第11页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 其次，以y轴为投影轴，列出力投影的平衡方程Fy=0：第12页，共32页，2022年，5月20日，19点0分，星期二6.1平面应

5、力问题及其基本方程式 由此，得到平面应力问题的平衡方程： 这两个微分方程中包含三个未知数x、 y、xy =yx：因此，决定应力分量的问题是超静定，还必须考虑变形情况才能解决问题。(6-4)第13页，共32页，2022年，5月20日，19点0分，星期二xy (2)几何方程式 下面从平面应力问题的几何学方面，导出应变分量与位移分量之间的关系式，即几何方程式。O第14页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 线段PA的正应变： 线段PB的正应变：(a)(b)第15页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式

6、剪应变xy： 线段PA的转角为： 线段PB的转角为：第16页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 剪应变xy： 式(a)、(b)、(c)是应变分量与位移分量之间的关系式，现归纳为：称为平面应力问题的几何方程式，又称柯西方程式(c)(6-5)第17页，共32页，2022年，5月20日，19点0分，星期二 由式（6-5）可见，当板内各点的位移分量u、v为已知函数时，就可确定各点的应变分量x、 y、 xy；反之，假定三个应变分量函数，那么按式（6-5）的前两个式子就可以求出位移函数u、v，若用此两个位移分量函数代入第三个方程式求xy，就会与假定的xy

7、不同。这样就出现了矛盾。这一矛盾是因为板内任一点的应变分量之间有相互联系所造成的，如果在假定各应变分量函数时不反映出这种联系，那就将使变形不连续，即板变形将发生空隙或裂缝。从数学上讲，式（6-5）的应变分量是三个，而位移分量只有两个，因此三个应变分量不能相互独立，而必然有一定的联系，这个关系叫应变协调方程或应变相容条件。现推导如下：6.1平面应力问题及其基本方程式第18页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 应变相容方程或叫应变协调方程： 称为变形连续方程或圣维南方程(6-6)第19页，共32页，2022年，5月20日，19点0分，星期二6.1

8、平面应力问题及其基本方程式 （3）物理方程式 上面建立了二个平衡微分方程式和三个几何方程式，五个方程式中共有八个未知函数x、 y 、xy 、u、v、x、 y、 xy，尚需补充三个方程式才可能求解。这三个方程式就是应变分量与应力分量之间的关系式，即物理方程式。物理方程式就是材料力学中广义虎克定律，平面应力状态下的广义虎克定律为：第20页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 或写成：(6-7)(6-8)第21页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 用矩阵形式表示为： 简记为：称为平面应力问题的弹性矩

9、阵(6-9)(6-10)第22页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 （4）边界条件 上述式（6-4）、（6-5）、（6-7）或（6-8）一共有八个基本方程式。这八个基本方程式中包含八个未知函数（坐标x,y的函数）；三个应力分量；三个应变分量；两个位移分量。方程的数目恰好等于未知函数的数目。因此，在适当的边界条件下，由上述方程式求解未知函数是可能的。 根据边界条件的不同，平面应力问题分为位移边界问题、应力边界问题和混合边界问题。第23页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 （4）边界条件 在位移

10、边界问题中，物体在全部边界上的位移分量是已知的，也就是在全部边界上有：式中，us和vs在边界上是坐标的已知函数。这就是平面应力问题的位移边界条件。(6-11)第24页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 （4）边界条件 在应力边界问题中，物体在全部边界上所受的面力是已知的，也就是说，面力分量X和Y在全部边界上是坐标的已知函数。根据面力分量与边界上应力分量之间的关系，可以把面力已知条件转换为应力已知的条件，这就是所谓应力边界条件，推导如下：xNXYYX第25页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式

11、在物体边界上取直角三角形的微块PAB，它的斜面AB与物体边界面重合，如图6-4所示。用N表示边界面的外法线法向，其方向余弦为：xNXYYX第26页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 设该边界面AB 的长度为ds，则截面PA和PB的长度分别为lds和mds，另设微块的厚度为1。由微块的平衡条件Fx=0，得第27页，共32页，2022年，5月20日，19点0分，星期二6.1平面应力问题及其基本方程式 同理由微块的平衡条件FY=0，得另外一方程式。这两个方程式为： 式（6-12）就是平面应力问题的应力边界条件。如果考虑第三个平衡条件FM=0，则可以

12、再写出一个方程式，但是在ds趋近于零时，设一方程式将成为剪应力互等方程式。(6-12)第28页，共32页，2022年，5月20日，19点0分，星期二6.2解题方法及有限单元法概念 1.解题方法简介 传统的弹性力学解题方法有三种：位移法，应立法和混合法。 当平面应力问题按位移法求解时，以位移分量u、v为基本未知函数。将几何方程式（6-5）代入物理方程式（6-8），得应力分量与位移分量之间的关系式，再将所得关系式代入平衡微分方程式（6-4）和应力边界条件（6-12），简化后得到：(6-13)第29页，共32页，2022年，5月20日，19点0分，星期二6.2解题方法及有限单元法概念 和：(6-14)第30页，共32页，2022年，5月20日，19点0分，星期二6.2解题方法及有限单元法概念 当平面应力问题按应力法求解时，以应力分量为基本未知函数。将几何方程式