椭球摆

文化 2024-08-28 17:12 吉林

上个月发了一个椭圆摆的微振动问题（见：物竞-曲率圆-微振动）：

它这个问题的振动方向是沿左右的，不是垂直纸面的，当然题干描述的确实不严谨，因为严格来讲这个模型确实可能是一个二维振动模型。

然后它这个左右微振动的等效摆长和周期结果是

那天有个老师问了这样一个问题，沿垂直纸面振动的周期是多少？怎么计算，我们也讨论了一下，后来由于计算太麻烦以及太晚了就没继续进行了，然后也忘了，今天想起来还是算一算比较好，他这个问题来源是24年6月浙江高考物理第九题：

题中提到垂直纸面微振动等效于小球重垂线于AB交点的单摆，为什么是又或者为什么这样等效？合理性在哪里？我第一反应就是做的近似，但我今天算完发现事情没有那么简单。

我一开始研究了一下有没有什么好算的技巧，比如找焦点长短轴，发现并不容易，直接上方程似乎是最简单粗暴的，他这个题干是倾角30°，我们把这个倾角设成是任意的。

很容易发现，在绳子的约束下，小球的运动轨迹是在一个椭球面上的，并且是两个轴相等的轴对称的椭球面。

当θ=0时，AB水平，沿左右方向振动轨迹是椭圆，等效摆长（曲率半径）和周期就是

沿垂直纸面方向的振动时，运动轨迹恰好是一个圆，摆长和周期是

当θ=π/2时，AB竖直，沿左右方向振动轨迹依然是一个椭圆，等效摆长和周期是

此时的椭球面是关于z轴对称的，那么垂直纸面的振动和左右振动是一样的：

这时小球重垂线与AB重合，如果取极限的话并不是A点（在后面写了），而小球到A的距离是 $OA=l-d$ ，这里 $OA\ne R_{12}$ ，有点差距。

最后再回到求解任意倾角时垂直纸面方向的微振动了，首先它是一个倾斜的椭球面，它的方程可以由坐标系旋转得到，在x’y’z’坐标系（S＇）中，该椭球面的方程是

$\frac{x\prime^2}{b^2}+\frac{y\prime^2}{a^2}+\frac{z\prime^2}{b^2}=1$

其中 $a=l\,\, , b=\sqrt{l^2-d^2}$ ，坐标系沿x轴旋转为xyz坐标系（S系），转动矩阵是

$K=\left( \begin{matrix} 1& & \\ & \cos \theta& \sin \theta\\ & -\sin \theta& \cos \theta\\ \end{matrix} \right)$

得到坐标变换

这样代入原方程就可以得到我们的椭球面方程：

看起来就很复杂了，我们先直接令x=0，得到

$\frac{\left( y\cos \theta +z\sin \theta \right) ^2}{a^2}+\frac{\left( y\sin \theta -z\cos \theta \right) ^2}{b^2}=1$ 这个就是沿水平方向振动的运动轨迹，z的最小值就是平衡位置,在ｚ最大最小处时，ｙ只有一个值，只需要判断对于y的判别式为零即可：

$\Delta =\sqrt{\left( \frac{z\sin 2\theta}{a^2}-\frac{z\sin 2\theta}{b^2} \right) ^2-4\left( \frac{\cos ^2\theta}{a^2}+\frac{\sin ^2\theta}{b^2} \right) \left( \frac{z^2\sin ^2\theta}{a}+\frac{z^2\cos ^2\theta}{b}-1 \right)}=0$

可解得

别看他复杂，结果是真简单，计算该点ｙ坐标

再计算该点的曲率半径，当然在这个坐标系下曲率半径不好求，那就换回去在S'系中计算

$\frac{y'^2}{a^2}+\frac{z'^2}{b^2}=1$ 在

处的曲率半径

最后结果是

$R=\frac{l^2\left( l^2-d^2 \right)}{\left( l^2-d^2\cos ^2\theta \right) ^{\frac{3}{2}}}$ 结果没有问题，这时我们得到了最低点的坐标是 $(0,y_0,z_0)$ 。

这时我们令椭球面方程中的 $y=y_0$ ，就得到了垂直于纸面方向的振动轨迹（上图紫色轨迹）,然后化简（过程不愿意看，省略）：

$\frac{x^2}{b^2}+\frac{1}{a^2b^2}\left( l^2-d^2\sin ^2\theta \right) \left[ z-\frac{d^2y_0\sin \theta \cos \theta}{\left( l^2-d^2\sin ^2\theta \right)} \right] ^2=1-\frac{y_{0}^{2}}{\left( l^2-d^2\sin ^2\theta \right)}$

轨迹还是一个椭圆，我们需要计算在 $x,z=z_0$ 处的轨迹，注意看这个椭圆方程并没有旋转只是平移，处理一下

计算曲率半径

最后就是

$R=\frac{l^2-d^2}{\sqrt{\left( l^2-d^2\cos ^2\theta \right)}}$ 这个个就是最后结果了，当θ=0时 $R=\sqrt{l^2-d^2}$ ，当θ=π/2时， $R=\frac{l^2-d^2}{l}$ ，与最开始分析的结果一致。

OK，我们再回到那个高考题上θ=π/6，OA和AB垂直，那么垂直纸面方向振动的等效摆长是 $R=\frac{l^2-d^2}{\sqrt{l^2-\frac{3}{4}d^2}}$ ，并且根据几何关系

可以得到

