理论力学10 质点运动微分方程
- 格式:ppt
- 大小:688.50 KB
- 文档页数:21
变质量动力学曾凡林哈尔滨工业大学理论力学教研组本讲主要内容1、变质量质点的运动微分方程2、变质量动力学在火箭发射中的应用3、变质量质点的动力学普遍定理1、变质量质点的运动微分方程(1) 变质量质点的运动微分方程m 在时刻t ,质点的质量为m ,速度为vv 1在时刻t+d t ,并入速度为v 1的微小质量d mm +d m v 并入后,系统质量变为m +d m ,速度变为v +质点系在t 瞬时的动量:11d m m =+×p v v t +d t 质点系在t+d t 瞬时的动量:2(d )(d )m m =++p v v 根据动量定理有:(e)21d d t=-=p p p F (e)1d d d d d d m m m m t+×+×-×=v v v v F 略去高阶微量d m ·d v ,并在等式两边同时除以d t , 得:(e)1d d ()d d m m t t --=v v v F 式中v 1-v=v r 为微小质量在并入前相对于质点m 的相对速度, 令d d r m t f =F v 则有:(e)d d m tf =+v F F —变质量质点的运动微分方程方程形式与常质量质点运动微分方程相似,仅在右端多了一项F ϕ,它具有力的量纲,常称为反推力。
当d m /d t >0 时,F ϕ与v r 同向;当d m /d t <0 时,F ϕ与v r 反向。
1、变质量质点的运动微分方程(2) 常用的几种质量变化规律i 质量按线性规律变化1)1(0<-=t t m m b b ,由知,其反推力为:b 0d d m t m-=r 0rd d mm t f b ==-F v v 当v r 为常量时,反推力也为常量,且与v r 方向相反。
ii 质量按指数规律变化tm m b -=e 0由知,其反推力为:0d d t m m e t b b -=-r 0rd d tmm e t b f b -==-F v v 令a ϕ表示仅在反推力F ϕ作用下变质量质点的加速度,则:0rrtt m e m m e b f f b b b ---===-F v a v 当v r 为常量时,a ϕ也为常量,即由反推力引起的加速度为常量。
【郑重说明】《理论力学》课程的习题及解答方面的参考书很多,学习者可以通过各种形式阅读与学习,按照学院对教学工作的要求,为了满足学习者使用不同媒体学习的实际需要,通过各种渠道收集、整理了部分习题及参考解答,仅供学习者学习时参考。
由于理论力学的题目解答比较灵活,技巧性也比较强,下面这些解答不一定是最好的方法,也可能会存在不够完善的地方,希望阅读时注意之。
学习理论力学课程更重要的是对物理概念的掌握与理解,学习处理问题的思想与方法,仅盲目的做题目或者阅读现成的答案,很难达到理想的结果。
质点动力学思考题与习题及参考解答思考题(1) 有一质量为m 的珠子, 沿一根置于水平面内的铁丝滑动, 采用自然坐标法描述. 珠子受重力g m W=, 铁丝施与的约束力b Nb n Nn t Nt Ne F e F e F F ++=.t Nt e F 即为滑动摩擦力f F, 设动摩擦因数为μ. 试判断下列各式正误: (1) mg F f μ=; (2) Nb f F F μ= (3)Nn f F F μ=;(4) 22Nb Nnf F F F +=μ(2) 用极坐标系描述单摆的运动. 某甲如思考题(2图(a)规定θ角正向, 得到动力学方程θθsin mg ml -= ; 某乙如思考题(2图(b)规定θ角正向, 则得到θθsin mg ml += . 你认为谁的做法正确?(a) (b)思考题(2图(3) 质量为m 的质点, 由静止开始自高处自由落下. 设空气阻力f F与速度成正比, 比例系数为k . 某甲建立竖直向上的坐标如思考题(3图(a), 得到方程为y k mg y m+-=. 某乙建立竖直向下的坐标如思考题(3图(b), 得到方程为y k mg y m-=.他们列出的方程对吗?(a) (b)思考题(3(4)有人认为: 用极坐标系讨论质点的平面运动时, 如果0≡r F , 则沿径向动量守恒,==rm p r 常量;若0≡θF , 则沿横向动量守恒. 这种看法对吗? (5) 试判断以下二论断是否正确:(1) 若质点对固定点O 的角动量守恒, 则对过O 点的任意固定轴的角动量守恒. (2) 若质点对固定轴的角动量守恒, 则对该轴上任一固定点的角动量守恒.(6) 一质点动量守恒, 它对空间任一固定点的角动量是否守恒? 如质点对空间某一固定点角动量守恒, 该质点动量是否守恒?(7) 当质点做匀速直线运动时, 其动量是否守恒? 角动量是否守恒?(8) 在固定的直角坐标系Oxyz 中, 质量为m 的质点的速度k v j v i v v z y x++=, 所受合力为k F j F i F F z y x ++=. 能否将质点的动能定理r F mv d )21(d 2⋅=向x 轴方向投影而得出分量方程x F mv x x d )21(d 2= 该方程是否正确?思考题解答(1) 仅(4)式正确.(2) 甲正确. 乙错在角度不可以定义为从动线指向定线.(3) 乙的方程正确. 甲错在空气阻力亦应为yk -,y 取负值,y k -取正值. (4) 仅对固定方向才有动量守恒的分量形式. 径向和横向均不是空间固定方向. (5) (1)对;(2)错. (6) 一质点动量守恒,则对空间任一固定点角动量守恒. 质点对空间某一固定点角动量守恒,其动量不一定守恒.(7) 质点作匀速直线运动时,其动量和角动量均守恒.(8) 动能定理是标量方程,不可能投影而得出分量方程. 但xF mv x d )21(d 2=是正确的. 仿照动能定理的导出,用x t v x d d =乘牛顿第二定律的x 分量方程x xF t v m=d d 即可证明.质点动力学习题及参考解答【1】研究自由电子在沿x 轴的振荡电场中的运动. 已知电场强度i t E E)cos(0ϕω+=,ϕω,,0E 为常量. 电子电量为e -, 质量为m . 初始时, 即当0=t 时i x r00=, i v v 00=. 忽略重力及阻力, 求电子的运动学方程.【解】力为时间的函数,积分两次可得)cos(200ϕωω+++=t m eE t V X x ,其中ϕωcos 2000m eE x X -=,ϕωsin 00m eE v V +=.【2】 以很大的初速度0v自地球表面竖直上抛一质点, 设地球无自转并忽略空气阻力, 求质点能达到的最大高度. 已知地球半径为R , 地球表面处重力加速度为g .【解】以地心O 为原点,建立x 轴经抛出点竖直向上. 质点受万有引力沿x 轴负方向. 所以2x GMm xm -= . 因为2R GMmmg =,故g R GM 2=. 故有22x g R x -= . 做变换)2(d d d d d d d d 2x x x x x t x x x x ===,则x x g R x d )2(d 222-= . 积分并用0=t 时R x =,0v x = 定积分常数,得到 )11()(212202R x g R v x -=- . 质点达最大高度时H R x +=,0=x,可求出 1220)21(2--=Rg v g v H .三点讨论:(1)令∞=H ,对应Rg v 20=为第二宇宙速度.(2)若Rg v 220<<,则回到重力场模型所得结果. (3)题中不考虑地球自转及空气阻力,均不大合理,试进一步讨论之.【3】 将质量为m 的质点竖直上抛, 设空气阻力与速度平方成正比, 其大小22gv mk F R =.如上抛初速度为0v , 试证该质点落回抛出点时的速率2201v k v v +=.【解】质点运动微分方程为(Oy 轴竖直向上);上升阶段22y g mk mg y m--=,下降阶段22y g mk mg ym +-=. 【4】向电场强度为E 、磁感应强度为B 的均匀稳定电磁场中入射一电子. 已知B E⊥, 电子初速0v 与E 和B 均垂直, 如题4图所示. 试求电子的运动规律. 设电子电量为e -.题4图【解】令m eB=ω,电子运动微分方程为y xω-=, (1) m eEx y-= ω, (2)0=z . (3)对(2)式求导,利用(1)式得02=+y yω,解出)sin(αω+=t A y . 0=t 时0=y 故0=α,由t A y ωωcos = ,且0=t 时m eBv Ee y0+-= ,故B Bv E A 0+-=,则t B Bv E y ωsin 0+-= . 积分得)cos 1()(20t m eB eB Bv E m y -+-=. 代入(1)式积分可得t m eB eB Bv E m t B E x sin )(20--=.【5】 旋轮线如题5图所示, 可理解为一半径为a 的圆轮在直线上做无滑滚动时轮缘上一点P 的轨迹, 其参数方程为)sin (ϕϕ+=a x , )cos 1(ϕ-=a y . 在重力场中, 设y 轴竖直向上, 一质点沿光滑旋轮线滑动, 试证质点运动具有等时性(绕O 点运动周期与振幅无关).题5图【解】(旋轮线是如图圆轮在直线AB 上作无滑滚动时P 点的轨迹,曲线上P 点切线方向即为轮上P 点速度方向. 因无滑,0P 为瞬心,故P 点切线与P P 0垂直,因此可知P 点切线与x 轴夹角为2ϕ. )以曲线最低点(0=ϕ)为自然坐标原点,弧长正方向与t e 一致. 质点运动微分方程为2sinϕmg s m -= .对曲线参数方程求微分,得ϕϕd )cos 1(d +=a x 和ϕϕd sin d a y =,所以ϕϕd 2cos 2d d d 22a y x s =+=,积分并用0=ϕ时0=s 定积分常数,得2sin 4ϕa s =. 代入质点运动微分方程消去ϕ,得到4=+s a gs ,s 作简谐振动而具有等时性. 其解为)cos(0αω+=t A s ,a g40=ω与振幅无关.【6】 一小球质量为m , 系在不可伸长的轻绳之一端, 可在光滑水平桌面上滑动. 绳的另一端穿过桌面上的小孔, 握在一个人的手中使它向下做匀速运动, 速率为a , 如题【6图所示. 设初始时绳是拉直的, 小球与小孔的距离为R , 其初速度在垂直绳方向上的投影为0v . 试求小球的运动规律及绳的张力.题6图【解】小球运动微分方程为T F r r m -=-)(2θ , (1) 0)2(=+θθr r m , (2)a r-= . (3) 由(3)式求出at R r -=,代入(2)式求出)/(0at R t v -=θ,再由(1)式求出3220)(at R R mv F T -=.【7】 一质量为m 的珠子串在一半径为R 的铁丝做成的圆环上, 圆环水平放置. 设珠子的初始速率为0v , 珠子与圆环间动摩擦因数为μ, 求珠子经过多少弧长后停止运动 (根据牛顿第二定律求解).【解】珠子的运动微分方程为2b 2n d d N N F F t v m+-=μ, (1)n 2/N F mv =ρ, (2)mg F N -=b 0, (3)R =ρ(约束方程). (4)把(2)、(3)、(4)式代入(1)式,作变换sv t v d /)21(d d d 2=,可求出]/)ln[()2/(224020Rg g R v v R s ++=μ.【8】 质量为m 的小球沿光滑的、半长轴为a 、半短轴为b 的椭圆弧滑下, 此椭圆弧在竖直平面内且短轴沿竖直方向. 设小球自长轴端点开始运动时其初速度为零. 求小球达到椭圆弧最低点时对椭圆弧的压力 (根据牛顿第二定律求解). 【解】以椭圆最低点为自然坐标原点O ,弧长正方向指向小球初始位置,θ为切向与水平方向的夹角,小球的运动微分方程为θsin mg vm -= , (1) θρcos /2mg F mv N -=. (2)Oy 竖直向上,将s y d /d sin =θ代入(1)式得s y g s v v d /d d /d -=,积分可求出小球达最低点时gb v 22=. 由轨道方程22x a a by --=求出当0=x 时0='y ,2/a b y ='',由公式可求出22/32)1(1a b y y ='+''=ρ. 再由(2)式求出0=θ时)/21(/cos 22a b mg mv mg F N +=+=ρθ.【9】 力1F 和2F分别作用在长方体的顶角A 和B 上, 长方体的尺寸和坐标系如题【9图所示. 试计算1F 和2F对原点O 及3个坐标轴的力矩.题9图【解】11bF M x =,11aF M y -=,01=z M ,2222/b a bcF M x +=,2222/b a acF M y +-=,02=z M .【10】 已知质量为0m 的质点做螺旋运动, 其运动学方程为t r x ωcos 0=, t r y ωsin 0=,kt z =,k r ,,0ω为常量. 试求: (1)t 时刻质点对坐标原点的角动量;(2) t 时刻质点对过),,(c b a P 点, 方向余弦为),,(n m l 的轴的角动量.【解】由运动学方程求出→v ,根据定义即可求出→→→→→→++--=⨯=k r m j t t t r km i t t t r km v r m L ωωωωωωω200000000)sin (cos )cos (sin ,)]cos ()sin )([(]cos )()sin ([000000),,(a t r k t r c kt m m t r c kt b t r k l m L n m l -+-----=ωωωωωω)sin cos (00200t br t ar r n m ωωωωω--+.【11】 如题【11图所示, 质量为m 的小球安装在长为l 的细轻杆的A 端, 杆的B 端与轴21O O 垂直地固连. 小球在液体中可绕21O O 轴做定轴转动, 轴承1O 和2O 是光滑的. 转动中小球所受液体阻力与角速度成正比, ωαm F R =,α为常量. 设初始角速度为0ω,试求经多少时间后, 角速度减小为初始值的一半,以及在这段时间内小球所转圈数.(忽略杆的质量及所受阻力.)题 11图【解】由对21O O 轴的角动量定理ωαωm l ml t -=)(d d2,积分可得lt /0e αωω-=,求出α/)2ln (l t =. 将角动量定理化为l /d d θαω-=,积分可以求得αωαωθπ4/)r a d (2/00l l ==(圈)【12】 质量为m 的质点沿椭圆轨道运动, 其运动学方程为kt a x cos =, kt b y sin = (k b a ,,为常量). 用两种方法计算质点所受合力在0=t 到k t 4π=时间内所做的功.【解】(1)由动能定理)(4121212222122b a mk mv mv W -=-=.(2)用曲线积分算⎰⎰+=⋅=→→2121)d d (y ym x x m r d F W ,把轨道参数方程kt b y kt a x sin ,cos ==代入,则曲线积分化为对t 的积分,可得同样结果.【13】 试用动能定理求解7题.【解】珠子的动能定理为sF F mv N N d )21(d 2b 2n 2--=μ,参见3.7提示【14】 有一小球质量为m , 沿如题【14图所示的光滑的水平的对数螺旋线轨道滑动. 螺旋线轨道方程为θa e r r -=0, a 为常数. 已知当极角0=θ时,小球初速为0v . 求轨道对小球的水平约束力N F 的大小. (用角动量及动能定理求解, 图中δ为θe 与v 方向间夹角,a =δtg.)题14图【解】因机械能守恒,小球动能不变,因此0v v =.过O 点作z 轴竖直向上(垂直纸面向外),质点对z 轴的角动量δcos rmv L z =. 质点所受对z 轴力矩δsin N z rF M -=. 由对z 轴的角动量定理得δδsin )cos (d d0N rF rmv t -=.由于θθθθθ ar ar t r r v a r -=-===-e d d d d 0,θθ r v =. 故a v v r =-=θδtan . 将它代入角动量定理方程,得到N N arF rF rmv -=-=δtan 0 . 而δδsin sin 0v v v r r -=-== ,所以θδδδa N a r mv a r mv ar mv ar mv F e 11tan 1tan sin 2020220222020+=+=+==.【15】 已知质点所受力F 的3个分量为z a y a x a F x 131211++=,z a y a x a F y232221++=, z a y a x a F z 333231++=,系数)3,2,1,(=j i a ij 都是常量. 这些ij a 满足什么条件时与力F相关的势能存在? 在这些条件被满足的条件下, 计算其势能.【解】当0=⨯∇→F 时势能存在,要求311332232112,,a a a a a a ===. 以原点为势能零点,则)222(21132312233222211xz a zy a xy a z a y a x a V +++++-=.【16】 一带有电荷q 的质点在电偶极子的场中所受的力为3c o s 2r pq F r θ=,3sin r pq F θθ=,p 为偶极距, r 为质点到偶极子中心的距离.试证此力场为有势场.【解】)/cos (d d d )d d (d 2r pq r F r F e r e r F r F r r θθθθθ-=+=+⋅=⋅→→→→→,故为有势场 【17】 如题17图所示, 自由质点在Oxy 平面内运动, 静止中心A 和B 均以与距离成正比的力吸引质点M , 比例系数为k . 试证明势能存在并求出质点的势能.v题【17图【解】y ky x kx y ky ky x b x k b x k r F d 2d 2d )(d )]()([d --=--+--+-=⋅→→)](d [22y x k +-=.故势能存在. 以O 为势能零点,则)(22y x k V +=.【18】 试用机械能守恒定律求解8题.【解】根据机械能守恒定律,以椭圆弧最低点为势能零点,mgbmv =221,可知gb v 2=,参见3.8提示.【20】 将质量为m 的质点竖直抛上于有阻力的媒质中。
4.1一质点受一与距离成反比的引力作用在一直线上运动,质点的质量为m ,比例系数为k ,如此质点从距原点O 为a 的地方由静止开始运动,求其到达O 点所需的时间。
解:质点受引力为:xk F -=,其运动微分方程为:xk tm-=d d v (1)即: x k xm -=d d v v分离变量积分:⎰⎰-=x axx k m d d 0v v vxa k m ln212=v)ln(2d d xa mk tx -==v (2)(v 与x 反向,取负值) )ln00ln ),0((∞→→>∴∈xa x xa a x令:y ayex aex xa y yyd 2d )ln(22---===,代入(2)式得;mk ty aey2d d 22-=-分离变量积分:)0:0:(∞→→y a x⎰⎰=-∞t yt mk y ea 0d 2d 22t mk a22π2=故到达O 点所需的时间为: km a t 2π=4.2一质点受力3K xa x F +-=作用,求势能)(x V 与运动微分方程的解。
解:C x a x x xa x x F x V ++=+--=-=⎰⎰2232K 21d )K (d )(适当选取势能零点,使0=C ,则222K 21)(xa x x V +=机械能 =++=2222K 2121xa x xm E 常量 (1)将(1)改写成2222K 242xa x E xm --= (2)质点运动微分方程:32K xa x xm +-= 22K 22xa x xmx +-=⇒ (3)(3)+(2)得22K 44)(2x E xx x m -=+ 即0)K(K 4d d 2222=-+E x mtx (4)(4)式通解:⎪⎪⎭⎫⎝⎛++=02 K2cos K θt m A Ex当0=x时,222K 21xa x E += 解得KK K)(2max 2a EE x -+=,KK 2aEA -=所以 ⎪⎪⎭⎫⎝⎛+-+=022K2cos KK Kθt m aE E x4.3若质点受有心力作用而在圆θcos 2a r =上运动时,则5228rh ma F -=,式中m 为质量,h 为速度矩。