机械系统的理论建模.docx
- 文档编号:1641148
- 上传时间:2022-10-23
- 格式:DOCX
- 页数:92
- 大小:1.42MB
机械系统的理论建模.docx
《机械系统的理论建模.docx》由会员分享,可在线阅读,更多相关《机械系统的理论建模.docx(92页珍藏版)》请在冰豆网上搜索。
机械系统的理论建模
第1章机械系统的理论建模
目前,对机械结构动态性能的研究主要有三种基本方法,即理论建模及分析方法、实验建模及分析方法和二者相结合的方法。
所谓理论建模及分析方法是基于结构动力学原理,根据结构的设计方案、图样、先验知识和资料等建立起模拟机械结构动力特性的动力学模型,而无需依赖于已有的机械设备。
通过对该动力学模型的分析计算,即可获得该机械结构各种模拟的动力特性。
这不仅可以检验其动力特性是否满足设计目标,是否需要对结构进行修改,还可以通过对理论模型的计算机仿真,预估结构设计及其改进后的动力特性或对其进行动态优化设计。
所以,理论建模及分析方法,可以在机械结构设计方案具体实施之前,建立其动力学模型,利用计算机进行模拟仿真,对各种设计方案反复进行比较、修改,使其动态特性逼近设计目标函数的要求。
从而可经济、迅速地达到优化设计的目标,把提高机械结构动态性能的问题解决在方案及图样设计阶段。
理论建模及其分析方法的不足,在于建立能够确切模拟机械结构动力特性的动力学模型较为困难,就目前的各种理论建模方法而言,都存在一些难以确定的因素,如难于对机械结构各种工况下的边界条件考虑和处理得与实际工况完全吻合,也难于把机械结构中各种结合部来的误差,影响了所建动力学模型的模拟精度。
所以,提高理论模型对机械结构动态性能的模拟精度,使之满足工程实际的需要,是机械结构动态性能理论分析方法必须解决的首要问题。
目前,对一般机械结构的理论建模,除应慎重使用集中参数法或传递矩阵法建模外,最主要的建模方法是有限元法。
因此,本章首先讨论有限元及其建模方法,包括静态有限元法和动力学问题官限单元法。
然后讨论在机械结构动力分析中常用的另外两种理论建模及分析方法:
集中参数法和传递矩阵法。
有限单元法是20世纪50年代初期根据变分原理发展起来的一种强有力的数值近似解法。
该方法以计算机为手段,采用分割近似,进而逼近整体的研究思想求解数学物理问题。
目前,有限元法已在许多领域成为分析、解决工程和数学物理问题的有力工具。
由于弹性力学问题静态有限元法的理论与方法均已成熟,有着丰富的书籍和资料,因此,本章只讨论静态有限元方程的建立,以其说明弹性力学和变分法等通用力学、数学原理与有限元法的关系,并简述静态有限元法的分析过程,直接给出常用单元的位移模式及计算公式。
若要对此部分内容作详细了解,可参见有关书籍。
本章还将讨论动力学问题有限单元建模方法。
顺便指出,本章所讨论的变分法及其近似解法,也是后续内容“机械系统动态性能的最优控制”和“动态优化”的理论基础。
1.1有限单元法的预备知识
1908年,瑞士科学家里兹首先提出用变分法处理弹性薄板问题,开创了弹性力学中直接求解泛函极值问题的近似解法。
后经人们的长期努力,形成了变分问题直接解法中最重要的里兹法。
里兹法不从微分方程出发,而是根据某个物理问题建立其泛函表达式,并根据某种泛函驻值条件,直接求解泛函极值的近似解,从而使变分解法具有了重要的实用意义。
因为人们发现并在数学上得到证明:
在连续介质问题中,许多物理、力学问题既可以转化为微分方程的定解问题,也可以归结为变分极值问题,它们的表达形式不同,但却是等价的。
例如,图1-1所示等截面梁的横向弯曲平衡方程可用材料力学方法建立,也可由变分法求得。
设梁的抗弯刚度为EJ,受分布载荷q(x)作用产生下弯变形w(x),两端固定的边界条件为
w(0)=w'(0)=0,w(l)=w'(l)=0(a)图1-l梁的横向弯曲简图
由材料力学知,等截面梁挠曲线的近似微分方程为
d2w(x)/dx2=M(x)/(EI)(b)
式中,M(x)为梁横截面上的弯矩。
若梁横截面上的剪力为Q(x),则有如下关系
dQ(x)/dx=q(x);dM(x)/dx=Q(x)
或d2M(x)/dx2=dQ(x)/dx=q(x)(c)
将式(c)代人式(b),即可得等截面梁的横向弯曲平衡方程
EId4w(x)/dx4-q(x)=0(d)
用变分原理导出梁的平衡方程如下:
在梁达到平衡时,梁和载荷作为整体的势能达到最小值。
梁的势能是梁在弯曲时所取的弯曲能,其值为
(e)
另外裁荷q(x)的势能在梁弯曲变形w(x)的位移中降低
(f)
所以,梁和载荷作为整体时的总势能为
(g)
平衡条件为总势能达到最小值,即δU=0(h)
式(g)的变分为(i)
通过分部积分,并利用固定端点条件d(δw)/dx│x=0,l=0,δw│x=0,l=0,得
(j)
根据变分法基本原理,即得梁的平衡方程为
EId4w(x)/dx4-q(x)=0(k)
可见,两种方法获得完全相同的结果。
这说明满足微分方程及其边界条件的函数将使泛函取极值或驻值,反过来,使泛函取极值或驻值的函数恰是满足微分方程及其边界条件的解答。
所以,就某个物理问题所建立的泛函极值问题的解和就该问题所建立的微分方程边值问题的解是一致的。
而且,从求近似解的角度出发,求泛函极值的近似解常常要比求微分方程的近似解更为方便[1,2]。
这就为我们实际求解工程及数学物理问题提供了一条重要途径:
直接从某个物理问题的泛函变分求其近似解,或者把微分方程的定解问题转化成相应的泛函变分求其近似解。
里兹法就属于这一类数值近似解法。
需要特别指出,我们在这里反复强调的是求解泛函极值的近似解,这是因为,从变分法出发,求解泛函极值归结为求解欧拉方程,这又回到了微分方程的求解问题。
变分法和欧拉方程代表同一个物理问题,从欧拉方程求解和从变分法求解具有相同的效果。
但欧拉方程的求解常常是困难的,甚至是不可能的,而从泛函变分求近似解往往并不困难,这也正是变分法备受重视的原因。
里兹法正是从泛函变分求极值近似解的有效方法之一。
目前常用的有限单元法均与泛函变分问题直接解法中的里兹法或伽辽金法有着紧密的联系。
里兹法和伽辽金法具有可以直接求出结构解析解的近似表达式,便于作理论分析等优点,但它们都是以结构整体为研究对象的近似计算方法,必须预先给定待求结构一个变形允许的坐标函数(基函数),而对于边界形状比较复杂的结构、要想找到合适的坐标函数却非常困难,即使找到,也需用到相当高次的多项式或某些函数的特殊结构。
因此,这种方法只能对很简单的结构求解,而对于由板梁组合而成且形状复杂的机械结构,是无能为力的[2,3]。
有限元法则与此不同。
虽然有限元法是在求解泛函变分问题的里兹法基础上发展起来的数值近似解法,但有限元法是建立在结构离散并使之单元理想化基础上的近似方法,具有结构上的分割近似性质。
因此,有限元法的基础是变分原理和剖分插值。
即首先把复杂连续结构离散为有限个形状简单的单元,如梁、三角形、四边形单元等,再对每一个单元选择一组简单的坐标函数,便可方便地将里兹法或伽辽金法应用于这些单元;然后,根据变形协调条件把这些单元重新组合,便可求得整体结构的数值解。
有限元法的数学基础是变分原理,而力学基础则是能量原理。
根据弹性力学最小势能原理,变形弹性体在受到外力作用而处于平衡状态时,在很多可能的变形允许曲线中,使总势能为最小的那条曲线是真正的变形曲线。
而弹性体的总势能是一个以位移为自变量函数,以坐标x,y,z为自变量的泛函表达式。
根据变分原理,这条真正的变形曲线就是泛函的极值曲线。
有限元法把复杂连续结构离散为有限个形状简单的单元,因此很容易将这些单元的假设变形曲线设置得很好,利用里兹法求解这些单元,便可获得这些单元真正的变形曲线,再根据各单元变形协调条件,把各单元变形曲线装配成整个弹性体的初始变形曲线,并用里兹法求解,便可获得整个弹性体的真实变形曲线。
这就是基于弹性力学能量原理,把连续弹性体的变形求解问题转化为泛函的变分极值问题求里兹解的有限元法的形成过程。
1.1.1变分法简介
有限元法的数学基础是变分法。
变分法是研究泛函极值的一种经典方法[4,5]。
有限元法是在弹性力学能量原理的基础上发展起来的一种数值计算方法,而弹性体的能量表达式是一个泛函,因此,求解弹性体的能量表达式问题可以归结为求泛函极值的问题。
1.函数与泛函
我们知道,对于变量x的某一变化域中的每一个x值,若y都有一值与之对应,则称y是x的函数,记作y=f(x)。
如果对于某一类函数{y(x)}中的每一个函数y(x),因变量J都有一个确定的值与之相对应,则称因变量J是函数y(x)的泛函数,简称泛函,记为J=J[y(x)]。
因此可以说,函数是变量与变量的关系,泛函是变量与函数的关系。
例如,图1-2所示的平面中有两固定点A(xa,ya)和B(xb,yb),连接这两点的曲线弧长l是曲线函数y=y(x)的泛函。
图l-2求弧长的变分问题
由弧长的微分知
(dl)2=(dx)2+(dy)2
dl/dx=[1+(dy/dx)2]1/2=[1+y′2]1/2
所以
A、B两点间的连线可以是各种曲线形式。
显然,对不同的曲线y(x),就有不同的长度l与之对应。
所以,A、B间连线的弧长l是曲线函数y(x)的函数,即泛函。
记为
式中,L(y')=(l+y')2。
一般情况下,L也是x、y的函数,因此可写成
(1-1)
两点间的最短弧长是直线y*(x),见图1-2。
即
lmin=J*=minJ[y(x)]=J[y*(x)]
我们把式(1-1)这种建立在函数L和变量J之间的关系叫做泛函关系。
满足此泛函关系式的任一函数y(x)均称为泛函J[y(x)]的容许函数。
2.泛函的极值与变分
求泛函的极大值或极小值问题均称为变分问题。
对于泛函J[y(x)],其自变函数y(x)的增量是指两个不同的y(x)值之差y(x)-y0(x)=Δy。
当这种增量很小时称为变分,用δy(x)或δy表示。
即
δy(x)=y(x)-y0(x)
对于自变函数y(x)的变分δy(x)所引起的泛函的增量可表示分[1,4]
ΔJ[y(x)+δy(x)]-J[y(x)]=L[y(x),δy(x)]+R[y(x),δy(x)]max│δy│(1-2)
式中,L[y(x),δy(x)]是δy的线性连续泛函;R[y(x),δy(x)]是δy(x)的高阶无穷小项。
当max│δy│→0时,R[y(x),δy(x)]→0,则L[y(x),δy(x)]就称为泛函J[y(x)的变分。
即
δJ=L[y(x),δy(x)](1-3)
由此可见,泛函的变分是泛函增量的线性主部。
若泛函有变分,而且增量ΔJ可用式(1-2)表达时,则称泛函是可微的。
泛函的变分也可定义为
δJ=∂J[y(x)+αδy(x)]/∂α│α=0(1-4)
式中,α为一正的小量。
根据式(1-4),利用函数的微分法则可方便地进行泛函变分的计算。
变分规则:
设L1和L2是x、y、y'的函数,则有如下变分规则[6]:
(1)δ(L1+L2)=δL1+δL2
(2)δ(L1L2)=L1δL2+δL1L2
(3)δ∫L(x,y,y')dx=∫δL(x,y,y')dx
(4)δy'=dδy/dx
如果泛函J[y(x)]在任何一条与y0(x)接近的曲线上所取的值不小于J[y0(x)],即
ΔJ=J[y(x)]-J[y0(x)]≥0(1-5)
则称泛函J[y(x)]在y0(x)曲线上达到极小值。
反之,若
ΔJ=J[y(x)]-J[y0(x)]≤0(1-6)
则称泛函J[y(x)]在y0(x)曲线上达到极大值。
事实上,所谓泛函极值问题就是在泛函的容许函数集中寻求一个容许函数y*(x),使泛函J[y(x)]在其上的值达到极大或极小[5,7,8]。
根据泛函变分原理,若可做泛函J[y(x)]在函数y0(x)上达到极值时,则泛函在y=y0(x)上的一阶变分等于零,即
δJ=0(1-7)
考虑(1-1)
在自变量的取值区间xa=a≤x≤xb=b
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 机械 系统 理论 建模
![提示](https://static.bdocx.com/images/bang_tan.gif)