$\begin{cases} OA+OB=2l\\ OB^2-OA^2=4d^2\\ \end{cases}\Rightarrow \left\{ \begin{array}{c} OA=l-\frac{d^2}{l}\\ OB=l+\frac{d^2}{l}\\ \end{array} \right.$

以及 $OA=OB\cos \left( 2\alpha \right)$ ，得到

又因为

$2l\cos \left( \frac{\pi}{2}-2\alpha +\frac{\pi}{6} \right) =2d\cos \frac{\pi}{6}$

那么

舍掉了两个不符合条件的解，那么曲率半径是 $R=\frac{2d^2}{\sqrt{3d^2-\frac{3}{4}d^2}}=\frac{4}{3}d$ ，然后再算一算小球悬挂点到小球重垂线与AB交点的距离：

$h=\frac{2}{\sqrt{3}}OA=\frac{4}{3}d$ 哦吼，破案了，竟然真的是一样的。

那这是特例吗？我们不再看题，依然研究一般情况的小球悬挂点到小球重垂线与AB交点的距离，由几何关系

得到

按相似三角形的性质还有

$\frac{h}{l\cos \alpha -d\sin \theta}=\frac{l_2}{l}$ 得到

哎，你别说，还真不是特例，这椭球摆的垂直纸面微振动的等效曲率圆中心就在椭球长轴上。

上述过程也可以用参数方程来计算，这个模型写拉氏量解EL方程也不容易。

完结撒花！欢迎讨论！

http://mp.weixin.qq.com/s?__biz=MzI1Mjk2MTE5OA==&mid=2247510055&idx=1&sn=1a63b441fd37dce916bdafa542792096

Cosmos and Us

物理是世界上至高无上的艺术。在这里，我会谈论这个世界中最基本的艺术：一种叫做物理学的艺术。这里有宇宙诞生、基本粒子与宇宙结构、天体塌缩到黑洞、黑洞蒸发的故事。我会从大爆炸说起，直到现在，到……未知的未来……以及谈论科学与我们的生活。

最新文章

浅浅讨论一个关于ln|1-z|的环路积分

没想到什么标题

网传一道极限题

新闻？搞点无聊的东西

3.7：由吉布斯熵公式推导系综

3.8：等温等压系综

数理表情包VI：学不完了

标准宇宙学模型的一些时期

知识永不完结！

数理表情包V：硝苯氮

部分思考问题的提示与解答241113

广义相对论24：水星近日点进动与验证

狭义相对论前传V：部分拖曳假说与裴索实验

English to English I

广义相对论23：史瓦西时空中的自由落体与圆周运动

中科大2019年617普通物理A考研题

认识物质变化之窗——广义热力学方程组

公共数学资料（高数线代概率论）

数理化自学参考资料（高中和大学）

别吵了！

结束了，但是还没结束

The Elements Song

AI元年

“苹果落地”的热力学诠释

它来了，它带着OCR来了！

第四类永动机为何不能实现

自然现象的热力学原理

代数几何中的曲面专题（第一篇）：曲面上的几何（上）

代数几何中的曲面专题（第二篇）：曲面上的几何（下）

是否应有热力学第四定律与第五定律

1.8Maxwell总结的电磁现象最初的20个方程

椭球摆

一种余热资源再利用的热力循环及实现循环的热气机

经典力学6.0角动量：角动量定理

经典力学8.0守恒律与对称性

经典力学8.1力学守恒律与几何

论工程热力学方程组在提升热气机效率上的作用

开放性物理思考问题v240827

希腊字母读音

关于工程热力学统一理论的思考

南京大学2023年802普通物理

李政道逝世（1926.11.24-2024.8.4）

中科大2024年828量子力学考研题

中国科学院大学2024年811量子力学考研题

经典力学5.0功能原理：动能定理与能量守恒

文章更新周期、公式显示bug等说明

西安交通大学2023年722量子力学考研题

中科大2024年617普通物理A考研题

中科大2020年617普通物理A考研题

公众号文章目录20240715版

分类

时事

民生

政务

教育

文化

科技

财富

体娱

健康

情感

旅行

百科

职场

楼市

企业

乐活

学术

汽车

时尚

创业

美食

幽默

美体

文摘

原创标签

时事社会财经军事教育体育科技汽车科学房产搞笑综艺明星音乐动漫游戏时尚健康旅游美食生活摄影宠物职场育儿情感小说曲艺文化历史三农文学娱乐电影视频图片新闻宗教电视剧纪录片广告创意壁纸头像心灵鸡汤星座命理教育培训艺术文化金融财经健康医疗美妆时尚餐饮美食母婴育儿社会新闻工业农业时事政治星座占卜幽默笑话独立短篇连载作品文化历史科技互联网

发布位置

广东北京山东江苏河南浙江山西福建河北上海四川陕西湖南安徽湖北内蒙古江西云南广西甘肃辽宁黑龙江贵州新疆重庆吉林天津海南青海宁夏西藏香港澳门台湾美国加拿大澳大利亚日本新加坡英国西班牙新西兰韩国泰国法国德国意大利缅甸菲律宾马来西亚越南荷兰柬埔寨俄罗斯巴西智利卢森堡芬兰瑞典比利时瑞士土耳其斐济挪威朝鲜尼日利亚阿根廷匈牙利爱尔兰印度老挝葡萄牙乌克兰印度尼西亚哈萨克斯坦塔吉克斯坦希腊南非蒙古奥地利肯尼亚加纳丹麦津巴布韦埃及坦桑尼亚捷克阿联酋安哥拉