力学中的泛函分析和变分原理第六讲
- 格式:pdf
- 大小:439.67 KB
- 文档页数:10
七、弹性力学参量的下标表示法前面给出的力分量、应力分量、应变分量和位移分量,其表示方法引用的是记号法;这是一种公认的弹性力学参量表示方法。
下标表示法书写简洁,便于力学问题的理论推导。
1. 下标符号具有相同性质的一组物理量,可用一个带下标的字母表示:如:位移分量u, v, w 表示为u 1, u 2, u 3,缩写为u i (i =1,2,3)坐标x, y, z 表示为x 1, x 2, x 3,缩写为x i (i =1,2,3)单位矢量i, j, k 表示e i (i =1,2,3)。
体力分量X, Y, Z 表示为X 1, X 2, X 3,缩写为X i (i =1,2,3)应力分量:z zy zx yz y yxxz xy x 可表示为:333231232221131211 缩写为:)3,2,1;3,2,1( j i ij4. 克罗内克(Kroneker)符号具有如下性质 )cos(j i ij e ej i e eji ji ij 01 100010001333231232221131211 ij ij (1)3ii j i ij A A ij 也称换名算子同理:ijkj ik A a (2)选取可能位移:十、利用位移变分原理的近似解法m mm m mm mm m w C w w v B v v u A u u 000其中系数是完全任意的m m m C B A 、、1、瑞雷—里兹法(1)是在边界上满足位移边界条件的设定函数000w v u 、、(2)是在边界上为零的设定函数m m m w v u 、、可见,由(1)、(2)选取出来的是可能位移w v u 、、。
第二章 泛函极值及变分法(补充内容)2.1 变分的基本概念2.1.1 泛函和变分泛函是一种广义的函数,是指对于某一类函数{y (x )}中的每一个函数y (x ),变量J 有一值与之对应,或者说数J 对应于函数y (x )的关系成立,则我们称变量J 是函数y (x )的泛函,记为J [y (x )]。
例1:如果表示两固定端点A (x A ,y A ),B (x B ,y B )间的曲线长度J (图2.1.1),则由微积分相关知识容易得到:dx dx dy J BAx x ⎰+=2)/(1 (2.1.1)显然,对于不同的曲线y (x ),对应于不同的长度J ,即J 是函数y (x )的函数,J =J [y (x )]。
图2.1.1 两点间任一曲线的长度例2:历史上著名的变分问题之一——最速降线问题,如果2.1.2所示。
设在不同铅垂线上的两点P 1与P 2连接成某一曲线,质点P 在重力作用下沿曲线由点P 1自由滑落到点P 2,这里不考虑摩擦作用影响,希望得到质点沿什么样的曲线滑落所需时间最短。
图2.1.2 最速降线问题选取一个表示曲线的函数y (x ),设质点从P 1到P 2沿曲线y =y (x )运动,则其运动速度为:dsv dt ==其中,S 表示曲线的弧长,t 表示时间,于是:dt =设重力加速度为g ,则gy v 2=。
因为P 1和P 2点的横坐标分别为x 1到x 2,那么质点从P 1到P 2所用时间便为:1[()]x x J y x =⎰211/2211[()]2[()()]x x y x dx g y x y x ⎧⎫'+=⎨⎬-⎩⎭⎰(2.1.2)则最速降线问题对应于泛函J [y (x )]取最小值。
回顾函数的微分:对于函数的微分有两种定义: 一种是通常的定义,即函数的增量:),()()()(x x x x A x y x x y y ∆+∆=-∆+=∆ρ (2.1.3) 其中A (x )与∆x 无关,且有∆x →0时ρ(x ,∆x )→0,于是就称函数y (x )是可微的,其线性部分称为函数的微分()()dy A x x y x x '=∆=∆,函数的微分就是函数增量的主部。
变分原理泛函是指某一个量,它的值依赖于其它一个或者几个函数。
因此泛函也称为函数的函数。
变分法的基本问题是求解泛函的极值。
对于弹性力学问题,根据能量关系可以使偏微分方程的边值问题转化为代数方程。
弹性体的应变能是基本未知量应力或者应变分量的函数,当然应力或者应变分量是坐标的函数。
因此,应变能就是泛函。
在数学分析中,讨论函数和函数的极值。
变分法讨论泛函的极值,是极值问题的推广。
下面简单介绍复变函数的定义和基本性质。
如果需要深入探讨复变函数问题,请查阅参考资料。
§1 泛函和泛函的极值首先引入泛函的概念。
泛函是指某一个量,它的值依赖于其它一个或者几个函数。
因此泛函也称为函数的函数。
变分法的基本问题是求解泛函的极值作为变分法的简单例题。
考察x,y平面上连接两个定点的所有曲线中,求满足边界条件的任意曲线y(x)中的最短曲线。
(补充图)设P1(x1,y1)和P2(x2,y2)为平面上给定的两点,y(x)为连接两点的任意曲线。
于是,这一曲线的长度为连接P1,P2两点的曲线有无数条,每一条曲线都有一个L值与其对应。
满足边界条件的y(x)称为容许函数,问题是要从这些曲线,容许函数中找出使得曲线长度L最小的一条。
根据上式,L [y]依赖于y(x),而y(x)是x的函数,因此称y(x)为自变函数;L [y]是倚赖于自变函数的函数,称为泛函。
求解最短程线问题,即在满足边界条件在x=x1时,y(x1)=y1,y'(x1)= y'1在x=x2时,y(x2)=y2,y'(x1)= y'2的函数y(x)中,求使得泛函L [y]为极值的特定函数。
因此y(x)称为容许函数。
上述问题应用变分法可以概括为求解泛函在边界条件y(x1)=y1,y(x2)=y2的极小值问题。
§2 泛函极值的必要条件-欧拉方程假设函数y(x)是使得泛函L [y]为最小的特定函数(真实的)。
变分法有兴趣研究的是邻近于y(x)的任意容许函数引起泛函L [ ]的改变。
112第六章 能量泛函的转换形式及其应用§6.1 总位能泛函转换形式及其应用由§4.1节中的(4-16)式,定义了总位能泛函,即⎰⎰σ--ε=∏S i i Vi i ij S u T V u F A d d ])([P (4-16)该泛函为单变量变分原理,其自变量要求满足位移应变关系及位移边界条件,即)(21,,i j j i ij u u +=ε 0=-i i u u所以,这种变分原理是有条件的,并可以进一步证明总位能原理是极小值原理,解的收敛性得到保证。
这种原理是目前广为流行的绝大部分有限元素模型的基础,比较理想的情形是“保续元”的建立,而放松某些边界协调条件则构成了有限元素法中的“非保续元”。
【例1】 梁元素的总位能泛函及其变换。
图6-1所示的一维梁,承受横向分布载荷)(x p ,简支端(L x =)作用一集中力矩M ,梁的另一端为固持。
显然,其边界条件为0=x :0)0()0(='=w w L x =:0)(=L w ,及M L M =)( 6-1)总位能泛函根据定义可写为 V U +=∏p (6-2)其中⎰''=L x w EJ U 02d )(21 (应变能) (6-3) )(d 0L w M x w p V L '+-=⎰ (外力位能) (6-4)上面各式中,w 表示挠度,它是坐标x 的函数,而w '与w ''分别代表x w d d 及22d d x w。
现在对总位能取一阶变分,)(δd δd δδδδ0p L w M x w p x w w EJ V U LL '--''''=+=∏⎰⎰ (6-5)当弯曲刚度EJ 沿长度不变时,可将它放在积分号之前,再利用Green 公式,可得[][]⎰⎰+'''-'''=''''LLLL x w EJw w w EJ w w EJ x w w EJ 0)4(000d δδδd δ (6-6)将(6-6)式代入(6-5)式中,利用条件(6-1)式,整理后可得图6-1 一维弯曲梁113w M w EJ x w p EJw L x L'-''+-=∏=⎰δ][d δ][δ0)4(p (6-7)现令(6-7)的0δp =H ,利用变分法中的预备定理,可得到0)4(=-p EJw (6-8)0=-''=M w EJ L x (6-9)(6-8)式即为平衡方程,与材料力学所导出的公式完全一致,(6-9)式为力的边界条件,即相当于(6-1)式中的最后一个公式。