杆梁结构的有限元分析.pptx

《杆梁结构的有限元分析.pptx》由会员分享，可在线阅读，更多相关《杆梁结构的有限元分析.pptx（50页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、【基本变量问题的基本变量基本变量问题的基本变量 由于该问题是沿由于该问题是沿x方向的一维问题，因此只有沿方向的一维问题，因此只有沿x方向的基本变量，即方向的基本变量，即 定义沿定义沿x方向移动为位移：方向移动为位移：u(x)定义沿定义沿x方向的相对伸长方向的相对伸长(或缩短或缩短)量为应变：量为应变：x(x)定义沿定义沿x方向的单位横截面上的受力为应力：方向的单位横截面上的受力为应力：第第1页页/共共50页页【基本方程问题的基本方程基本方程问题的基本方程 该问题的三大类基本方程和边界条件如下该问题的三大类基本方程和边界条件如下:取出杆件的任意一个截面，可得到平衡方程取出杆件的任意一个截面，可得

2、到平衡方程(无体力无体力)为为 取出杆件取出杆件x位置处的一段长度位置处的一段长度dx，设它的伸长为，设它的伸长为du，则它，则它的相对伸长量为的相对伸长量为 该方程称为几何方程。该方程称为几何方程。第第2页页/共共50页页 由该材料的拉伸试验，可得到该材料的虎克定律为由该材料的拉伸试验，可得到该材料的虎克定律为:该方程也称为物理方程。该方程也称为物理方程。边界条件边界条件(BC,boundary condition)第第3页页/共共50页页从求解思路来说，可以有两类方法来对该问题进行求从求解思路来说，可以有两类方法来对该问题进行求解，即解，即(1)直接求解方法：该问题比较简单，因此，可以由直

3、接求解方法：该问题比较简单，因此，可以由3个个方程来直接求解方程来直接求解3个变量。个变量。(2)基于试函数的间接方法：可以先选取一个变量基于试函数的间接方法：可以先选取一个变量(如位如位移移)作为最基本的待求变量，将其它变量都用它来表达，并作为最基本的待求变量，将其它变量都用它来表达，并采用间接的近似求解方法；具体的做法为：先对待求的位采用间接的近似求解方法；具体的做法为：先对待求的位移变量假设一种事先满足位移边界条件的可能解移变量假设一种事先满足位移边界条件的可能解(其中有一其中有一些待定的系数些待定的系数)，称为，称为试函数试函数(trail function)，让该受力系统，让该受力系

4、统的势能取最小值来最后确定出可能解的势能取最小值来最后确定出可能解(试函数试函数)中的那些中的那些待定待定系数系数(unknown constant)；也可以让该受力系统的内部变形；也可以让该受力系统的内部变形虚功等于外部施加力的虚功，来求出试函数中的那些待定虚功等于外部施加力的虚功，来求出试函数中的那些待定系数。系数。下面针对图下面针对图3-2所示的一端固定的拉杆问题，分别讨论所示的一端固定的拉杆问题，分别讨论基于直接求解方法以及基于试函数的间接方法的求解过程。基于直接求解方法以及基于试函数的间接方法的求解过程。第第4页页/共共50页页【求解原理问题的直接求解求解原理问题的直接求解 第第5页

5、页/共共50页页【求解原理问题的虚功原理求解求解原理问题的虚功原理求解 先以一个简单的结构静力平衡问题来描述虚功原理的基本思先以一个简单的结构静力平衡问题来描述虚功原理的基本思想，然后再具体求解一端固定的拉杆问题。想，然后再具体求解一端固定的拉杆问题。第第6页页/共共50页页第第7页页/共共50页页弹性力学中的虚功原理可表述为：在外力作用下处于平衡状弹性力学中的虚功原理可表述为：在外力作用下处于平衡状态的变形体，当给物体以微小虚位移时，外力所做的总虚功等于态的变形体，当给物体以微小虚位移时，外力所做的总虚功等于物体的总虚应变能物体的总虚应变能(即应力在由虚位移所产生虚应变上所作的功即应力在由虚

6、位移所产生虚应变上所作的功)。注意这里的虚位移是指仅满足位移边界条件。注意这里的虚位移是指仅满足位移边界条件BC(u)的许可位移。的许可位移。下面应用虚功应力来具体求解如图下面应用虚功应力来具体求解如图3-2所示的一端固定的拉杆所示的一端固定的拉杆问题，设有满足位移边界条件的位移场问题，设有满足位移边界条件的位移场 第第8页页/共共50页页第第9页页/共共50页页【求解原理问题的最小势能原理求解求解原理问题的最小势能原理求解 第第10页页/共共50页页第第11页页/共共50页页局部坐标系中的杆单元描述局部坐标系中的杆单元描述 最简单的标准单元就是杆单元，下面就要研究它的试函数最简单的标准单元就

