一种鞋底振动传递特性分析方法技术

本发明专利技术公开了一种鞋底振动传递特性分析方法，包括：步骤S1，建立鞋底三维实体模型；步骤S2，建立所述步骤S1中鞋底模型的有限元模型并进行模态分析，获得其固有振动特性；步骤S3，基于所述步骤S2中的有限元模态分析，对鞋底进行稳态动力学分析，获得其在三种工况下对振动的响应情况；步骤S4，利用Python编程语言计算并绘制鞋底等效机械导纳分布云图，及获得鞋底不同区域等效机械导纳的分布情况；步骤S5，基于鞋底不同区域等效机械导纳的分布情况，利用Python编程语言计算并获得鞋底不同区域的等效振动传递率分布情况。本发明专利技术可以利用有限元分析的方法获得鞋底对于振动的响应情况及振动在鞋底不同区域的传递情况。动在鞋底不同区域的传递情况。动在鞋底不同区域的传递情况。

【技术实现步骤摘要】
一种鞋底振动传递特性分析方法


[0001]本专利技术涉及振动分析领域，具体涉及一种鞋底振动传递特性分析方法。

技术介绍

[0002]大量研究表明，人体足部长期暴露于振动环境中，能够使他们患振动诱发白足的风险增加。而足部在接触振动时，鞋底是与来自地面的冲击振动最直接接触的部分，同时也会起到减振、隔振的作用。对于足部和鞋底振动的研究主要采用运动学测量的方法，有着诸多的缺点，如试验周期长、成本高等。

技术实现思路

