第5章 临界力的近似解法

第5章临界力的近似解法任课教师：强士中卫星第5章临界力的近似解法Timoshenko能量法瑞利－里兹法勃布诺夫－伽辽金法有限差分法5.1Timoshenko能量法在临界力作用下，体系从一种平衡状态过渡到另一种平衡状态时，应变能的增量与外力功的改变应该相等。ΔU＝ΔW(1)根据这一条件即可确定临界荷载。ΔUΔW：体系是稳定的ΔUΔW：体系是不稳定的如图所示受外力N作用的两端铰支杆，应变能增量和外力功的增量都是从压杆的直线平衡状态算起。应变能增量：(2)外力功的增量：(3)将(2)、(3)代入(1)，得：∴2220201212lldvUUEIdxdxMUUdxEI202lNdvWPdxdx222200122lldvNdvEIdxdxdxdx222020lcrldvEIdxdxNdvdxdx(4)按Timoshenko能量法计算临界荷载时，需要假定变形曲线；按Timoshenko能量法计算得到的临界荷载通常是近似的临界荷载；临界荷载的近似程度取决于假定的变形曲线与实际变形曲线的接近程度；通常由Timoshenko能量法计算得到的临界荷载都大于精确值。(1)设利用边界条件x=0，x=l时，v=0，得：c=0，b=-al∴v=a(x2-lx)因而：代入式(4)，求得临界荷载：精确值：误差：2vaxbxc2dvaxldx222dvadx220232220441232lTcrlEIadxEIalEINallaxldx2229.87crEIEINll%21.58%TcrcrcrNNeN(2)设满足边界条件x=0，x=l时，v=0因而：代入式(4)，积分后求得临界荷载：精确值：误差：323224qxvxlxlEI25272229.88241201720160TcrqlEIEINqllEI2229.87crEIEINll%0.13%TcrcrcrNNeN3234624dvqxlxldxEI2222dvqxlxdxEI(3)设满足边界条件x=0，x=l时，v=0因而：代入式(4)，积分后求得临界荷载：精确值：误差：sinxvAl242322222TcrAEIEIlNAll22crEINl%0.00%TcrcrcrNNeNcosdvAxdxll22222sindvAxdxllTimoshenko能量法求解临界荷载，可避免建立和求解联立方程，但需设定屈曲时变形曲线；Timoshenko能量法求解临界荷载的关键是假定的变形曲线必须合适，尽可能接近实际屈曲形式且便于计算；假定的变形曲线应尽可能地满足较多的边界条件；假定的变形曲线形状必须合理；假定的变形曲线方程必须便于积分运算；ΔU选用不同的表达形式，导致不同的近似程度(1)若则因而：临界荷载：精确值：误差：2vaxlx222222501260lNaxlxNalUdxEIEI23220226lNNalWaxldx2325210660TcrEIalalNEIl2229.87crEIEINll%1.32%TcrcrcrNNeN2MNvNaxlx(2)若则因而：临界荷载：精确值：误差：323224qxvxxllEI2225229232323014.2714102224lNqxNqlUxxlldxEIEIEI2427232322204.21631046224lNqNqlWxxlldxEIEI29.871TcrEINl2229.87crEIEINll%0.014%TcrcrcrNNeN323224NqxMNvxxllEI变形曲线的一般形式(5)式中：—满足给定边界条件的函数；ai—未知常数参变数(广义坐标)。代入式(4)：(6)为求得临界荷载，应选择ai使N为最小，即N取极小值的条件：(i=1,2,3,…,n)所以：(7)11221nnniiivaxaxaxax2012121201,,,,nliininlniiiEIaxdxAaaaNBaaaaxdx2210iiiiiABBAaaNAABaBBaBa0iiAABaBa0iiABNaaix(7)式中：同理：代入式(7)，可得：(i=1,2,3,…,n)令：(i=1,2,3,…,n)则：(i=1,2,3,…,n)(8)201nliiiiiAEIaxdxaa012nliiiiiiEIaxaxdxa012nliiiiEIaxxdx012nliiiiEIaxxdx200112nnlliiiiiiiiiBaxdxaxxdxaa010nliiiiiiaNxxEIxxdx0lijijijcNxxEIxxdx10nijiica(8)式可展开为：参数ai不能全等于零，故系数行列式应等于零，即(9)(9)式为稳定方程微分方程n个齐次线性代数方程111122121122221122000nnnnnnnnncacacacacacacacaca1112121222120nnnnnnccccccDccc5.2瑞利－里兹法瑞利－里兹法是建立在势能驻值原理基础上的近似方法。瑞利－里兹法是采用具有n个广义坐标的位移函数近似代替真实的位移曲线。瑞利－里兹法是将具有无限多个变量的泛函变分问题转化为有限多个变量的函数极值问题。设失稳时的变形曲线如下：(10)临界状态的变形状态就取决于参数a1、a2、…、an。11221nnniiivaxaxaxax总势能的变分：由稳定的能量准则，有为微小的任意值，所以：(i=1,2,3,…,n)(11)在线性问题中，势能可表示为：将式(10)代入上式，可得：(12)12112nniiniaaaaaaaa10niiiaa0iaia2222001122lldvdvEIdxNdxdxdx2200111122nnlliiiiiiEIaxdxNaxdx将式(12)代入式(11)，可得：(i=1,2,…,n)进而可得，稳定方程：(13)式中：由稳定方程，可以确定临界荷载。由解变分问题获得稳定方程的方法，称为瑞利－里兹(Rayleigh-Ritz)法010nliiiiiiaNxxEIxxdx1112121222120nnnnnnccccccDccc0lijijijcNxxEIxxdx举例：一端固定，一端自由的压杆，受轴力N作用，用瑞利－里兹法确定其临界力。设变形曲线代入势能表达式：由可得：误差：21vax22222231111001122223llEIaxdxNaxdxEIalNal31114403EIalNala23RcrEINl2222234%%21.6%4RcrcrcrEIEINNlleEINl设变形曲线代入势能表达式：令及可得：可得稳定方程：式中：展开行列式：，可解得：∴误差：2412vaxax2223121200464222llEINaaxdxaxaxdx10a2*NlNEI2222432241112212236442425357aEIlaaalalNlaalal20a224122241212203511820557EINlaEIlNlaEINlaEIlNla****2123501821557NNNN*2*451050NN*min2.4688N22.4688RcrEINl222222.46884%%0.057%4RcrcrcrEIEINNlleEINl若势能表达式取:仍取变形曲线任意截面弯矩：，其中：，故：总势能：由可得：误差：21vax25311840153NlNlaaEI22.5RcrEINl222222.54%%1.32%4RcrcrcrEIEINNlleEINl2200122llMNdvdxdxEIdxMNfv21fal221MNalx222222532121100142222153llNalxNNlNlaxdxaEIEI5.3勃布诺夫－伽辽金法勃布诺夫－伽辽金法是建立在虚位移原理基础上的近似方法；近似变形曲线应同时满足几何边界条件及静力边界条件；根据虚位移原理：(14)设位移函数：(15)则：将(15)代入(14)，则：424200ldvdvEINvdxdxdx11221nnniiivaxaxaxax11221nnniiivxaxaxaxa4242010nliiidvdvEINadxdxdx为任意值，可得伽辽金方程组：(i=1,2,…,n)(16)将式(15)代入式(16)，经积分，可得到含有n个参数ai的线性齐次代数方程；令ai的系数行列式D=0，从而确定临界荷载。伽辽金方程与里兹法方程，都满足同一能量关系，两种方法的结果一致;伽辽金方法的优点是经过积分运算直接得到线性方程组，而瑞利－里兹法须进行微分运算；424200lidvdvEINdxdxdxia伽辽金方程组的其他形式：(i=1,2,…,n)(17)或：(i=1,2,…,n)(18)举例：用伽辽金法确定一端固定一端铰支压杆临界荷载。取试解位移函数：满足几何边界条件：x=0及x=l时，v=0;x=l时，v`=0;满足静力边界条件：M0=0，即x=0时，v``=0;2200lidvEIpvxdxdx2200lidvEIMxdxdx33432valxlxx由位移函数，可得：代入伽辽金方程：积分运算后，可得：∴误差：323198vallxx21643vaxlx1683vaxl(4)148va233411048643320lEIaaNxlxlxlxxdx57620535EIlNl221GcrEINl22222210.7%%4.26%0.7GcrcrcrEIEIllNNeEINl5.4有限差分法差分法是求解微