我之前分析了上海交通大学工程力学系刘延柱教授走路力学论文错误,老虎又找了篇刘教授写的科普(也可能是我上次没看见)链接在:https://166.111.121.20:9080/mathjournal/LXYS702/lxys702030.caj.pdf (主要内容见文末附件)
这又是基本错误,而且在兜圈子。
错在哪? 不会运用F=ma,没有物理直觉,数学也不行(注一)。
你必须考虑需要摩擦力最大的时候,才能得出摩擦力对步长的限制。
显然,我直觉认为这个时候需要的摩擦力最大。
那么为什么一脚离地时摩擦力更小?
F = ma !
单脚着地时,地面与鞋底的摩擦力必然等于人体质量乘以其水平方向的加速度。但我们根据经验就知道,你尽量匀速行进时,水平加速度是很小的。所以,你单脚着地时摩擦力必然较小。具体计算如下
假设腿长为L,腿与垂直线夹角为[ix]\theta[/ix],那么[ix]x = L \sin\theta, dx/dt = L\cos\theta \frac{d\theta}{dt}[/ix]
水平加速度为 [ix]d^2x/dt^2 = L\cos\theta\frac{d^2\theta}{dt^2} - L \sin\theta (\frac{d\theta}{dt})^2 [/ix]。
上面的结果其实不要这样计算中学生也能看出来:向心加速度的水平分量加上切向加速度的水平分量,不过有些人喜欢看符号,我也正好测试一下我们的这个写符号的功能。
第一项这个[ix]\frac{d^2\theta}{dt^2}[/ix]是多少?以脚为转轴,力矩正比于 [ix]\sin\theta[/ix],因此[ix]\frac{d^2\theta}{dt^2}[/ix]也正比与[ix]\sin\theta[/ix] . 如果用刘教授的这个模型示意图,
[ix]J\frac{d^2\theta}{dt^2} = mg L \sin\theta ===> \frac{d^2\theta}{dt^2} = mg L \sin\theta/J \approx \frac{g\sin\theta}{L}[/ix]
有了这个,第二项也好办,[tx]\frac{d\theta}{dt} \frac{d}{dt} ({d\theta/dt} ) = g/L \sin\theta \frac{d\theta}{dt} ==> \frac{d}{dt} (\frac{d\theta}{dt})^2 = 2g/L \sin\theta \frac{d\theta}{dt} [/tx]
[ix](Ld\theta/dt)^2 = v^2 + 2gL (1-\cos\theta) ==> L (d\theta/dt)^2 = v^2/L + 2g(1-\cos\theta)[/ix] (其实这就是机械能守恒)
那么水平加速度为
[tx]d^2x/dt^2 = -L \sin\theta (\frac{d\theta}{dt})^2 + g\cos\theta\sin\theta = g\sin\theta(3\cos\theta-2) - \sin\theta\frac{v^2}{L}[/tx]
因此,摩擦力 [ix]F = ma = mg\sin\theta(3\cos\theta-2) - \sin\theta\frac{mv^2}{L} [/ix]
因此,摩擦力 [ix]F= m\sin\theta[g(3\cos\theta-2) - v^2/L][/ix] (注二)
(地面支撑力[ix]N = m \cos\theta [g(3\cos\theta- 2) - v^2/L][/ix])
其中 v为步行速度。而且我们还发现一个有趣的事情,走得快,水平加速度还小些,因为后面那项是负的,这也许说明为什么滑路要快走。
人是怎么走路的?通过后脚脚掌的转动,等于后退伸长,使重心向前腿推进,然后后脚才离地,因此后脚离地时,着地的那条腿与地面的角度是很小的,因此上面的[ix]\sin\theta[/ix]项很小。
因此,我们转了这么半天明白了一件事,人“匀速”行走时,水平加速度较小。f=ma,因此单脚着地时,摩擦力也小。
所以,你如果想计算摩擦系数对步长的限制,不能看单脚着地的时候,而要看双脚着地的时候。为什么这个时候摩擦力大?我们之前进行了计算,摩擦力正比于 [ix] \tan\theta/2[/ix],但这个角度大 。
注一:刘教授的文中套用了小角度单摆的方程,实际是一个[ix]\sin\theta = \theta[/ix]的近似。这就无法研究角度稍大的情况,而且是毫无必要的。单板的小角度近似是在需要计算角度-时间关系时的简化(否则要用到椭圆积分),但在我们的问题里根本不需要得出时间关系,完全没有必要只考虑小角度。
注二: 刘教授的文中假定 [ix] F= K \theta[/ix],其中有个未知常数K,也不知是怎么拍脑袋编出来的
附:刘延柱教授走路力学