为什么小姐姐能摇一晚上不倒?
The following article is from 数学模型 Author 大模头
推荐阅读：建议收藏：高中各科思维导图汇总
西安大唐不夜城“不倒翁”女孩街头表演的视频曾一夜走红网络。在大唐不夜城步行街，“不倒翁”小姐姐身姿轻盈眼神妩媚令人梦回大唐，一颦一笑将中国唐朝美人的妩媚娇羞演绎得淋漓尽致。
https://v.qq.com/txp/iframe/player.html?vid=i0959y2darg&auto=0
这勾人的小眼神，我觉得画面太美还能看 100 遍。尽管很美，但其实小姐姐表演得很辛苦。别只看到小姐姐在演出时体态轻盈、游刃有余，这背后可是下足了苦功夫。
大唐不夜城的真人不倒翁演员都是经过层层选拔的，女性演员身高要在 163 厘米左右，体重不超过 50 公斤。扮演不倒翁的演员冯佳晨[2]，今年 23岁，身高 163 厘米，体重 45 公斤。虽有近10年的舞蹈功底，但这个表演刚开始的时候对她来说真的不容易。因表演需要将下半身完全固定在如图 3 所示的铁架子上，仅靠腰部发力运作不倒翁的底座，冯佳晨每天表演完，胯和膝盖都会被磨青[3]。
不倒翁，是一种最常见的玩具。通常形状像一个蛋型、上轻而下重，扳倒后还能自动竖立起来。历史最早记载唐代的捕醉仙[4]就是一种不倒翁。不倒翁不倒的原理并不难以理解：上轻下重的物体比较稳定，重心越低越稳定。当不倒翁在竖立状态处于平衡时，重心和接触点的距离最小，即重心最低。偏离平衡位置后，重心总是升高的。因此，这种状态的平衡是稳定平衡。所以不倒翁无论如何摇摆，总是不倒的。
但不倒翁究竟是如何运动的呢？小姐姐又是如何操控不倒翁，让其自由的摇摆呢？ 本文将建立不倒翁小姐姐的数学模型，研究不倒翁的运动方程，并通过计算给出小姐姐摇晃的最佳策略。
模型
先不考虑小姐姐的主动发力和地面滚动阻力，建立一个简单的不倒翁模型。在此基础上，再建立一个考虑小姐姐的主动发力和地面滚动阻力的模型。
简单的不倒翁模型
我们把不倒翁简化为仅在二维平面（图 5 中的 平面）内运动的半球。半球和地面均为刚性，两者无滑移。如图 5 所示的半球，其质量为 ，质心为 ，半球球心为 。
容易求得半球质心位置 和对 轴的转动惯量[5]
若该半球某时刻的位置如图 6 所示， 为半球重力，地面对半球的支撑力和摩擦力分别为 和 。
为了确定运动方程，需要知道半球是如何运动的。我们用 与垂直方向的角度变量 来描述半球的状态。考虑点 和 的坐标
因此，质心 的速度为
其中 和 是 和 对时间的一阶导数。将 点速度对时间求导数，得到质心的加速度
现在让我们通过考虑作用在物体上的所有力来研究运动方程。运动方程不仅要考虑刚体的平动，还必须考虑到关于质心转动：
有三个方程，却有四个未知数（，，，）。要求解这个问题，还需要一个描述接触条件的表达式。对于不倒翁，假设半球与地面接触点没有滑移，则有
从而可得： 上式微分方程也可由欧拉-拉格朗日方程[6]导出，感兴趣的读者可自行尝试。上式微分方程可由 MATLAB 函数 ode45 函数给出数值解，并结合式 可完全确定不同时刻半球的位置和角度。计算结果如图 7 所示。
模拟结果表明，如果地面没有滚动阻力（机械能守恒），不倒翁将不停的摇摆。角度 为 0 时速度，角速度 最大，重心 最低。
不倒翁小姐姐模型
在前文中，本文仅给出了一个不倒翁玩具在理想地面上的运动状态。而实际上不倒翁小姐姐的底座在摇晃的过程中是受到滚动阻力的，小姐姐的腰部在不停地扭动。为了了解小姐姐到底是如何将不倒翁摇晃起来的，接下来将建立一个考虑滚动阻力和小姐姐摇晃的模型。
根据网络上的视频和图片，假设整个不倒翁底座重 kg，是由质量为 kg，半径为 m，厚为 m 的铁制半球形圆壳和质量为 kg，半径为 的球缺（）配重以及支架（忽略支架的质量）组成。小姐姐身高 ，体重 。如图 8 左图所示。
实际摇晃过程中，小姐姐只能通过摇晃上半身带动底座摇晃。因此可将小姐姐简化为两个由铰链连接的圆柱（图 8 中图）。假设小姐姐的密度[7]为 ，则容易求得圆柱半径为 m。假设小姐姐腰以下长 m，腰以上长 m。小姐姐的腰部以下固定在 T 字型支架上，因此我们将小姐姐的下半身与底座看成一体（图 8 右图），质量为 kg，转动惯量为 。小姐姐上半身质量为 kg，转动惯量为 。
图 8 右图是对两部分的受力分析，将小姐姐上半身对下半身的作用力简化为力 和力偶矩 。在底座实际滚动过程中，由于地面和底座的变形，地面对底座的支撑力 的作用点并不在 点，偏离距离 可由滚动阻力系数给出。与简单的不倒翁模型类似，我们可以定义出 、、、 和 的坐标、速度和加速度，然后分别给出系统两部分的动力学方程，对于底座和小姐姐的下半身组成的整体有 对于小姐姐的上半身有 根据以上六式以及 可以求解出 ，由于表达式比较复杂，这里不再列出。小姐姐可以通过 来控制 ，即小姐姐可以通过摇晃上半身来使整个不倒翁底座摇晃起来。假设小姐姐摇晃的 随时间的变化满足： 其中 为摇晃幅度， 为摇晃周期。
求解
经过我们的分析，终于得到了小姐姐的摇晃方程。可有人说看不懂方程，也不想看方程，就要看小姐姐摇。没办法，为了满足读者的要求，我们只好写段 MATLAB 程序模拟一下。不倒翁底座的摇晃方程可以由 ode45 函数直接求解，但为了加快程序运行速度并实时显示，使用了 Euler 法求解，时间步长取 s。模拟的主要步骤如下：
根据时间确定上半身相对于底座的角度 、角速度 和角加速度 根据角速度更新角度： 根据角加速度更新角速度：，其中角加速度 由运动方程给出。 确定底座球心 的水平位置：。 根据底座球心 的水平位置，以及底座的转动角度，确定半圆位置。 更新时间步 。如果 s，结束；否则跳到第 1 步。
为了模拟的视觉效果，我们还需要动态展示出结果。为此专门写了一个函数用来画底座半圆和性感的小姐姐。绘制结果如图 9 所示。很多人看了我画的小姐姐，表示震惊：实在是太像了，将小姐姐妙曼的身姿展现得淋漓尽致，这大长腿，这伸出的纤纤玉手。
根据运动方程，算出某时刻半圆心的位移以及角度后，需要对半圆和小姐姐进行平移和旋转，平移容易。旋转可根据以下公式进行： 以上公式将点 绕 旋转 角。将以上公式写成 MATLAB 函数 rotxy。
结果
为了研究小姐姐究竟是如何通过摆动上半身将底座摇晃起来的。我们研究了不同情况下，底座摇晃的规律。
若小姐姐不摇晃，即摇晃幅度，并且没有滚动阻力，即，则模型退化为简单的不倒翁模型，对于，结果如图 10 所示。
结果表明，即使小姐姐不摇晃，小姐姐的腰也受到周期性的作用力。若加入滚动阻力的考虑并假设 ，则结果如图 11 所示。
结果表明，在滚动阻力的作用下，不倒翁底座会逐渐停止摇晃。也就是说，在实际情况中，想要保持不倒翁底座的摇晃，小姐姐必须持续不断的摇晃身体。
接下来我们讨论小姐姐摇晃身体的周期对不倒翁底座运动的影响。假设小姐姐摇晃身体的幅度 。研究表明，小姐姐身体的摇晃与不倒翁底座的摇晃之间会产生共振效应，当小姐姐摇晃身体的周期 s 时，不倒翁底座摇晃的幅度最大，结果如图 12 所示。
而更小的周期（如 s，结果见图 13）或更大的周期（如 s，结果见图 14）都不能让不倒翁底座大幅度地摇晃。
结论
通过运动和受力分析可以求解出不倒翁小姐姐的运动方程，并可以确定出小姐姐与底座发生共振的最佳摇晃周期。 小姐姐能大幅度摇晃的最佳周期为 2.2 s 左右。若小姐姐照此周期不停的摇，一小时大约要晃 下。 在本文的最佳周期模拟中，小姐姐晃动幅度约为 ，这相当于要晃出 1 m 左右的距离，一个周期内腰移动的距离约为 4 m。小姐姐腰间平均水平发力在 130 N 左右。因此一个周期小姐姐做功约 520 J，摇 1 小时做功约 850 kJ。因为仅考虑了水平方向的力，实际上这低估了小姐姐实际做功。 根据上面两条结论，可知“不倒翁”小姐姐摇一晚上真的很累，为“不倒翁”小姐姐点赞！ 听说不少群众模仿“不倒翁”小姐姐并闪了腰，请大家小心模仿。此外，骨科专家也告诫大家谨慎模仿。 需要此本书的朋友，可关注公众号后，在对话框回复：命题研究，四个字获取下载链接。不会操作的朋友，可文末添加微信好友获取，仅供学习，禁止商用。如侵权，请联系删除
家里种植的猕猴桃，最近开花了，给大家分享下美丽的猕猴桃花
文章文章选自：《 数学模型》如若转载，文末需注明：文章来源于公众号：妙解之慧（ID:WanZhuanShuXue1）妙解之慧由陕西西安孙冰钰老师创建，专注分享初，高中数学优质资源，旨在服务于全国师生，让更多朋友受益。如侵权，请联系删除
2020年二轮备考线上直播：导数双变量 不等式证明（含PPT）
2020届湖北高三数学二轮复习的认识与思考（含PDF)
导数的应用：含参数的不等式恒成立问题（含PPT）
基于核心素养，突出数学理解：2020届概率与统计复习备考定位（含PPT）
近5年全国1卷试题分析，近3年文理8大专题分析及2020二轮复习备考建议（含256页PPT)
2011—2019年新课标全国1卷理科数学高考分析及2020年高考预测（含Word）
李尚志老师谈教与学 | 第一课：黑发要知勤学早，线上课堂效果好
2020高考命题组已经确定，披露2020高考数学命题步骤和方向的变化
基于核心素养发展的教学创新:从有效教学到优质课堂(含PPT）
学习方法：高中数学靠“悟”不靠“练”，在于“走心”但不能“心走”！
没方向、没兴趣、没自信？我的任教经历, 告诉你什么叫“空心病”
往期精选文章回顾
一：高中数学
二：初中数学
中考作文 | 7-9年级写作文的100个怎么办，初中作文所有问题都在里面！
三：小学数学
四：学习方法
新学期先抓习惯，再谈成绩！1-9年级学生好习惯培养指南，请老师转给家长
五：数学之美
六：家乡特产
家里父亲种植的猕猴桃，猕猴桃被誉为“维C之冠，水果之王”富含多种维生素，目前正在热销中，欢迎各位朋友们订购品尝，第一次购买，如果担心桃味，吃完了给钱都行。早上摘，中午装箱，下午快递运走，上市时间每年8月份---10月中旬。
目前家里的猕猴桃树正在发芽，过上一两星期就会开花，开花前需要疏花骨朵，去掉两侧小的，留下中间大的，将来接猕猴桃，开花之后需要人工授粉
文章来源于网络，转发只为让更多朋友学习，尊重原创，如有侵权，请联系删除
点上方蓝字即可阅读相关文章
长按并识别下方二维码，即可关注此公众号，主要分享一些初，高中数学好题，学习方法，由于时间有限，初高中文章交替更新，关注后可在对话框回复数字：1 有小惊喜等着您呢
长按下方二维码即可识别，添加作者“微信”添加时请备注：地区身份，欢迎各位朋友，相识即是缘，感恩遇见