CN114330040A - 一种鞋底振动传递特性分析方法 - Google Patents
一种鞋底振动传递特性分析方法 Download PDFInfo
- Publication number
- CN114330040A CN114330040A CN202111145559.XA CN202111145559A CN114330040A CN 114330040 A CN114330040 A CN 114330040A CN 202111145559 A CN202111145559 A CN 202111145559A CN 114330040 A CN114330040 A CN 114330040A
- Authority
- CN
- China
- Prior art keywords
- sole
- excitation
- equivalent
- vibration
- equivalent mechanical
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 28
- 230000005540 biological transmission Effects 0.000 title claims abstract description 10
- 238000009826 distribution Methods 0.000 claims abstract description 46
- 239000007787 solid Substances 0.000 claims abstract description 23
- 230000004044 response Effects 0.000 claims abstract description 17
- 238000012546 transfer Methods 0.000 claims abstract description 9
- 230000005284 excitation Effects 0.000 claims description 59
- 210000001872 metatarsal bone Anatomy 0.000 claims description 19
- 238000000034 method Methods 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000010586 diagram Methods 0.000 claims description 9
- 210000004744 fore-foot Anatomy 0.000 claims description 8
- 239000000463 material Substances 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 210000002683 foot Anatomy 0.000 description 7
- 238000013016 damping Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000000638 stimulation Effects 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000002087 whitening effect Effects 0.000 description 1
Images
Landscapes
- Footwear And Its Accessory, Manufacturing Method And Apparatuses (AREA)
Abstract
本发明公开了一种鞋底振动传递特性分析方法,包括:步骤S1,建立鞋底三维实体模型;步骤S2,建立所述步骤S1中鞋底模型的有限元模型并进行模态分析,获得其固有振动特性;步骤S3,基于所述步骤S2中的有限元模态分析,对鞋底进行稳态动力学分析,获得其在三种工况下对振动的响应情况;步骤S4,利用Python编程语言计算并绘制鞋底等效机械导纳分布云图,及获得鞋底不同区域等效机械导纳的分布情况;步骤S5,基于鞋底不同区域等效机械导纳的分布情况,利用Python编程语言计算并获得鞋底不同区域的等效振动传递率分布情况。本发明可以利用有限元分析的方法获得鞋底对于振动的响应情况及振动在鞋底不同区域的传递情况。
Description
技术领域
本发明涉及振动分析领域,具体涉及一种鞋底振动传递特性分析方法。
背景技术
大量研究表明,人体足部长期暴露于振动环境中,能够使他们患振动诱发白足的风险增加。而足部在接触振动时,鞋底是与来自地面的冲击振动最直接接触的部分,同时也会起到减振、隔振的作用。对于足部和鞋底振动的研究主要采用运动学测量的方法,有着诸多的缺点,如试验周期长、成本高等。
发明内容
本发明的主要目的在于克服现有技术中的上述缺陷,提出一种鞋底振动传递特性分析方法,借助于数值建模和有限元分析的方法,不仅可以研究振动在鞋底中的响应及传递情况,还能为制鞋业与脚用个人防护设备的开发提供评估标准与指导。
为了解决上述的技术问题,本发明采用如下技术方案:一种鞋底振动传递特性分析方法,包括如下步骤:
步骤S1,建立鞋底三维实体模型;
步骤S2,建立鞋底三维实体模型的有限元模型并进行模态分析,获得鞋底三维实体模型的固有振动特性;
步骤S3,基于有限元模型的模态分析,对鞋底进行稳态动力学分析,获得鞋底在三种工况下对振动的响应情况;所述三种工况包括前脚掌激励、全脚掌激励及脚后跟激励;
步骤S4,利用Python编程语言计算并绘制鞋底等效机械导纳分布云图,及获得鞋底不同区域等效机械导纳的分布情况;
步骤S5,基于鞋底不同区域等效机械导纳的分布情况,利用Python编程语言计算并获得鞋底不同区域的等效振动传递率分布情况。
在一较佳实施例中:所述步骤S2具体包括:
步骤S21:将鞋底三维实体模型导入Abaqus有限元分析软件中;
步骤S22:对鞋底三维实体模型按照解剖学原理划分成不同的区域,具体包括4个大的区域,再将这4个大的区域细分为9个小的区域;4个大的区域包括:足跟区、足弓区、跖骨区、趾骨区;其中足跟区包括足跟内侧区HM、足跟外侧区HL,足弓区包括足弓内侧区MF1、足弓外侧区MF2,跖骨区包括第一跖骨区M1、第二、三跖骨区M23、第四、五跖骨区M45,趾骨区包括第一趾骨区T1、第二-五趾骨区T25共9个区域;
步骤S23:对所述步骤S22中的鞋底三维实体模型赋予材料属性、施加边界条件、划分网格单元得到鞋底三维实体模型的有限元模型;
步骤S24:对鞋底三维实体模型的有限元模型在一定频率范围内进行模态分析,得到鞋底的固有振动特性。
在一较佳实施例中:所述步骤S3具体包括:
步骤S31:分别对三种不同工况下的鞋底施加一定频率范围内的简谐激励力;
步骤S32:对所述步骤S31中的鞋底进行稳态动力学分析,得到鞋底对于振动的响应情况;
步骤S33:将鞋底下表面各节点的速度作为历史输出变量,输出各节点速度到Abaqus结果文件.Odb文件中。
在一较佳实施例中:所述步骤S4具体包括:
步骤S41:利用Python编程语言从所述步骤S33中的.Odb结果文件中读取激励面及鞋底下表面各节点的响应速度到矩阵A,读取激励面各节点的节点激励力到矩阵B;
步骤S42:根据响应速度矩阵A与激励力矩阵B,并结合等效机械导纳计算公式得到鞋底激励面及鞋底下表面各节点在全频率范围内的等效机械导纳数据,并写入.csv表格文件中;
步骤S43:利用Python编程语言求取鞋底激励面及鞋底下表面各节点等效机械导纳数据在全频率范围内的平均值,并绘制其分布云图;
步骤S44:求取鞋底激励面及鞋底下表面各区域等效机械导纳的平均值,得到鞋底激励面及鞋底下表面不同区域等效机械导纳的数据,并绘制柱形图得到其等效机械导纳的分布情况。
在一较佳实施例中:所述步骤S5具体包括:
步骤S51:根据步骤S44所述的鞋底激励等效机械导纳与鞋底下表面不同区域等效机械导纳,并结合等效传递率计算公式计算得到鞋底下表面不同区域的等效传递率;
步骤S52:利用Python编程软件绘制鞋底下表面不同区域等效传递率的柱形图,得到鞋底下表面不同区域等效传递率分布情况。
由上述对本发明的描述可知,与现有技术相比,本发明具有如下有益效果:
1)本发明方法中,对模型的质量密度、杨氏模量和泊松比分别进行赋值,使模型与实物的材质相同,从而达到仿真的目的;
2)本发明方法中,采用等效机械导纳与等效振动传递率来对进行分析,能够清晰的反映鞋底对于振动的响应情况,及鞋底振动传递情况;
3)本发明方法中,利用Python编程语言绘制等效机械导纳分布云图及柱形图和等效振动传递率分布柱形图,能够很清晰的反映振动能量在鞋底各区域的分布情况,以及振动能量在鞋底各区域的放大衰减情况;
4)通过本发明对鞋底进行振动响应及传递特性分析,能够为制鞋业与脚用个人防护设备的开发提供评估标准与指导。
附图说明
图1为本发明方法的主要步骤流程示意图;
图2为鞋底三维实体模型图;
图3为鞋底区域划分图;
图4(a)为前脚掌激励工况下鞋底有限元模型图;
图4(b)为全脚掌激励工况下鞋底有限元模型图;
图4(c)为脚后跟激励工况下鞋底有限元模型图;
图5(a)为前脚掌激励工况下鞋底下表面等效机械导纳分布云图;
图5(b)为全脚掌激励工况下鞋底下表面等效机械导纳分布云图;
图5(c)为脚后跟激励工况下鞋底下表面等效机械导纳分布云图;
图6(a)为前脚掌激励工况下鞋底各区域等效机械导纳分布的柱状图;
图6(b)为全脚掌激励工况下鞋底各区域等效机械导纳分布的柱状图;
图6(c)为脚后跟激励工况下鞋底各区域等效机械导纳分布的柱状图;
图7(a)为前脚掌激励工况下鞋底下表面等效振动传递率分布的柱状图;
图7(b)为全脚掌激励工况下鞋底下表面等效振动传递率分布的柱状图;
图7(c)为脚后跟激励工况下鞋底下表面等效振动传递率分布的柱状图;
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述;显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例,基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要说明的是,术语“上”、“下”、“内”、“外”、“顶/底端”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“设置有”、“套设/接”、“连接”等,应做广义理解,例如“连接”,可以是壁挂连接,也可以是可拆卸连接,或一体地连接,可以是机械连接,也可以是电连接,可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通,对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
一种鞋底振动传递特性分析方法,如图1所示,包括:
步骤S1,在UG中建立简化的鞋底三维实体模型,如图2所示;
步骤S2,建立所述步骤S1中鞋底三维实体模型的有限元模型并进行模态分析,获得其固有振动特性,具体包括:
步骤S21:将鞋底三维实体模型导入Abaqus有限元分析软件中;
步骤S22:对鞋底三维实体模型按照解剖学原理划分成不同的区域,具体包括4个大的区域,再将这4个大的区域细分为9个小的区域;
具体的,4个大的区域包括:足跟区、足弓区、跖骨区、趾骨区。其中足跟区包括足跟内侧区HM、足跟外侧区HL,足弓区包括足弓内侧区MF1、足弓外侧区MF2,跖骨区包括第一跖骨区M1、第二、三跖骨区M23、第四、五跖骨区M45,趾骨区包括第一趾骨区T1、第二-五趾骨区T25共9个区域,如图3所示。
步骤S23:对所述步骤S22中的鞋底三维实体模型赋予材料属性、施加边界条件、划分网格单元得到鞋底有限元模型,其中材料属性包括质量密度、杨氏模量、泊松比和材料阻尼;
具体的,鞋底密度设置为1230kg/m3,弹性模量设置为4MPa,泊松比设置为0.4,材料阻尼设置为0.1;鞋底边界条件为无约束的自由边界条件;网格单元尺寸为4mm。
步骤S24:对所述步骤S23中的鞋底有限元模型在0-360Hz频率范围内进行模态分析,得到鞋底的固有振动特性。
步骤S3,基于所述步骤S2中的有限元模态分析,对鞋底进行稳态动力学分析,获得其在三种工况下对振动的响应情况,具体包括:
步骤S31:分别对三种不同工况下的鞋底施加频率范围为1-360Hz的简谐激励力,包括前脚掌激励、全脚掌激励及脚后跟激励三种工况;
具体的,继续沿用所述步骤S23中的鞋底材料属性,鞋底无约束的自由边界条件,网格单元尺寸。另外,分别对在鞋底上表面不同区域施加z方向向下,频率在1-360Hz内逐一增加、幅值为1N的单位简谐激励力,以模拟前脚掌激励、全脚掌激励及脚后跟激励三种工况。其中前脚掌激励施加在第二、三跖骨区域M23,如图4(a)所示;全脚掌激励施加在整个鞋底上表面,如图4(b)所示;脚后跟激励施加在足后跟圆形区域,如图4(c)所示。
步骤S32:对所述步骤S31中的模型进行稳态动力学分析,得到鞋底对于振动的响应情况;
步骤S33:将鞋底下表面各节点z方向的速度V3作为历史输出变量,输出各节点速度到Abaqus结果文件.Odb文件中。
步骤S4,利用Python编程语言计算并绘制鞋底等效机械导纳分布云图,及获得鞋底不同区域等效机械导纳的分布情况,具体包括:
步骤S41:利用Python编程语言从所述步骤S33中的.Odb结果文件中读取激励面及鞋底下表面各节点的响应速度到矩阵A,读取激励面各节点的节点激励力到矩阵B;
步骤S42:根据所述步骤S41中的响应速度矩阵A与激励力矩阵B,并结合等效机械导纳计算公式得到鞋底激励面及鞋底下表面各节点在全频率范围内的等效机械导纳数据,并写入.csv表格文件中;
步骤S43:利用Python编程语言求取鞋底激励面及鞋底下表面各节点等效机械导纳数据在全频率范围内的平均值,并绘制其分布云图;
具体的,前脚掌激励工况下鞋底下表面等效机械导纳分布云图如图5(a)所示,全脚掌激励工况下鞋底下表面等效机械导纳分布云图如图5(b)所示,脚后跟激励工况下鞋底下表面等效机械导纳分布云图如图5(c)所示。
步骤S44:求取鞋底激励面及鞋底下表面各区域等效机械导纳的平均值,得到鞋底激励面及鞋底下表面不同区域等效机械导纳的数据,并绘制柱形图得到其等效机械导纳的分布情况。
具体的,前脚掌激励工况下鞋底不同区域等效机械导纳分布的柱状图如图6(a)所示,全脚掌激励工况下鞋底不同区域等效机械导纳分布的柱状图如图6(b)所示,脚后跟激励工况下鞋底不同区域等效机械导纳分布的柱状图如图6(c)所示。
步骤S5,基于鞋底不同区域等效机械导纳的分布情况,利用Python编程语言计算并获得鞋底不同区域的等效振动传递率分布情况,具体包括:
步骤S51:根据步骤S44所述的鞋底激励等效机械导纳(即等效输入导纳)与鞋底下表面不同区域等效机械导纳(即等效传递导纳),并结合等效传递率计算公式计算得到鞋底下表面不同区域的等效传递率;
步骤S52:利用Python编程软件绘制鞋底下表面不同区域等效传递率的柱形图,得到鞋底下表面不同区域等效传递率分布情况。
具体的,前脚掌激励工况下鞋底不同区域等效振动传递率分布的柱状图如图7(a)所示,全脚掌激励工况下鞋底不同区域等效振动传递率分布的柱状图如图7(b)所示,脚后跟激励工况下鞋底不同区域等效振动传递率分布的柱状图如图7(c)所示。
以上所述,仅为本发明较佳的具体实施方式,但本发明的设计构思并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,利用此构思对本发明进行非实质性的改动,均属于侵犯本发明保护范围的行为。
Claims (5)
1.一种鞋底振动传递特性分析方法,其特征在于包括如下步骤:
步骤S1,建立鞋底三维实体模型;
步骤S2,建立鞋底三维实体模型的有限元模型并进行模态分析,获得鞋底三维实体模型的固有振动特性;
步骤S3,基于有限元模型的模态分析,对鞋底进行稳态动力学分析,获得鞋底在三种工况下对振动的响应情况;所述三种工况包括前脚掌激励、全脚掌激励及脚后跟激励;
步骤S4,利用Python编程语言计算并绘制鞋底等效机械导纳分布云图,及获得鞋底不同区域等效机械导纳的分布情况;
步骤S5,基于鞋底不同区域等效机械导纳的分布情况,利用Python编程语言计算并获得鞋底不同区域的等效振动传递率分布情况。
2.如权利要求1所述的一种鞋底振动传递特性分析方法,其特征在于,所述步骤S2具体包括:
步骤S21:将鞋底三维实体模型导入Abaqus有限元分析软件中;
步骤S22:对鞋底三维实体模型按照解剖学原理划分成不同的区域,具体包括4个大的区域,再将这4个大的区域细分为9个小的区域;4个大的区域包括:足跟区、足弓区、跖骨区、趾骨区;其中足跟区包括足跟内侧区HM、足跟外侧区HL,足弓区包括足弓内侧区MF1、足弓外侧区MF2,跖骨区包括第一跖骨区M1、第二、三跖骨区M23、第四、五跖骨区M45,趾骨区包括第一趾骨区T1、第二-五趾骨区T25共9个区域;
步骤S23:对所述步骤S22中的鞋底三维实体模型赋予材料属性、施加边界条件、划分网格单元得到鞋底三维实体模型的有限元模型;
步骤S24:对鞋底三维实体模型的有限元模型在一定频率范围内进行模态分析,得到鞋底的固有振动特性。
3.如权利要求2所述的一种鞋底振动传递特性分析方法,其特征在于,所述步骤S3具体包括:
步骤S31:分别对三种不同工况下的鞋底施加一定频率范围内的简谐激励力;
步骤S32:对所述步骤S31中的鞋底进行稳态动力学分析,得到鞋底对于振动的响应情况;
步骤S33:将鞋底下表面各节点的速度作为历史输出变量,输出各节点速度到Abaqus结果文件.Odb文件中。
4.如权利要求3所述的一种鞋底振动传递特性分析方法,其特征在于,所述步骤S4具体包括:
步骤S41:利用Python编程语言从所述步骤S33中的.Odb结果文件中读取激励面及鞋底下表面各节点的响应速度到矩阵A,读取激励面各节点的节点激励力到矩阵B;
步骤S42:根据响应速度矩阵A与激励力矩阵B,并结合等效机械导纳计算公式得到鞋底激励面及鞋底下表面各节点在全频率范围内的等效机械导纳数据,并写入.csv表格文件中;
步骤S43:利用Python编程语言求取鞋底激励面及鞋底下表面各节点等效机械导纳数据在全频率范围内的平均值,并绘制其分布云图;
步骤S44:求取鞋底激励面及鞋底下表面各区域等效机械导纳的平均值,得到鞋底激励面及鞋底下表面不同区域等效机械导纳的数据,并绘制柱形图得到其等效机械导纳的分布情况。
5.如权利要求4所述的一种鞋底振动传递特性分析方法,其特征在于,所述步骤S5具体包括:
步骤S51:根据步骤S44所述的鞋底激励等效机械导纳与鞋底下表面不同区域等效机械导纳,并结合等效传递率计算公式计算得到鞋底下表面不同区域的等效传递率;
步骤S52:利用Python编程软件绘制鞋底下表面不同区域等效传递率的柱形图,得到鞋底下表面不同区域等效传递率分布情况。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111145559.XA CN114330040A (zh) | 2021-09-28 | 2021-09-28 | 一种鞋底振动传递特性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111145559.XA CN114330040A (zh) | 2021-09-28 | 2021-09-28 | 一种鞋底振动传递特性分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114330040A true CN114330040A (zh) | 2022-04-12 |
Family
ID=81045545
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111145559.XA Pending CN114330040A (zh) | 2021-09-28 | 2021-09-28 | 一种鞋底振动传递特性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114330040A (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150101134A1 (en) * | 2013-10-16 | 2015-04-16 | Adidas Ag | Three-Dimensional Shoe Manufacturing |
CN111222217A (zh) * | 2018-11-23 | 2020-06-02 | 财团法人工业技术研究院 | 鞋垫设计方法及鞋垫设计系统 |
CN111274702A (zh) * | 2020-01-20 | 2020-06-12 | 华侨大学 | 一种鞋底足跟区缓冲结构鞋底模型设计方法 |
CN112084606A (zh) * | 2020-09-29 | 2020-12-15 | 华侨大学 | 一种球体多孔足跟区填充结构鞋底及其优化设计方法 |
CN112100853A (zh) * | 2020-09-16 | 2020-12-18 | 河海大学常州校区 | 一种面向支撑减压功能的鞋垫个性化设计方法 |
CN112199790A (zh) * | 2020-09-29 | 2021-01-08 | 华侨大学 | 一种正多面体多孔足跟区填充结构鞋底及其设计方法 |
CN112989669A (zh) * | 2021-04-01 | 2021-06-18 | 南通大学 | 一种针对糖尿病足垫的个性化设计方法 |
-
2021
- 2021-09-28 CN CN202111145559.XA patent/CN114330040A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150101134A1 (en) * | 2013-10-16 | 2015-04-16 | Adidas Ag | Three-Dimensional Shoe Manufacturing |
CN111222217A (zh) * | 2018-11-23 | 2020-06-02 | 财团法人工业技术研究院 | 鞋垫设计方法及鞋垫设计系统 |
CN111274702A (zh) * | 2020-01-20 | 2020-06-12 | 华侨大学 | 一种鞋底足跟区缓冲结构鞋底模型设计方法 |
CN112100853A (zh) * | 2020-09-16 | 2020-12-18 | 河海大学常州校区 | 一种面向支撑减压功能的鞋垫个性化设计方法 |
CN112084606A (zh) * | 2020-09-29 | 2020-12-15 | 华侨大学 | 一种球体多孔足跟区填充结构鞋底及其优化设计方法 |
CN112199790A (zh) * | 2020-09-29 | 2021-01-08 | 华侨大学 | 一种正多面体多孔足跟区填充结构鞋底及其设计方法 |
CN112989669A (zh) * | 2021-04-01 | 2021-06-18 | 南通大学 | 一种针对糖尿病足垫的个性化设计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113836771B (zh) | 一种多胞结构鞋底振动能量传递评估方法 | |
CN113657008A (zh) | 一种基于机械导纳的鞋底足跟区减振结构优化设计方法 | |
Shorten et al. | The ‘heel impact’force peak during running is neither ‘heel’nor ‘impact’and does not quantify shoe cushioning effects | |
CN102262699B (zh) | 基于无网格伽辽金与质点弹簧耦合的软组织形变仿真方法 | |
CA2783977A1 (en) | Method and system for partitioning parallel simulation models | |
GB2352036A (en) | Near wellbore modelling method and apparatus | |
CN105188530A (zh) | 步态仪和程序 | |
CN103838963A (zh) | 一种文胸压力舒适性评价方法 | |
CN105653786B (zh) | 一种考虑封闭框架剪切变形和刚性节点的内力计算方法 | |
CN113435082B (zh) | 一种用于骨折手术机器人的力学调控仿真系统 | |
US20200265943A1 (en) | Systems and methods for approximating musculoskeletal dynamics | |
CN102955881A (zh) | 一种计算集成电路芯片焊点热疲劳失效概率的方法 | |
CN109662392A (zh) | 一种基于落地运动的足部和鞋底有限元建模的测试方法 | |
CN112084606B (zh) | 一种球体多孔足跟区填充结构鞋底及其优化设计方法 | |
CN113033044A (zh) | 预测脚踝护具防护效果的建模方法 | |
CN112199790A (zh) | 一种正多面体多孔足跟区填充结构鞋底及其设计方法 | |
CN103622709A (zh) | 一种鞋垫式步行输入传感装置 | |
CN114330040A (zh) | 一种鞋底振动传递特性分析方法 | |
Ibbitt et al. | Re-scaling the topographic index to improve the representation of physical processes in catchment models | |
DeBerardinis et al. | Enhancing the accuracy of vertical ground reaction force measurement during walking using pressure-measuring insoles | |
JP2012093907A (ja) | 音響解析プログラム、音響解析方法および音響解析装置 | |
Wang et al. | A hybrid shallow water approach with unstructured triangular grids for urban flood modeling | |
Hamidon et al. | Prediction of future climate change for rainfall in the upper Kurau River Basin, Perak using statistical downscaling model (SDSM) | |
CN112715322A (zh) | 一种农业灌溉用水获取方法和装置 | |
CN213582170U (zh) | 一种正多面体多孔足跟区填充结构鞋底 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |