当前课程知识点:Finite Element Method (FEM) Analysis and Applications >  4、Mechanical description of deformed bodies with complex geometry 2 >  4.3 The description and solution of plane beam under pure bending >  Video 4.3

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

Video 4.3在线视频

下一节:Video 4.4

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

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

下面讨论平面纯弯梁的描述

我们首先考虑一个简支梁

这个简支梁是一个长度为l,高度为h,厚度为b

那么上面作用有均布载荷,叫p一横

左边和右边,简支

那么我们可以用两种建模方法来描述它

第一种是通用建模方法

也就是前面我们介绍的内部加表面

内部呢可以取出一个微元体

也就是dxdy_t

这个地方呢,厚度就是b

来进行微元体的分析

来建立最一般的三大类变量、三大类方程

这就是我们通用的建模方法

另外呢,由于我们这个对象

假定是一个细长梁,又是一个平面梁

那这种情况就是等截面

那么我们可以针对它的特征

又是细长、又是等截面

这样的话,来进行一个等效的建模

我们把这种等效的建模叫特征建模方法

它是一种简化方法

那么这种特征呢就是细长

那么我们只用沿着x坐标来描述来刻画它

这个自变量就变成一个只变x

那同样我们的这个变形

它主要是在承受载荷以后呢它就是一个弯曲

弯曲呢就是有一个挠度

那么这个挠度呢,它就是关于y方向的一个分量

所以我们就完全可以只用挠度来描述它的位移场

当然我们要进行这样的特征就要进行简化

这个简化呢就是小变形、还有直法线假定

那么下面我们用特征建模的方法来进行一个分析

那么对于这个简支梁

我们从中间,x坐标的一个位置处,切出dx这一段

那么我们把这一段切出来

我们通透的切出来,也就是沿高度方向完整的给它切下来

就是上面还带着有分布力,也就是px一横

那么在切下来的两个侧面上

我们把力合成一下

就是说它有一个弯矩,还有一个剪力

同样随着dx的变化它有一个弯矩

和剪力它有一个增量

当然我们要假定

它的应力呢沿着高度方向是一个线性分布

实际上这就是我们后面要讲的直法线假定

那么我们把这个特征建模叫做

具有全高度梁的dx_b_h

什么意思呢

dx就是取出这么一个宽度

_b呢就是厚度方向完全把它取下来

这个_h呢就是沿高度方向完全把它取下来

那这样我们就用三大类变量来描述它

这样就比较简单了

我们在位移变量里面我们就用中性层

也就是说中间没有拉伸也没有压缩的这么一层

作为它的位移的,y方向挠度的这么位移的描述变量

那么我们把这个v写成x的函数

同时要指明一下它关于y呢是y带三角这个是等于零

它是表明是中性层

同样应力呢,我们也是σx,xy的函数

εx,是xy的函数

那么其他的应力分量很小,我们一般不考虑

还有呢就是要满足直法线假定

我们后面要给大家一个具体的推导

那同样我们取出来的这个微段

来推导三大类方程

同样它也是一个微元体

这个微元体呢也有一个x方向的合力平衡

y方向的合力平衡,还有力矩平衡

那么在x方向的合力平衡的话

我们就可以写出来

因为我们的应力σx和弯矩呢要达到一个平衡

所以我们就取出一个应力针对中性层求矩

得到的弯矩的平衡,这是x方向的平衡

在y方向,我们同样,Q,剪力

还有在剪力右边还有一个剪力的增量

与上面的压力,上面的这个分布力进行一个平衡

我们就可以得到dq加上一个px,乘上一个dx等于零

那这样的话我们就可以得到

这样的一个关于剪力的平衡方程

那同样我们对任意一点求力矩的平衡

那同样我们可以得到dm减去Q乘dx要等于零

那这样的话我们得到一个Q等于dmdx

这是力矩的平衡

在几何方程这一块

那么我们取出一个微元体

这个微元体我们假定一个中性层

把坐标取在中性层

基于中性层的这个坐标我们叫作y一个三角

假定它的变形呢是沿着有一个圆心做了一个弯曲的变化

也就是说做了一个圆的这么一个变形

这个圆的变形呢,在里面的纤维我们取出一层

也就是说在y这个位置呢取出一层

取出这么一层纤维

那么变形前的这一层呢是中性层的这一层

它是叫ds

这一层呢是变形以后

那么我们来求一下相对于y这一层的相对伸长量是多少

那么它就等于这一层的长度减去变形前的长度

因为变形前是直的,就是ds

