CN113836771B - Multi-cell structure sole vibration energy transmission evaluation method - Google Patents
Multi-cell structure sole vibration energy transmission evaluation method Download PDFInfo
- Publication number
- CN113836771B CN113836771B CN202111145580.XA CN202111145580A CN113836771B CN 113836771 B CN113836771 B CN 113836771B CN 202111145580 A CN202111145580 A CN 202111145580A CN 113836771 B CN113836771 B CN 113836771B
- Authority
- CN
- China
- Prior art keywords
- sole
- cell structure
- equivalent
- cell
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
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
The invention discloses a multi-cell structure sole vibration energy transmission evaluation method, which comprises the following steps: step S1, an original sole three-dimensional solid model is established; s2, constructing multi-cell structures with different lattice types, and respectively filling the multi-cell structures into the heel area of the sole to obtain soles with multi-cell structures with different lattice types; step S3, constructing a finite element model of the multi-cell structure sole in the step S2; s4, steady-state dynamics analysis is carried out on the multi-cell structure sole, and the response condition of the multi-cell structure sole to vibration is obtained; s5, calculating and drawing a multi-cell structure sole equivalent mechanical admittance distribution cloud chart by using a Python programming language, and obtaining the distribution conditions of the multi-cell structure sole equivalent mechanical admittances in different areas; and S6, calculating and obtaining the distribution situation of the equivalent vibration transmissivities of different areas of the sole with the multicellular structure by using a Python programming language based on the distribution situation of the equivalent mechanical admittances of the different areas of the sole with the multicellular structure.
Description
Technical Field
The invention relates to the field of vibration analysis, in particular to a multi-cell structure sole vibration energy transmission evaluation method.
Background
The multi-cell structure is used for the design of the buffering structure of the sole because of the excellent buffering and energy absorbing effects. Research shows that the existence of the multi-cell structure can reduce the impact of the foot from the ground, and the effect of buffering and absorbing energy is achieved. While the studies of the multi-cell structure sole in terms of vibration have been rarely conducted. In addition, the current research on foot and sole vibration mainly adopts a kinematic measurement method, and has a plurality of defects such as long test period, high cost and the like.
Disclosure of Invention
The invention aims to overcome the defects in the prior art, and provides a multi-cell structure sole vibration energy transmission evaluation method which can not only study the response and transmission conditions of vibration in a multi-cell structure sole by means of numerical modeling and finite element analysis, but also provide guidance for the shoe industry and the development of foot personal protection equipment.
The invention adopts the following technical scheme: a method of evaluating vibrational energy transfer of a multi-cellular structure sole comprising:
step S1, an original sole three-dimensional solid model is established;
s2, constructing multi-cell structures with different lattice types, and respectively filling the multi-cell structures into the heel area of the sole to obtain soles with multi-cell structures with different lattice types;
s3, constructing a finite element model of the multi-cell structure sole;
s4, carrying out steady-state dynamics analysis on the finite element model of the multi-cell structure sole to obtain the response condition of the finite element model to vibration;
s5, calculating and drawing an equivalent mechanical admittance distribution cloud picture of the multi-cell structure sole by using a Python programming language, and obtaining the distribution situation of the equivalent mechanical admittances of different areas of the multi-cell structure sole;
step S6, calculating and obtaining the distribution situation of equivalent vibration transmissivities of different areas of the sole with the multicellular structure by using a Python programming language based on the distribution situation of equivalent mechanical admittances of the different areas of the sole with the multicellular structure;
and S7, repeatedly executing the steps S3-S6 on the multi-cell structure soles with different lattice types in the step S2 to obtain equivalent mechanical admittance distribution cloud patterns and distribution conditions of the multi-cell structure soles with different lattice types and distribution conditions of equivalent vibration transmissibility.
In a preferred embodiment: the step S2 specifically includes:
step S21: constructing different types of multicellular structures, including Grid-type multicellular structures, star-type multicellular structures and X-type multicellular structures;
step S22: and filling different types of multi-cell structures into the heel area of the sole by taking the heel area of the original sole three-dimensional solid model as a multi-cell structure filling area to obtain soles with different types of multi-cell structures.
In a preferred embodiment: the step S3 specifically includes:
step S31: introducing the multi-cell structured sole into Abaqus finite element analysis software;
step S32: the multi-cell structure sole is divided into different areas according to an anatomical principle, and particularly comprises 4 large areas of a heel area, an arch area, a metatarsal area and a phalangeal area; wherein the heel region comprises a medial heel region HM, a lateral heel region HL, the arch region comprises a medial arch region MF1, a lateral arch region MF2, the metatarsal region comprises a first metatarsal region M1, a second, a third metatarsal region M23, a fourth, a fifth metatarsal region M45, and the phalangeal region comprises 9 total regions of a first phalangeal region T1, a second-fifth phalangeal region T25;
step S33: and (3) endowing the multi-cell structure sole with material properties, applying boundary conditions and dividing grid cells in the step (S32) to obtain a finite element model of the multi-cell structure sole.
In a preferred embodiment: the step S4 specifically includes:
step S41: performing modal analysis on the finite element model of the multi-cell structure sole within a certain frequency range to obtain the inherent vibration characteristic of the multi-cell structure sole;
step S42: applying unit simple harmonic excitation force within a certain frequency range to an excitation surface of a finite element model of the multi-cell structure sole;
step S43: performing modal analysis based on the finite element model of the multi-cell structure sole in the step S41, and performing steady-state dynamics analysis on the multi-cell structure sole model in the step S42 to obtain the response condition of the multi-cell structure sole to vibration;
step S44: and taking the speed of each node on the surface of the multi-cell sole shoe as a history output variable, and outputting the speed of each node to an Abaqus result file and an odb file.
In a preferred embodiment: the step S5 specifically includes:
step S51: reading response speeds of all nodes on the excitation surface of the multi-cell structure sole and the lower surface of the sole from the odb file to a matrix A, and reading node excitation forces of all nodes on the excitation surface to a matrix B;
step S52: according to the response speed matrix A and the excitation force matrix B, combining an equivalent mechanical admittance calculation formula to obtain equivalent mechanical admittance data of each node of the excitation surface of the multi-cell structure sole and the lower surface of the sole in a full frequency range, and writing the equivalent mechanical admittance data into a csv table file;
step S53: calculating the average value of equivalent mechanical admittance data of each node of the excitation surface of the multi-cell structure sole and the lower surface of the sole in a full frequency range by using a Python programming language, and drawing a distribution cloud chart of the average value;
step S54: and (3) obtaining the average value of the equivalent mechanical admittances of the areas of the sole excitation surface and the lower surface of the sole of the multi-cell structure, and obtaining the data of the equivalent mechanical admittances of the areas of the sole excitation surface and the lower surface of the sole of the multi-cell structure.
In a preferred embodiment: the step S6 specifically includes:
step S61: equivalent input admittance data are obtained from equivalent mechanical admittance of the excitation surface of the multi-cell structure sole, and equivalent transfer admittance data are obtained from equivalent mechanical admittance of different areas of the lower surface of the sole;
step S62: and calculating the equivalent vibration transmissivities of different areas on the lower surface of the multi-cell structure sole according to the equivalent input admittance data and the equivalent transmission admittance data of the multi-cell structure sole and by combining an equivalent vibration transmissibility calculation formula.
In a preferred embodiment: the step S7 specifically includes:
step S71: repeatedly executing the steps S3-S6 on the multi-cell structure soles with different lattice types to obtain equivalent mechanical admittance distribution cloud charts of the multi-cell structure soles with different lattice types and equivalent mechanical admittance and equivalent vibration transmissibility data of different areas of the excitation surfaces of the multi-cell structure soles with different lattice types and the lower surfaces of the soles;
step S72: drawing a histogram of equivalent mechanical admittance and equivalent vibration transmissibility distribution of the soles with different types of multicellular structures in different areas of the soles according to the equivalent mechanical admittance and equivalent vibration transmissibility data of the excitation surfaces of the soles with different multicellular structures with different lattice types and different areas of the lower surfaces of the soles;
step S73: and according to the equivalent mechanical admittance and equivalent vibration transmissibility data of different areas of the lower surface of the multi-cell structure sole with different lattice types, calculating the average value of the equivalent mechanical admittance and the equivalent vibration transmissibility of the whole lower surface of the multi-cell structure sole with different lattice types, and drawing a histogram for comparison.
As can be seen from the above description of the present invention, compared with the prior art, the present invention has the following advantages:
1) In the method, the mass density, young modulus and Poisson's ratio of the model are respectively assigned, so that the model and a physical material are the same, and the purpose of simulation is achieved;
2) In the method, the equivalent mechanical admittance and the equivalent vibration transmissibility are adopted to analyze the vibration energy distribution and the transmission condition, so that the response condition of the multi-cell structure sole to vibration and the vibration energy transmission condition of the sole can be clearly reflected;
3) In the method, the Python programming language is utilized to draw the equivalent mechanical admittance distribution cloud chart, the column chart and the equivalent vibration transmissibility distribution column chart, so that the distribution condition of vibration energy in different areas of soles with different multicellular structures and the amplification and attenuation condition of vibration energy in different areas of soles with multicellular structures can be clearly reflected;
4) The invention can evaluate the vibration energy distribution and transmission conditions of soles with different types of multicellular structures, and can provide guidance for the shoe industry and the development of foot personal protection equipment.
Drawings
FIG. 1 is a schematic flow chart of main steps of the method of the invention;
FIG. 2 is a three-dimensional solid model diagram of a sole;
FIG. 3 is a cutaway view of a sole region;
FIG. 4 is a finite element model diagram of the sole under heel excitation conditions;
FIG. 5 (a) is a cloud plot of equivalent mechanical admittance distribution of the lower surface of the original sole;
FIG. 5 (b) is a cloud chart of equivalent mechanical admittance distribution of the lower surface of the Grid-type multicellular structure shoe sole;
FIG. 5 (c) is a cloud chart of equivalent mechanical admittance distribution of the lower surface of the Star multicellular structure sole;
FIG. 5 (d) is a cloud chart of equivalent mechanical admittance distribution of the lower surface of the X-shaped multicellular structure sole;
FIG. 6 (a) is a bar graph of equivalent mechanical admittance distribution for different regions of different types of multi-cellular structure soles;
FIG. 6 (b) is a bar graph of equivalent vibration transmissivity distribution at different areas of the underside of a different type of multi-cellular structure shoe sole;
FIG. 7 (a) is a bar graph comparing the equivalent mechanical admittance of the excitation surface and the lower surface of a different type of multi-cellular structure sole;
figure 7 (b) is a bar graph of equivalent vibration transmissibility versus lower surface of different types of multicellular structures.
Detailed Description
The technical solutions 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 apparent that the described embodiments are only some embodiments of the present invention, not all embodiments, and that all other embodiments obtained by persons of ordinary skill in the art without making creative efforts based on the embodiments in the present invention are within the protection scope of the present invention.
In the description of the present invention, it should be noted that the positional or positional relationship indicated by the terms such as "upper", "lower", "inner", "outer", "top/bottom", etc. are based on the positional or positional relationship shown in the drawings, are merely for convenience of describing the present invention and simplifying the description, and do not indicate or imply that the apparatus or elements referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus should not be construed as limiting the present invention. Furthermore, the terms "first," "second," and the like, 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, unless explicitly specified and limited otherwise, the terms "mounted," configured to, "" engaged with, "" connected to, "and the like are to be construed broadly, and may be, for example," connected to, "wall-mounted," connected to, removably connected to, or integrally connected to, mechanically connected to, electrically connected to, directly connected to, or indirectly connected to, through an intermediary, and may be in communication with each other between two elements, as will be apparent to those of ordinary skill in the art, in view of the detailed description of the terms herein.
A method for evaluating vibration energy transfer of a multi-cellular structure sole, as shown in fig. 1, comprising:
step S1, a simplified original sole three-dimensional solid model is built in UG, as shown in FIG. 2 (a);
step S2, constructing multi-cell structures with different lattice types, and respectively filling the multi-cell structures into the heel area of the sole to obtain soles with the multi-cell structures with different lattice types, wherein the steps specifically comprise:
step S21: constructing different types of multicellular structures, including Grid-type multicellular structures, star-type multicellular structures and X-type multicellular structures;
specifically, three lattice type multicellular structures with cell size of 8mm, rod diameter of 2mm, length width of 48mm×64mm×16mm are constructed based on secondary development of Grasshopper and UG, and the lattice types comprise Grid type, star type and X type.
Step S22: and filling different types of multi-cell structures into the heel area of the sole by taking the heel area of the original sole as a multi-cell structure filling area to obtain soles with different types of multi-cell structures.
Specifically, a cuboid with the length, width and height consistent with those of the multicellular structure is dug out in the heel area of the original sole, so that the original sole with the cavity is obtained; secondly, filling different types of multicellular structures into cuboid cavities in the heel area of the sole respectively; and finally, carrying out Boolean summation operation on the multi-cell structure and the original sole containing the cavity to obtain soles with different types of multi-cell structures. The Grid-type multi-cell structure sole is shown in fig. 2 (b), the Star-type multi-cell structure sole is shown in fig. 2 (c), and the X-type multi-cell structure sole is shown in fig. 2 (d).
Step S3, constructing a finite element model of the multi-cell structure sole in step S2, which specifically comprises the following steps:
step S31: importing the three-dimensional solid model of the multi-cell structure sole in the step S22 into Abaqus finite element analysis software;
step S32: the multi-cell structure sole is divided into different areas according to an anatomical principle, and specifically comprises 4 large areas such as a heel area, an arch area, a metatarsal area, a phalange area and the like, wherein the 4 large areas are further subdivided into 9 small areas;
specifically, the heel region includes a medial heel region HM, a lateral heel region HL, the arch region includes a medial arch region MF1, a lateral arch region MF2, the metatarsal region includes a first metatarsal region M1, a second, third metatarsal region M23, a fourth, fifth metatarsal region M45, and the phalangeal region includes 9 total regions of a first phalangeal region T1, a second-fifth phalangeal region T25, as shown in fig. 3.
Step S33: and (3) endowing the multi-cell structure sole with material properties, applying boundary conditions and dividing grid cells in the step (S32) to obtain a finite element model of the multi-cell structure sole.
Specifically, the sole density is set to 1230kg/m3, the elastic modulus is set to 4MPa, the Poisson's ratio is set to 0.4, and the material damping is set to 0.1; the sole boundary condition is an unconstrained free boundary condition; the mesh cell size was 4mm.
Step S4, steady-state dynamics analysis is carried out on the multi-cell structure sole to obtain the response condition of the multi-cell structure sole to vibration, and the method specifically comprises the following steps:
step S41: performing modal analysis on the finite element model of the multi-cell structure sole in the step S33 within the frequency range of 0-360Hz to obtain the inherent vibration characteristic of the multi-cell structure sole;
step S42: applying unit simple harmonic excitation force within the frequency range of 1-360Hz to the excitation surface of the multi-cell structure sole finite element model in the step S33;
specifically, a z-direction downward force is applied to the circular area of the heel of the multi-cell sole, the frequency is increased one by one within 1-360Hz, and the unit simple harmonic excitation force with the amplitude of 1N is applied to simulate the heel excitation working condition, as shown in fig. 4.
Step S43: based on the finite element modal analysis of the multi-cell structure sole in the step S41, steady-state dynamic analysis is carried out on the multi-cell structure sole model in the step S42, so that the response condition of the multi-cell structure sole to vibration is obtained;
step S44: and taking the speed of each node on the surface of the multi-cell sole shoe as a history output variable, and outputting the speed of each node to an Abaqus result file and an odb file.
Step S5, calculating and drawing a multi-cell structure sole equivalent mechanical admittance distribution cloud chart by using a Python programming language, and obtaining the distribution situation of the multi-cell structure sole equivalent mechanical admittance in different areas, wherein the method specifically comprises the following steps:
step S51: reading response speeds of all nodes on the excitation surface and the lower surface of the sole of the multi-cell structure from the odb result file in the step S44 to a matrix A by using a Python programming language, and reading node excitation forces of all nodes on the excitation surface to a matrix B;
step S52: according to the response speed matrix A and the excitation force matrix B in the step S51, combining an equivalent mechanical admittance calculation formula to obtain equivalent mechanical admittance data of each node of the excitation surface of the multi-cell structure sole and the lower surface of the sole in a full frequency range, and writing the equivalent mechanical admittance data into a csv table file;
step S53: calculating the average value of equivalent mechanical admittance data of each node of the excitation surface of the multi-cell structure sole and the lower surface of the sole in a full frequency range by using a Python programming language, and drawing a distribution cloud chart of the average value;
step S54: and (3) obtaining the average value of the equivalent mechanical admittances of the areas of the sole excitation surface and the lower surface of the sole of the multi-cell structure, and obtaining the data of the equivalent mechanical admittances of the areas of the sole excitation surface and the lower surface of the sole of the multi-cell structure.
Step S6, calculating and obtaining the distribution situation of the equivalent vibration transmissivities of different areas of the sole with the multicellular structure by using a Python programming language based on the distribution situation of the equivalent mechanical admittances of the different areas of the sole with the multicellular structure, wherein the distribution situation comprises the following specific steps:
step S61: obtaining equivalent input admittance data from equivalent mechanical admittances of the excitation surface of the multi-cell structure sole in the step S54, and obtaining equivalent transfer admittance data from equivalent mechanical admittances of different areas of the lower surface of the sole;
step S62: and according to the equivalent input admittance data and the equivalent transmission admittance data of the multi-cell structure sole in the step S61, and combining an equivalent vibration transmission rate calculation formula, calculating to obtain the equivalent vibration transmission rates of different areas of the lower surface of the multi-cell structure sole.
Step S7, repeatedly executing the steps S3-S6 on the soles with different types of multi-cell structures in the step S2 to obtain equivalent mechanical admittance distribution cloud patterns and distribution conditions of equivalent vibration transmissivities of the soles with different types of multi-cell structures, wherein the method specifically comprises the following steps:
step S71: repeatedly executing the steps S3-S6 on the soles with different types of multi-cell structures in the step S2 to obtain equivalent mechanical admittance distribution cloud charts of the soles with different types of multi-cell structures and equivalent mechanical admittance and equivalent vibration transmissibility data of different areas of the excitation surfaces of the soles with different types of multi-cell structures and the lower surfaces of the soles;
specifically, the original sole lower surface equivalent mechanical admittance distribution cloud chart is shown in fig. 5 (a), the Grid type multi-cell structure sole lower surface equivalent mechanical admittance distribution cloud chart is shown in fig. 5 (b), the Star type multi-cell structure sole lower surface equivalent mechanical admittance distribution cloud chart is shown in fig. 5 (c), and the X type multi-cell structure sole lower surface equivalent mechanical admittance distribution cloud chart is shown in fig. 5 (d).
Step S72: drawing a histogram of the equivalent mechanical admittance and the equivalent vibration transmissibility distribution of the soles with different types of multi-cell structures in different areas of the soles according to the equivalent mechanical admittance and the equivalent vibration transmissibility data of the sole excitation surfaces with different types of multi-cell structures and the different areas of the lower surface of the sole in the step S71;
specifically, a histogram of equivalent mechanical admittance distribution of different areas of the sole with different types of multi-cellular structures is shown in fig. 6 (a), and a histogram of equivalent vibration transmissivity distribution of different areas of the lower surface of the sole with different types of multi-cellular structures is shown in fig. 6 (b).
Step S73: according to the data of equivalent mechanical admittance and equivalent vibration transmissibility of different areas of the lower surface of the sole with different types of multi-cell structures in the step S71, the average value of the equivalent mechanical admittance and equivalent vibration transmissibility of the whole lower surface of the sole with different types of multi-cell structures is obtained, and a histogram is drawn for comparison.
Specifically, a histogram of comparison of equivalent mechanical admittances of the lower surfaces of the different types of multi-cellular structure shoes is shown in fig. 7 (a), and a histogram of comparison of equivalent vibration transmissivities of the lower surfaces of the different types of multi-cellular structure shoes is shown in fig. 7 (b).
The foregoing 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 will be able to make insubstantial modifications of the present invention within the scope of the present invention disclosed herein by this concept, which falls within the actions of invading the protection scope of the present invention.
Claims (4)
1. The multi-cell structure sole vibration energy transmission evaluation method is characterized by comprising the following steps of:
step S1, an original sole three-dimensional solid model is established;
s2, constructing multi-cell structures with different lattice types, and respectively filling the multi-cell structures into the heel area of the sole to obtain soles with multi-cell structures with different lattice types;
s3, constructing a finite element model of the multi-cell structure sole;
s4, carrying out steady-state dynamics analysis on the finite element model of the multi-cell structure sole to obtain the response condition of the finite element model to vibration;
s5, calculating and drawing an equivalent mechanical admittance distribution cloud picture of the multi-cell structure sole by using a Python programming language, and obtaining the distribution situation of the equivalent mechanical admittances of different areas of the multi-cell structure sole;
step S6, calculating and obtaining the distribution situation of equivalent vibration transmissivities of different areas of the sole with the multicellular structure by using a Python programming language based on the distribution situation of equivalent mechanical admittances of the different areas of the sole with the multicellular structure;
step S7, repeatedly executing the steps S3-S6 on the multi-cell structure soles with different lattice types in the step S2 to obtain equivalent mechanical admittance distribution cloud patterns and distribution conditions of equivalent vibration transmissibility of the multi-cell structure soles with different lattice types;
the step S2 specifically includes:
step S21: constructing multi-cell structures with different lattice types, including Grid type multi-cell structures, star type multi-cell structures and X type multi-cell structures;
step S22: taking the heel area of the original sole three-dimensional solid model as a multi-cell structure filling area, and filling multi-cell structures with different lattice types into the sole heel area to obtain soles with different multi-cell structures;
the step S3 specifically includes:
step S31: introducing the multi-cell structured sole into Abaqus finite element analysis software;
step S32: the multi-cellular structure sole is divided into 4 large areas according to an anatomical principle, namely a heel area, an arch area, a metatarsal area and a phalange area, wherein the heel area comprises a heel inner side area (HM) and a heel outer side area (HL), the arch area comprises an arch inner side area (MF 1) and an arch outer side area (MF 2), the metatarsal area comprises a first metatarsal area (M1), a second metatarsal area (M23), a third metatarsal area (M45) and a fourth metatarsal area (M45), and the phalange area comprises 9 areas in total of a first phalange area (T1) and a second-five phalange area (T25);
step S33: endowing the multi-cell structure sole in the step S32 with material properties, applying boundary conditions and dividing grid cells to obtain a finite element model of the multi-cell structure sole;
the step S4 specifically includes:
step S41: performing modal analysis on the finite element model of the multi-cell structure sole within a certain frequency range to obtain the inherent vibration characteristic of the multi-cell structure sole;
step S42: applying unit simple harmonic excitation force within a certain frequency range to an excitation surface of a finite element model of the multi-cell structure sole;
step S43: performing modal analysis based on the finite element model of the multi-cell structure sole in the step S41, and performing steady-state dynamics analysis on the multi-cell structure sole model in the step S42 to obtain the response condition of the multi-cell structure sole to vibration;
step S44: and taking the speed of each node on the surface of the multi-cell sole shoe as a history output variable, and outputting the speed of each node to an Abaqus result file and an odb file.
2. A multi-cell structure sole vibration energy transfer evaluation method according to claim 1, wherein: the step S5 specifically includes:
step S51: reading response speeds of all nodes on the excitation surface and the lower surface of the sole of the multi-cell structure from the odb file to a matrix A, and reading node excitation forces of all nodes on the excitation surface to a matrix B by using a Python programming language;
step S52: according to the response speed matrix A and the excitation force matrix B, combining an equivalent mechanical admittance calculation formula to obtain equivalent mechanical admittance data of each node of the excitation surface of the multi-cell structure sole and the lower surface of the sole in a full frequency range, and writing the equivalent mechanical admittance data into a csv table file;
step S53: calculating the average value of equivalent mechanical admittance data of each node of the excitation surface of the multi-cell structure sole and the lower surface of the sole in a full frequency range by using a Python programming language, and drawing a distribution cloud chart of the average value;
step S54: and (3) obtaining the average value of the equivalent mechanical admittances of the areas of the sole excitation surface and the lower surface of the sole of the multi-cell structure, and obtaining the data of the equivalent mechanical admittances of the areas of the sole excitation surface and the lower surface of the sole of the multi-cell structure.
3. A multi-cell structure sole vibration energy transfer evaluation method according to claim 1, wherein: the step S6 specifically includes:
step S61: equivalent input admittance data are obtained from equivalent mechanical admittances of the excitation surface of the multi-cell structure sole, and equivalent transfer admittance data are obtained from equivalent mechanical admittances of different areas of the lower surface of the sole;
step S62: and calculating the equivalent vibration transmissivities of different areas on the lower surface of the multi-cell structure sole according to the equivalent input admittance data and the equivalent transmission admittance data of the multi-cell structure sole and by combining an equivalent vibration transmissibility calculation formula.
4. A multi-cell structure sole vibration energy transfer evaluation method according to claim 1, wherein: the step S7 specifically includes:
step S71: repeatedly executing the steps S3-S6 on the multi-cell structure soles with different lattice types to obtain equivalent mechanical admittance distribution cloud charts of the multi-cell structure soles with different lattice types and equivalent mechanical admittance and equivalent vibration transmissibility data of different areas of the excitation surfaces of the multi-cell structure soles with different lattice types and the lower surfaces of the soles;
step S72: drawing a histogram of equivalent mechanical admittance and equivalent vibration transmissibility distribution of the soles with different types of multicellular structures in different areas of the soles according to the equivalent mechanical admittance and equivalent vibration transmissibility data of the excitation surfaces of the soles with different multicellular structures with different lattice types and different areas of the lower surfaces of the soles;
step S73: and according to the equivalent mechanical admittance and equivalent vibration transmissibility data of different areas of the lower surface of the multi-cell structure sole with different lattice types, calculating the average value of the equivalent mechanical admittance and the equivalent vibration transmissibility of the whole lower surface of the multi-cell structure sole with different lattice types, and drawing a histogram for comparison.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111145580.XA CN113836771B (en) | 2021-09-28 | 2021-09-28 | Multi-cell structure sole vibration energy transmission evaluation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111145580.XA CN113836771B (en) | 2021-09-28 | 2021-09-28 | Multi-cell structure sole vibration energy transmission evaluation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113836771A CN113836771A (en) | 2021-12-24 |
| CN113836771B true CN113836771B (en) | 2023-05-26 |
Family
ID=78967180
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111145580.XA Active CN113836771B (en) | 2021-09-28 | 2021-09-28 | Multi-cell structure sole vibration energy transmission evaluation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113836771B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114580041B (en) * | 2022-02-25 | 2024-08-16 | 华侨大学 | Vibration reduction sole structure optimization design method based on variable-size lattice filling |
| CN114722679A (en) * | 2022-04-28 | 2022-07-08 | 华侨大学 | Toe cap and toe cap optimization method based on multi-cell structure |
| CN114970167B (en) * | 2022-05-31 | 2024-07-26 | 华侨大学 | Design method and device for lattice variable density topology optimization structure of inner cavity of sole |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104033530B (en) * | 2014-04-24 | 2017-11-28 | 华北电力大学（保定） | A kind of method for improving and controlling Multi-halfwaves Filters energy absorption efficiency |
| CN208053276U (en) * | 2018-02-08 | 2018-11-06 | 长沙理工大学 | A kind of more born of the same parents' thin-wall tubes of functionally gradient |
| US20200170336A1 (en) * | 2018-11-30 | 2020-06-04 | Footwear Unlimited Inc. | Three layer shoe construction with improved cushioning, flexibility and traction |
| CN112084606B (en) * | 2020-09-29 | 2024-06-07 | 华侨大学 | Ball porous heel area filling structure sole and optimal design method thereof |
| CN113190917A (en) * | 2021-04-22 | 2021-07-30 | 武汉理工大学 | Multi-cell energy absorption structure combination optimization design method |
-
2021
- 2021-09-28 CN CN202111145580.XA patent/CN113836771B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN113836771A (en) | 2021-12-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113836771B (en) | Multi-cell structure sole vibration energy transmission evaluation method | |
| CN113657008A (en) | Shoe sole heel area vibration reduction structure optimization design method based on mechanical admittance | |
| CN102262699B (en) | Soft tissue deformation simulation method based on coupling of mesh-free Galerkin and mass spring | |
| CN111274702B (en) | Design method of sole model with buffer structure in sole heel area | |
| CN112199790B (en) | Regular polyhedron porous heel area filling structure sole and design method thereof | |
| CN113821958A (en) | Optimized design method for buffering multi-cellular sole structure | |
| CN113435082B (en) | Mechanical regulation and control simulation system for fracture surgery robot | |
| CN102955881A (en) | Method for calculating thermal fatigue failure probability of welding point of integrated circuit chip | |
| CN106021977B (en) | Subcutaneus adipose tissue biomethanics modeling method based on linear elasticity and hyperelastic model | |
| CN109388886A (en) | A kind of calculation method of porous material acoustic absorptivity | |
| Barforooshi et al. | Study neo-Hookean and Yeoh hyper-elastic models in dielectric elastomer-based micro-beam resonators | |
| CN112084606B (en) | Ball porous heel area filling structure sole and optimal design method thereof | |
| CN109101752B (en) | Method for calculating natural vibration frequency of local structure of complex hydraulic structure | |
| CN112199875A (en) | Component welding point random vibration fatigue life distribution prediction method based on rain flow method | |
| CN114330040B (en) | Sole vibration transfer characteristic analysis method | |
| Sarkar et al. | Performance of 2D-spectral finite element method in dynamic analysis of concrete gravity dams | |
| CN117195648A (en) | Wind-proof cultivation parameter analysis method and device based on dynamic growth monitoring | |
| CN112904604B (en) | Display panel support column performance evaluation method and device and display panel | |
| Wilke et al. | A quadratically convergent unstructured remeshing strategy for shape optimization | |
| CN213582170U (en) | Sole with regular polyhedron porous heel area filling structure | |
| CN114282403B (en) | Efficient high-precision habitat simulation method for coupled habitat suitable model | |
| Edwards et al. | Identification of tipping points in supply chain dynamics using effective dimension and resilience index | |
| Ahmad | High-resolution wave propagation method for stratified flows | |
| Turnea et al. | Numerical simulations of blunt-prosthesis-rod complex | |
| Zawra et al. | Stress Recovery in the Reduced Space for Parametric Reduced Models in Microelectronics |
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 | ||
| GR01 | Patent grant |