当前课程知识点:Finite Element Method (FEM) Analysis and Applications > 9、Finite element analysis of continuum structure (2) > 9.1 Three-dimensional 4-node tetrahedral element and MATLAB programming > Video 9.1
返回《Finite Element Method (FEM) Analysis and Applications》慕课在线视频课程列表
返回《Finite Element Method (FEM) Analysis and Applications》慕课在线视频列表
同学们好,首先我们回顾一下上一讲的主要内容
上一讲我们学习了平面3节点三角形单元及MATLAB编程
平面4节点矩形单元及MATLAB编程
轴对称单元
分布力的处理
平面矩形薄板分析的MATLAB编程
平面矩形薄板的ANSYS分析实例
这一讲我们将讲解
空间4节点四面体单元及MATLAB编程
空间8节点六面体单元及MATLAB编程
参数单元的原理
数值积分
典型空间问题的MATLAB编程
最后我们将介绍典型空间问题的ANSYS分析实例
空间4节点四面体单元及MATLAB编程
首先进行节点的描述
我们看一下这个单元
这个单元有4个节点,节点1,2,3,4
每个节点有3个坐标(xi,yi,zi)
我们几何坐标就是1号节点到4号节点
分别给出它对应的3个坐标
节点的位移呢
每一个节点它有3个方向的自由度
也就是u,v,w
一共有4个节点,所以一共有12个自由度
也就是说它是一个节点位移的列阵12X1
那么和它对应的是节点力的列阵
同样也是,每个节点有3个力的分量
基于节点的描述,我们要进行单元的场的描述
要把三大类变量,也就是位移场、应力场、应变场
分别用节点位移来进行描述
首先进行位移场的描述
位移场我们就是要取位移模式
我们把所有节点的x方向的位移来对
x方向的位移场进行插值
也就是u(x,y,z)这个位移场
我们有4个
有u1,u2,u3,u4来进行插值
我们可以取4个待定系数
还是按照唯一确定性原则
从低阶到高阶这么一个位移模式的两个原则
我们可以设
对于y方向的位移场,也就是v(x,y,z)
同样我们也设定4个待定系数,分别是
对于z方向的位移场也是这样的
我们分别为
我们要用节点条件来求这12个待定系数
也就是说这个位移场在1号节点要等于ui,vi,wi
2、3、4号节点也要分别等于相应的节点位移分量
这样我们就可以得到一组方程组
由这个方程组我们把
a0,a1,a2,a3一直到c0,c1,c2,c3
把这12个系数全部求解出来
求解出来过后就把它表达成节点位移的描述
也就是说将12个待定系数用12个节点位移
也就是u1,v1,w1一直到u4,v4,w4来代替
这样我们重新来写位移场的函数
我们对x方向的位移也就是u(x,y,z)
把它可以写成
对于y方向位移分量的场,也就是v(x,y,z)
同样我们也可以写出来
对于z方向的位移分量,这个场的表达,就是w(x,y,z)
我们同样可以得到
其中这个Ni就是形状函数
它可以表达成坐标的函数
也就是说我们可以写成
这里面的ai,bi,ci,di这些系数都可以表达成
1、2、3、4号节点坐标的关系
相应的系数可以查阅相应的资料
我们把位移场函数写成矩阵形式
也就是说把3个分量写成一个列向量
这样我们就写出N1,N2,N3,N4这么一个组合的
3X12的N矩阵,我们叫形状函数矩阵
再乘上qe,也就是节点位移的列阵
对于应变场,首先我们由几何方程
几何方程,我们三维问题有6个独立的分量
也就是3个主方向的应变还有3个剪切的应变
我们把它写成一个算子矩阵乘上相应的3个分量的关系
这是一个偏导的6X3的矩阵
作用在位移上面
这个位移场我们前面已经得到了
它是N矩阵乘上qe列阵
我们把它代进去
也就是说对于这个几何算子矩阵,这个偏导矩阵
它就作用在N矩阵上面
这样我们就得到B矩阵
这个B矩阵我们同样可以进行一个计算
这样我们也通过几何方程把应变场表达成
一个几何矩阵乘上节点的位移列阵qe
具体对于几何矩阵的表达
它具体是作用在N矩阵上
对于每一个分块矩阵Bi
每一个分块矩阵都是6X3的
6X3具体的表达
它就是算子矩阵作用在3X3的这么一个对角线的
形状函数矩阵上面
具体也可以得到各个系数的表达
每个元素的表达
同样这些元素和前面形状函数里面的元素是一样的
在得到应变场的表达以后
我们再来看看应力场的表达
那我们只要把弹性系数矩阵作用在应变上
也就是说物理方程, 我们就可以得到应力的表达
这个应力的表达是6X1的
也就是说是有6个分量
那么我们把B矩阵乘D矩阵合到一块
相乘以后我们把它叫做应力函数矩阵
在得到了三大类变量的表达
也就是说把一个单元里面三大类变量的这个场
分别表达成我们所对应的节点位移的关系
也就是说位移等于
应变等于
应力等于
基于这些表达,我们再来看看这个单元的势能
我们都知道单元的势能是由应力乘应变
在整个域上进行积分
这是应变能减外力功
外力功是体积力乘上相应的位移
加上面积力乘相应的位移
我们把应力、应变还有相应的位移
当然这都是单元里的场的函数
把前面已经得到的场的表达代进去
这样我们就可以得到单元的应变能
表达成基于节点位移qe的表达式
具体表达出来就是
我们把中间这个积分叫单元刚度矩阵
同样后面也是,减去外力功这一块
这后面这个积分,我们叫单元节点的等效载荷
我们看看,单元的刚度矩阵等于
同样,如果有体积力和面积力的作用
通过形状函数矩阵对它进行作用以后
作出相应的积分
我们就可以得到相应的等效节点载荷
那么我们也是把单元的势能表达成基于节点的描述
那节点的位移就是我们的待定系数了
我们用最小势能原理对势能取极值
也就是对这个待定系数求偏导,让它等于0
这样我们就得到
这就是我们的刚度方程
这个是12X12,节点位移是12X1
节点载荷是12X1的
对于空间4节点四面体单元的性质我们进行一个讨论
首先看一下位移
我们看看位移场
前面我们基于节点位移进行了一个插值
也就是说唯一确定性、从低阶到高阶
我们分别对x,y,z方向的位移进行了一个插值
我们看一下,它是一个完全的线性函数
也就是它的描述能力是一个完全线性的描述能力
我们把这个关系代到应变场里去
应变场实际上就是一个偏导算子对它进行作用
前面完全线性的位移场代进去以后
进行这么一个算子矩阵的作用
我们就分别得到了一些系数
我们可以看看
这个应变场得到的各个分量全部都是系数
也就是常数
同样,应力场
应力场是在应变场的基础上
前面再乘上一个弹性系数矩阵
所以它也是一个常系数的各个分量
那么简单总结一下空间4节点四面体单元的性质
第一个就是单元内的应力及应变都为常数
所以也称为空间问题的常应力及常应变体单元
它和平面3节点三角形单元的性质是完全一样的
也是对应的
第二个性质就是单元的节点位移是在整体坐标系下建的
也就是说我们每一个节点都有3个分量
都是在整体坐标系下建立的
因此,没有针对节点位移的坐标变换的问题
第三个就是对于应变梯度较大
也就是变化比较剧烈的区域
那么我们这个单元由于是常应力或者常应变单元
它的计算误差是比较大的
对于变化比较剧烈的地方就一定要划更多的单元进行加密
这样才能得到比较满意的结果
对于空间4节点四面体单元的MATLAB编程
同样我们要定义3个函数
也就是说进行刚度矩阵计算的Stiffness
进行单元组装的Assembly
以及单元应力计算的Stress
对于Stiffness这个函数
要求输入弹性模量、泊松比
还需要输入4个节点的x,y,z坐标
输出是单元刚度矩阵k,它是12X12的
对于Assembly这个函数
它是输入单元刚度矩阵k
还要求输入单元的节点编号,也就是i,j,m,n
输出组装过后的整体单元刚度矩阵KK
对于计算单元应力的函数Stress
它要求输入弹性模量、泊松比
以及4个节点的坐标
还要求输入单元的位移列阵,u(12X1)的
它输出单元的应力Stress,它是6X1的
那这个6X1也就对应着这个体单元的6个应力分量
-Finite element, infinite capabilities
--Video
-1.1 Classification of mechanics:particle、rigid body、deformed body mechanics
--1.1 Test
-1.2 Main points for deformed body mechanics
--1.2 Test
-1.3 Methods to solve differential equation solving method
--1.3 Test
-1.4 Function approximation
--1.4 Test
-1.5 Function approximation defined on complex domains
--1.5 Test
-1.6 The core of finite element: subdomain function approximation for complex domains
--1.6 Test
-1.7 History and software of FEM development
--1.7 Test
-Discussion
-Homework
-2.1 Principles of mechanic analysis of springs
--2.1 Test
-2.2 Comparison between spring element and bar element
--2.2 Test
-2.3 Coordinate transformation of bar element
--2.3 Test
-2.4 An example of a four-bar structure
--2.4 Test
-2.5 ANSYS case analysis of four-bar structure
--ANSYS
-Discussion
-3.1 Mechanical description and basic assumptions for deformed body
--3.1 Test
-3.2 Index notation
--3.2 Test
-3.3 Thoughts on three major variables and three major equations
--3.3 Test
-3.4 Test
-3.4 Construction of equilibrium Equation of Plane Problem
-3.5 Test
-3.5 Construction of strain-displacement relations for plane problems
-3.6 Test
-3.6 Construction of constitutive relations for plane problems
-3.7 Test
-3.7 Two kinds of boundary conditions
- Discussion
-- Discussion
-4.1 Test
-4.1 Discussion of several special cases
-4.2 Test
-4.2 A complete solution of a simple bar under uniaxial tension based on elastic mechanics
-4.3 Test
-4.3 The description and solution of plane beam under pure bending
-4.4 Test
-4.4 Complete description of 3D elastic problem
-4.5 Test
-4.5 Description and understanding of tensor
-Discussion
-5.1 Test
-5.1Main method classification and trial function method for solving deformed body mechanics equation
-5.2 Test
-5.2 Trial function method for solving pure bending beam: residual value method
-5.3 Test
-5.3How to reduce the order of the derivative of trial function
-5.4 Test
-5.4 The principle of virtual work for solving plane bending beam
-5.5 Test
-5.5 The variational basis of the principle of minimum potential energy for solving the plane bending
-5.6 Test
-5.6 The general energy principle of elastic problem
-Discussion
-6.1Test
-6.1 Classic method and finite element method based on trial function
-6.2 Test
-6.2 Natural discretization and approximated discretization in finite element method
-6.3 Test
-6.3 Basic steps in the finite element method
-6.4 Test
-6.4 Comparison of classic method and finite element method
-Discussion
-7.1 Test
-7.1 Construction and MATLAB programming of bar element in local coordinate system
-7.2 Test
-7.2 Construction and MATLAB programming of plane pure bending beam element in local coordinate syste
-7.3 Construction of three-dimensional beam element in local coordinate system
-7.4 Test
-7.4 Beam element coordinate transformation
-7.5 Test
-7.5 Treatment of distributed force
-7.6 Case Analysis and MATLAB programming of portal frame structure
-7.7 ANSYS case analysis of portal frame structure
-8.1 Test
-8.1 Two-dimensional 3-node triangular element and MATLAB programming
-8.2 Test
-8.2 Two-dimensional 4-node rectangular element and MATLAB programming
-8.3 Test
-8.3 Axisymmetric element
-8.4 Test
-8.4 Treatment of distributed force
-8.5 MATLAB programming of 2D plane rectangular thin plate
-8.6 Finite element GUI operation and command flow of a plane rectangular thin plate on ANSYS softwar
-Discussion
-9.1 Three-dimensional 4-node tetrahedral element and MATLAB programming
-9.2 Three-dimensional 8-node hexahedral element and MATLAB programming
-9.3 Principle of the isoparametric element
-9.4Test
-9.4Numerical integration
-9.5 MATLAB programming for typical 2D problems
-9.6 ANSYS analysis case of typical 3Dl problem
-Discussion
-10.1Test
-10.1Node number and storage bandwidth
-10.2Test
-10.2 Properties of shape function matrix and stiffness matrix
-10.3Test
-10.3 Treatment of boundary conditions and calculation of reaction forces
-10.4Test
-10.4 Requirements for construction and convergence of displacement function
-10.5Test
-10.5C0 element and C1 element
-10.6 Test
-10.6 Patch test of element
-10.7 Test
-10.7 Accuracy and property of numerical solutions of finite element analysis
-10.8Test
-10.8 Error and average processing of element stress calculation result
-10.9 Test
-10.9 Error control and the accuracy improving method of h method and p method
-Discussion
-11.1 Test
-11.1 1D high-order element
-11.2 Test
-11.2 2D high-order element
-11.3 Test
-11.3 3D high-order element
-11.4 Test
-11.4 Bending plate element based on thin plate theory
-11.5 Test
-11.5 Sub-structure and super-element
-12.1Test
-12.1 Finite element analysis for structural vibration: basic principle
-12.2 Test
-12.2 Case of finite element analysis for structural vibration
-12.3 Test
-12.3 Finite element analysis for elastic-plastic problems: basic principle
-12.4 Test
-12.4 Finite element analysis for elastic-plastic problems: solving non-linear equations
-Discussion
-13.1 Test
-13.1 Finite element analysis for heat transfer: basic principle
-13.2 Test
-13.2 Case of finite element analysis for heat transfer
-13.3 Test
-13.3 Finite element analysis for thermal stress problems: basic principle
-13.4 Test
-13.4 Finite element analysis for thermal stress problems: solving non-linear equation
-Discussion
-2D problem: finite element analysis of a 2D perforated plate
-3D problem: meshing control of a flower-shaped chuck
-Modal analysis of vibration: Modal analysis of a cable-stayed bridge
-Elastic-plastic analysis: elastic-plastic analysis of a thick-walled cylinder under internal pressur
-Heat transfer analysis: transient problem of temperature field during steel cylinder cooling process
-Thermal stress analysis: temperature and assembly stress analysis of truss structure
-Probability of structure: Probabilistic design analysis of large hydraulic press frame
-Modeling and application of methods: Modeling and analysis of p-type elements for plane problem