CN113657008A - Shoe sole heel area vibration reduction structure optimization design method based on mechanical admittance - Google Patents

Shoe sole heel area vibration reduction structure optimization design method based on mechanical admittance Download PDF

Info

Publication number
CN113657008A
CN113657008A CN202111143694.0A CN202111143694A CN113657008A CN 113657008 A CN113657008 A CN 113657008A CN 202111143694 A CN202111143694 A CN 202111143694A CN 113657008 A CN113657008 A CN 113657008A
Authority
CN
China
Prior art keywords
sole
cellular structure
mechanical admittance
area
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.)
Granted
Application number
CN202111143694.0A
Other languages
Chinese (zh)
Other versions
CN113657008B (en
Inventor
刘晓颖
郝艳华
刘斌
吴旭阳
王思成
叶成涛
江志城
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huaqiao University
Original Assignee
Huaqiao University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Huaqiao University filed Critical Huaqiao University
Priority to CN202111143694.0A priority Critical patent/CN113657008B/en
Publication of CN113657008A publication Critical patent/CN113657008A/en
Application granted granted Critical
Publication of CN113657008B publication Critical patent/CN113657008B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/16Customisation or personalisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/26Composites
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Footwear And Its Accessory, Manufacturing Method And Apparatuses (AREA)

Abstract

本发明公开了一种基于机械导纳的鞋底足跟区减振结构优化设计方法,包括:步骤S1,建立优化前的鞋底模型;步骤S2,以鞋底足跟区作为优化设计区域,将不同类型多胞结构填充入鞋底足跟区,得到不同类型多胞结构鞋底;步骤S3,构建所述步骤S2中多胞结构鞋底的有限元模型;步骤S4,对多胞结构鞋底进行稳态动力学分析,获得其对振动的响应情况;步骤S5,利用Python编程语言绘制多胞结构鞋底下表面等效机械导纳分布云图及多胞结构下表面等效机械导纳的平均值;步骤S7,对比不同类型多胞结构鞋底的等效机械导纳分布云图并绘制不同类型多胞结构鞋底下表面等效机械导纳平均值对比的柱状图。

Figure 202111143694

The invention discloses a method for optimizing the design of a vibration-damping structure in the heel area of a sole based on mechanical admittance, comprising: step S1, establishing a sole model before optimization; Filling the multi-cellular structure into the heel area of the sole to obtain different types of multi-cellular structure soles; step S3, constructing the finite element model of the multi-cellular structure sole in the step S2; step S4, performing steady-state dynamic analysis on the multi-cellular structure sole , obtain its response to vibration; step S5, use Python programming language to draw the distribution cloud map of the equivalent mechanical admittance on the lower surface of the multicellular structure sole and the average value of the equivalent mechanical admittance on the lower surface of the multicellular structure; step S7, compare different The distribution cloud map of the equivalent mechanical admittance of the soles with different types of multi-cellular structures and the histogram of the average value of the equivalent mechanical admittance on the lower surface of the soles with different types of multi-cellular structures is drawn.

Figure 202111143694

Description

Shoe sole heel area vibration reduction structure optimization design method based on mechanical admittance
Technical Field
The invention relates to the field of vibration analysis, in particular to a vibration energy transfer evaluation method for a multi-cell structure sole.
Background
The vibration transmitted from the foot to the human body can cause serious damage to the human body, and the sole can play a role in damping as a damping structure for isolating the impact vibration of the ground from the foot of the human body. The cellular structure is widely applied to the sole due to excellent vibration damping performance, and most sole-based researches adopt a test method, so that the cost is high, the test period is long, and time and labor are wasted.
Disclosure of Invention
The invention mainly aims to overcome the defects in the prior art and provides a method for optimally designing a vibration reduction structure of a sole heel area based on mechanical admittance, which can not only study the response and transmission conditions of vibration in a multi-cellular structure sole, but also provide guidance for the shoe industry and the development of personal protection equipment for feet by means of a numerical modeling and finite element analysis method.
The technical scheme adopted by the invention for solving the technical problems is as follows: a method for optimally designing a damping structure of a sole heel area based on mechanical admittance comprises the following steps:
step S1, establishing a simplified sole three-dimensional solid model before optimization in UG;
step S2, filling the multi-cell structures with different lattice types into the sole heel area by taking the sole heel area as an optimized design area to obtain the multi-cell structure soles with different lattice types;
Step S3, constructing a finite element model of the multi-cellular structure sole in finite element analysis software Abaqus;
step S4, carrying out steady-state dynamic analysis on the multi-cellular structure sole to obtain the response condition of the multi-cellular structure sole to vibration;
step S5, drawing an equivalent mechanical admittance distribution cloud picture of the lower surface of the sole with the multi-cell structure and an average value of equivalent mechanical admittance of the lower surface of the multi-cell structure by using a Python programming language;
step S6, repeating the steps S3-S5 to obtain equivalent mechanical admittance distribution cloud charts of the soles with the multi-cell structures of different lattice types and an average value of equivalent mechanical admittance on the lower surfaces of the soles;
and step S7, comparing the equivalent mechanical admittance distribution cloud charts of the soles with the different lattice types and drawing a histogram of the comparison of the equivalent mechanical admittance average values of the lower surfaces of the soles with the different lattice types and the multi-cellular structures to obtain an optimal sole model with the vibration reduction structure in the heel area of the sole.
In a preferred embodiment: the step S2 specifically includes:
step S21: constructing multi-cell structures of different lattice types, including Cross type, Diamond type, Grid type, Star type and X type;
and step S22, filling the multi-cell structures with different lattice types into the heel area of the sole by taking the heel area of the sole as an optimized design area to obtain the multi-cell structure soles with different lattice types.
In a preferred embodiment: the step S3 specifically includes:
step S31: importing the multi-cellular structure sole model into Abaqus finite element software;
step S32: the sole with the multi-cellular structure is divided into regions according to the anatomical principle: dividing the sole of the multi-cellular structure into four large areas, namely a heel area, an arch area, a metatarsal area and a phalange area, wherein the heel area comprises a heel inner area HM and a heel outer area HL, the arch area comprises an arch inner area MF1 and an arch outer area MF2, the metatarsal area comprises a first metatarsal area M1, a second metatarsal area M23, a fourth metatarsal area M45 and a fifth metatarsal area M45, and the phalange area comprises 9 areas, namely a first phalange area T1 and a second phalange area T25;
step S33: and giving material properties to the multi-cell structure sole in the step S32, applying boundary conditions, and dividing grid units to obtain a finite element model of the multi-cell structure sole.
In a preferred embodiment: the step S4 specifically includes:
step S41: carrying out finite element modal analysis on the finite element model of the multi-cellular structure sole within the frequency range of 0-360Hz to obtain the inherent vibration characteristic of the finite element model;
step S42, based on the result of finite element modal analysis, applying unit simple harmonic excitation force to the finite element model of the multi-cell structure, and carrying out modal-based steady-state dynamics analysis to obtain the response result of the finite element model to vibration;
Step S43, setting the historical output variable in Abaqus as the z-direction speed V3 of each node on the lower surface of the multi-cellular structure sole, and outputting each node speed to the result file, odb file.
In a preferred embodiment: the step S5 specifically includes:
step S51: reading speed data of each node in a frequency range of 1-360Hz from an Abaqus result file and an odb file by using a Python programming language;
step S52: calculating and obtaining equivalent mechanical admittance data of each node on the lower surface of the multi-cellular structure sole within the range of 1-360Hz by an equivalent mechanical admittance calculation formula according to the speed data of each node on the lower surface of the multi-cellular structure sole and the unit simple harmonic excitation force data applied to the multi-cellular structure sole;
step S53: calculating the average value of equivalent mechanical admittance of each node on the lower surface of the multi-cell structure in the frequency range of 1-360Hz, and storing the average value in a table file csv file;
step S54: reading data of the equivalent mechanical admittance average value of the lower surface of the sole with the multi-cell structure in the csv file by using a Python programming language, and drawing an equivalent mechanical admittance distribution cloud chart of the lower surface of the sole;
step S55: and calculating the average value of the equivalent mechanical admittance of the lower surface of the sole with the multi-cellular structure, and taking the average value as an evaluation index for measuring the vibration reduction performance of the sole with the multi-cellular structure.
In a preferred embodiment: the step S7 specifically includes:
step S71: drawing a histogram of comparison of the equivalent mechanical admittance average values of the multi-cellular structure soles with different lattice types in Python according to the average value of the equivalent mechanical admittance of the lower surface of the multi-cellular structure soles with different lattice types;
step S72: and preferably selecting a sole model with the damping structure in the sole heel area with the optimal damping performance according to equivalent mechanical admittance distribution cloud charts of soles with multi-cellular structures of different lattice types and a histogram of comparison of equivalent mechanical admittance average values.
Compared with the prior art, the technical scheme of the invention has the following beneficial effects:
1) in the method, the mass density, the Young modulus and the Poisson ratio of the model are respectively assigned, so that the material of the model is the same as that of a real object, and the aim of simulation is fulfilled;
2) in the method, the equivalent mechanical admittance nephogram and the histogram for comparing the equivalent mechanical admittance are adopted to jointly evaluate the damping effect of the sole, and the evaluation is carried out from two aspects of qualitative and quantitative, so that the evaluation effect is clearer and more effective;
3) in the method, the sole model with the damping structure in the sole heel area with excellent damping performance is optimized.
Drawings
FIG. 1 is a schematic flow chart of the main steps of the method of the present invention;
FIG. 2 is a diagram of a three-dimensional solid model of a sole before optimization;
FIG. 3(a) is a Cross type multi-cellular structure sole;
FIG. 3(b) is a sole of a Diamond type multi-cellular structure;
FIG. 3(c) is a Grid type multi-cellular structure sole;
FIG. 3(d) is a Star-shaped multi-cellular sole;
FIG. 3(e) is a view of an X-shaped multi-cellular sole;
FIG. 4 is a sole region partition view;
FIG. 5 is a diagram of a finite element model of a sole with a multi-cell structure under a heel excitation condition;
FIG. 6(a) is a cloud view of the equivalent mechanical admittance distribution of the lower surface of a sole of a Cross-type multi-cellular structure;
FIG. 6(b) is a cloud view of equivalent mechanical admittance distribution of the lower surface of a sole of a Diamond-type multi-cellular structure;
FIG. 6(c) is a cloud view of the equivalent mechanical admittance distribution of the lower surface of the sole of Grid-type multi-cellular structure;
FIG. 6(d) is a cloud of equivalent mechanical admittance distribution of the lower surface of the Star-type multi-cellular structure sole;
FIG. 6(e) is a cloud view of the equivalent mechanical admittance distribution of the lower surface of the sole of the X-shaped multi-cellular structure;
FIG. 7 is a histogram of equivalent mechanical admittance averages for the lower surface of a multi-cellular sole of different lattice types;
FIG. 8 is a diagram of an optimal model in a preferred embodiment of the present invention.
Detailed Description
The technical solution in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention; it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all embodiments, and all other embodiments obtained by those skilled in the art without any inventive work are within the scope of the present invention.
In the description of the present invention, it should be noted that the terms "upper", "lower", "inner", "outer", "top/bottom", and the like indicate orientations or positional relationships based on those shown in the drawings, and are only for convenience of description and simplification of description, but do not indicate or imply that the referred device or element must have a specific orientation, be constructed in a specific orientation, and be operated, and thus should not be construed as limiting the present invention. Furthermore, the terms "first" and "second" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
In the description of the present invention, it should be noted that, unless otherwise explicitly specified or limited, the terms "mounted," "disposed," "sleeved/connected," "connected," and the like, are used in a broad sense, and for example, "connected" may be a wall-mounted connection, a detachable connection, an integral connection, a mechanical connection, an electrical connection, a direct connection, an indirect connection through an intermediate medium, and a communication between two elements, and those skilled in the art will understand the specific meaning of the terms in the present invention specifically.
A method for optimally designing a damping structure of a heel area of a shoe sole based on mechanical admittance, as shown in fig. 1, comprises:
step S1, establishing a simplified sole three-dimensional solid model before optimization in UG, as shown in FIG. 2;
step S2, using the sole heel area as the optimization design area, filling different types of multi-cell structures into the sole heel area to obtain multi-cell structure soles with different lattice types, which specifically comprises:
step S21: constructing multi-cell structures of different lattice types, including Cross type, Diamond type, Grid type, Star type and X type;
specifically, five lattice types of multi-cell structures with the cell size of 8mm, the rod diameter of 2mm and the length, width and height of 48mm multiplied by 64mm multiplied by 16mm are constructed based on Grasshopper and UG secondary development, and the lattice types comprise Cross type, Diamond type, Grid type, Star type and X type.
Step S22, filling the multi-cell structures with different lattice types into the heel area of the sole by taking the heel area of the sole as an optimized design area to obtain the multi-cell structure soles with different lattice types;
specifically, a Cross type multi-cellular structure sole is shown in fig. 3(a), a Diamond type multi-cellular structure sole is shown in fig. 3(b), a Grid type multi-cellular structure sole is shown in fig. 3(c), a Star type multi-cellular structure sole is shown in fig. 3(d), and an X type multi-cellular structure sole is shown in fig. 3 (e).
Step S3, constructing a finite element model of the multi-cellular sole in step S2 in a finite element analysis software Abaqus, specifically including:
step S31: importing the multi-cellular structure sole model in the step S22 into Abaqus finite element software;
step S32: dividing the multi-cellular structure sole into regions according to the anatomical principle;
specifically, the sole with the multi-cellular structure is divided into four large areas, namely a heel area, an arch area, a metatarsal area and a phalange area, wherein the heel area comprises a medial heel area HM and a lateral heel area HL, the arch area comprises a medial arch area MF1 and a lateral arch area MF2, the metatarsal area comprises a first metatarsal area M1, a second metatarsal area M23, a third metatarsal area M23, a fourth metatarsal area M45 and the phalange area comprises 9 areas, namely a first phalange area T1 and a second phalange area T25, as shown in fig. 4.
Step S33: and giving material properties to the multi-cell structure sole in the step S32, applying boundary conditions, and dividing grid units to obtain a finite element model of the multi-cell structure sole.
Specifically, the material properties include: the density of the sole is 1230kg/m3, the elastic modulus is 4MPa, the Poisson ratio is 0.4, and the material damping is 0.1; the sole boundary condition is an unconstrained free boundary condition; the cell size when dividing the grid cells was 4 mm.
Step S4, carrying out steady state dynamic analysis on the multi-cellular structure sole to obtain the response condition of the multi-cellular structure sole to vibration, which specifically comprises the following steps:
step S41: carrying out finite element modal analysis on the finite element model of the multi-cell structure in the step S33 within the frequency range of 0-360Hz to obtain the inherent vibration characteristics of the finite element model;
step S42, based on the finite element modal analysis result of the step S41, unit simple harmonic excitation force is applied to the finite element model of the multi-cell structure, and steady state dynamic analysis based on the mode is carried out to obtain the response result of the finite element model to the vibration;
specifically, the excitation force applied to the sole with the cellular structure is a unit simple harmonic excitation force with the frequency gradually increased in the frequency range of 1-360Hz, the amplitude of the unit simple harmonic excitation force is 1N, and the direction of the unit simple harmonic excitation force is downward in the z direction, and the unit simple harmonic excitation force is applied to a circular area of the heel of the sole with the cellular structure, as shown in fig. 5.
Step S43, setting the historical output variable in Abaqus as the z-direction speed V3 of each node on the lower surface of the multi-cellular structure sole, and outputting each node speed to the result file, odb file.
Step S5, drawing a multi-cell structure sole lower surface equivalent mechanical admittance distribution cloud chart and an average value of the multi-cell structure lower surface equivalent mechanical admittance by using Python programming language, which specifically comprises the following steps:
Step S51: reading speed data of each node in a frequency range of 1-360Hz from an Abaqus result file and an odb file by using a Python programming language;
step S52: calculating and obtaining equivalent mechanical admittance data of each node on the lower surface of the multi-cellular structure sole within the range of 1-360Hz by an equivalent mechanical admittance calculation formula according to the speed data of each node on the lower surface of the multi-cellular structure sole and the unit simple harmonic excitation force data applied to the multi-cellular structure sole;
step S53: calculating the average value of equivalent mechanical admittance of each node on the lower surface of the multi-cell structure in the frequency range of 1-360Hz, and storing the average value in a table file csv file;
step S54: reading data of the equivalent mechanical admittance average value of the lower surface of the sole with the multi-cell structure in the csv file by using a Python programming language, and drawing an equivalent mechanical admittance distribution cloud chart of the lower surface of the sole;
step S55: and calculating the average value of the equivalent mechanical admittance of the lower surface of the sole with the multi-cellular structure, and taking the average value as an evaluation index for measuring the vibration reduction performance of the sole with the multi-cellular structure.
And S6, repeating the steps S3-S5 to obtain equivalent mechanical admittance distribution cloud charts of the soles with the multi-cell structures of different lattice types and an average value of equivalent mechanical admittance of the lower surfaces of the soles.
Specifically, a Cross type multi-cellular structure sole lower surface equivalent mechanical admittance distribution cloud chart is shown in fig. 6(a), a Diamond type multi-cellular structure sole lower surface equivalent mechanical admittance distribution cloud chart is shown in fig. 6(b), a Grid type multi-cellular structure sole lower surface equivalent mechanical admittance distribution cloud chart is shown in fig. 6(c), a Star type multi-cellular structure sole lower surface equivalent mechanical admittance distribution cloud chart is shown in fig. 6(d), and an X type multi-cellular structure sole lower surface equivalent mechanical admittance distribution cloud chart is shown in fig. 6 (e).
Step S7, comparing the equivalent mechanical admittance distribution cloud charts of the multi-cellular structure soles with different lattice types, drawing a histogram of the equivalent mechanical admittance average value comparison of the lower surfaces of the multi-cellular structure soles with different lattice types, selecting one of the optimal models as the optimal sole model with the damping structure in the heel area of the sole, and specifically comprising the following steps:
step S71: according to the average value of the equivalent mechanical admittances of the lower surfaces of the soles with the multi-cell structures of different lattice types in the step S6, drawing a histogram of comparison of the average values of the equivalent mechanical admittances of the soles with the multi-cell structures of different lattice types in Python;
step S72: according to the equivalent mechanical admittance distribution cloud charts of the soles with the cellular structures of different lattice types and the histogram of comparison of the equivalent mechanical admittance average values, a sole-heel-area vibration-reduction-structure sole model with the optimal vibration-reduction performance is preferably selected, as shown in fig. 7. In this example, the final selected sole model of the damping structure in the heel area of the sole has a Star-type lattice structure, a cell size of 8mm, a stem diameter of 3mm, and a length, width and height of 48mm x 64mm x 16mm, and is filled into the heel area of the sole as shown in fig. 8.
The above description is only a preferred embodiment of the present invention, but the design concept of the present invention is not limited thereto, and any person skilled in the art can make insubstantial changes in the technical scope of the present invention within the technical scope of the present invention, and the actions infringe the protection scope of the present invention are included in the present invention.

