当前课程知识点:Finite Element Method (FEM) Analysis and Applications >  11、High-order and complex element >  11.4 Bending plate element based on thin plate theory >  Video 11.4

返回《Finite Element Method (FEM) Analysis and Applications》慕课在线视频课程列表

Video 11.4在线视频

下一节:Video 11.5

返回《Finite Element Method (FEM) Analysis and Applications》慕课在线视频列表

Video 11.4课程教案、知识点、字幕

以前我们学了细长梁的单元,我们叫梁单元

下面我们来学习一下

基于薄板理论的弯曲板单元

首先我们看看基本变量和方程

这是一个薄板

这个薄板z方向是厚度方向

x方向和y方向是它的面内的方向

这个薄板有几个特征

它的厚度远小于其它两个方向的尺寸

这就是薄板的一个特征

那么我们要引入一定的假设

对厚度方向的特征进行简化

这个简化我们叫Kirchhoff假定

由于薄板中要保持转角的连续,所以它要承受弯矩

因而薄板问题是一个C1问题,它要保持转角的连续性

我们看看,对薄板的受力问题进行简化的

Kirchhoff假定

第一条就是薄板的中面法线变形后仍然为直线

并且厚度方向的正应变很小,可以近似为0

这对应着我们细长梁的直法线假定

也就是说

第二个假定是薄板中面没有横向位移,则会导出

这是在中面上

那么离开中面的位置上z,它有u和v的位移

这就相当于是我们的直法线假定了

第三个假定就是应力σzz引起的变形很小,可以忽略

我们看看,薄板的各个不同的横截面上的应力情况

这个相应的应力有σxx,σyy,τzx,τzy

我们把每一个面上的应力进行合成以后就可以得到

Mx,My,Mxy,Myx,另外还有Qx,Qy

它分别对应着各个侧面上的弯矩、扭矩和剪力

我们把三大类变量归纳一下

位移就是中面上的位移,就是w(x,y)

它是在z=0这个中面上,它作为一个基本位移

u和v满足直法线假定,也就是分别等于

应变在前面的基础上求一个一阶导数,这样就可以得到

应力是在应变的基础上乘上一个相应的弹性系数矩阵

我们再看看三大方程

平衡方程

平衡方程我们也是取出中间的一个微元体

就是dxdy,连厚度取完就是h

取完以后我们建立一个垂直方向的平衡

我们就可以得到

也可以把它写成弯矩的平衡关系

其中的D0等于

这个μ是泊松比

它为薄板弯曲问题里的弯曲刚度

物理方程,这就是广义虎克定律

写出来就是这样的

那么我们把物理方程

分别沿着面上把它求一个合力

这样我们可以把物理方程

写成一个广义力与广义位移的关系

这个M就是

D是弹性系数矩阵,它等于

几何方程是和应变的定义是一样的

实际上就是应变的表达式

边界条件我们有三种情况

第一种是S1

边界上是要等于给定的挠度和转角

这个n是边界法线的方向

在S2上面要等于指定的挠度和指定的力矩

n也是边界的法线方向

在S3边界上,它也是要等于给定的力矩和横向载荷

这个n为边界的法向,m为边界的切向

Qn为边界截面上单元长度的横向剪力

我们定义薄板的总势能

这是对dxdy面上进行积分

同样对Qn,这是给定的边界力

乘上相应的挠度

力矩Mn乘上相应的转角

这个地方我们看看

к实际上由三项组成

就是εx,εy,γxy

这三项分别又是等于

这个地方的w就是我们的基本变量

所有都可以表达成薄板中面的挠度w

也就是z=0这个位置的挠度

它作为我们薄板问题的基本变量

在基于三大类变量和三大类方程描述的基础上

我们来构造薄板的板单元

我们来看看,4节点矩形的薄板单元

首先4个节点

每个节点我们有3个自由度

也就是说挠度w,另外绕着x轴的转角