7、是杆单元，下面就要研究它的试函数描述，以及计算它的应变能和外力功。描述，以及计算它的应变能和外力功。【单元构造杆单元的描述单元构造杆单元的描述 单元的描述包括单元的几何及节点描述、位移场、应变场、单元的描述包括单元的几何及节点描述、位移场、应变场、应力场、势能，也就是要充分利用描述问题的三大类变量以及应力场、势能，也就是要充分利用描述问题的三大类变量以及三大类方程来计算单元的势能，然后，由最小势能原理三大类方程来计算单元的势能，然后，由最小势能原理(或虚功或虚功原理原理)来得到单元的方程。实际上，单元内位移场的描述就是它来得到单元的方程。实际上，单元内位移场的描述就是它的试函数的选取。的试函数

8、的选取。第第12页页/共共50页页第第13页页/共共50页页第第14页页/共共50页页第第15页页/共共50页页第第16页页/共共50页页第第17页页/共共50页页【典型例题变截面杆单元的推导典型例题变截面杆单元的推导 如图如图3-5所示，有一受轴载荷的线性变截面杆件，两端的截所示，有一受轴载荷的线性变截面杆件，两端的截面积为面积为A1和和A2，长度为，长度为l，材料的弹性模量为，材料的弹性模量为E，试建立描述该，试建立描述该杆件的一个杆单元。杆件的一个杆单元。第第18页页/共共50页页第第19页页/共共50页页第第20页页/共共50页页杆单元的坐标变换杆单元的坐标变换 1.平面杆单元的坐标变

9、换平面杆单元的坐标变换 在工程实际中，杆单元可能处于在工程实际中，杆单元可能处于整体坐标系整体坐标系(global coordinate system)中的任意一个位置，如图中的任意一个位置，如图3-6所示，这需要所示，这需要将原来在将原来在局部坐标系局部坐标系(local coordinate system)中所得到的单元中所得到的单元表达等价地变换到整体坐标系中，这样，不同位置的单元才有表达等价地变换到整体坐标系中，这样，不同位置的单元才有公共的坐标基准，以便对各个单元进行集成公共的坐标基准，以便对各个单元进行集成(即组装即组装)。图。图3-6中中的整体坐标系为的整体坐标系为()，杆单元的

10、局部坐标系为，杆单元的局部坐标系为(ox)。第第21页页/共共50页页第第22页页/共共50页页第第23页页/共共50页页第第24页页/共共50页页第第25页页/共共50页页2.空间杆单元的坐标变换空间杆单元的坐标变换 第第26页页/共共50页页第第27页页/共共50页页第第28页页/共共50页页杆单元分析的杆单元分析的MATLAB程序程序 学习有限元方法的一个最佳途径，就是在充分掌学习有限元方法的一个最佳途径，就是在充分掌握基本概念的基础上亲自编写有限元分析程序，这就握基本概念的基础上亲自编写有限元分析程序，这就需要一个良好的编程环境或平台，需要一个良好的编程环境或平台，MATLAB软件就是

11、软件就是这样一个平台，它以功能强大、编程逻辑直观、使用这样一个平台，它以功能强大、编程逻辑直观、使用方便见长，从本节开始，将自主学习方便见长，从本节开始，将自主学习MATLAB程序，程序，并编写以下程序：并编写以下程序：1D杆单元的有限元分析程序杆单元的有限元分析程序(Bar1D2Node)2D杆单元的有限元分析程序杆单元的有限元分析程序(Bar2D2Node)第第29页页/共共50页页杆结构分析的算例杆结构分析的算例【典型例题四杆桁架结构的有限元分析典型例题四杆桁架结构的有限元分析 如图如图3-8所示的结构，各杆的弹性模量和横截面积所示的结构，各杆的弹性模量和横截面积都为都为E=29.510

12、4N/mm2,A=100mm2，试求解该结构的，试求解该结构的节点位移、单元应力以及支反力。节点位移、单元应力以及支反力。第第30页页/共共50页页第第31页页/共共50页页第第32页页/共共50页页第第33页页/共共50页页第第34页页/共共50页页第第35页页/共共50页页第第36页页/共共50页页3.2 梁件有限元分析的标准化表征与算例梁件有限元分析的标准化表征与算例 第第37页页/共共50页页第第38页页/共共50页页第第39页页/共共50页页第第40页页/共共50页页第第41页页/共共50页页第第42页页/共共50页页第第43页页/共共50页页第第44页页/共共50页页第第45页页/共共50页页第第46页页/共共50页页第第47页页/共共50页页第第48页页/共共50页页第第49页页/共共50页页感谢您的观看。感谢您的观看。第第50页页/共共50页页