Claims (6)

1.一种基于机械导纳的鞋底足跟区减振结构优化设计方法,其特征在于包括如下步骤:1. a shoe sole heel area vibration-damping structure optimization design method based on mechanical admittance is characterized in that comprising the steps: 步骤S1,在UG中建立优化前简化的鞋底三维实体模型;Step S1, establishing a simplified three-dimensional solid model of the sole before optimization in UG; 步骤S2,以鞋底足跟区作为优化设计区域,将不同晶格类型的多胞结构填充入鞋底足跟区,得到不同晶格类型的多胞结构鞋底;Step S2, using the sole and heel area as the optimized design area, filling the multi-cellular structure of different lattice types into the sole and heel area, to obtain the multi-cellular structure sole of different lattice types; 步骤S3,在有限元分析软件Abaqus中构建所述多胞结构鞋底的有限元模型;Step S3, building the finite element model of the multi-cellular structure sole in the finite element analysis software Abaqus; 步骤S4,对多胞结构鞋底进行稳态动力学分析,获得其对振动的响应情况;In step S4, steady-state dynamics analysis is performed on the sole of the multi-cellular structure to obtain its response to vibration; 步骤S5,利用Python编程语言绘制多胞结构鞋底下表面等效机械导纳分布云图及多胞结构下表面等效机械导纳的平均值;Step S5, using the Python programming language to draw the distribution cloud map of the equivalent mechanical admittance on the lower surface of the multi-cellular structure sole and the average value of the equivalent mechanical admittance on the lower surface of the multi-cellular structure; 步骤S6,重复所述步骤S3-S5,得到不同晶格类型的多胞结构鞋底的等效机械导纳分布云图及鞋底下表面等效机械导纳的平均值;Step S6, repeating the steps S3-S5, to obtain the distribution nephogram of the equivalent mechanical admittance of the multicellular structure sole with different lattice types and the average value of the equivalent mechanical admittance of the lower surface of the sole; 步骤S7,对比不同晶格类型多胞结构鞋底的等效机械导纳分布云图并绘制不同晶格类型多胞结构鞋底下表面等效机械导纳平均值对比的柱状图,得到最优的鞋底足跟区减振结构鞋底模型。Step S7, comparing the distribution cloud maps of the equivalent mechanical admittance of the soles with different lattice types of multi-cellular structures and drawing a histogram of the average value of the equivalent mechanical admittances of the lower surfaces of the soles of the multi-cellular structures of different lattice types, to obtain the optimal sole foot. Sole model with vibration-absorbing structure in the heel area. 2.根据权利要求1所述的一种基于机械导纳的鞋底足跟区减振结构优化设计方法,其特征在于:所述步骤S2具体包括:2. a kind of based on mechanical admittance according to claim 1, it is characterized in that: described step S2 specifically comprises: 步骤S21:构建不同晶格类型的多胞结构,包括Cross型、Diamond型、Grid型、Star型、X型;Step S21: constructing multicellular structures of different lattice types, including Cross-type, Diamond-type, Grid-type, Star-type, and X-type; 步骤S22,以鞋底足跟区作为优化设计区域,将不同晶格类型的多胞结构填充入鞋底足跟区域,得到不同晶格类型的多胞结构鞋底。Step S22 , using the sole and heel area as an optimized design area, filling the multi-cellular structure of different lattice types into the sole and heel area to obtain a multi-cellular structure sole of different lattice types. 3.根据权利要求1所述的一种基于机械导纳的鞋底足跟区减振结构优化设计方法,其特征在于:所述步骤S3具体包括:3. a kind of based on mechanical admittance according to claim 1, it is characterized in that: described step S3 specifically comprises: 步骤S31:将所述多胞结构鞋底模型导入Abaqus有限元软件中;Step S31: import the multi-cellular structure sole model into Abaqus finite element software; 步骤S32:将多胞结构鞋底按照解剖学原理划分区域:将多胞结构鞋底分为足跟区、足弓区、跖骨区、趾骨区四个大的区域,其中足跟区包括足跟内侧区HM、足跟外侧区HL,足弓区包括足弓内侧区MF1、足弓外侧区MF2,跖骨区包括第一跖骨区M1、第二、三跖骨区M23、第四、五跖骨区M45,趾骨区包括第一趾骨区T1、第二-五趾骨区T25共9个区域;Step S32: Divide the multi-cellular structure sole into regions according to anatomical principles: divide the multi-cellular structure sole into four large regions: the heel region, the arch region, the metatarsal region, and the phalangeal region, wherein the heel region includes the inner heel region HM, lateral heel area HL, arch area includes medial arch area MF1, lateral arch area MF2, metatarsal area includes first metatarsal area M1, second and third metatarsal area M23, fourth and fifth metatarsal area M45, phalanges The area includes the first phalangeal area T1, the second-fifth phalangeal area T25, a total of 9 areas; 步骤S33:对所述步骤S32中的多胞结构鞋底赋予材料属性、施加边界条件、划分网格单元,得到多胞结构鞋底的有限元模型。Step S33: Endow material properties, apply boundary conditions, and divide mesh elements to the multi-cellular structure sole in step S32 to obtain a finite element model of the multi-cellular structure sole. 4.根据权利要求1所述的一种基于机械导纳的鞋底足跟区减振结构优化设计方法,其特征在于:所述步骤S4具体包括:4. a kind of based on mechanical admittance according to claim 1, it is characterized in that: described step S4 specifically comprises: 步骤S41:对所述多胞结构鞋底的有限元模型在0-360Hz频率范围内进行有限元模态分析,获得其固有振动特性;Step S41: performing a finite element modal analysis on the finite element model of the multi-cellular structure sole in the frequency range of 0-360 Hz to obtain its natural vibration characteristics; 步骤S42,基于有限元模态分析的结果,对多胞结构有限元模型施加单位简谐激励力,并进行基于模态的稳态动力学分析,获得其对振动的响应结果;Step S42, based on the result of the finite element modal analysis, apply a unit simple harmonic excitation force to the finite element model of the multi-cellular structure, and perform a modal-based steady-state dynamic analysis to obtain its response result to vibration; 步骤S43,在Abaqus中设置历史输出变量为多胞结构鞋底下表面各节点的z方向速度V3,输出各节点速度到结果文件.odb文件中。Step S43, the historical output variable is set in Abaqus as the z-direction velocity V3 of each node on the lower surface of the multi-cellular structure sole, and the velocity of each node is output to the result file .odb file. 5.根据权利要求1所述的一种基于机械导纳的鞋底足跟区减振结构优化设计方法,其特征在于:所述步骤S5具体包括:5. a kind of mechanical admittance-based optimal design method for vibration-damping structure of sole and heel region according to claim 1, is characterized in that: described step S5 specifically comprises: 步骤S51:利用Python编程语言从Abaqus结果文件.odb文件中读取各节点在1-360Hz频率范围内的速度数据;Step S51: using the Python programming language to read the speed data of each node in the frequency range of 1-360 Hz from the Abaqus result file .odb file; 步骤S52:根据多胞结构鞋底下表面各节点的速度数据,及对多胞结构鞋底施加的单位简谐激励力数据,由等效机械导纳计算公式计算并获得多胞结构鞋底下表面各节点在1-360Hz范围内的等效机械导纳数据;Step S52: According to the speed data of each node on the lower surface of the multi-cellular structure sole, and the unit simple harmonic excitation force data applied to the multi-cellular structure sole, calculate and obtain each node on the lower surface of the multi-cellular structure sole by calculating the equivalent mechanical admittance formula Equivalent mechanical admittance data in the range 1-360Hz; 步骤S53:求取多胞结构下表面各节点在1-360Hz频率范围内等效机械导纳的平均值,并将其存入表格文件.csv文件中;Step S53: Obtain the average value of the equivalent mechanical admittance of each node on the lower surface of the multicellular structure in the frequency range of 1-360 Hz, and store it in the table file .csv file; 步骤S54:利用Python编程语言读取表格文件.csv文件中的多胞结构鞋底下表面等效机械导纳平均值的数据,并绘制鞋底下表面等效机械导纳分布云图;Step S54: use the Python programming language to read the data of the average value of the equivalent mechanical admittance on the lower surface of the multi-cellular structure sole in the table file .csv file, and draw the distribution cloud diagram of the equivalent mechanical admittance on the lower surface of the sole; 步骤S55:求取多胞结构鞋底下表面等效机械导纳的平均值,作为衡量多胞结构鞋底减振性能的评价指标。Step S55: Obtain the average value of the equivalent mechanical admittance of the lower surface of the sole with the multi-cellular structure, and use it as an evaluation index for measuring the vibration reduction performance of the sole with the multi-cellular structure. 6.根据权利要求1所述的一种基于机械导纳的鞋底足跟区减振结构优化设计方法,其特征在于:所述步骤S7具体包括:6. a kind of mechanical admittance-based optimal design method for vibration-damping structure of sole heel area according to claim 1, is characterized in that: described step S7 specifically comprises: 步骤S71:根据所述不同晶格类型的多胞结构鞋底下表面等效机械导纳的平均值,在Python中绘制不同晶格类型多胞结构鞋底等效机械导纳平均值对比的柱状图;Step S71 : according to the average value of the equivalent mechanical admittance of the lower surface of the multicellular structure sole of the different lattice types, draw a bar chart comparing the average value of the equivalent mechanical admittance of the multicellular structure sole with different lattice types in Python; 步骤S72:根据不同晶格类型的多胞结构鞋底的等效机械导纳分布云图,并结合其等效机械导纳平均值对比的柱状图,优选出一种减振性能最优的鞋底足跟区减振结构鞋底模型。Step S72: According to the distribution cloud map of the equivalent mechanical admittance of the multi-cellular structure sole of different lattice types, and combined with the histogram of the average value of the equivalent mechanical admittance, a sole and heel with the best vibration damping performance are selected. Zone damping structure sole model.
CN202111143694.0A 2021-09-28 2021-09-28 An optimization design method for vibration reduction structure in the heel area of the sole based on mechanical admittance Active CN113657008B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111143694.0A CN113657008B (en) 2021-09-28 2021-09-28 An optimization design method for vibration reduction structure in the heel area of the sole based on mechanical admittance

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111143694.0A CN113657008B (en) 2021-09-28 2021-09-28 An optimization design method for vibration reduction structure in the heel area of the sole based on mechanical admittance