和绕着y轴的转角

我们把每个节点的自由度表达成qi

它有3项,就是wi,θxi,θyi

由于一共有4个节点,所以一共有12个自由度

我们的基本变量就是w,就是中面的挠度

12个节点我们可以确定12个待定系数

我们分别取为,从a1一直取到a12

前面有一次项,完整的二次项,完整的三次项

另外还有两项是四次项,是不完全的

当然我们取这两项是考虑了一个对称性

所以分别取

按照收敛性准则一,完备性要求

就是说要有刚体位移,要有常应变

那么我们从a1到a3的完全一次项就是描述刚体位移的

这个是完全的二次项

它可以描述薄板弯曲问题的常应变

所以它是满足完备性的要求

我们来考查一下收敛性准则二

也就是连续性的要求

假定有两个单元,①号单元②单元

我们分别要考查这两个单元交界

这个交界有沿着y轴的交界,还有沿着x轴的交界

这两个交界,我们看看它的挠度的连续性的问题

我们首先考查一下沿着y方向的交界

也就是沿着L1和L2这条线

我们看看这条线的坐标就是x=0,就是沿着y轴了

我们把这个x=0代入w函数里面去

这样我们就得到了这么4项

同样我们看看,沿着这个方向的转角是

在x=0上,这个时候的转角是

这个里面的待定系数我们看看是a1,a3,a6,a10

这4个系数完全要用1号节点的函数值

也就是(w1,θx1)和2号节点的函数值(w2,θx2)唯一确定

我们看看,它确实能唯一确定

所以沿着L1和L2这个接线上应该没有问题

我们再考查一下沿着x方向,也就是沿着L3,L4这两条线

它的连续性的问题

同样我们对这两条线的交点A点

考查一下它的一阶导数,也就是

我们看看它是怎么过渡

我们一阶导数求了过后就得到这么4项

同样这4个系数也要由节点1和节点2的节点值来确定

对于节点1和节点2的自由度

每个节点有3个自由度,就是(w1,θx1,θy1)

2号节点也是3个

那么在前面的L1和L2的交线处我们已经用了4个条件

分别已经用了(w1,θx1)和(w2,θx2)

这个时候节点条件里面只剩下θy1和θy2

由θy1和θy2不能确定刚才我们考查L3和L4在交线处

在A点处的那4个待定系数

不能唯一确定

所以我们就可以简单归纳一下

在单元交界面上w是连续的

但单元之间的法向导数连续性一般不能满足

所以这种单元是非协调的

但由于这种单元能通过拼片试验

所以当单元尺寸,或者是网格不断缩小的时候

计算结果还是可以收敛于精确解的

我们再研究一下3节点三角形的薄板单元

3个节点,每个节点3个自由度

w,沿着x和y方向的偏导

这样我们3个节点一共是9个自由度

3节点三角形单元

我们由9个自由度可以确定9个待定系数

而一个完全的三次多项式要包含10项

也就是说如果取完全三次多项式

就是从a1一直取到a10

我们唯一性确定原则里面只能确定9个

所以中间有一个问题

我们考虑一个对称性

我们把这个完全三次多项式的10项里面的两项

这两项我们把它合到一块,也就是令a8=a9

这样我们把这两个自由度就合成一个

插值的基底函数就分别为一次项、二次项,三次项是

这样我们就把它写成这么一个插值的关系

其中这个a是9X1的

这就是从a1到a9的待定系数

我们把节点条件代入到插值关系里面去

就可以得到这么一个节点的位移要等于这么一个插值关系

这个表达是通过H矩阵来转换的

H矩阵就是由单元3个节点的几何位置来确定的系数矩阵

我们由上面这个方程要求出a

再代入的形状函数矩阵那个方程里面去

求a的时候,要对H矩阵求逆

H的行列式是这么一个表达

从行列式可以看出来

当三角形的高满足一个关系,也就是说

这个行列式就为0,它就不能求逆

