文档文库
手机版
投诉建议
热门搜索： 心得体会 演讲稿 思想汇报
首页 > 数学建模人体膝关节受力分析西南财经大学校赛2015

数学建模人体膝关节受力分析西南财经大学校赛2015

发布时间：   来源：文档文库   
小 中 大
字号：
手机查看
2015年西南财经大学数学建模校赛




人体膝关节运动学问题

参赛队员信息
姓名
专业年级
学号
电话邮箱
2015/5/4



人体膝关节运动学问题
摘要：
对于问题一，分析对峰力矩的影响因素及其影响大小。首先处理了数据中的异常和缺失数据，用excel做出了各变量相对于峰力矩的散点图，发现速度对其影响不受其他变量干扰，求出其它变量在不同角速度下与峰力矩的相关系数。结果是年龄，左右脚与峰力矩相关系数极小；身高，峰力矩角度与峰力矩相关性不大；峰力矩角度与峰力矩是测试系统同时生成的观测指标，不能作为自变量。性别x1，体重x4，屈伸膝x6，角速度w与峰力矩y的相关系数分别为0.571075,0.5195775，0.48212，-0.49915。对相关性较强的变量建立多元回归分析模型，用matlab软件解得方程为
y0.655123.0532x11.2939x431.5904x60.1314
对于问题二，分析人胫股关节接触力与屈膝角度、身体各部位倾斜度的关系。我们将身体简化为以胯和膝两处为转折点，躯干、大腿、小腿为三段均匀杆的模型，杆的宽度是膝盖的宽度。假设人体质量均匀分布与长度正相关。（1）人的重心落在脚的地面的接触点上，过重心垂线两侧质量相等（两侧杆长度之和相等）。（2）胫股关节接触力力矩等于重力力矩。（3）查找资料确定人体各关节活动角度范围，完成模型建立。用LINGO软件求解得胫股关节最大接触力是体重的7.1倍。此时人体姿势为大腿与地面水平，与小腿夹角为45º，小腿与地面夹角为45º，腰部与水平面夹角为80.12º。对于问题中的说法验证结果为人体屈膝30º，膝关节承受压力是体重的1.54641倍；屈膝60º，膝关节压力为体重的4.0926倍；屈膝90º，所承受的压力是体重的6.44204倍。在一定误差范围下说法是正确的。
对于问题三，分析人体在上下台阶时胫股关节接触力与腿部动作速度的关系。（1）沿用问题二中均匀杆模型，建立力矩平衡方程。（2）假设始终单腿承重，重心落在承重腿与地面接触点上。取上（下）一级台阶时间为一周期，完成承重腿由弯曲到直立（直立到弯曲）的动作。根据非承重腿刚离开（接触）地面时与地面和垂线构成直角三角形建立方程。（3）首先研究小腿与垂线角度，列出其与时间的关系式,再找到力与该角度关系,用LINGO解出力最大时的角度,确定当时测试者的姿势。LINGO的结果显示上下楼梯胫股关节接触力最大时姿势相同，在小腿与竖直面的夹角为41.38º膝盖受力最大，力为5505.44N，是体重的7.86倍，平均受力为2752.72N。
对于问题四，定性分析了在举重过程中胫股关节接触力与其产生的对上半身支撑力的关系。建模后并进行两次修正。模型一中沿用问题三中均匀杆模型，根据大小腿在水平方向上分量相同建立联系，又根据虚功原理（膝盖水平方向做功等于支撑力竖直方向做功）列出支撑力G'与胫股关节接触力N关系式：，从中看出大腿与Ntan
G'
1+
bsin1c2b2cos21
1