Publications (2)

Publication Number Publication Date
CN113657008A true CN113657008A (en) 2021-11-16
CN113657008B CN113657008B (en) 2025-02-07

Family

ID=78494459

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111143694.0A Active CN113657008B (en) 2021-09-28 2021-09-28 An optimization design method for vibration reduction structure in the heel area of the sole based on mechanical admittance

Country Status (1)

Country Link
CN (1) CN113657008B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114580041A (en) * 2022-02-25 2022-06-03 华侨大学 An optimization design method for vibration-damping sole structure based on variable-size lattice filling
CN114722679A (en) * 2022-04-28 2022-07-08 华侨大学 An optimization method for the toe of a pointe shoe based on a multicellular structure
CN114970167A (en) * 2022-05-31 2022-08-30 华侨大学 Design method and device for lattice variable density topological optimization structure of inner cavity of sole

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180271211A1 (en) * 2017-03-27 2018-09-27 Adidas Ag Footwear midsole with warped lattice structure and method of making the same
CN111274702A (en) * 2020-01-20 2020-06-12 华侨大学 A design method of a sole model for a buffer structure in the heel area of the sole
CN112084606A (en) * 2020-09-29 2020-12-15 华侨大学 Sole with spherical porous heel area filling structure and optimal design method thereof
CN112199790A (en) * 2020-09-29 2021-01-08 华侨大学 A regular polyhedron porous heel area filling structure sole and design method thereof
CN212782026U (en) * 2020-09-29 2021-03-23 华侨大学 Sole with spherical porous heel area filling structure
CN112861410A (en) * 2021-03-15 2021-05-28 上海理工大学 Medical image-based orthopedic insole manufacturing method
CN213582170U (en) * 2020-09-29 2021-06-29 华侨大学 A regular polyhedron porous heel area filling structure sole
CN113190917A (en) * 2021-04-22 2021-07-30 武汉理工大学 Multi-cell energy absorption structure combination optimization design method

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20180271211A1 (en) * 2017-03-27 2018-09-27 Adidas Ag Footwear midsole with warped lattice structure and method of making the same
CN111274702A (en) * 2020-01-20 2020-06-12 华侨大学 A design method of a sole model for a buffer structure in the heel area of the sole
CN112084606A (en) * 2020-09-29 2020-12-15 华侨大学 Sole with spherical porous heel area filling structure and optimal design method thereof
CN112199790A (en) * 2020-09-29 2021-01-08 华侨大学 A regular polyhedron porous heel area filling structure sole and design method thereof
CN212782026U (en) * 2020-09-29 2021-03-23 华侨大学 Sole with spherical porous heel area filling structure
CN213582170U (en) * 2020-09-29 2021-06-29 华侨大学 A regular polyhedron porous heel area filling structure sole
CN112861410A (en) * 2021-03-15 2021-05-28 上海理工大学 Medical image-based orthopedic insole manufacturing method
CN113190917A (en) * 2021-04-22 2021-07-30 武汉理工大学 Multi-cell energy absorption structure combination optimization design method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
凌静秀 等: "多胞鞋底结构性能分析及胞元参数影响", 《福建工程学院院报》, vol. 17, no. 6, 25 December 2019 (2019-12-25), pages 561 - 563 *
柯思成 等: "有限元建模技术在预测鞋底减震性能中的应用研究", 《皮革科学与工程》, vol. 30, no. 1, 4 February 2020 (2020-02-04), pages 68 - 73 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114580041A (en) * 2022-02-25 2022-06-03 华侨大学 An optimization design method for vibration-damping sole structure based on variable-size lattice filling
CN114580041B (en) * 2022-02-25 2024-08-16 华侨大学 An optimization design method for vibration-damping sole structure based on variable-size lattice filling
CN114722679A (en) * 2022-04-28 2022-07-08 华侨大学 An optimization method for the toe of a pointe shoe based on a multicellular structure
CN114970167A (en) * 2022-05-31 2022-08-30 华侨大学 Design method and device for lattice variable density topological optimization structure of inner cavity of sole