这种情况也包含等腰直角三角形

不能求逆就意味着这种类型的单元有缺陷

其原因就是前面强行令一个三次项

令它们的系数相等

就是硬让它们相等

所以这种具有“零函数”的单元是一个“亏损”单元

这个困难可以通过引入面积坐标的方法加以克服

将薄板的挠度用我们面积坐标

就是说L1,L2,L3来进行表达

那么我们可以取出9项

那么这是一次项,这是三次项

其中在三次项中引入了一个系数c

由于这个函数不是完全的三次式

所以不能保证w满足常应变的要求

但是当c取1/2的时候

前面所取的表达式正好就满足常应变

所以我们就把c取为1/2

这样,我们经过推导可以得到

3节点三角形薄板单元的形状函数矩阵是这样的

这个是1X9的,里面分别是N1,N2,N3

N1,N2,N3具体的表达是3项

N1基于面积坐标的表达是这么一个关系

对于N2和N3也可以通过上面N1表达的下标的轮换来得到

上面这个表达式的bi,ci

是等于三个节点坐标的这么一个关系

我们把所构造的常用的其它情况下的薄板单元

我们把它列到这

我们看看,这是一个ACM非协调单元

它是4节点的,每个节点取3个自由度

它一共有12项,有12个待定系数,从a1到a12

它的位移模式是这样一个关系

这是BFS-16协调单元

它是4节点,每个节点取的自由度是

另外它还加了一个

这个里面的插值函数关系要取为Hermite插值关系

BFS-24高精度的协调元

它每个节点的自由度取为

然后把两次偏导数取全了

也就是说刚才只取了一项

现在把两阶偏导也取上

每一个节点一共有6个自由度

所以4个节点一共24个自由度

同样也选取Hermite插值多项式

3节点三角形薄板单元

有Tocher单元,我们叫T-9单元

这就是我们前面介绍过的

每个节点有3个自由度

同样我们取9项

它不是一个完整的三次项

也就是说这两项是合到一块,把它耦合起来

Tocher(T-10)非协调元

它取10项,每个节点是3个自由度

实际上它只有9个自由度,但是取了10项

它利用里兹法把a9压缩掉

Adini(A-9)非协调元

每个节点3个自由度

它在完整的三次多项式里面

把交叉项xy忽略掉

原来是10项,现在就取为9项

cowper(C-18)协调元

每个节点自由度是取6个

3个节点,所以一共是18个自由度

但是它取了21项

它完整地取了一次项、二次项、三次项

四次项还有五次项

一共取了21项

那么它通过施加3个约束方程

把21个待定系数压缩为18个

这样就可以满足唯一确定性的原则

Finite Element Method (FEM) Analysis and Applications课程列表:

0、Course summary

-Finite element, infinite capabilities

--Video

1、Introduction

-1.1 Classification of mechanics:particle、rigid body、deformed body mechanics

--Video 1.1

--1.1 Test

-1.2 Main points for deformed body mechanics

--Video 1.2

--1.2 Test

-1.3 Methods to solve differential equation solving method

--Video 1.3

--1.3 Test

-1.4 Function approximation

--Video 1.4

--1.4 Test

-1.5 Function approximation defined on complex domains

--Video 1.5

--1.5 Test

-1.6 The core of finite element: subdomain function approximation for complex domains

--Video 1.6

--1.6 Test

-1.7 History and software of FEM development

--Video 1.7

--1.7 Test

-Discussion

--Discussion

-Homework

2、Finite element method of bar system based on direct stiffness method

-2.1 Principles of mechanic analysis of springs

--Video 2.1

--2.1 Test

-2.2 Comparison between spring element and bar element

--Video 2.2

--2.2 Test

-2.3 Coordinate transformation of bar element

--Video 2.3

--2.3 Test

-2.4 An example of a four-bar structure

--Video 2.4

--2.4 Test

