CN111274702A - Design method of sole model with buffer structure in sole heel area - Google Patents
Design method of sole model with buffer structure in sole heel area Download PDFInfo
- Publication number
- CN111274702A CN111274702A CN202010065537.1A CN202010065537A CN111274702A CN 111274702 A CN111274702 A CN 111274702A CN 202010065537 A CN202010065537 A CN 202010065537A CN 111274702 A CN111274702 A CN 111274702A
- Authority
- CN
- China
- Prior art keywords
- sole
- model
- heel
- area
- heel area
- 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
Links
Images
Landscapes
- Footwear And Its Accessory, Manufacturing Method And Apparatuses (AREA)
Abstract
The invention discloses a design method of a sole model with a buffer structure in a sole heel area, which comprises the following steps: step S1, establishing an optimized front sole model: step S2, establishing a finite element model of the sole model in the step S1 and carrying out static analysis; step S3, selecting an area in the heel area of the sole as the optimized area of the sole buffer structure, and respectively establishing four different porous structure models in the heel area to obtain four heel area buffer porous structure models respectively of M1, M2, M3 and M4. And step S4, setting boundary conditions for the four sole structures, performing power flow visualization programming in a phython language, and performing dynamic analysis correspondingly to obtain power flow visualization cloud pictures of the four sole structures. And step S5, comparing the power flow cloud charts of the four models, selecting one optimal model, and using the optimal model as an optimal heel buffer structure. The invention can ensure that the sole is in motion and improve the buffer performance.
Description
Technical Field
The invention relates to the field of optimization design, in particular to a design method of a sole model of a buffer structure in a sole heel area.
Background
With the rapid development of the footwear industry, the cushioning structure of the sole plays a crucial role in the movement, with the primary purpose of reducing impact to reduce foot injuries. For the evaluation of the shock absorption performance of the shoe sole, besides subjective evaluation, it is most common to perform kinematic measurements, such as comparing the shock absorption effect of different shoe soles by measuring the pressure distribution of the sole and the shoe sole during exercise.
In the traditional method, a real object is required to be made first, and then a tester wears the method to carry out test comparison, so that the method for evaluating the shock absorption performance of the sole needs a large number of real objects to carry out experimental analysis, and is high in cost and low in efficiency.
Disclosure of Invention
The invention mainly aims to overcome the defects in the prior art and provides a design method of a sole model of a buffer structure in a heel area of a sole.
The invention adopts the following technical scheme:
a design method of a sole model of a buffer structure in a sole heel area is characterized by comprising the following steps:
step S1, establishing a sole model before optimization;
step S2, establishing a finite element model of the sole model in the step S1 and carrying out static analysis;
s3, selecting an area in the heel area of the sole as an optimized area of the sole buffer structure, and respectively establishing four different porous structure models in the optimized area to obtain four heel area buffer porous structure models which are respectively M1, M2, M3 and M4;
step S4, setting boundary conditions of the four heel buffer porous structure models, performing phython language power flow visualization programming in abaqus software, and performing dynamic analysis correspondingly to obtain power flow visual cloud charts of four different porous structure models;
and step S5, comparing the power flow cloud charts of the four different porous structure models, selecting an optimal model from the four different porous structure models, and taking the optimal model as an optimal heel buffer structure.
Preferably, the step S2 specifically includes:
step S21: importing a sole model and a ground model into abaqus;
step S22: respectively endowing the sole model and the ground model with material properties, wherein the material properties comprise mass density, Young modulus and Poisson ratio;
step S23: establishing a contact pair of the outer surface of the sole and the ground surface;
step S24: applying a load and a boundary condition;
step S25: the integral model formed by the sole model and the ground model is divided into surface meshes, and the integral model is converted into a finite element model with moderate meshes and node quantity and high simulation degree.
Preferably, the step S3 specifically includes:
step S31: segmenting the heel region of the sole in UG;
step S32: performing circular Boolean reduction operation on the heel area to obtain an optimized model M1;
step S33: performing triangular Boolean reduction operation on the heel area to obtain an optimized model M2;
step S34: performing quadrilateral Boolean reduction operation on the heel area to obtain an optimized model M3;
step S35: and performing hexagonal Boolean reduction operation on the heel area to obtain an optimized model M4.
Preferably, the step S4 specifically includes:
step S41: the abaqus is introduced into four heel region buffer porous structure models M1, M2, M3 and M4;
step S42: carrying out meshing, material assignment and contact pair setting;
step S42: performing dynamic analysis on the heel landing condition of the model;
step S43: and writing a script by utilizing a python language to obtain a power flow cloud picture.
Preferably, the step S5 specifically includes:
step S51: comparing and analyzing a stress cloud picture and a dynamic picture of the whole process under the working condition that the heel of the sole touches the ground;
step S52: comparing and analyzing the power flow cloud pictures;
step S53: and providing an optimal optimized sole structure based on the analysis result.
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, 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, various geometric structures are provided, and the energy aspect of the sole is analyzed by using a power flow method, so that the advantages and the disadvantages of the soles with different structures are effectively analyzed.
3) The method can guide the design and production of the sports shoes, and saves a large amount of cost compared with the traditional method of making a real object and then wearing the sports shoes on a human body for testing.
4) The invention uses numerical model and finite element method to analyze, can directly test the damping performance of the sole in the model, arbitrarily change the structure of the heel area of the sole, and has high stress condition and efficiency of the sole of the model.
Drawings
FIG. 1 is a schematic flow chart of the main steps of the method of the present invention;
FIG. 2 is a view of a sole model;
FIG. 3 is an assembly view of the sole and the ground;
FIG. 4 is a four sole optimized configuration;
FIG. 5 is a model of the heel of the sole when it lands;
FIG. 6(a) is a cloud view of sole displacement;
FIG. 6(b) is a sole stress cloud;
fig. 7 is a power flow vector diagram of the sole.
The invention is described in further detail below with reference to the figures and specific examples.
Detailed Description
The invention is further described below by means of specific embodiments.
A method for designing a sole model of a cushioning structure in a heel region of a sole, as shown in fig. 1 and 2, comprising:
step S1: and establishing a sole model before optimization in UG.
Step S2: establishing a finite element model of the sole model in the step S1 and performing static analysis, as shown in fig. 3, specifically including:
step S21: the sole model and the bottom surface model are imported into abaqus. abaqus is a set of powerful finite element software for engineering simulation, and comprises a rich unit library capable of simulating any geometric shape and various material model libraries capable of simulating the performance of typical engineering materials.
Step S22: and respectively endowing the sole model and the ground model with material properties, wherein the material properties comprise mass density, Young modulus and Poisson ratio, so that the model and the real object are made of the same material, and the aim of simulation is fulfilled.
Specifically, the Young model of the sole is set to be 4MPa, and the Poisson ratio is set to be 0.4; the Young's modulus of the ground was set at 17000MPa, and the Poisson's ratio was 0.3.
Step S23: establishing a contact pair of the outer surface of the sole and the ground surface;
specifically, a self-contact algorithm was used, and the coefficient of friction was set to 0.6.
Step S24: applying a load and a boundary condition;
step S25: the integral model formed by the sole model and the bottom model is divided into surface meshes, so that the integral model of three-dimensional simulation is converted into a finite element entity unit model with moderate mesh and node quantity and high simulation degree.
Step S3: selecting an area in a heel area of a sole as an optimized area of a sole buffering structure, respectively establishing four different porous structure models in the heel area, and obtaining four heel area buffering porous structure models which are respectively M1, M2, M3 and M4, as shown in FIG. 4, specifically comprising:
step S31: segmenting the heel region of the sole in UG;
step S32: performing circular Boolean reduction operation on the heel area to obtain an optimized model M1;
step S33: performing triangular Boolean reduction operation on the heel area to obtain an optimized model M2;
step S34: performing quadrilateral Boolean reduction operation on the heel area to obtain an optimized model M3;
step S35: and performing hexagonal Boolean reduction operation on the heel area to obtain an optimized model M4.
Step S4: setting boundary conditions for the four types, performing phython language for power flow visualization programming in abaqus software, and performing dynamic analysis correspondingly to obtain power flow visualization cloud charts of the four sole structures, specifically comprising:
step S41: introducing abaqus into models M1, M2, M3 and M4, and adjusting the position of the sole to simulate the heel-strike condition as shown in FIG. 5;
step S42: carrying out meshing, material assignment and contact pair setting;
specifically, a contact pair of the outer surface of the sole and the upper surface of the ground is established, wherein the outer surface of the sole is a first surface, the upper surface of the ground is a second surface, the friction coefficient is 0.6, the sole is set to initial speeds of 1650mm/s in the x direction and 230mm/s in the z direction, and 300N pressure is applied to the upper surface of the heel of the sole.
Step S42: performing dynamic analysis on the heel landing condition of the model;
step S43: and deducing and calculating the power flow, and writing a script by utilizing a python language to obtain a power flow cloud diagram.
Step S5: comparing the power flow cloud charts of the four models, selecting one optimal model, and using the optimal model as an optimal heel buffer structure.
Step S51: comparing and analyzing the stress cloud pictures and dynamic pictures of the whole process under the working condition that the heel of the sole touches the ground, obtaining the stress condition of the sole, the maximum stress and the maximum displacement of the sole, and preliminarily judging the quality of the sole structure by comparing the maximum stress and the maximum displacement, as shown in fig. 6a and 6 b;
step S52: comparing and analyzing the power flow cloud pictures, as shown in fig. 7, the power flow cloud pictures are vector pictures, the size of arrows in the pictures can reflect the size of power flow, the direction of the arrows reflects the direction of the power flow, the outward transmission condition of sole energy at a stressed position can be obtained, and a better sole structure can be found by comparing the distribution condition of the power flow;
step S53: and an optimal sole optimization structure is provided based on the analysis result, the stress cloud picture and the power flow cloud picture are comprehensively compared, and the sole structure with smaller maximum stress and smaller power flow arrow is the optimal structure.
The above description is only an embodiment of the present invention, but the design concept of the present invention is not limited thereto, and any insubstantial modifications made by using the design concept should fall within the scope of infringing the present invention.
Claims (5)
1. A design method of a sole model of a buffer structure in a sole heel area is characterized by comprising the following steps:
step S1, establishing a sole model before optimization;
step S2, establishing a finite element model of the sole model in the step S1 and carrying out static analysis;
s3, selecting an area in the heel area of the sole as an optimized area of the sole buffer structure, and respectively establishing four different porous structure models in the optimized area to obtain four heel area buffer porous structure models which are respectively M1, M2, M3 and M4;
step S4, setting boundary conditions of the four heel buffer porous structure models, performing phython language power flow visualization programming in abaqus software, and performing dynamic analysis correspondingly to obtain power flow visual cloud charts of four different porous structure models;
and step S5, comparing the power flow cloud charts of the four different porous structure models, selecting one optimal model, and using the optimal model as an optimal heel buffer structure.
2. The method according to claim 1, wherein said step S2 specifically includes:
step S21: importing a sole model and a ground model into abaqus;
step S22: respectively endowing the sole model and the ground model with material properties, wherein the material properties comprise mass density, Young modulus and Poisson ratio;
step S23: establishing a contact pair of the outer surface of the sole and the ground surface;
step S24: applying a load and a boundary condition;
step S25: the integral model formed by the sole model and the ground model is divided into surface meshes, and the integral model is converted into a finite element model with moderate meshes and node quantity and high simulation degree.
3. The method according to claim 1, wherein said step S3 specifically includes:
step S31: segmenting the heel region of the sole in UG;
step S32: performing circular Boolean reduction operation on the heel area to obtain an optimized model M1;
step S33: performing triangular Boolean reduction operation on the heel area to obtain an optimized model M2;
step S34: performing quadrilateral Boolean reduction operation on the heel area to obtain an optimized model M3;
step S35: and performing hexagonal Boolean reduction operation on the heel area to obtain an optimized model M4.
4. The method according to claim 1, wherein said step S4 specifically includes:
step S41: the abaqus is introduced into four heel region buffer porous structure models M1, M2, M3 and M4;
step S42: carrying out meshing, material assignment and contact pair setting;
step S42: performing dynamic analysis on the heel landing condition of the model;
step S43: and writing a script by utilizing a python language to obtain a power flow cloud picture.
5. The method according to claim 1, wherein said step S5 specifically includes:
step S51: comparing and analyzing a stress cloud picture and a dynamic picture of the whole process under the working condition that the heel of the sole touches the ground;
step S52: comparing and analyzing the power flow cloud pictures;
step S53: and providing an optimal optimized sole structure based on the analysis result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010065537.1A CN111274702B (en) | 2020-01-20 | 2020-01-20 | Design method of sole model with buffer structure in sole heel area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010065537.1A CN111274702B (en) | 2020-01-20 | 2020-01-20 | Design method of sole model with buffer structure in sole heel area |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111274702A true CN111274702A (en) | 2020-06-12 |
| CN111274702B CN111274702B (en) | 2022-06-07 |
Family
ID=71002238
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010065537.1A Active CN111274702B (en) | 2020-01-20 | 2020-01-20 | Design method of sole model with buffer structure in sole heel area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111274702B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 | 华侨大学 | Sole with regular polyhedron porous heel area filling structure and design method thereof |
| CN112206553A (en) * | 2020-09-24 | 2021-01-12 | 天津大学 | Method for evaluating competitive adsorption of polystyrene microspheres and humic acid in sand filtration |
| CN113558347A (en) * | 2021-06-16 | 2021-10-29 | 福建师范大学 | Design method for middle structure of outsole of endurance running shoe |
| CN113806988A (en) * | 2021-09-28 | 2021-12-17 | 华侨大学 | Visual analysis method for vibration energy transfer of sole with lattice structure |
| CN114330040A (en) * | 2021-09-28 | 2022-04-12 | 华侨大学 | Sole vibration transmission characteristic analysis method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101327111B1 (en) * | 2012-05-31 | 2013-11-08 | 서울시립대학교 산학협력단 | Topology optimization method of linking finite element analysis program and numerical analysis program and design method for ultrasonic transducer using the same |
| CN105956349A (en) * | 2016-07-04 | 2016-09-21 | 上海理工大学 | Insole production system based on finite element mechanical analysis |
| CN106372374A (en) * | 2016-11-03 | 2017-02-01 | 西安交通大学 | Individualized design method for diabetic foot shoe pad |
| CN109662392A (en) * | 2018-12-20 | 2019-04-23 | 上海工程技术大学 | A kind of test method of foot and sole finite element modeling based on landing movement |
-
2020
- 2020-01-20 CN CN202010065537.1A patent/CN111274702B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101327111B1 (en) * | 2012-05-31 | 2013-11-08 | 서울시립대학교 산학협력단 | Topology optimization method of linking finite element analysis program and numerical analysis program and design method for ultrasonic transducer using the same |
| CN105956349A (en) * | 2016-07-04 | 2016-09-21 | 上海理工大学 | Insole production system based on finite element mechanical analysis |
| CN106372374A (en) * | 2016-11-03 | 2017-02-01 | 西安交通大学 | Individualized design method for diabetic foot shoe pad |
| CN109662392A (en) * | 2018-12-20 | 2019-04-23 | 上海工程技术大学 | A kind of test method of foot and sole finite element modeling based on landing movement |
Non-Patent Citations (1)
| Title |
|---|
| 章浩伟 等: "基于三维膝-踝-足有限元模型的足跟痛足底压力生物力学分析", 《医用生物力学》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112206553A (en) * | 2020-09-24 | 2021-01-12 | 天津大学 | Method for evaluating competitive adsorption of polystyrene microspheres and humic acid in sand filtration |
| 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 | 华侨大学 | Sole with regular polyhedron porous heel area filling structure and design method thereof |
| CN112084606B (en) * | 2020-09-29 | 2024-06-07 | 华侨大学 | Ball porous heel area filling structure sole and optimal design method thereof |
| CN112199790B (en) * | 2020-09-29 | 2024-06-07 | 华侨大学 | Regular polyhedron porous heel area filling structure sole and design method thereof |
| CN113558347A (en) * | 2021-06-16 | 2021-10-29 | 福建师范大学 | Design method for middle structure of outsole of endurance running shoe |
| CN113558347B (en) * | 2021-06-16 | 2022-08-16 | 福建师范大学 | Design method for middle structure of outsole of endurance running shoe |
| CN113806988A (en) * | 2021-09-28 | 2021-12-17 | 华侨大学 | Visual analysis method for vibration energy transfer of sole with lattice structure |
| CN114330040A (en) * | 2021-09-28 | 2022-04-12 | 华侨大学 | Sole vibration transmission characteristic analysis method |
| CN114330040B (en) * | 2021-09-28 | 2024-07-26 | 华侨大学 | Sole vibration transfer characteristic analysis method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111274702B (en) | 2022-06-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111274702B (en) | Design method of sole model with buffer structure in sole heel area | |
| US7415398B2 (en) | Method of simulating viscoelastic material | |
| Asadi et al. | A bonded particle model simulation of shear strength and asperity degradation for rough rock fractures | |
| CN110473597B (en) | Conglomerate mechanical property evaluation and analysis method and system | |
| JP3958666B2 (en) | Method for calculating energy loss in viscoelastic material, and method for evaluating energy loss of golf ball using the method | |
| CN108959686A (en) | A kind of correction method for finite element model based on sensitivity analysis | |
| CN110398385A (en) | A kind of evaluation method for Evaluations for Ride Comfort of Vehicle Seats | |
| CN109662392A (en) | A kind of test method of foot and sole finite element modeling based on landing movement | |
| Zeller | Cloth simulation on the gpu | |
| Tang et al. | A multilayer nodes update method in FEM simulation of large depth fretting wear | |
| CN109086542A (en) | A kind of simulation optimization method, system and the device of machine box for server pedestal convex closure | |
| Bakroon et al. | Arbitrary lagrangian-eulerian finite element formulations applied to geotechnical problems | |
| CN115455754A (en) | Mine hydraulic support design method based on digital twinning | |
| CN109766663A (en) | A kind of Programs of Large Pressurized Vessel overlong welded seam welding residual stress and deform efficient computation processing method | |
| CN104166747B (en) | A kind of bolted joint Dynamic Characteristics Analysis Method for considering suface processing quality | |
| CN113033044A (en) | Modeling method for predicting protective effect of ankle protector | |
| JP5432549B2 (en) | Rubber material simulation method | |
| Li et al. | A rapid prototyping method for crop models using the discrete element method | |
| CN115114823A (en) | Part finite element analysis method, equipment and medium based on actual boundary constraint | |
| CN108256246B (en) | Design method and device of porous material based on Unity3D | |
| Wijnhoven et al. | Testing archaeological mail armour in a virtual environment: 3rd century BC to 10th century AD | |
| Volino et al. | From measured physical parameters to the haptic feeling of fabric | |
| CN104156512B (en) | Method for setting up finite element model simulating corn ear power contact | |
| US8229711B2 (en) | Method of analysis for kinetic properties of golf club head and golf club therefor | |
| JP2021133814A (en) | Tire simulation method, program and tire simulation device |
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 |