建模一个傅科摆
这个例子展示了如何建模一个傅科摆。傅科摆是法国物理学家福柯里昂。这是为了证明地球绕着它的轴旋转。傅科摆振动平面的旋转在一天由于地球的自转。振荡完成整个圆的飞机在一个时间间隔,这取决于地理纬度。
福柯的最著名的钟摆是安装在巴黎先贤祠。这是一个28公斤的金属球体附加到一个67米长的电线。这个例子模拟一个67米长的钟摆在巴黎的地理纬度。
万博1manbetx仿真软件®模型
最简单的方法解决傅科摆的问题在仿真软件®是建立一个模型,解决了耦合系统的微分方程。万博1manbetx这个模型如图1所示。下面的方程描述给出了傅科摆。有关的物理模型和这些方程的推导过程,明白了分析和物理。
打开模型
类型sldemo_foucault
在MATLAB®命令窗口打开这个模型。该模型仿真数据日志到变量sldemo_foucault_output
。记录信号有一个蓝色的指示器。有关更多信息,请参见马克信号记录。
图1:傅科摆模型
初始条件
这个模型加载的常量和初始条件sldemo_foucault_data.m
文件。这个文件的内容如表1所示。你可以直接修改仿真参数在MATLAB工作区。摆的初始振幅必须小而摆的长度,因为只有小振动微分方程是有效的。
表1:初始条件
g = 9.83;%重力加速度(m /秒^ 2)L = 67;%摆长度(m) initial_x = L / 100;%初始x坐标(m) initial_y = 0;%初始y坐标(m) initial_xdot = 0;%初始x速度(米/秒)initial_ydot = 0;%初始y速度(米/秒)ω= 2 *π/ 86400;%地球的旋转角速度对其轴(rad /秒)λ= 49/180 *π;%的纬度(rad)
运行仿真
按下“播放”按钮在工具栏窗口运行仿真模型。仿真步骤硬解算器将使用一个变量,ode23t。将模拟一个傅科摆3600秒(你可以改变仿真时间)。模型使用一个默认的相对宽容RelTol = 1 e-6
。
图2:傅科摆仿真结果(模拟时间为3600秒)
结果
上面的仿真结果如图2所示。仿真计算钟摆x和y坐标,x和y速度分量的钟摆。
钟摆摆动平面完成360度扫描超过24小时。地理纬度的扫描周期是一个函数λ
(见推导分析和物理)。
图3:动画块显示钟摆摆动平面旋转多少一个小时
运行仿真后,双击动画块将结果。
注:“动画效果”部分的示例需要信号处理工具箱™。双击动画块会导致一个错误如果没有安装。所有其他部分的示例将正常工作,而信号处理工具箱。
的sldemo_foucault_animate.m
文件块的位置摆鲍勃在不同的时间点。你可以清楚地看到钟摆摆动平面旋转。
注意:如果您正在运行的模拟大相对宽容,结果将是数值不稳定在很长一段时间。确保您使用的是硬变步求解器进行求解。阅读更多关于僵硬的数值不稳定和解决性能问题的“探索解决变步使用的模式”例子。
关闭模式
关闭模式。生成的数据。
分析和物理
本节分析了傅科摆并描述了其背后的物理学。钟摆可以建模为一个质点悬线的长度l
。钟摆位于地理纬度λ
。方便使用的坐标系如图4所示:惯性坐标系I(相对于地球的中心),和惯性坐标系N(相对于一个观察者在地球表面)。惯性框架加速旋转的结果。
图4:的惯性和惯性框架问题
点O惯性坐标系的原点n是地球表面上的点在悬挂点的钟摆。选择非惯性坐标系,z轴点远离地球的中心和垂直于地球表面。轴点南部和轴点西方。
在介绍中提到的,傅科摆振动平面的旋转。振荡平面完成一个完整的旋转刚学步的小孩
由以下公式给出,Tday
是一天的时间(即时间地球围绕它的轴一次)。
sin因素需要进一步讨论。常常错误地认为,钟摆的摆动平面固定在惯性坐标系相对于地球的中心。这仅仅是真正的在北极和南极。为了消除这种混淆,思考点(见图4),摆的暂停。在惯性坐标系,S点转一圈。钟摆鲍勃是悬浮在一线的长度。为简单起见忽略空气摩擦。在惯性坐标系,只有两种力量作用于鲍勃-线张力T
和重力成品
。
向量r
给的位置摆鲍勃,B(参见图4)。牛顿第二定律指出,力量作用于身体的总和等于质量乘以加速度的身体。
在这个证明,点表示时间衍生品,箭头表示向量,大写表示单一向量(i, j, k在x, y,和z轴)。一个点以上向量箭头表示的时间导数向量。箭头上方的点表示向量的时间导数。见下面的加速度和径向加速度之间的区别。
总加速度:
径向加速度:
重力加速度指向地球中心(- z)。
分解加速条件:
单位向量的时间衍生品出现,因为惯性参考系N是在太空中旋转。这意味着单一向量i, j, k在太空中旋转。下面给出时间衍生品。ω是地球的旋转角速度绕着它的轴。标量ω是角速度的值。角速度矢量ω是向量。它的方向是由右手定则来确定。
重写的时间导数向量r相对于ω。
同样,表达的第二次导数向量r。
为了简化方程,假设ω为地球很小。这让我们忽略上面的方程的第三个任期。事实上,第二项(已远小于第一项)是四个数量级大于第三项。这减少了方程如下形式:
牛顿第二定律可以编写和分解为x, y,和z组件如下:
角振荡幅度很小。因此,我们可以忽略垂直速度和垂直加速度(z-dot和z-double-dot)。弦张力组件可以使用小角度近似表达,这也大大简化了问题,使其二维(见下文)。
特性微分方程
最后的物理问题可以耦合方程组所描述的下面。x和y坐标指定的位置摆鲍勃所看到的地球上的观察者。
分析解决方案(近似)
下面是一个傅科摆解析解的问题。不幸的是,这并不确切。如果你试图解析解代入微分方程,uncanceledω的平方的顺序将保持。然而,由于ω是很小的,可以忽略uncanceled条款实际用途。
实际微分方程系统是不对称的
在推导,条款涉及ω的平方都被忽略了。这导致了xy对称性在微分方程。如果考虑ω的平方项,微分方程系统变得不对称(见下文)。
你可以很容易地修改当前傅科摆模型考虑到不对称微分方程。只需编辑相应的增益模块包含g / L
和添加必要的表达式。这种变化将引入一个很小的整体校正的数值结果。