-2.5 ANSYS case analysis of four-bar structure

--Video 2.5

--ANSYS

-Discussion

--Discussion

3、Mechanical description of deformed bodies with complex geometry 1

-3.1 Mechanical description and basic assumptions for deformed body

--Video 3.1

--3.1 Test

-3.2 Index notation

--Video 3.2

--3.2 Test

-3.3 Thoughts on three major variables and three major equations

--Video 3.3

--3.3 Test

-3.4 Test

-3.4 Construction of equilibrium Equation of Plane Problem

--Video 3.4

-3.5 Test

-3.5 Construction of strain-displacement relations for plane problems

--Video 3.5

-3.6 Test

-3.6 Construction of constitutive relations for plane problems

--Video 3.6

-3.7 Test

-3.7 Two kinds of boundary conditions

--Video 3.7

- Discussion

-- Discussion

4、Mechanical description of deformed bodies with complex geometry 2

-4.1 Test

-4.1 Discussion of several special cases

--Video 4.1

-4.2 Test

-4.2 A complete solution of a simple bar under uniaxial tension based on elastic mechanics

--Video 4.2

-4.3 Test

-4.3 The description and solution of plane beam under pure bending

--Video 4.3

-4.4 Test

-4.4 Complete description of 3D elastic problem

--Video 4.4

-4.5 Test

-4.5 Description and understanding of tensor

--Video 4.5

-Discussion

--Discussion

5、Principle of trial function method for solving mechanical equations of deformed body

-5.1 Test

-5.1Main method classification and trial function method for solving deformed body mechanics equation

--Video 5.1

-5.2 Test

-5.2 Trial function method for solving pure bending beam: residual value method

--Video 5.2

-5.3 Test

-5.3How to reduce the order of the derivative of trial function

--Video 5.3

-5.4 Test

-5.4 The principle of virtual work for solving plane bending beam

--Video 5.4

-5.5 Test

-5.5 The variational basis of the principle of minimum potential energy for solving the plane bending

--Video 5.5

-5.6 Test

-5.6 The general energy principle of elastic problem

--Video 5.6

-Discussion

--Discussion

6、Classic implementation and finite element implementation based on trial function method

-6.1Test

-6.1 Classic method and finite element method based on trial function

--Video 6.1

-6.2 Test

-6.2 Natural discretization and approximated discretization in finite element method

--Video 6.2

-6.3 Test

-6.3 Basic steps in the finite element method

--Video 6.3

-6.4 Test

-6.4 Comparison of classic method and finite element method

--VIDEO 6.4

-Discussion

--Discussion

7、Finite element analysis of bar and beam structures

-7.1 Test

-7.1 Construction and MATLAB programming of bar element in local coordinate system

--Video 7.1

-7.2 Test

-7.2 Construction and MATLAB programming of plane pure bending beam element in local coordinate syste

--Video 7.2

-7.3 Construction of three-dimensional beam element in local coordinate system

--Video 7.3

-7.4 Test

-7.4 Beam element coordinate transformation

--Video 7.4

-7.5 Test

-7.5 Treatment of distributed force

--Video 7.5

-7.6 Case Analysis and MATLAB programming of portal frame structure

--Video 7.6

-7.7 ANSYS case analysis of portal frame structure

--Video 7.7

8、Finite element analysis of continuum structure (1)

-8.1 Test

-8.1 Two-dimensional 3-node triangular element and MATLAB programming

--Video 8.1

-8.2 Test

-8.2 Two-dimensional 4-node rectangular element and MATLAB programming

--Video 8.2

-8.3 Test

-8.3 Axisymmetric element

--Video 8.3

-8.4 Test

-8.4 Treatment of distributed force

--Video 8.4

-8.5 MATLAB programming of 2D plane rectangular thin plate

--Video 8.5

-8.6 Finite element GUI operation and command flow of a plane rectangular thin plate on ANSYS softwar

--Video 8.6

-Discussion

