地球引力场的悬链线方程
By 苏剑林 | 2011-05-15 | 68081位读者 |之前曾在《自然极值》系列文章中提到过均匀重力场下的悬链线形状问题,并且在那文章中向读者提出:在一个质点(地球)引力场中的悬链线形状会是怎么样的。说实话,提出这个问题的时候,我还不懂怎么解答这个问题,不过现在会了,回头一看,已经几个月了,时间过得真快...
与之前的思路一样,我们依旧采用的是“平衡态公理”,即总势能最小。从天体力学中我们知道,任意两个质点间的势能为−Gm1m2r。对于本题的悬链线问题,我们可以把地球放到坐标原点位置,而悬链的两个固定点分别为(x1,y1)和(x2,y2),链的总长度为l。即
∫x2x1√dx2+dy2=l
设悬链的线密度为1(无量纲处理,长度即质量)。那么悬链上的点(x,y)到(x+dx,y+dy)一小段的引力势能可以记为
dE=−GM√dx2+dy2√x2+y2
总势能显然是E=−GM∫x2x1√dx2+dy2√x2+y2
选取适当的坐标对被积函数进行化简是很必要的,我们采取极坐标系:
√dx2+dy2√x2+y2=√r2dθ2+dr2r=√r2˙θ2+1rdr
其中˙θ=dθdr
于是E=−GM∫r2r1√r2˙θ2+1rdr=−GM∫r2r1F(˙θ,r)dr
代入拉格朗日方程就得到(将r看成自变量,θ作为函数):
0=∂F∂θ=ddr(∂F∂˙θ)=ddr(r˙θ√r2˙θ2+1)
于是积分得到r˙θ√r2˙θ2+1=Const
并且化简得到:r2˙θ2=C21(C1>0)
要求C1>0只是为了描述上的方便,即±˙θ=C1r⇒C2±θ=C1lnr
我们约定,在θ=C1lnr−C2<0时,上面去负号,反之取正号,则可以写成:
C2+|θ|=C1lnr
可以进一步化简成r=AeB|θ|的形式。其形状如下图:
转载到请包括本文地址:https://spaces.ac.cn/archives/1361
更详细的转载事宜请参考:《科学空间FAQ》
如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。
如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!
如果您需要引用本文,请参考:
苏剑林. (May. 15, 2011). 《地球引力场的悬链线方程 》[Blog post]. Retrieved from https://spaces.ac.cn/archives/1361
@online{kexuefm-1361,
title={地球引力场的悬链线方程},
author={苏剑林},
year={2011},
month={May},
url={\url{https://spaces.ac.cn/archives/1361}},
}
May 15th, 2011
这个图像好奇怪啊
April 25th, 2014
"代入拉格朗日方程就得到"下面的公式的第二个等号应该是减号,后面计算就都有问题了。。。感觉上,悬链的形状不应该有那个尖点,应该是光滑曲线(至少有一阶连续导数)。
我当初也困惑了好久,为此我仔细检验过,但是没有错呀,拉格朗日方程是
∂F∂θ=ddr(∂F∂˙θ),但是我们有∂F∂θ=0,所以上面是没有计算错误的。
April 26th, 2014
我也没看出推导中的错误。图像还有一个明显的问题:悬链应该凸向引力中心才对。我英语不太好,你可以看看这篇文献:http://www.mathematik.uni-muenchen.de/~hinz/catenaria.pdf
好,谢谢。我仔细阅读后再与您讨论。
关于评论问题,是因为带有链接的评论(如含有http字样的)将会进入审核状态,所以推迟了显示时间,万望谅解。现在已经取消了这个限制。^_^
那篇文章的过程我看不太懂,不过结论和你的不太一样。。。。
我已经大概知道问题所在,但是我没有想到解决问题的思路。这个问题比较奇怪。过一段时间我会专门写写它
请问题主找到问题的根源了没
无
August 1st, 2014
这个曲线很像双曲线,也有两渐近线的 ,不是方程错了就是图画错了,用软件数值求解微分方程画图就行
这个问题我还在困惑着,极坐标下的测地线问题~~所以暂时先不做修改,等我理清了再说...
March 12th, 2020
可能这是极大值
August 23rd, 2023
是因为势能零点默认为无穷远处,而你又没有考虑绳长限制,求出来的积分是错的。应该加一个拉格朗日乘子到绳长积分上,再与势能相加,带人拉格朗日方程