再除上一个变形前的长度

这样我们叫相对伸长量

这样的话我们就做一个

这一层的相对于转动的这个圆心的它的转动的半径

是R减上y

然后再乘上一个变形的角度呢

是弯曲的转动的角度呢是dθ

那这是变形以后的长度

变形以前的长度呢是因为它的半径它就是R

R乘上一个dθ

那再除上一个原长,就是R乘上dθ

这样我们就得到了相对伸长量就是εx它的一个表达

那同样,刚才说的,它的变化

实际上这个R呢就是一个曲率半径

那这个曲率K呢与曲率半径的关系

那么它是一个倒数的关系

那么我们由高等数学

当这个梁在中性层它有一个挠度函数的时候

那它的曲率的计算公式

那么我们可以表达成这样一个公式

那么由于我们是小变形

把高阶量去掉以后

我们就可以得到一个曲率就是关于挠度的二阶导数

那么我们把这个曲率代到几何方程里面去

也就是说它的倒数,就是曲率半径

我们放进去,就得到了几何方程的这么一个表达

那么物理方程,同样,我们也可以得到

这是一个一维的虎克定理的物理方程

总结一下前面建立的三大类方程

那么我们进行化简归纳以后

我们就可以得到y方向的平衡方程

就是[方程]

那x方向的平衡方程呢就是那个力矩和应力的平衡

我们把关于挠度的函数代进去以后

就是关于挠度的二阶导数

那几何方程同样我们也可以得到[方程]

那么物理方程呢,前面再乘上一个E

注意一下这个里面,这个EI,这个I呢

它实际上是关于梁的截面的惯性矩

关于中性轴,我们既不伸长、又不缩短的

那个中性层的垂直的那个轴的惯性矩

那么可以看出来,这些方程变量也很多

但实际上我们从第一个方程可以看出来

那px这是我们给的外部的条件,给的一个分布载荷

那么实际上它的变量就是v

那么我们可以看出来

从这个方程里面就可以直接独立的针对v进行求解

当v求解以后,就可以求出

同样这个截面的弯矩、应变还有应力

所以这个方程呢应该说是一个解耦的方程

好,下面我们对简支梁给出边界条件

那么看起来,也有位移边界条件,也有力的边界条件

那么位移的边界条件,左边固定

那也就是说x等于零的时候,v就等于零

右边呢也是垂直方向是固定的

所以说它x等于l的时候,v也是等于零

关于力的边界条件,看起来是上面的分布载荷

看起来好像是力的边界

但是注意一下,我们刚才前面建的模啊是进行特征建模

也就是说,中间我们取的微段,是把上下取透的

也就是说把px也取到那个微段里面去了

也就是说这个px已经在那个平衡方程里面已经体现

所以这个时候px就不是边界上的力了

那么我们的边界在简支梁里面就是左端点和右端点是边界

中间都不叫边界

那么对于力的边界,左边,由于承受的外加的弯矩为零

对于右边也是,外加的施加的弯矩为零

所以说它的真正的力的边界条件是M也就是弯矩

在x等于零的时候它也等于零

M弯矩在x等于l的时候它也等于零

那么我们看看前面在建立方程的时候

其中有一个关于弯矩的,也就是x方向平衡方程的时候

有一个弯矩和挠度的这么一个关系

也就是弯矩等于零,完全可以写成这么一个挠度的表达

那这样就是它等于零

就是挠度二阶导数在x等于零的时候要等于零

挠度的二阶导数在x等于l的时候它要等于零

好我们再把简支梁所有的方程和边界条件

我们再简单的归纳一下

那么我们的平衡方程,也叫控制方程

我们解耦以后,我们就是一个

就是关于挠度的四阶导数乘上EI要等于分布力

假定我们要求解的话这个分布力是一个均布载荷

那么这个分布力就是一个常数

那么边界条件是左端点固定,挠度等于零

右端点也是固定,它挠度等于零

左端点的弯矩等于零,那么它的二阶导数为零

右端点的弯矩为零,也是它的二阶导数为零

那么要求这么一个一维的四次方程

它是一个常微分方程,所以说也比较简单

可以看出来它就是一个四次多项式

那么四次多项式呢,我们再加上四个边界条件

这样的话我们就可以把四个待定系数把它确定下来

最后呢,我们得到的这个结果是这样的

它的挠度,简支梁中点的挠度

也就说是在x等于二分之l的时候

它的挠度的计算结果是这样的

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 4.3笔记与讨论

也许你还感兴趣的课程:

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