--Discussion

9、Finite element analysis of continuum structure (2)

-9.1 Three-dimensional 4-node tetrahedral element and MATLAB programming

--Video 9.1

-9.2 Three-dimensional 8-node hexahedral element and MATLAB programming

--Video 9.2

-9.3 Principle of the isoparametric element

--Video 9.3

-9.4Test

-9.4Numerical integration

--Video 9.4

-9.5 MATLAB programming for typical 2D problems

--Video 9.5

-9.6 ANSYS analysis case of typical 3Dl problem

--Video 9.6

-Discussion

--Discussion

10、Basic properties in finite element method

-10.1Test

-10.1Node number and storage bandwidth

--Video 10.1

-10.2Test

-10.2 Properties of shape function matrix and stiffness matrix

--Video 10.2

-10.3Test

-10.3 Treatment of boundary conditions and calculation of reaction forces

--Video 10.3

-10.4Test

-10.4 Requirements for construction and convergence of displacement function

--Video 10.4

-10.5Test

-10.5C0 element and C1 element

--Video 10.5

-10.6 Test

-10.6 Patch test of element

--Video 10.6

-10.7 Test

-10.7 Accuracy and property of numerical solutions of finite element analysis

--Video 10.7

-10.8Test

-10.8 Error and average processing of element stress calculation result

--Video 10.8

-10.9 Test

-10.9 Error control and the accuracy improving method of h method and p method

--Video 10.9

-Discussion

--Discussion

11、High-order and complex element

-11.1 Test

-11.1 1D high-order element

--Video 11.1

-11.2 Test

-11.2 2D high-order element

--Video 11.2

-11.3 Test

-11.3 3D high-order element

--Video 11.3

-11.4 Test

-11.4 Bending plate element based on thin plate theory

--Video 11.4

-11.5 Test

-11.5 Sub-structure and super-element

--Video 11.5

12、Introduction to the application of finite element analysis (1)

-12.1Test

-12.1 Finite element analysis for structural vibration: basic principle

--Video 12.1

-12.2 Test

-12.2 Case of finite element analysis for structural vibration

--Video 12.2

-12.3 Test

-12.3 Finite element analysis for elastic-plastic problems: basic principle

--Video 12.3

-12.4 Test

-12.4 Finite element analysis for elastic-plastic problems: solving non-linear equations

--Video 12.4

-Discussion

--Discussion

13、Introduction to the application field of finite element analysis (2)

-13.1 Test

-13.1 Finite element analysis for heat transfer: basic principle

--Video 13.1

-13.2 Test

-13.2 Case of finite element analysis for heat transfer

--Video 13.2

-13.3 Test

-13.3 Finite element analysis for thermal stress problems: basic principle

--Video 13.3

-13.4 Test

-13.4 Finite element analysis for thermal stress problems: solving non-linear equation

--Video 13.4

-Discussion

--Discussion

14、Project

-2D problem: finite element analysis of a 2D perforated plate

--Video I-1

-3D problem: meshing control of a flower-shaped chuck

--Video I-2

-Modal analysis of vibration: Modal analysis of a cable-stayed bridge

--Video I-3

-Elastic-plastic analysis: elastic-plastic analysis of a thick-walled cylinder under internal pressur

--Video I-4

-Heat transfer analysis: transient problem of temperature field during steel cylinder cooling process

--Video I-5

-Thermal stress analysis: temperature and assembly stress analysis of truss structure

--Video I-6

-Probability of structure: Probabilistic design analysis of large hydraulic press frame

--Video I-7

-Modeling and application of methods: Modeling and analysis of p-type elements for plane problem

--Video I-8

Video 11.4笔记与讨论

也许你还感兴趣的课程:

© 柠檬大学-慕课导航 课程版权归原始院校所有,
本网站仅通过互联网进行慕课课程索引,不提供在线课程学习和视频,请同学们点击报名到课程提供网站进行学习。