找回密码
 欢迎注册
查看: 39360|回复: 57

[原创] 算算这个摆的运动方程

[复制链接]
发表于 2014-10-15 10:45:33 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?欢迎注册

×
两个物体当成质点质量都是1  杆长度1,系统摩擦不计,重力加速度为10.就这么放手。求上面滑块随时间的运动方程,求下面摆的运动轨迹……
QQ图片20141015104145.jpg
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2014-10-15 11:00:08 | 显示全部楼层
杆没有质量   开始滑块位置在(0,0) 那个初始角度就取60度了
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-10-15 11:22:00 | 显示全部楼层
理论力学中很普通的题目,使用牛顿力学或者分析力学的方法都可以列出运动方程,但后者更方便。因此直接考虑设定虚位移,建立动力学普遍方程 或者 设定独立的广义坐标,列出拉格朗日方程即可。
最后结果是一致的,这个微分方程组没有解析解,需要数值求解。

点评

轨迹是椭圆!系统不会水平运动!系统会振动!(ps:不是“震动”!)  发表于 2014-10-19 10:38
大致轨迹是什么 整体会不会水平运动下去 还是呆在原地震动  发表于 2014-10-15 17:19
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-10-15 11:25:51 | 显示全部楼层
感觉上面太轻会跳起来。要改条件成卡在类似过山车式轨道上吗?

点评

各个物体的属性都标明了 都是国际制单位 只需要数值计算轨迹就行了  发表于 2014-10-15 17:17
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-10-15 18:03:00 来自手机 | 显示全部楼层
由于棍子的约束,球相对方块的运动方向应该和棍子垂直。由于水平方向没有外力,水平冲量为零,两者水平速度和为零,然后再根据能量守恒,可以求出给定球的高度的速度情况
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-10-15 18:32:34 来自手机 | 显示全部楼层
不考虑方块和地面间弹力是不会跳起的
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-10-15 18:50:57 | 显示全部楼层
下面小球的摆动轨迹很容易求得。
设球的坐标B(x,y),则根据水平方向上的动量守恒,可以知道AB质点系的中心横坐标为定值,不妨设为0,则A的坐标为(-x,0)
从而由|AB|=1得:
4x^2+y^2=1
是个椭圆。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2014-10-15 21:13:39 来自手机 | 显示全部楼层
还有滑块随时间的运动函数
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-10-15 21:34:52 | 显示全部楼层
系统约束为完整约束,故可用拉格朗日方程。
沿用一楼图,显然系统只有两个自由度,选取 `x` 和 `\varphi` 为广义坐标。于是系统动能为$$T=\frac{1}{2}m \dot x^2+\frac{1}{2}m [(\dot x+\dot \varphi l \cos \varphi)^2+(\dot \varphi l \sin \varphi)^2]$$以方块所在水平面为零势面,则有系统势能$$V=-mg l\cos \varphi$$拉格朗日函数$$L=T-V$$利用第二类拉格朗日方程$$\begin{cases}\D\frac{\dif}{\dif t}\frac{\partial L}{\partial \dot x}-\frac{\partial L}{\partial  x}=0 &&&&& (1) \\\\
\D\frac{\dif}{\dif t}\frac{\partial L}{\partial \dot \varphi}-\frac{\partial L}{\partial  \varphi}=0
&&&&& (2)\end{cases}$$
于是得到运动方程$$\begin{cases}2m\ddot x+ml\ddot \varphi \cos \varphi-ml\dot \varphi^2\sin\varphi=0 \\
\D \ddot \varphi+\frac{1}{l}\ddot x\cos \varphi+\frac{g}{l}=0
\end{cases}$$
初始条件 `x(0)=x^{\ast},\dot x(0)=0\;\mathrm{m/s}`,`\varphi(0)=\varphi^{\ast},\dot \varphi(0)=0 \;\mathrm{rad/s}`

将 `g=10 \;\mathrm{m/s^2}`,`m=1\;\mathrm{kg}`,`l=1\;\mathrm m` 代入上述方程即可。

补充内容 (2014-10-16 19:34):
上面运动方程中的最后一项 `\frac{g}{l}`应该是 `\frac{g}{l}\sin \varphi`

点评

@倪举鹏,分析力学是有严格数学基础的,你这有数学依据吗?摩擦力可以当成主动力,求出对应广义位移上的广义力放在方程右边即可。你真的可以看看书了,跟你推荐了好多次,你就是不听,这样学也只能学个半吊子。  发表于 2015-4-11 23:51
如果滑块有摩擦力 这个力乘以走过的路程的耗散功,负值可不可以添加到动能项里  发表于 2015-4-11 21:46
方程会得到不一样的x值与φ值,是我的方程有问题么  发表于 2014-10-17 20:40
用牛顿加速度方程φ''(t)=-10sin(φ(t)) φ(0)=pi/3 φ'(0)=0 x"(t)=1/2*sin(φ(t))(10cos(φ(t))+(φ'(t))^2) x(0)=-3^0.5/4 x'(0)=0也可以得到质心水平不动,可是用相同的时间数值t代入,你的方程跟我这   发表于 2014-10-17 20:39
@倪举鹏,是时间。这个例子中的主动力(重力)是保守力,故方程右端是零。L是时间的函数,只不过不显含时间。  发表于 2014-10-17 18:30
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2014-10-16 09:40:17 | 显示全部楼层
分析力学的知识,
忽略了杆的质量

点评

回复倪举鹏:轨迹就是真的椭圆!此系统的质心就是杆的中点横坐标一直就没有变化!因为系统水平动量守恒!@倪举鹏  发表于 2014-10-19 10:33
有个问题,下面摆的轨迹只是形状像椭圆的数值方程,还是轨迹根本就是真的椭圆?此系统的质心就是杆的中点横坐标一直就没有变化吗?  发表于 2014-10-16 15:03
杆质量不忽略一样可以做的,只要知道或者求出杆的转动惯量即可。  发表于 2014-10-16 09:57
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-5-8 12:45 , Processed in 0.094164 second(s), 20 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表