地面垂直时即使N很小G'也趋于无限大。模型二修正了模型一中N为恒力的假设，得
b
(NNsin112出cG'
bsin1
1
c2b2cos21
，可见G'有一极限值，同时G'不仅与大小腿拉力有关还与大腿弯
曲角度成正弦而非正切关系。模型三保持踝关节位置不变，修正了模型二中关于膝盖位移沿水平方向的假设，修正后结果G'
sin1
(N1cos1N2cos2sin2。
sin(12
关键词：回归模型人体均匀杆模型优化问题胫股关节受力分析
1

1.问题重述
1、膝关节力量的测试分析
采用CON-TREX等速测力系统采集实验数据：选择膝屈/伸两个实验项目，进行四种方案测试：静止130°用力、运动60º/s、180º/s、300º/s，分别进行5次。测试者上身进行固定，要求双手握住两侧扶手，测试时必须用尽全力。
测试数据见文件：data1：数据项包括：测试者编号、性别（1男2女）、年龄、身高、体重、左/右腿（1左2右）、屈/伸（1伸2屈）、静止130º峰力矩、60º/s峰力矩、60º/s峰力矩角度、180º/s峰力矩、180º/s峰力矩角度、300º/s峰力矩、300º/s峰力矩角度。
试分析测试数据有那些特征，即：峰力矩的值与那些因素有关，以及关系的强弱。2、膝关节承重分析
体重负荷下，胫股关节接触力随屈膝角度增大而增加。有资料显示，人体屈膝30º，膝关节承受压力和体重相等，屈膝60º，膝关节压力为体重的4倍，屈膝90º，所承受的压力是体重的6倍。
试建立数学模型，分析在体重负荷、静止、双脚支撑状况下，胫股关节接触力与屈膝角度、身体各部位倾斜度的关系，确定最大胫股关节接触力及对应的屈膝角度、小腿等的倾斜度。并说明上段说法是否正确（可在一定误差下）。3、台阶运动对膝关节的影响
爬楼梯属于负重运动，上下台阶时下肢各关节的运动幅度、关节负荷以及肌肉活动等均与在平地上静止、行走有差异，膝关节起主要承重和缓冲作用。
有资料显示，正常人在爬楼梯时膝关节承受的压力会在瞬间增加3倍。即，一位体重为70公斤的人在爬楼梯时其两侧膝关节所承受的压力则高达280公斤。同时，爬楼梯速度越快，膝关节承受的压力就越大。
考察台阶：长90cm、宽28cm、高18cm，测试者：170cm、70kg，速度：96步/分。试建立数学模型，分析上下台阶时，胫股关节接触力与上下楼梯时腿部动作、速度等的关系。分析上下楼梯是否有差异、上下楼梯最大膝关节压力各是多少、平均膝关节压力各是多少。并说明上段说法是否正确。4、运动对膝关节的影响
若时间容许的话，请选取步行（例如快步走）、武术（例如太极拳）、球类（例如篮球）、田径（例如跳远）等一个或多个运动项目，对运动对膝关节的影响进行进一步讨论。
2.问题分析
2.1问题1
分析对峰力矩的影响因素及其影响大小。这可以看成是多元回归模型。我们先用插值法处理了数据中的异常数据，求出每个变量相对于峰力矩的相关性。发现年龄，身高，左右腿，峰力矩角度与峰力矩相关性不强，将其从回归中剔除。同时用Excel散点图发现在不同的角速度下峰力矩的变化趋势几乎一致，说明角速度对峰力矩的影响不受其他因素的干扰，故先分析一种速度下其他因素对峰力矩的影响。最后加入角速度因素并对模型进行优化。2.2问题2
求出膝关节最大受力的情况。这可以看成是优化模型，并且运用了力学受力原理。可以将人体简化为三段轻杆（小腿，大腿，躯干）和2个节点(膝关节，腰的受力模型。
2

通过受力分析建立方程，用LINGO求解。2.3问题3
分析膝盖在上楼下楼时所受的压力。可以看成是优化模型，运用力学受力原理，延用问题2的假设，以一步为周期，建立有关力学模型，用LINGO求解。2.4问题4
分析人体在举重时胫骨关节的受力和人体产生的支撑力的关系。做定性分析，延用问题2的假设，建立有关力学模型，求出表达式。
3模型假设
1假设统计的数据真实有效，与现实无偏差；
2假设实验对象除了给出的变量以外其他情况完全相同；
3）人体在力学研究中简化为大腿，躯干，小腿三部分，股，膝为两处折点；4）人体在力学研究中质量均匀分布，重心在经过脚的与地面垂直的线上；5）人体在上下楼运动中完成一个周期后的姿势不变；6）人体重心在上下阶梯换承力腿时瞬间转移到承力腿上。
4符号说明
x1,x2,x3,x4,x,5x6,x7：性别，年龄，身高，体重，左右腿，伸屈膝，峰力矩角度
：运动角速度
y：峰力矩值
a,b,c：分别表示躯干，大腿，小腿的长度
a1：躯干在重心线左边的长度a2：躯干在重心线右边的长度
：小腿与地面夹角
：大腿与水平面的夹角：躯干与水平面夹角
F：膝盖所受到的压力
F1：肌力F2：重力l1：肌力的力臂l2：重力的力臂
N1：大腿对膝关节的拉力
3

N2：小腿对膝关节的拉力
d：膝关节的受力宽度
m：人体质量
g：重力加速度
：上楼运动中年小腿与竖直面夹角
0：上楼初始状态时小腿与竖直面夹角
：下楼运动中小腿与竖直面夹角
t：下楼末状态时小腿与竖直面夹角
h：台阶的高度k：台阶的宽度
v：人体运动速度
G'：支撑力
m0g：举重时的物重
dx：举重时膝关节水平方向上的位移
dy：举重时主动力G'的虚位移
dr：举重时主动力N1,N2的虚位移
N：经股关节所受的横向力恒力
Nx：经股关节所受的横向力变力，N1,N2在X方向上的合力s0：人体脚掌长度
s：改进后重心线与后跟接触点距离
5模型建立与求解
5.1问题一模型的建立和求解5.11数据处理
根据分析，文件data1中出现三个异常值，分别为M8缺失（第一位测试者在右膝屈膝情况下第二次测得300º/s峰力矩），K29缺失（第二位测试者在左膝伸膝情况下第三次测得180º/s峰力矩），H74数值异常（第四位测试者右膝伸膝情况下的第四次静止
4

130º峰力矩）。对于缺失数据，用插值法进行修正，取相同情况下测试的其他几组数据的平均值作为修正数据，修正后的结果分别为66.47（M8）,112.70(K29；对于异常数据（7777.72）直接修正为（77.72）。
为求结果的精确，将每位测试者4种情况（左腿伸膝，左腿屈膝，右腿伸膝，右腿屈膝）下测试的五组数据取其平均值作为最终数据。5.12模型建立和求解
峰力矩受到多个变量的影响，适用于多元回归模型。
首先用Excle软件绘制不同角速度下的峰力矩折线图（图1）。

图1不同角速度下的峰力矩折线图
分析发现角速度对于峰力矩的影响不受其他变量的干扰，曲线变化趋势基本一致。故在分析其他变量对峰力矩的影响时可以逐个分析每一个角速度下的影响。
然后用Excel2013版中correl(arrange1,arrange2函数求出每个变量对应的相关系数，分析每个变量在不同角速度下与峰力矩的相联性，结果如下表（表1）所示。相关性别x1年龄x2身高x3体重x4左右脚x5伸屈膝x6角度x7系数
静止0.5904-0.06020.42170.5012-0.0400.5309600.5437-0.06780.37870.51340.01010.5687-0.39361800.5525-0.07250.37470.54860.09210.52150.22393000.5941-0.06750.42180.514920.07390.307170.3438
表1变量与峰力矩的相关性表
分析发现年龄，身高，左右脚与峰力矩相关性过弱，不将其作为模型的有效变量。而峰力矩力矩角度是指力矩曲线中，峰力矩所对应的角度，是作为测试系统自带的观测指标，不能作为独立的自变量进行处理，且相关性分析结果偏小，故也不作为有效变量。
剩下的有效变量x1性别，x3身高，x4体重，x5伸屈膝在不同的速度下与峰力矩的相关性大致相同，证明原假设（角速度对于峰力矩的影响不受其他变量的干扰）成立。对其也做相关性分析，R值为-0.499155736，呈现较强的负相关，也是有效变量。
对剩下的变量性别，体重，伸屈膝和角速度建立多元回归模型如下

y01x12x43x64
5

其中分别为性别，体重，伸屈膝和角速度。通过MATLAB求解,程序见
附录1.
算法如下：
（1）利用MATLAB统计工具箱中的命令regress求解回归方程；（2）利用MATLAB统计工具箱中的命令rcoplot作出残差图；
（3）利用MATLAB编程剔除残差异常点重复（1）（2）步骤，直到残差正常，此时得到回归方程；
求解结果见下表（表2）
02413值-19.6304
置信区间[-34.4241
-4.8368]
-25.58021.742928.0193[22.4559[1.4833[25.065928.7046]2.0026]30.9727]表2第一次回归结果表
R2=0.7335F=370.8149P=0.0000
-0.1470[-0.1598-0.1342]
残差图见下图（图2）

图2第一次回归残差图
图中红色数据点为残差异常点，将其剔除后再回归，结果见下表（表3）
01
-22.7857
2
1.2957
3
1.2957
4
-0.1279[-0.1363-0.1196]
-0.6125值
置信区间[-4.2375
-3.0124]
[20.8613[1.2247[28.795024.7100]1.3667]32.4773]表3第二次回归结果表
R2=0.75F=105.70P=0.0000
残差图见下图（图3）
6


图3第二次回归残差图
再次剔除残差异常点，重新回归，结果见下表（表4）
01
-23.0532
2
1.2939
3
31.5904
4
-0.1314
-0.6551值
置信区间[-4.1596
-2.8493]
[21.1701[1.2251[29.790324.9364]1.3626]33.3905]表4第三次回归结果表
R2=0.79F=111.50P=0.0000
[-0.1395-0.1232]
残差图见下图（图4）

图4第三次回归残差图
故回归方程为
y0.655123.0532x11.2939x431.5904x60.1314
7

5.2问题二的模型建立和求解5.21模型建立(1受力结构
理想的成年人的身体一般常用8个头长作为身高比例。从脚部向上算起，小腿约占两个头长，大腿两个头长，躯干三个头长。脚与地面的接触为一点，头竖直于躯干，忽略双臂的影响，我们将身体模型简化为一胯和膝两处为转折点的三段轻杆，杆的宽度是膝盖的受力宽度如下图（图5）。
假设人体质量均匀分布，则在分析受力结构时，人体可以再次简化为轻杆重心线的连线要使人站立不倒，人体重心作地面的垂线必须落在脚与地面接触的点上。受力结构见下图（图6）。

图5人体简化模型图图6人体重心模型图
要使人站立不倒，人体重心作地面的垂线必须落在脚与地面接触的点上。人各部分长度与质量成正比，为保证在做肢体运动时人体保证稳定状态，不发生偏倒，重心需在经过脚的与地面垂直的线上，即重心线左右两部分质量应该相等，也就是两部分长度相同，且各部分长度应该满足图中集合关系，保证人的各部分连接在一起。
即满足：
cccosseca1bccosseca2
a2cosbccosseccos
a1a2a
a为躯干的长度，b为大腿的长度，c为小腿的长度，a1为躯干在重心线左边的长度，
a2为躯干在重心线右边的长度，为小腿与水平方向夹角，为大腿与水平方向的夹角（正时表示大腿在水平面以上，为负则表示大腿在水平面以下），为躯干
8


与水平方向夹角。
资料显示人体的踝关节运动角度满足正常背伸至70度，跖屈至140度，约有70度的活动范围，即活动角度为45度~90度，因为人体不能完成骨干的完全重合，故膝关节活动角度为15度~180度，腰部活动角度为0度~180度。即满足：
11421
12
0
（2）受力分析
经股关节受力是来自大腿和小腿肌肉拉力的合力，是以膝关节侧面图中心位置为支点的一个杠杆模型。膝关节受力与力矩有关，其宽度在受力分析时不能忽略不计。一力作用在距支点一定距离的地方，就会引起角运动，这个角运动的趋势就是力矩，用以表示旋转力的大小，力矩是作用力和力臂的乘积。力臂是力的作用线到旋转轴之间的垂直距离。膝盖只有在所受力矩和为0时才能保持静止状态，不发生旋转。则满足肌力对膝关节的力矩=重力对膝关节的力矩。受力图如下（图7），膝关节支点为宽度的中点。


图7力矩受力分析图
即满足：
Fl11F2l2l1ccos1l2d
2F1mg

9

F2表示肌力，F1表示重力，l1表示重力的力臂，l2表示肌力的力臂，a表示躯干长度,c表示小腿的长度，为小腿与地面的夹角，d表示膝关节的受力宽度，m表示人体质量，g是重力加速度。
膝关节受到的压力为大腿对其的拉力和小腿对其的拉力的合力。如下图（图8）所示。

图8膝关节受压力示意图
即满足：
2
F2N12N22N1N2cos(
F为膝关节所受的压力，N1为大腿对膝关节的拉力，N2为小腿对膝关节的拉力，为
小腿与水平面的夹角，为大腿与水平面的夹角。
人体膝关节的运动是由关节周围的结缔组织带动，由于结缔组织在单一肌群里，各点及结缔组织的作用力皆相等，即大腿肌肉对膝关节的拉力N1与小腿肌肉对膝关节的拉力N2相等。又因为其受力处处相等，即人体膝盖很类比于滑轮模型，各方向力是相等的，即肌力F2与拉力相等。如下图（图9）所示。


图9人体结缔组织受力图
10

即满足：
N1=N2=F2
N1为大腿对膝关节的拉力，N2为小腿对膝关节的拉力，F2为肌力。
综上模型为：
maxF
cccosseca1bccosseca2
a2cosbccosseccosaaa12Fl11F2l2Fmg1
l1ccos
lds..t22222
FN1N22N1N2cos(NNF
22
1
11
42
112
0

5.22模型求解
利用lingo求解。程序见附录2.结果见下表（表5）

45
0



0
F

0
F/mg

7.1
0

80.12

5685.07N
表5膝盖最大受力结果表
一个标准成年人（体重80kg，身高180cm，膝盖受力宽度5cm），膝盖受力最大为5685.07N是体重的7.1倍，此时人体姿势为膝盖与地面水平，与小腿夹角为45º，小腿与地面夹角为45º，腰部与水平面夹角为80.12º。
对于问题中的说法：人体屈膝30º，膝关节承受压力和体重相等，屈膝60º，膝关节压力为体重的4倍，屈膝90º，所承受的压力是体重的6倍。其屈膝角度相当于模型
中的。检验结果见下表（表6）。
F（受力）F/mg（倍数）
屈膝30º1237.1281.5464
屈膝60º3274.1414.0926
表6结果检验表
11
屈膝90º5153.6396.4420

因此，在一定误差下，题目说法是正确的。5.3问题三模型的建立和求解5,31模型建立（1）上楼梯
在做上楼运动时，人体一直处于单脚支撑状态，重心在换脚瞬间转移，始终落在支
撑脚与地面的接触点上。
对于双腿在爬楼梯时的初始状态，其分析图见下图（图10）

图10上楼初始状态分析图
在初始状态时，人的承重腿大小腿形成一定角度与台阶面接触，非承重腿，大小腿拉伸成一条直线与地面接触，双腿交叉处与承重腿和地面的接触点在同一条垂直于店面的线上，且非承重腿与地面接触点到该线的距离为一个步长，等于台阶的宽度k。非承重腿，垂直线和步长形成一个直角三角形，满足勾股定理。
即满足：
[(bccos0h]2k2(bc2
b,c分别为大腿，小腿的长度，bc为全部下肢的长度，h为台阶高度，k为台阶宽
度，0为初始状态下小腿与垂直面的夹角。
运动过程中，每一次爬楼梯都是做重复运动。即可以研究一个周期的运动状态来表示整个过程的运动状态。承重腿的小腿在刚接触台阶时保持同一倾斜角度0，与大腿的夹角固定。运动过程中即承重腿交换过程中，小腿与大腿匀速拉伸成一条直线，在拉伸过程中小腿与竖直面的夹角会减少到0，这时非承重腿变为承重腿，以相同形状与台阶接触。其运动中受力分析见下图（11）

12

图11上楼力矩分析图
沿用问题二中对于膝关节轻杆—滑轮模型的构建。膝盖处必须保持重力与肌力的力矩相等，使得膝盖不发生左右旋转，支点为膝盖受力宽度的中点。理想的成年人的身体小腿约占2个头长，大腿占2个头长，故大腿和小腿的长度可以近似看作一样。
即满足：
bcF1l1F2l2
l1bsinl2dF1mg
b,c分别为大腿和小腿的长度，F2为膝盖所受肌力，F1表示重力，l1表示重力的力臂，l2
表示肌力的力臂，表示运动过程中小腿与竖直面的夹角，d为膝盖受力力臂。与第二问一样，结缔组织在单一肌群里，各点及结缔组织的作用力皆相等，即大腿肌肉对膝关节的拉力N1与小腿肌肉对膝关节的拉力N2相等。又因为其受力处处相等，即人体膝盖类似于滑轮模型，各方向力是相等的，即肌力F1与拉力相等。膝关节受到压力即为大腿和小腿对其施加力的合力。受力分析见下图（图12）

图12上楼膝盖压力分析图
即满足：
2
F2N12N22N1N2cos(2


2


N1N2F2
F为膝盖所受压力，N1为大腿对膝关节的拉力，N2为小腿对膝关节的拉力，F2为肌力，
13


为运动中小腿与竖直面的夹角，为运动中大腿与水平面的夹角。
上楼梯的周期为小腿与竖直面夹角从变为0的时间，即完成一个台阶上楼运动的时间。同时人体用速度v（步/分钟）上楼梯，一步则可以完成上一个阶梯的动作，所以每步所用时间也可表示周期。
即满足：
T10
1T1=
v
T1为上楼的周期，0为初始状态小腿与竖直面的夹角，为大小腿角度变化的角速度，v为人上楼速度
综上所述模型为：
maxF
[(bccos0h]2k2(bc2bcFlF2l211
l1bsinld2

s..tF1mg
2
F2N12N22N1N2cos(2
2
NNF
221
01T=1
v
（2）下楼梯
下楼时末状态，即完成一次下楼动作时的状态如下图（图13）
图13下楼末状态分析图

14

与上楼时初状态类似。承重腿大小腿形成一定角度与台阶面接触，非承重腿大小腿拉伸成一条直线与台阶面接触，双腿交叉处与承重腿和地面的接触点在同一条垂直于地面的垂直线上，非承重腿，步长和垂直线构成直角三角形，满足勾股定理。
即满足：
[(bcostccosth]2k2(bc2
b,c分别为大腿，小腿的长度，bc为全部下肢的长度，h为台阶高度，k为台阶宽度，
t为最终状态下小腿与垂直面的夹角。
运动过程中，承重腿在初始时保持大小腿拉伸呈直线。运动过程中即承重腿交换过程中，小腿与大腿匀速弯曲成一定角度，在拉伸过程中小腿与竖直面的夹角会从0增加到t，这时非承重腿变为承重腿，以相同形状与台阶接触。
下楼过程的受力分析与上楼类似，其分析图见图14和图15。


图14下楼力矩分析图图15下楼膝盖合力分析图即满足：
bc
Fl11F2l2l1bsin
l2dF1mg
2
F2N12N22N1N2cos(2
N1N2F2
b,c分别为大腿和小腿的长度，F1为重力，F2表示肌力，l1表示重力的力臂，l2表示肌
15


力的力臂，表示运动过程中小腿与竖直面的夹角，d为膝盖受力宽度。F为膝盖所受压力，N1为大腿对膝关节的拉力，N2为小腿对膝关节的拉力，为运动中小腿与竖直面的夹角。
下楼依旧取一步为一个周期，周期大小是小腿与竖直面夹角从0变为的时间，即完成一个台阶上楼运动的时间。
即满足：
T2t
1T2=
v
T2为上楼的周期，t为末状态小腿与竖直面的夹角，为大小腿角度变化的角速度，v为人上楼速度。
综上所述模型为：
maxF
[(bccosth]2k2(bc2bcF1l1F2l2
l1bsinld2s..tF1mg
2
F2N12N22N1N2cos(2
N1N2F2T2t
1T2=
v
5.32模型求解
用lingo对上下楼梯的模型进行求解，程序见附录3.结果如下。
上楼梯时当为41.380时，即小腿与竖直面夹角为41.380时膝盖受力最大，力为5505.44N，平均受力为2752.72N是体重的7.86倍。下楼梯时结果遇上楼一致。与人运动速度无关。
5.4问题四模型的建立和求解5.41模型建立和求解（1）模型1
构建人体做举重运动时的膝关节受力模型。沿用第二问中均匀杆模型。在举重时，人体大腿与小腿弯曲成一定角度，杠铃通过躯干将重量施加于下肢。假设大小腿长度不一样，大小腿与竖直面的夹角也不一样，即大小腿的倾斜程度不同。根据几何图形，大小腿应满足在水平方向的投影长度相同。其受力结构如下图（图16）。
16




图16下肢复杠杆示意图
即满足：
bsin1csin2，bcos1ccos2
b,c分别为人大腿和小腿的长度，1,2分别为大腿和小腿与水平面的夹角。
举重时人体对杠铃的支撑力应是重力和杠铃重力的合力的反作用力。即满足：
G'mgm0g
G'为胫股关节向上的支撑力，m为人体质量，m0为杠铃质量，g为重力加速度。
在举重过程中大小腿从弯曲状态逐渐伸直，即弯曲角度逐渐变小。其运动过程见下图（图17）。


图17模型1的复杠杆作用图
以过膝关节平行于地面的直线为X轴建立坐标系，膝关节原始位置在X轴左边，运动过程中逐渐向右移动，在X轴上形成一段位移，大小腿与水平面的夹角也逐渐变大，最终人体成站直状态，即大小腿拉伸成直线且与地面垂直，在Y轴上形成两段位移
17

dy1,dy2。
即满足：
dy1bbsin1
dy2ccsin2dxbcos1ccos2
其中dy1为大腿在竖直方向上的位移，dy2为小腿在竖直方向上的位移，dx为膝关节在水平方向的位移，b,c分别为人大腿和小腿的长度，1,2分别为大腿，小腿与水平方向的夹角。
根据虚功原理，横向力所做的功应与支撑力的反作用力做的功相同。即满足：
NdxG'(dy1dy2
综上模型一如下：
bcos1ccos2dybbsin
11
dy2ccsin2

dxbcosccos12
NdxG'(dy1dy2G'mgm0g
求解得：
G'
1+
结果分析：当1
Ntan1
bsin1cbcos1
2
2
2
m0
(1+
Ntan1bsin1cbcos1
2
2
2
mg

2
的作用力N，也能产生相当大的支撑力G'，复杠杆的机械力是相当大的。因此，人体
时，即膝关节接近伸直时，tan1的值相当大，因此即使很小
下肢能产生极大的支撑力使人举起质量相当大的m0。
（2）模型一修正（模型二）
模型缺陷：
模型一中，由于大腿和小腿肌肉对胫股关节所合成的横向力N被当作恒力来处理,事实上力N不仅随时间变化而变化,而且随着膝关节角度的变化,其合力的大小也将随之变化。
模型改进：
18

人体在举重时，其重心线始终通过支撑点，即大小腿对胫股关节的拉力一起作用使得人体保持平衡。故将横向力N用N1,N2在水平面的分量表示。
即满足：
NNxN1cos1N2cos2
Nx为N1,N2在X方向上的合力,N1为大腿对膝盖的拉力，N2为小腿对膝盖的拉力。
则新模型为：
bcos1ccos2
dy1bbsin1dy2ccsin2
dxbcos1ccos2
NdxG'(dydy
12
x
NxN1cos1N2cos2
G'mgm0g
求解得:
bb
(N1N2sin1(N1N2sin1
ccmG'm0
bsin1bsin1
(1g1
222222
cbcos1cbcos1

结果分析:当190时,即膝关节接近伸直时,则sin11,cos10,支撑力
G'有一极限值,且G'max(N1N2。同时表明支撑力G'不仅与N1,N2有关，并且与
膝关节角的关系为正弦关系,而并非正切关系,因此其实际的机械效益并非很大。固有
(m0max
(N1N2
m，人体上举m0不仅与N1,N2有关，并且与膝关节角的关系为正弦g
关系,不能无限大。在模型一中，膝关节在接近伸直时,很小的肌力就能产生很大的支撑力。改进的模型首次将肌力看成是变化的,结果表明支撑力有最大值,较符合人体运动的实际情况。下肢伸膝时复杠杆产生的对上体的支撑力有极限值，不可能超过构成复杠杆的两肌群肌张力总和。(2模型二修正（模型三）
模型缺陷：
上诉两个模型都是假设膝关节运动方向水平。但是实际生活中的约束应该是踝关节为一定点，膝关节运功并非水平。
模型改进：
由于模型2在建立的过程中破坏了原有的踝关节约束，且未加说明，在力学原理的应用中容易使人产生误解和错误的概念。虽然，虚位移在完整、双面和定常的理想约束
19

中不是唯一的，但应在不破坏约束的条件下，质点为约束所允许的无限小位移。所以膝关节在大小腿拉力的作用下的虚位移方向与小腿垂直，则在几何关系上满足勾股定理。受力分析见下图（图18）。

图18主动力示意图和虚位移示意图
即满足：
dxdrcos(2
2dydrsin(2
2


其中dx为膝关节在水平方向的位移，dy为膝关节在竖直方向的位移，dr为膝关节在垂直于小腿方向上的虚位移，2为小腿与水平方向的夹角。
由于将人简化为刚体，刚体上任意两点的虚位移在两点连线上的投影必须相等。即满足：
drcos(12dycos(1
22

1为大腿与水平方向的夹角。
根据虚位移原理大小腿横向力所做的功应与支撑力的反作用力做的功相同。
即满足：
NxdxG'dy0
则修正模型为：
20


dxdrcos(22
dydrsin(
2
2
drcos(dycos(
1212
22NdxG'dy0x
NxN1cos1N2cos2
G'mgm0g1
求解得：
G'
sin1
(N1cos1N2cos2sin2
sin(12
sin1
(N1cos1N2cos2sin2
sin(12m0m
g
结果分析：若将bcos1ccos2带入模型化简后可得模型二的表达式。模型三为模型二的另一种表达。
7模型的评价和改进
7.1模型的评价（1）回归模型评价
回归分析法在分析多因素模型时，更加简单和方便。运用回归模型，只要采用的模型和数据相同，通过标准的统计方法可以计算出唯一结果。
回归分析可以准确地计量各个因素之间的相关程度的高低，提高预测方程式的效果。
多元回归对推广变量之间的函数关系很有帮助，多元回归的许多应用都涉及到主要变量之间的非线性关系。（2）膝盖受力模型评价
将人体简化为轻杆模型能够在展示人体下肢受力情况的同时将实际问题大大简化。在建模过程中考虑了人体受力的生物学特点，使得模型不单单是数学模型，还能够与现实更好的拟合。（3）举重模型评价
分析膝关节胫股关节在伸屈膝过程中的作用时,简化模型具有十分重要的地位。由于关节结构和肌肉力学的复杂性,简化模型仅是为了解释和说明膝关节胫股关节屈肌群在伸膝过程中的作用,属于解释性模型,在实际中只能作为定性的说明。以上三个模型都属于解释性模型，模型二、三则表示上举力有极值。在应用虚位移原理求解人体(多刚体系统平衡问题时，没有必要将系统拆开,忽略约束力的限制可将问题大为简化。7.2模型的改进
对于膝关节受力模型，原模型中将人的脚部简化为一点，人的重心需在脚与地面的接触点上，即重心线是固定的，见下图（图19）。而现实生活中，人的脚掌有一定的长
21

度，重心线可以在以脚长为宽度的范围内垂直于地面移动，使腰部，膝关节和踝关节有更大的活动范围。其移动情况见下图（图20）。


图19原重心假设图图20修正后重心假设图
此时重心线可以在以脚后跟为端点，以人体脚掌长为长度的水平射线范围内移动。即满足：
0ss0
此时为保持身体平衡，人体姿势任需满足在重心线左右两边身体各部分长度相等。即满足：
a1b2c1a2b1c2
在几何关系上满足：
c2ssecc1cc2b2c1cossec
b1bb2a2b1cosseca1aa2
将新的约束条件加入模型，结合重心线移动规律将能更好的拟合人体膝关节的受力情况。

22

8参考文献
[1]姜启源.数学模型[M].北京：高等教育出版社,2011。
[2]卓金武.MATLAB在数学建模中的应用[M].北京：北京航空航天大学的出版社,
2010.
[3]潘慧炬.人体膝关节复杠杆模型的修正[J].体育科学,2000.0320卷第2期.[4]RW斯达西,刘普和等译.生物物理和医学物理学原理[M].北京:科学出版
社,1962.
[5]体院运动生物力学教材编写组编.运动生物力学[M].北京:人民体育出版
社,1990.
[6]李良标,吕秋平等.运动生物力学[M].北京:北京体育学院出版社,1991.[7]高师运动生物力学教材编写组编.运动生物力学[M].北京:高等教育出版
社,1999.
[8]周衍柏.理论力学教程[M].北京:人民教育出版社,1979.[9]赵景员.力学[M].北京:人民教育出版社,1979.[10]J·菲利普,等,胡勋力等译.运动学和应用解剖学[M].北京:人民体育出版社,
1985.
[11]刘建智.膝关节复杠杆简化模型的研究与教学[J].陕西:西安体育学院学报,
2002.
9附录

附录1.
Matlab编程求解回归方程：
>>x1=[];x4=[];x6=[];x8=[];%双击x1x4x6x8后赋值>>y=[];%用同样方法给y赋值
>>x=[ones(544,1x1x4x6x8];%给x赋值
>>[b,bint,r,rint,stats]=regress(y,x%用regress命令求解>>rcoplot(r,rint%画出残差图
%由于残差图中有较多值，所以需要将其剔除>>id=find((rint(:,1.*rint(:,2>0;>>x1(id=11111;>>x4(id=11111;>>x6(id=11111;>>x8(id=11111;>>y(id=11111;
>>x1=x1(find(x1<11111;>>x4=x4(find(x4<11111;>>x6=x6(find(x6<11111;>>x8=x8(find(x8<11111;
>>y=y(find(y<11111;%将异常值剔除>>x=[ones(521,1x1x4x6x8];
>>[b,bint,r,rint,stats]=regress(y,x;%再次用regress做回归分析>>rcoplot(r,rint%画出剔除异常值后的残差图
23


附录2
求解静止状态膝盖受力最大的情况sets:
a/1..6/:l;b/1..3/:m,x,y;c/1..2/:F;endsetsdata:
m=402020;enddata
pi=3.1415926;x(1>=pi/4;x(1<=pi/2;
x(1+x(2>=pi/12;x(1+x(2<=pi;x(2+x(3>=0;
x(2+x(3<=pi;!根据所查到的资料，限制踝关节膝关节的最大弯曲程度;l(3=0.45;l(4=0.45;!为大腿和小腿的长度赋值;
l(4+l(4*@cos(x(1/@cos(x(2+l(1=l(3-l(4*@cos(x(1/@cos(x(2+l(2;!重心所在垂线左右两端质量相等;
l(2*@cos(x(3=(l(3-l(4*@cos(x(1/@cos(x(2*@cos(x(2;!列出角x(3和x(2的关系;
l(1+l(2=0.9;!上半身的长度定为0.9;
F(1*l(5=F(2*l(6;!重力端的力矩与膝盖端的力矩相等;F(1=(m(1+m(2+m(3*10;!重力等于质量乘以g重力加速度;l(5=l(4*@cos(x(1;!算出重力端的力臂长度;l(6=0.05;!查询资料定膝盖端的力臂为0.05米;
@for(b(i:y(i=x(i*180/pi;!将角度由弧度制转换为角度制，便于观察;max=@sqrt(F(1^2+F(2^2-2*F(1*F(2*@cos(pi-x(1-x(2;!根据合力公式算出膝盖所受合力，求出最大值;
附录3
上楼梯模型sets:
a/1..2/:F,l,m;endsets
pi=3.1415926;
b=0.45;c=0.45;h=0.18;k=0.28;!分别给bchk赋初值;T=5/8;!根据题目给出的每分钟96步计算出周期;F(2*l(2=F(1*l(1;!力矩平衡;
l(1=b*@sin(m(1;!从受力分析图可求出重力端力臂与角度的关系;l(2=0.05;F(1=700;
24

((b+c*@cos(m(2+h^2+k^2=(b+c^2;!从受力分析图可算出上楼梯末状态的角度m(2;
w=m(2/T;!由得出的m(2和周期T算出角速度w;m(1=m(2-w*x;!得出角度m(1与时间x的关系;m(2<=pi;!限定最大角度m(2在180度以内;m(1<=m(2;
max=@sqrt(2*F(2^2-2*F(2^2*@cos(2*m(1;!根据合力公式算出膝盖所受合力，求得其最大值;得其最大值;下楼梯模型sets:
a/1..2/:F,l,m;endsets
pi=3.1415926;
b=0.45;c=0.45;h=0.18;k=0.28;!分别给bchk赋初值;T=5/8;!根据题目给出的每分钟96步计算出周期;F(2*l(2=F(1*l(1;!力矩平衡;
l(1=b*@sin(m(1;!从受力分析图可求出重力端力臂与角度的关系;l(2=0.05;F(1=700;
((b+c*@cos(m(2+h^2+k^2=(b+c^2;!从受力分析图可算出上楼梯初始状态的角度m(2;
w=m(2/T;!由得出的m(2和周期T算出角速度w;m(1=w*x;!得出角度m(1与时间x的关系;m(2<=pi;!限定最大角度m(2在180度以内;m(1<=m(2;
max=@sqrt(2*F(2^2-2*F(2^2*@cos(2*m(1;!根据合力公式算出膝盖所受合力，求得其最大值;
25

本文来源：https://www.2haoxitong.net/k/doc/6034b653332b3169a45177232f60ddccda38e6be.html

《数学建模人体膝关节受力分析西南财经大学校赛2015.doc》
将本文的Word文档下载到电脑，方便收藏和打印
推荐度：
点击下载文档

文档为doc格式
相
关
案
例

相关推荐

推荐内容

短语结构类型 长度、宽度、高度、深度等表示方法 老年节目策划文案 年会表彰会策划方案(2篇) 三年级数学《年月日》评课稿 关于青春奋斗的朗诵词 软件流程图 古代成语故事 常见古代成语故事 《菱形》教案[1] 《邓稼先》精神品质三个典型事例
网站内容来自网络，如有侵权请联系我们，立即删除！ Copyright © 范文写作网京ICP备16605803号