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.
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.