[0003]本专利技术的主要目的在于克服现有技术中的上述缺陷，提出一种鞋底振动传递特性分析方法，借助于数值建模和有限元分析的方法，不仅可以研究振动在鞋底中的响应及传递情况，还能为制鞋业与脚用个人防护设备的开发提供评估标准与指导。
[0004]为了解决上述的技术问题，本专利技术采用如下技术方案：一种鞋底振动传递特性分析方法，包括如下步骤：
[0005]步骤S1,建立鞋底三维实体模型；
[0006]步骤S2,建立鞋底三维实体模型的有限元模型并进行模态分析，获得鞋底三维实体模型的固有振动特性；
[0007]步骤S3，基于有限元模型的模态分析，对鞋底进行稳态动力学分析，获得鞋底在三种工况下对振动的响应情况；所述三种工况包括前脚掌激励、全脚掌激励及脚后跟激励；
[0008]步骤S4，利用Python编程语言计算并绘制鞋底等效机械导纳分布云图，及获得鞋底不同区域等效机械导纳的分布情况；
[0009]步骤S5，基于鞋底不同区域等效机械导纳的分布情况，利用Python编程语言计算并获得鞋底不同区域的等效振动传递率分布情况。
[0010]在一较佳实施例中：所述步骤S2具体包括：
[0011]步骤S21：将鞋底三维实体模型导入Abaqus有限元分析软件中；
[0012]步骤S22：对鞋底三维实体模型按照解剖学原理划分成不同的区域，具体包括4个大的区域，再将这4个大的区域细分为9个小的区域；4个大的区域包括：足跟区、足弓区、跖骨区、趾骨区；其中足跟区包括足跟内侧区HM、足跟外侧区HL，足弓区包括足弓内侧区MF1、足弓外侧区MF2，跖骨区包括第一跖骨区M1、第二、三跖骨区M23、第四、五跖骨区M45，趾骨区包括第一趾骨区T1、第二
‑
五趾骨区T25共9个区域；
[0013]步骤S23：对所述步骤S22中的鞋底三维实体模型赋予材料属性、施加边界条件、划分网格单元得到鞋底三维实体模型的有限元模型；
[0014]步骤S24：对鞋底三维实体模型的有限元模型在一定频率范围内进行模态分析，得到鞋底的固有振动特性。
[0015]在一较佳实施例中：所述步骤S3具体包括：
[0016]步骤S31：分别对三种不同工况下的鞋底施加一定频率范围内的简谐激励力；
[0017]步骤S32：对所述步骤S31中的鞋底进行稳态动力学分析，得到鞋底对于振动的响应情况；
[0018]步骤S33：将鞋底下表面各节点的速度作为历史输出变量，输出各节点速度到Abaqus结果文件.Odb文件中。
[0019]在一较佳实施例中：所述步骤S4具体包括：
[0020]步骤S41：利用Python编程语言从所述步骤S33中的.Odb结果文件中读取激励面及鞋底下表面各节点的响应速度到矩阵A，读取激励面各节点的节点激励力到矩阵B；
[0021]步骤S42：根据响应速度矩阵A与激励力矩阵B，并结合等效机械导纳计算公式得到鞋底激励面及鞋底下表面各节点在全频率范围内的等效机械导纳数据，并写入.csv表格文件中；
[0022]步骤S43：利用Python编程语言求取鞋底激励面及鞋底下表面各节点等效机械导纳数据在全频率范围内的平均值，并绘制其分布云图；
[0023]步骤S44：求取鞋底激励面及鞋底下表面各区域等效机械导纳的平均值，得到鞋底激励面及鞋底下表面不同区域等效机械导纳的数据，并绘制柱形图得到其等效机械导纳的分布情况。
[0024]在一较佳实施例中：所述步骤S5具体包括：
[0025]步骤S51：根据步骤S44所述的鞋底激励等效机械导纳与鞋底下表面不同区域等效机械导纳，并结合等效传递率计算公式计算得到鞋底下表面不同区域的等效传递率；
[0026]步骤S52：利用Python编程软件绘制鞋底下表面不同区域等效传递率的柱形图，得到鞋底下表面不同区域等效传递率分布情况。
[0027]由上述对本专利技术的描述可知，与现有技术相比，本专利技术具有如下有益效果：
[0028]1)本专利技术方法中，对模型的质量密度、杨氏模量和泊松比分别进行赋值，使模型与实物的材质相同，从而达到仿真的目的；
[0029]2)本专利技术方法中，采用等效机械导纳与等效振动传递率来对进行分析，能够清晰的反映鞋底对于振动的响应情况，及鞋底振动传递情况；
[0030]3)本专利技术方法中，利用Python编程语言绘制等效机械导纳分布云图及柱形图和等效振动传递率分布柱形图，能够很清晰的反映振动能量在鞋底各区域的分布情况，以及振动能量在鞋底各区域的放大衰减情况；
[0031]4)通过本专利技术对鞋底进行振动响应及传递特性分析，能够为制鞋业与脚用个人防护设备的开发提供评估标准与指导。
附图说明
[0032]图1为本专利技术方法的主要步骤流程示意图；
[0033]图2为鞋底三维实体模型图；
[0034]图3为鞋底区域划分图；
[0035]图4(a)为前脚掌激励工况下鞋底有限元模型图；
[0036]图4(b)为全脚掌激励工况下鞋底有限元模型图；
[0037]图4(c)为脚后跟激励工况下鞋底有限元模型图；
[0038]图5(a)为前脚掌激励工况下鞋底下表面等效机械导纳分布云图；
[0039]图5(b)为全脚掌激励工况下鞋底下表面等效机械导纳分布云图；
[0040]图5(c)为脚后跟激励工况下鞋底下表面等效机械导纳分布云图；
[0041]图6(a)为前脚掌激励工况下鞋底各区域等效机械导纳分布的柱状图；
[0042]图6(b)为全脚掌激励工况下鞋底各区域等效机械导纳分布的柱状图；
[0043]图6(c)为脚后跟激励工况下鞋底各区域等效机械导纳分布的柱状图；
[0044]图7(a)为前脚掌激励工况下鞋底下表面等效振动传递率分布的柱状图；
[0045]图7(b)为全脚掌激励工况下鞋底下表面等效振动传递率分布的柱状图；
[0046]图7(c)为脚后跟激励工况下鞋底下表面等效振动传递率分布的柱状图；
具体实施方式
[0047]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述；显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例，基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0048]在本专利技术的描述中，需要说明的是，术语本文档来自技高网...

【技术保护点】

【技术特征摘要】
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所述的一种鞋底振动...

【专利技术属性】
技术研发人员：刘晓颖，吴旭阳，黄贤伟，
申请(专利权)人：华侨大学，
类型：发明
国别省市：