Also Published As

Publication number Publication date
CN113657008B (en) 2025-02-07

Similar Documents

Publication Publication Date Title
CN113657008A (en) Shoe sole heel area vibration reduction structure optimization design method based on mechanical admittance
CN113836771B (en) Multi-cell structure sole vibration energy transmission evaluation method
CN102262699A (en) Soft tissue deformation simulation method based on coupling of mesh-free Galerkin and mass spring
CN109583088B (en) Method for Determination of Resilience Modulus of Non-uniform Spatial Distribution in Embankment Varying with Humidity
CN113821958A (en) A kind of optimal design method for cushioning multicellular sole structure
CN110289103B (en) Soft tissue local compression deformation simulation method
CN112199790B (en) Regular polyhedron porous heel area filling structure sole and design method thereof
CN105893716B (en) A kind of structural break Multidisciplinary systems analysis method based on fractal theory
CN109388886A (en) A kind of calculation method of porous material acoustic absorptivity
CN108661089B (en) A numerical analysis method for ultimate expansion shear force of pile foundation in expansive soil area
CN110457785A (en) A material information mapping method based on material point method for large deformation response of structure
CN112084606A (en) Sole with spherical porous heel area filling structure and optimal design method thereof
CN112715322A (en) Method and device for obtaining agricultural irrigation water
CN109101752B (en) Method for calculating natural vibration frequency of local structure of complex hydraulic structure
CN111341449B (en) Simulation method for virtual vascular intervention operation training
CN103065015B (en) A kind of bearing structure low-carbon (LC) material-saving method for designing based on internal force path geometry form
CN106503472B (en) A kind of equivalent time domain model building method considering soil with blower fan system dynamic interaction
CN112075989A (en) Regular polyhedron porous filling structure calcaneus prosthesis and optimal design method thereof
CN213582170U (en) A regular polyhedron porous heel area filling structure sole
CN102163263B (en) Fan Blade Vibration Displacement and Its Weibull Distribution Fitting Method
CN114330040B (en) Sole vibration transfer characteristic analysis method
CN108877944B (en) Virtual cutting method based on grid model including Kelvin viscoelastic model
CN107092730A (en) Suitable for the three-dimensional infinite element Artificial Boundaries method for building up of Explicit Analysis
CN116992796A (en) An adaptive low-dissipation SPH-HLLC Riemann solver coupling algorithm
CN114970239B (en) Multi-type monitoring data measuring point arrangement method, equipment and medium based on Bayesian system identification and heuristic deep reinforcement learning

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
GR01 Patent grant