新版本公告
~~空空如也
加载中
加载中
表情图片
评为精选
鼓励
加载中...
分享
加载中...
文件下载
加载中...

十分抱歉,这几天有事,没更……我们继续。

由于f(x)除上述情况外可以看作有3个实根,不妨设:

f(x)=2mglI1(xx1)(xx2)(xx3)(5.9)

由上文,我们知道这些实根的范围

不妨设x11x2xx31(5.10)

那么,由(5.3):

x˙=nmgl2I14(xx1)(xx2)(xx3)(5.11)

这里提出来一个4是为了后面好换元。

这个n的值是±1,用处同前。

这里说一下,对一般的拉格朗日陀螺,实际问题中只需要代入数值求解三个根即可,理论计算中,设出来就行,这几个根可以用三次方程求根公式求解,但是本身f(x)的系数就够复杂了,再带进三次方程中……反正我是不想这么做。不过保留这三个根对实际问题没有任何影响,只是多了一步求解罢了。

接下来,对(5.11)分离变量:

dx4(xx1)(xx2)(xx3)=nmgl2I1dt(5.12)

两边定积分:

x(0)x(t)dx4(xx1)(xx2)(xx3)=nmgl2I1t(5.13)

对上式左端的积分,理论力学上说,可以化成椭圆积分,然而,它并没有说如何做。然后,我去查了特殊函数概论这本“字典”,找到了类似情况的换元方法。我尝试使用,然而第一次发现积分变量的取值范围不对,故自己换了一次元,然后又发现k的取值范围不对,然后又换了一次元……而且每次都是快写到最后的时候才发现,真的折磨。不过最后试出来了如何换元,以下就不一一说明之前的过程了,直接给出换元的方法。(我去,我那个时候才知道有些书上巧妙的方法到底是怎么来的,不一定是什么灵光一闪,很可能是推导出来的人用头发献祭的……)

我们换元,令ξ=xx3x2x3(1145.14)

显然ξ[1,1]

代入(5.13)得(这一步就直接写结果了啊,就是硬算没啥物理含义):

x(0)x3x2x3x(t)x3x2x3dξ(1ξ2)(1(x2x3x1x3)2ξ2)=n(x2x3)mgl2I1t(5.15)

由定义,显然x2x3x1x3[0,1]

那么,(5.15)可以直接积分,然后把(5.14)代入:

F(arcsin(x(t)x3x2x3),x2x3x1x3)F(arcsin(x(0)x3x2x3),x2x3x1x3)=n(x2x3)mgl2I1t(5.16)

那么,便很容易解出x(t):

cos(θ)=x(t)=(x2x3)sn2(n(x2x3)mgl2I1t+F(arcsin(x(0)x3x2x3),x2x3x1x3),x2x3x1x3)+x3(5.17)

在进一步设元化简之前,我们需要讨论一下这个n。

从f(x)的图像可以看出,横坐标为x,纵坐标为f(x)的点在图像上函数值大于零,且在[-1,1]之间的部分运动,所以遇到函数值为零的点,或者横坐标为-1,1的点,n的值必须改变,以改变点运动的方向。由上文的图片可以看出,三实根情况下,一定不会遇到横坐标为-1的点,而若遇到横坐标为1的点,一定同时遇到函数值为零的点。也就是说,n在f(x)=0处改变符号,即n在

cos(θ)=x3cos(θ)=x2处改变符号。由(5.17)易知,当sn2=0时cos(θ)=x3,当sn2=1时cos(θ)=x2那么,同上文21楼的讨论,n是否改变符号对函数知无影响,所以可以当作常量取初始值。

接下来,为了方便,同上文一样定下参量(这里就和从上文取一样的了,反正是不同情况,不会搞混,反而可以一一对应,更加美观)

λ=n(x2x3)mgl2I1(5.18)

μ=F(arcsin(x(0)x3x2x3),x2x3x1x3)=F(arcsin(cos(θ0)x3x2x3),x2x3x1x3)(5.19)

k=x2x3x1x3(5.20)

那么,(5.17)便可以化成:

cos(θ)=(x2x3)sn2(λt+μ,k)+x3(5.21)

这样,结合(5.2)便可以得出拉格朗日-泊松情况在f(x)有三个实根情况下的一般解。

游客没有发表内容的权限。想参与大家的讨论?现在就 登录注册
文号 / 929155

名片发私信
学术分 0
总主题 10 帖总回复 121 楼拥有证书:进士 机友
注册于 2020-08-28 12:52最后登录 2025-03-19 23:47
主体类型:个人
所属领域:无
认证方式:手机号
IP归属地:广东

个人简介

火箭爱好者

文件下载
加载中...
视频暂不能访问,请登录试试
仅供内部学术交流或培训使用,请先保存到本地。本内容不代表科创观点,未经原作者同意,请勿转载。
音频暂不能访问,请登录试试
投诉或举报
加载中...
{{tip}}
请选择违规类型:
{{reason.type}}

空空如也

插入资源
全部
图片
视频
音频
附件
全部
未使用
已使用
正在上传
空空如也~
上传中..{{f.progress}}%
处理中..
上传失败,点击重试
等待中...
{{f.name}}
空空如也~
(视频){{r.oname}}
{{selectedResourcesId.indexOf(r.rid) + 1}}
处理中..
处理失败
插入表情
我的表情
共享表情
Emoji
上传
注意事项
最大尺寸100px,超过会被压缩。为保证效果,建议上传前自行处理。
建议上传自己DIY的表情,严禁上传侵权内容。
点击重试等待上传{{s.progress}}%处理中...已上传,正在处理中
空空如也~
处理中...
处理失败
加载中...
草稿箱
加载中...
此处只插入正文,如果要使用草稿中的其余内容,请点击继续创作。
{{fromNow(d.toc)}}
{{getDraftInfo(d)}}
标题:{{d.t}}
内容:{{d.c}}
继续创作
删除插入插入
插入公式
评论控制
加载中...
文号:{{pid}}
笔记
{{note.content}}
{{n.user.username}}
{{fromNow(n.toc)}} {{n.status === noteStatus.disabled ? "已屏蔽" : ""}} {{n.status === noteStatus.unknown ? "正在审核" : ""}} {{n.status === noteStatus.deleted ? '已删除' : ''}}
  • 编辑
  • 删除
  • {{n.status === 'disabled' ? "解除屏蔽" : "屏蔽" }}
我也是有底线的