CN113284569A - Method and device for optimizing plants in orthopedics department, electronic equipment and storage medium - Google Patents
Method and device for optimizing plants in orthopedics department, electronic equipment and storage medium Download PDFInfo
- Publication number
- CN113284569A CN113284569A CN202110406442.6A CN202110406442A CN113284569A CN 113284569 A CN113284569 A CN 113284569A CN 202110406442 A CN202110406442 A CN 202110406442A CN 113284569 A CN113284569 A CN 113284569A
- Authority
- CN
- China
- Prior art keywords
- dimensional model
- femur
- model
- implant
- femoral
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 230000000399 orthopedic effect Effects 0.000 title claims abstract description 20
- 238000003860 storage Methods 0.000 title claims abstract description 12
- 210000000689 upper leg Anatomy 0.000 claims abstract description 117
- 239000007943 implant Substances 0.000 claims abstract description 96
- 238000005457 optimization Methods 0.000 claims abstract description 54
- 238000012545 processing Methods 0.000 claims abstract description 22
- 239000000463 material Substances 0.000 claims abstract description 14
- 230000004044 response Effects 0.000 claims description 28
- 238000013461 design Methods 0.000 claims description 20
- 210000000988 bone and bone Anatomy 0.000 claims description 19
- 239000002639 bone cement Substances 0.000 claims description 15
- 238000012360 testing method Methods 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 9
- 230000001054 cortical effect Effects 0.000 claims description 7
- 230000002068 genetic effect Effects 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 238000009826 distribution Methods 0.000 claims description 5
- 238000009499 grossing Methods 0.000 claims description 4
- 238000013401 experimental design Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- 238000005192 partition Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 13
- 238000002513 implantation Methods 0.000 abstract description 5
- 230000035882 stress Effects 0.000 description 27
- 229910001069 Ti alloy Inorganic materials 0.000 description 22
- 238000013400 design of experiment Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 235000001674 Agaricus brunnescens Nutrition 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 208000001132 Osteoporosis Diseases 0.000 description 2
- WAIPAZQMEIHHTJ-UHFFFAOYSA-N [Cr].[Co] Chemical class [Cr].[Co] WAIPAZQMEIHHTJ-UHFFFAOYSA-N 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 208000020089 femoral neck fracture Diseases 0.000 description 2
- 210000004394 hip joint Anatomy 0.000 description 2
- 230000002980 postoperative effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 239000010935 stainless steel Substances 0.000 description 2
- 229910001220 stainless steel Inorganic materials 0.000 description 2
- 229910000883 Ti6Al4V Inorganic materials 0.000 description 1
- RTAQQCXQSZGOHL-UHFFFAOYSA-N Titanium Chemical compound [Ti] RTAQQCXQSZGOHL-UHFFFAOYSA-N 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 239000000316 bone substitute Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 210000002414 leg Anatomy 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000002271 resection Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 210000001694 thigh bone Anatomy 0.000 description 1
- 239000010936 titanium Substances 0.000 description 1
- 229910052719 titanium Inorganic materials 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- 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
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Prostheses (AREA)
Abstract
The invention provides a method, a device, electronic equipment and a storage medium for optimizing plants in orthopedics, wherein the method comprises the following steps: acquiring a human femur three-dimensional model; establishing a corresponding three-dimensional model of the internal implant according to the three-dimensional model of the human femur; embedding the human femur three-dimensional model and the implant three-dimensional model to obtain a complete femur implant three-dimensional model; parameter setting and optimization processing are carried out on the three-dimensional model of the femoral internal implant, the optimal solution of each target section of the three-dimensional model of the femoral internal implant is obtained, and the optimal solution is verified; and optimizing the size of the porous structure of each target segment according to the optimal solution to obtain the porous structure conforming to the optimal solution. According to the invention, the porous structure with material performance gradient is obtained by optimizing the endophyte, so that the femur stress shielding effect can be weakened as much as possible while the strength is ensured, and the success rate of the endophyte implantation is improved.
Description
Technical Field
The invention relates to the technical field of medical instruments, in particular to a method and a device for optimizing plants in orthopedics department, electronic equipment and a storage medium.
Background
With the gradual development of our society towards aging, the incidence of femoral neck fracture diseases of the aged gradually increases year by year. According to the relevant clinical report, the artificial hip joint replacement is applied to the clinical treatment of femoral neck fracture diseases of the old, so that the satisfactory clinical curative effect can be obtained, the occurrence of adverse events can be reduced, and the treatment safety can be guaranteed.
At present, the application in the field of artificial hip joints is important for the biomechanical compatibility of an implant material, the biomechanical compatibility has important influence on the effectiveness and reliability of the material implanted into a human body, the excellent biomechanical compatibility can ensure the friendly fusion of the implant and human bones, and the generation of rejection is avoided, wherein the elastic modulus is one of the most important physical properties for judging the biocompatibility of biomedical metal materials, particularly bone substitute materials. In the prior art, a more widely applied and mature endoprosthesis material is a titanium alloy material, although the elastic modulus of the titanium alloy is much lower than that of cobalt-chromium alloy and stainless steel, the elastic modulus of the titanium alloy is about ten times that of human bones, and too high metal elastic modulus can cause the interface between an endoprosthesis and human bones to loosen, affect the functions of an implantation device, or cause a stress shielding effect to cause the function degradation or absorption of bone tissues, and easily cause secondary injuries such as postoperative osteoporosis; and the elastic modulus of the inner plant is too low, so that the inner plant is easy to deform greatly under the action of stress and cannot play a role in fixing and supporting.
Disclosure of Invention
The invention provides a method and a device for optimizing orthopedic plants, electronic equipment and a storage medium, which are used for solving the technical problem that secondary damage such as postoperative osteoporosis and the like is easily caused by stress shielding generated by an elastic modulus of an internal implant higher than that of human bones so as to achieve the purposes of weakening the stress shielding effect and improving the success rate of implantation of the internal implant.
In a first aspect, the present invention provides a method for optimizing plants in the orthopedics department, comprising:
acquiring a human femur three-dimensional model;
establishing a corresponding three-dimensional model of the internal implant according to the three-dimensional model of the human femur;
embedding the human femur three-dimensional model and the implant three-dimensional model to obtain a complete femur implant three-dimensional model;
parameter setting and optimization processing are carried out on the three-dimensional model of the femoral internal implant, the optimal solution of the equivalent elastic modulus of each target section of the three-dimensional model of the femoral internal implant is obtained, and the optimal solution is verified;
and optimizing the size of the initial porous structure of each target segment according to the optimal solution to obtain a porous structure with elastic modulus gradient which accords with the optimal solution.
In a second aspect, the present invention provides an apparatus for optimizing plants in orthopedics, comprising:
the acquisition module is used for acquiring a human femur three-dimensional model;
the establishing module is used for establishing a corresponding three-dimensional model of the endoprosthesis according to the three-dimensional model of the human femur;
the embedding module is used for embedding the human femur three-dimensional model and the endoprosthesis three-dimensional model to obtain a complete femur endoprosthesis three-dimensional model;
the processing module is used for carrying out parameter setting and optimization processing on the three-dimensional model of the femoral internal implant, obtaining an optimal solution of the equivalent elastic modulus of each target section of the three-dimensional model of the femoral internal implant and verifying the optimal solution;
and the optimization module is used for optimizing the size of the initial porous structure of each target segment according to the optimal solution to obtain a porous structure with elastic modulus gradient which accords with the optimal solution.
In a third aspect, the present invention provides an electronic device comprising: a processor, a memory, and a bus, wherein,
the processor and the memory are communicated with each other through the bus;
the memory stores program instructions executable by the processor, the processor invoking the program instructions to perform a method as described in any of the above.
In a fourth aspect, the invention provides a non-transitory computer readable storage medium storing computer instructions that cause the computer to perform the method of any one of the above.
The invention provides a method, a device, electronic equipment and a storage medium for optimizing plants in orthopedics department, wherein the method comprises the steps of obtaining a human femur three-dimensional model, then establishing a corresponding internal implant three-dimensional model according to the human femur three-dimensional model, embedding the human femur three-dimensional model and the internal implant three-dimensional model to obtain a complete femur internal plant three-dimensional model, setting parameters and optimizing the complete femur internal plant three-dimensional model to obtain an optimal solution of equivalent elastic modulus of each target section, and optimizing the size of an initial porous structure of each target section according to the obtained optimal solution to obtain a porous structure meeting requirements. According to the invention, the porous structure of each target section of the femur implant three-dimensional model is optimized, so that the femur stress shielding effect is weakened while the strength of the implant is ensured, and the success rate of the implant of the femur is improved.
Drawings
In order to more clearly illustrate the technical solutions of the present invention or the prior art, the drawings needed for the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings according to the drawings without creative efforts.
FIG. 1 is a schematic diagram of the overall flow of the method for optimizing orthopedic endophytes provided by the present invention;
FIG. 2 is a schematic structural diagram of an orthopedic implant segment according to the present invention;
FIG. 3 is a schematic flow chart of the method for optimizing orthopedic endophytes provided by the present invention;
FIG. 4 is a schematic structural diagram of a three-dimensional model of a human femur for use in an orthopaedic implant optimization method according to the present invention;
FIG. 5 is a schematic structural diagram of an orthopedic implant optimization device provided by the present invention;
fig. 6 is a schematic structural diagram of an electronic device provided in the present invention.
Description of reference numerals:
1-femoral cortical bone; 2-femoral cancellous bone (bone cement);
3-titanium alloy inner plant.
Detailed Description
In order to more clearly illustrate the technical solutions of the present invention or the prior art, the drawings needed for the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and those skilled in the art can also obtain other drawings according to the drawings without creative efforts.
Fig. 1 is a general flow chart of the orthopedic endophyte optimization method provided by the invention. As shown in fig. 1, the method for optimizing orthopedic endophytes provided by the invention comprises the following steps:
step 101: acquiring a human femur three-dimensional model;
step 102: establishing a corresponding three-dimensional model of the internal implant according to the three-dimensional model of the human femur;
step 103: embedding the human femur three-dimensional model and the implant three-dimensional model to obtain a complete femur implant three-dimensional model;
step 104: parameter setting and optimization processing are carried out on the three-dimensional model of the femoral internal implant, the optimal solution of the equivalent elastic modulus of each target section of the three-dimensional model of the femoral internal implant is obtained, and the optimal solution is verified;
step 105: and optimizing the size of the initial porous structure of each target segment according to the optimal solution to obtain a porous structure with elastic modulus gradient which accords with the optimal solution.
Specifically, the endoprosthesis refers to a device for implanting into human femur, the material of the endoprosthesis can be stainless steel, cobalt-chromium alloy, titanium or titanium alloy, etc., in this embodiment, the preferable endoprosthesis is titanium alloy TC4 (Ti-6 Al-4V). Wherein each target segment refers to each portion of the model that needs to be processed.
In the implementation, a human femur three-dimensional model is obtained firstly, an internal plant three-dimensional model is established according to the human femur three-dimensional model in SolidWorks software, the two models are embedded and combined to obtain a complete femur internal plant three-dimensional model, and the model is exported to be an X _ T file; then, the complete three-dimensional model of the femoral internal implant is led into finite element analysis software WorkBench, parameter setting and optimization processing are carried out on the three-dimensional model of the femoral internal implant, the optimal solution of each target section of the three-dimensional model of the femoral internal implant is obtained, and the optimal solution of the equivalent elastic modulus of each target section is verified; and finally, optimizing the size of the initial porous structure of each target segment according to the optimal solution to obtain the porous structure with elastic modulus gradient which accords with the optimal solution.
In the embodiment of the invention, a three-dimensional model of a human femur is obtained, a corresponding three-dimensional model of an internal implant is established, the three-dimensional model of the internal implant is embedded into the three-dimensional model of the human femur to obtain a complete three-dimensional model of an internal plant of the femur, the three-dimensional model of the internal plant of the femur is subjected to parameter design and optimization treatment to obtain an optimal solution of each target section, and the size of a porous structure of each target section is optimized according to the optimal solution to obtain a porous structure meeting requirements. According to the invention, the femur stress shielding effect can be weakened while the internal plant strength is ensured by optimizing the three-dimensional model and the porous structure of the femur implant, and the success rate of the implantation of the femur implant is improved.
In another embodiment of the present invention, before embedding the three-dimensional model of human femur and the three-dimensional model of endoprosthesis and obtaining the complete three-dimensional model of femoral endoprosthesis, the embedding process includes:
performing segmentation processing on the three-dimensional model of the endoprosthesis;
wherein the endoprosthesis is cut into 8 segments.
Specifically, the three-dimensional model of the internal implant is a model obtained by three-dimensional software design according to a femur three-dimensional model conforming to an actual human body, and is more realistic.
In this embodiment, as shown in fig. 2, the preferable endoprosthesis is a titanium alloy endoprosthesis, the titanium alloy endoprosthesis model is divided into a-H eight segments by a segmented design, the divided endoprosthesis model is implanted into a cavity of a three-dimensional model of a human femur, and the distal femur is resected to form a smooth bottom surface for constraint.
In the embodiment of the invention, the three-dimensional model of the internal plant is processed in a segmented manner, so that the independent processing of the plants in each segment is realized, the stress shielding effect can be reduced, the gradient characteristic is adopted to ensure that the internal plant still has higher strength at a low sensitivity place, and the situation of insufficient strength caused by the integral specified setting of the elastic modulus is avoided.
In an embodiment of the present invention, as shown in fig. 3, the performing parameter setting on the three-dimensional model of the plant in femur comprises:
and carrying out material setting, grid division, contact setting, constraint setting, load loading and solving setting on the three-dimensional model of the plant in the femur to obtain initial stress distribution of the three-dimensional model of the plant in the femur.
Specifically, the material of the endophyte is titanium alloy; the parameter setting is operated in the finite element analysis software WorkBench, and may also be operated in other analysis software, which is not specifically limited herein.
In the embodiment, the material properties of each part of the three-dimensional model of the femoral internal implant are set, such as the outer cortex bone of the femur, the elastic modulus is set to be 13.7GPa, and the Poisson ratio is 0.33; the elastic modulus of the bone cement is set to be 2.2GPa, and the Poisson ratio is 0.3; the material of the plants in the titanium alloy is Ti-6Al-4V, the elastic modulus is set to be 114GPa, and the Poisson ratio is 0.36.
In this embodiment, mesh division of a three-dimensional model of the femoral implant is set, and a preferred global size in this embodiment is 2mm, wherein, the titanium alloy implant adopts a MultiZone and Hex Dominant method to divide the mesh, and the mesh quality is controlled by using the surface size in a narrow part; the bone cement adopts Hex Dominant and tetrahedron methods to divide grids, and the grid quality of a narrow area is adjusted by using the surface size and the edge size; the femoral cortical bone is divided into grids by adopting a Hex Dominant method, and other grids are set globally.
In the embodiment, the contact between each part in the three-dimensional model of the femoral implant is set, wherein the three-dimensional model of the titanium alloy implant needs to be set as an integral part, and each target segment is mutually influenced; the plants in the titanium alloy and the bone cement, the plants in the titanium alloy and the femoral cortical bone, and the bone cement and the femoral cortical bone are all set to be in binding contact. And setting constraint conditions of the distal end of the plant three-dimensional model in the femur, and setting a plane for distal end resection as fixed constraint.
In this embodiment, the load loading of the three-dimensional model of the femoral implant is set, considering that the weight of a male in asia is 70KG, and the downward pressure on the femur is 1400N when the male stands on one leg, and the pressure acts on the top of the titanium alloy plant instead of the mushroom head of the femur and is directed downward along the femur.
In this embodiment, the stress solution of the three-dimensional model of the femoral implant is set, and the set solution form is Von-Mises equivalent stress, wherein the equivalent stress of 7 Gruen regions on the outer surface of the femur is particularly required to be solved, so as to obtain the initial stress distribution condition of the titanium alloy implant before optimization.
In the embodiment of the invention, the initial stress distribution of the three-dimensional model of the femur implant is obtained by setting parameters of the three-dimensional model of the femur implant, which mainly relate to the parameter settings of six aspects of material setting, grid division, contact setting, constraint setting, load loading and solving setting. The obtained initial stress distribution of the plant three-dimensional model in the femur can be used for better optimizing the plant three-dimensional model in the femur.
In another embodiment of the present invention, as shown in fig. 3, the optimizing the plant three-dimensional model in femur includes:
acquiring initial stress of seven gurn partitions on the outer surface of the femur in the three-dimensional femoral plant model and equivalent elastic modulus corresponding to seven internal plant target sections;
carrying out parametric modeling on the obtained initial stress and the equivalent elastic modulus to obtain a parametric model;
importing the parameterized model into a response surface optimization module, and carrying out test design on the parameterized model to obtain a response surface model;
and carrying out optimization design on the response surface model based on a multi-objective genetic optimization algorithm to obtain the optimal solution of the equivalent elastic modulus of each target segment.
And determining the response surface model by acquiring the sample points and the related function of the experimental design.
In particular, DOE (DESIGN OF EXPERIMENT) plays a very important role in the whole process OF quality control, and is an important guarantee for improving product quality and process flow.
In the embodiment, the initial stress and the equivalent elastic modulus of each target segment of the three-dimensional model of the femur implant are obtained, then the obtained initial stress and the obtained equivalent elastic modulus are subjected to parameterization to generate a parameterized model, the parameterized model is led into a response surface module, and DOE (design of object) test is performed on the parameterized model to obtain a response surface model; and optimizing the response surface model based on a multi-objective genetic optimization algorithm to obtain the optimal solution of each target segment.
Firstly, parameterizing parts for Optimization, including equivalent stress values of 7 Gruen regions and equivalent elastic modulus in material properties of the later 7 sections of plants in the titanium alloy, and importing parameterized data into a Response Surface Optimization module Response Surface Optimization for further operation.
In this embodiment, a DOE test design is performed on an input parameterized model, a response surface model is obtained according to the DOE test design, and optimization processing is performed on the response surface model based on a multi-objective genetic optimization algorithm.
In the response surface model optimization setting, DOE test design is firstly needed to be used for sampling, a Latin Hypercube (Latin Hypercube) test design method is adopted, the method is suitable for the condition of influencing multiple factors, the test points are uniform, the test times are equal to the horizontal number, the test times can use any numerical value and are uniformly covered, the scale of the test can be effectively reduced, and the specific implementation steps are as follows:
step 1: dividing each dimension into m intervals which are not overlapped with each other, and enabling each interval to have the same probability;
step 2: randomly drawing one point in each interval in each dimension to form the whole Latin hypercube.
And step 3: randomly extracting points in each dimension to form a vector, and once each layer of samples is extracted, the layer is not extracted.
And calculating a response surface model based on the response surface type of the Kriging model according to the response surface design point obtained by the experimental sampling result. Based on the Kriging method, the output parameter is equal to the global design space plus the local deviation, and the expression is as follows:
y(x)＝F(β,x)+z(x)＝fT(x)β+z(x)
where β is a regression coefficient of the basis function, f (x) is a polynomial function of variable x representing the global model of the design space, and z (x) is a mean of 0 and a variance of σ2Gaussian random function of (2).
The design points obtained based on the Latin hypercube experiment design method form a correlation matrix:
the method comprises the following steps of obtaining a Kriging response surface model, calculating a fitting result of the response surface model, and determining a target Kriging response surface model according to the fitting result, wherein n is the total number of data points, the final Kriging response surface model is determined according to the sample points and a related function of experimental design, and the response surface model has a good effect on solving the nonlinear engineering optimization problem and has the characteristics of high calculation efficiency, short time, good response surface fitting effect, high result accuracy and the like.
After a response surface model is established, Optimization design of equivalent elastic modulus is carried out based on Optimization, the used genetic algorithm is a multi-objective genetic algorithm (MOGA), the maximum stress value of 7 Gruen areas on the outer surface of the femur is taken as target stress, the equivalent elastic modulus of plants in 7 sections of titanium alloy is taken as an input variable, iteration times are set, and final target elastic modulus Optimization is carried out to obtain the optimal reference point equivalent elastic modulus.
And importing the obtained optimal solution into the three-dimensional femoral implant model after the parameter setting is finished for verification, importing the obtained optimal elastic modulus into an initial titanium alloy implant model for stress analysis, resetting the equivalent elastic modulus of the lower end 7 section of the titanium alloy implant, keeping the other conditions unchanged, and performing stress analysis again to verify the effect of the optimization model.
In the embodiment of the invention, the initial stress and the equivalent elastic model of the model are obtained, the parameterized model is obtained, then the parameterized model is led into the response surface module, DOE test design is carried out, the response surface model is obtained, and the optimal solution of each target segment is obtained based on a multi-target genetic optimization algorithm. According to the method, the optimal solution of the equivalent elastic modulus of each target section is obtained through optimization processing of the model, so that the elastic modulus of the endophyte can be reduced, and the stress shielding effect of the femur is weakened.
In another embodiment of the present invention, the optimizing the size of the porous structure of each target segment according to the optimal solution comprises:
acquiring an initial porous structure of a plant in each target section in the three-dimensional model of the plant in the femur;
parameterizing the size and equivalent elastic modulus of the initial porous structure to be adjusted in the plant within each target segment;
and optimizing the size of the initial porous structure after the parameterization treatment based on a sequence quadratic programming method to obtain the porous structure with the elastic modulus gradient which accords with the optimal solution.
In particular, the sequential quadratic programming algorithm is one of the most effective methods for solving the constrained nonlinear optimization problem, and the specific implementation manner is not specifically stated herein.
In this embodiment, the initial porous structure of each target segment is subjected to size optimization according to the obtained optimal solution, and it should be noted that the size and the equivalent elastic modulus of the initial porous structure that needs to be adjusted in the plant in each target segment need to be parameterized. If a 4mm cube cell is taken as a unit, a porous structure is designed in the cell, and the porous structure is characterized in that the porosity of the whole cube is influenced by changing one dimension, and the equivalent elastic modulus of plants in the titanium alloy is reduced by changing the porosity.
The size and the equivalent elastic modulus of the internal cell element porous structure to be adjusted of each section of titanium alloy internal implant are parameterized, and WorkBench is also used for carrying out size optimization processing on the porous structure so as to achieve the aim of taking the equivalent elastic modulus required by design as an optimal target value. The method comprises the steps of taking the internal size of the porous structure of each target segment of the internal plant as an input variable, using a Direct Optimization method (Direct Optimization), and using a sequential quadratic programming method (NLPQL) as an Optimization algorithm for optimizing the size of the porous structure to obtain the optimal size of the porous structure and design the porous structure.
In the embodiment of the invention, the initial porous structure of each target segment is subjected to size optimization treatment according to the obtained optimal solution, and the porous structure with the elastic modulus gradient, which accords with the optimal solution, is obtained on the basis of a sequence quadratic programming method. According to the invention, by optimizing the size of the initial porous structure, the elastic modulus of the implant can be reduced, the stress shielding effect of the femur is weakened, and the success rate of the implant is improved.
In another embodiment of the present invention, before the parameter setting and optimizing process for the three-dimensional model of the femoral implant, the method further comprises:
and performing curved surface smoothing treatment on the outer surface of the femur of the three-dimensional model of the femur implant to obtain the three-dimensional model of the femur implant with a smooth surface.
Specifically, the curved surface smoothing treatment is to enable the femoral implant model to better perform finite element analysis and meet the quality requirement of finite element meshing.
In this embodiment, before the parameter setting and optimization process, the outer surface of the femur of the three-dimensional plant model in the femur needs to be smoothed. In this embodiment, preferably, the model is introduced into a finite element analysis software WorkBench, the spacecollaim component is used to open the three-dimensional model of the plant in the femur, and the irregular curved surface and the narrow curved surface of the outer surface of the femur are merged by using the curved surface merging function of the model, so as to obtain the three-dimensional model of the plant in the femur with a smooth surface.
In the embodiment of the invention, the femoral outer surface smoothing treatment is carried out on the three-dimensional model of the femoral implant before the parameter setting and optimization treatment is carried out on the model, so that the quality of the subsequent treatment of the model is ensured.
In another embodiment of the present invention, as shown in fig. 4, the obtaining a three-dimensional model of a human femur includes:
obtaining a femur model, wherein the femur model is obtained by scanning a human body;
and establishing a bone cement model with the same property as the femoral cancellous bone 2 in the femoral cortical bone 1 of the femoral model, and combining the femoral model and the bone cement model with each other to form the human femoral three-dimensional model.
Specifically, the femur model refers to a real human body model obtained by CT scanning.
In this embodiment, a femur model is obtained by scanning an actual human body, the interior of the solid femur model is hollowed, only a part of the exterior of the femur model corresponding to the femoral cortical bone 1 is reserved, a three-dimensional model of the cancellous bone of the femur is established according to the hollowed shape, the cancellous bone and the bone cement are combined into the bone cement to form a bone cement model, and the hollowed femur model and the bone cement model are combined to form the human femur three-dimensional model. It should be noted that the material properties of the femoral cancellous bone are similar to those of the cement bone cement used in the replacement operation.
The method comprises the steps of cutting off a femoral mushroom head in SolidWorks software according to a mode required by an operation, and drawing out a cavity for implanting an internal plant three-dimensional model in the femoral three-dimensional model according to the shape of the internal plant three-dimensional model along the cut-off mushroom head.
In the embodiment of the invention, the obtained human femur model is scanned, then the model is cut to be empty, the model and the generated bone cement model are combined to form the femur three-dimensional model which accords with the actual condition of a human body, and a cavity for implanting the plant three-dimensional model is cut to be empty in the model, so that the construction of the femur three-dimensional model is more accordant with the actual condition of the human body, and the femur three-dimensional model can adapt to crowds with different body characteristics by changing the parameter size, and has universality.
Fig. 5 is a schematic structural diagram of an orthopedic implant optimization device provided by the present invention, and as shown in fig. 5, the orthopedic implant optimization device provided by the present invention comprises:
an obtaining module 501, configured to obtain a human femur three-dimensional model;
an establishing module 502, configured to establish a corresponding three-dimensional model of an endoprosthesis according to the three-dimensional model of the human femur;
the embedding module 503 is configured to perform embedding processing on the human femur three-dimensional model and the endoprosthesis three-dimensional model to obtain a complete femur endoprosthesis three-dimensional model;
a processing module 504, configured to perform parameter setting and optimization processing on the three-dimensional femoral implant model, obtain an optimal solution of an equivalent elastic modulus of each target segment of the three-dimensional femoral implant model, and verify the optimal solution;
and an optimizing module 505, configured to perform size optimization on the initial porous structures of the target segments according to the optimal solution, and obtain a porous structure with an elastic modulus grade that meets the optimal solution.
Specifically, the three-dimensional model of the endoprosthesis is subjected to a segmentation process.
According to the device for optimizing the orthopedic implant, the three-dimensional model of the human femur is obtained through the obtaining module, the corresponding three-dimensional model of the implant is established according to the model, the two models are embedded through the embedding module to obtain a complete three-dimensional model of the implant in the femur, the processing module performs parameter setting and optimization processing on the three-dimensional model of the implant in the femur to obtain the optimal solution of the equivalent elastic modulus of each target section, and the optimizing module performs size optimization on the initial porous structure of each target section according to the optimal solution to obtain the porous structure which is in accordance with the optimal solution and has the elastic modulus grade. The invention can weaken the stress shielding effect of the thighbone while ensuring the strength, and improves the implantation success rate of the endophyte.
Since the principle of the apparatus according to the embodiment of the present invention is the same as that of the method according to the above embodiment, further details are not described herein for further explanation.
Fig. 6 is a schematic structural diagram of an electronic device according to an embodiment of the present invention, and as shown in fig. 6, the present invention provides an electronic device, including: a processor (processor)601, a memory (memory)602, and a bus 603;
the processor 601 and the memory 602 complete communication with each other through the bus 603;
An embodiment of the present invention provides a non-transitory computer-readable storage medium, where the non-transitory computer-readable storage medium stores computer instructions, and the computer instructions cause the computer to execute the method provided by the foregoing method embodiments, for example, the method includes: acquiring a human femur three-dimensional model; establishing a corresponding three-dimensional model of the internal implant according to the three-dimensional model of the human femur; embedding the human femur three-dimensional model and the implant three-dimensional model to obtain a complete femur implant three-dimensional model; parameter setting and optimization processing are carried out on the three-dimensional model of the femoral internal implant, the optimal solution of the equivalent elastic modulus of each target section of the three-dimensional model of the femoral internal implant is obtained, and the optimal solution is verified; and optimizing the size of the initial porous structure of each target segment according to the optimal solution to obtain a porous structure which accords with the optimal solution and has an elastic modulus grade.
Those of ordinary skill in the art will understand that: all or part of the steps for implementing the method embodiments may be implemented by hardware related to program instructions, and the program may be stored in a computer readable storage medium, and when executed, the program performs the steps including the method embodiments; and the aforementioned storage medium includes: various media that can store program codes, such as ROM, RAM, magnetic or optical disks.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.
Claims (10)
1. A method of orthopedic plant optimization, comprising:
acquiring a human femur three-dimensional model;
establishing a corresponding three-dimensional model of the internal implant according to the three-dimensional model of the human femur;
embedding the human femur three-dimensional model and the implant three-dimensional model to obtain a complete femur implant three-dimensional model;
parameter setting and optimization processing are carried out on the three-dimensional model of the femoral internal implant, the optimal solution of the equivalent elastic modulus of each target section of the three-dimensional model of the femoral internal implant is obtained, and the optimal solution is verified;
and optimizing the size of the porous structure of each target segment according to the optimal solution to obtain the porous structure with elastic modulus gradient which accords with the optimal solution.
2. The method for optimizing the orthopaedic implant according to claim 1, wherein before the embedding the three-dimensional model of the human femur and the three-dimensional model of the implant to obtain the complete three-dimensional model of the femoral implant, the method comprises:
performing segmentation processing on the three-dimensional model of the endoprosthesis;
wherein the endoprosthesis is cut into 8 segments.
3. The method for optimizing an orthopaedic implant according to claim 1, wherein the parameter setting of the three-dimensional model of the femoral implant comprises:
and carrying out material setting, grid division, contact setting, constraint setting, load loading and solving setting on the three-dimensional model of the plant in the femur to obtain initial stress distribution of the three-dimensional model of the plant in the femur.
4. The method for optimizing the orthopaedic implant according to claim 1, wherein the optimizing the three-dimensional model of the femoral implant comprises:
acquiring initial stress of seven gurn partitions on the outer surface of the femur in the three-dimensional femoral plant model and equivalent elastic modulus corresponding to seven internal plant target sections;
carrying out parametric modeling on the obtained initial stress and the equivalent elastic modulus to obtain a parametric model;
importing the parameterized model into a response surface optimization module, and carrying out test design on the parameterized model to obtain a response surface model;
and carrying out optimization design on the response surface model based on a multi-objective genetic optimization algorithm to obtain an optimal solution of the equivalent elastic modulus of each target segment.
And determining the response surface model by acquiring the sample points and the related function of the experimental design.
5. The method of orthopedic endoprosthesis optimization according to claim 4, wherein the optimizing the initial porous structure of each target segment in size according to the optimal solution comprises:
acquiring an initial porous structure of a plant in each target section in the three-dimensional model of the plant in the femur;
parameterizing the size and equivalent elastic modulus of the initial porous structure to be adjusted in the plant within each target segment;
and optimizing the size of the initial porous structure after the parameterization treatment based on a sequence quadratic programming method to obtain the porous structure with the elastic modulus gradient which accords with the optimal solution.
6. The method for optimizing orthopaedic implants according to claim 1, wherein before said parameter setting and optimizing process of said three-dimensional model of femoral implant, the method further comprises:
and performing curved surface smoothing treatment on the outer surface of the femur of the three-dimensional model of the femur implant to obtain the three-dimensional model of the femur implant with a smooth surface.
7. The method for optimizing an orthopaedic implant according to claim 1, wherein said obtaining a three-dimensional model of a human femur comprises:
obtaining a femur model, wherein the femur model is obtained by scanning a human body;
and establishing a bone cement model with the same property as the femoral cancellous bone in the femoral cortical bone of the femoral model, and combining the femoral model and the bone cement model to form the human femoral three-dimensional model.
8. An apparatus for orthopedic plant optimization, comprising:
the acquisition module is used for acquiring a human femur three-dimensional model;
the establishing module is used for establishing a corresponding three-dimensional model of the endoprosthesis according to the three-dimensional model of the human femur;
the embedding module is used for embedding the human femur three-dimensional model and the endoprosthesis three-dimensional model to obtain a complete femur endoprosthesis three-dimensional model;
the processing module is used for carrying out parameter setting and optimization processing on the three-dimensional model of the femoral internal implant, obtaining an optimal solution of the equivalent elastic modulus of each target section of the three-dimensional model of the femoral internal implant and verifying the optimal solution;
and the optimization module is used for optimizing the size of the initial porous structure of each target segment according to the optimal solution to obtain a porous structure with elastic modulus gradient which accords with the optimal solution.
9. An electronic device, comprising: a processor, a memory, and a bus, wherein,
the processor and the memory are communicated with each other through the bus;
the memory stores program instructions executable by the processor, the processor invoking the program instructions to perform the method of any of claims 1 to 7.
10. A non-transitory computer-readable storage medium storing computer instructions that cause a computer to perform the method of any one of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110406442.6A CN113284569A (en) | 2021-04-15 | 2021-04-15 | Method and device for optimizing plants in orthopedics department, electronic equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110406442.6A CN113284569A (en) | 2021-04-15 | 2021-04-15 | Method and device for optimizing plants in orthopedics department, electronic equipment and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113284569A true CN113284569A (en) | 2021-08-20 |
Family
ID=77276749
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110406442.6A Pending CN113284569A (en) | 2021-04-15 | 2021-04-15 | Method and device for optimizing plants in orthopedics department, electronic equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113284569A (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105740523A (en) * | 2016-01-26 | 2016-07-06 | 西安交通大学 | Host bone stress environment based custom prosthesis optimization design method |
| CN105997306A (en) * | 2016-04-25 | 2016-10-12 | 北京工业大学 | Design method for filling porous grid structure in bone implantation body |
| CN107391784A (en) * | 2017-06-03 | 2017-11-24 | 大连理工大学 | A kind of cancellous bone loose structure modeling method based on topological optimization technology |
| CN109977590A (en) * | 2019-04-09 | 2019-07-05 | 哈尔滨理工大学 | A kind of titanium implant structural optimization method based on age factor |
| CN110840626A (en) * | 2019-11-05 | 2020-02-28 | 吉林大学 | Femoral stem prosthesis design method |
| CN110866354A (en) * | 2019-11-08 | 2020-03-06 | 大连理工大学 | Optimized design method of polymer vascular stent structure considering scale effect |
| CN111967172A (en) * | 2020-07-06 | 2020-11-20 | 华南理工大学 | Optimal design method of internal high-pressure forming process based on kriging model |
-
2021
- 2021-04-15 CN CN202110406442.6A patent/CN113284569A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105740523A (en) * | 2016-01-26 | 2016-07-06 | 西安交通大学 | Host bone stress environment based custom prosthesis optimization design method |
| CN105997306A (en) * | 2016-04-25 | 2016-10-12 | 北京工业大学 | Design method for filling porous grid structure in bone implantation body |
| CN107391784A (en) * | 2017-06-03 | 2017-11-24 | 大连理工大学 | A kind of cancellous bone loose structure modeling method based on topological optimization technology |
| CN109977590A (en) * | 2019-04-09 | 2019-07-05 | 哈尔滨理工大学 | A kind of titanium implant structural optimization method based on age factor |
| CN110840626A (en) * | 2019-11-05 | 2020-02-28 | 吉林大学 | Femoral stem prosthesis design method |
| CN110866354A (en) * | 2019-11-08 | 2020-03-06 | 大连理工大学 | Optimized design method of polymer vascular stent structure considering scale effect |
| CN111967172A (en) * | 2020-07-06 | 2020-11-20 | 华南理工大学 | Optimal design method of internal high-pressure forming process based on kriging model |
Non-Patent Citations (1)
| Title |
|---|
| SAMI EMAD ALKHATIB 等: "Finite Element Analysis of Porous Titanium Alloy Hip Stem to Evaluate the Biomechanical Performance During Walking and Stair Climbing", JOURNAL OF BIONIC ENGINEERING, 6 December 2019 (2019-12-06), pages 1103 - 1115 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wieding et al. | Numerical optimization of open-porous bone scaffold structures to match the elastic properties of human cortical bone | |
| US20160256279A1 (en) | Patient-Specific Implant for Bone Defects and Methods for Designing and Fabricating Such Implants | |
| San Cheong et al. | Novel adaptive finite element algorithms to predict bone ingrowth in additive manufactured porous implants | |
| Helgason et al. | Risk of failure during gait for direct skeletal attachment of a femoral prosthesis: a finite element study | |
| Lacroix et al. | Three-dimensional simulation of fracture repair in the human tibia | |
| CN105748177B (en) | A kind of personalized implantable spinal prosthesis and its manufacture method with bionic micropore | |
| CN109091273A (en) | A kind of personalized pelvis endoprosthesis designs method based on topological structure optimization | |
| Cheong et al. | Partial bone formation in additive manufactured porous implants reduces predicted stress and danger of fatigue failure | |
| KR20160091322A (en) | Bone reconstruction and orthopedic implants | |
| Kluess | Finite element analysis in orthopaedic biomechanics | |
| US20190240029A1 (en) | Method for adjusting mechanical properties of implant and patient specific surgical implants | |
| CN112233801B (en) | Construction method and topological optimization design method of topological optimization mathematical model of endoprosthesis | |
| Ait Moussa et al. | Minimizing stress shielding and cement damage in cemented femoral component of a hip prosthesis through computational design optimization | |
| Rahmati et al. | An improved methodology for design of custom‐made hip prostheses to be fabricated using additive manufacturing technologies | |
| Van Cleynenbreugel et al. | Trabecular bone scaffolding using a biomimetic approach | |
| Liu et al. | Design of porous metal block augmentation to treat tibial bone defects in total knee arthroplasty based on topology optimization | |
| CN112220588B (en) | Method and system for generating controllable gradient bone tissue engineering scaffold | |
| KR20190100837A (en) | An elastic modulus adjustment method for implants having lattice scaffolds structure and patient specific surgical implants which use the method | |
| Wojnicz et al. | Mathematical approach to design 3D scaffolds for the 3D printable bone implant | |
| CN113284569A (en) | Method and device for optimizing plants in orthopedics department, electronic equipment and storage medium | |
| Gkaintes et al. | Computational modelling of human femur after total hip arthroplasty | |
| Polgar et al. | Development and numerical validation of a finite element model of the muscle standardized femur | |
| CN111297524B (en) | Knee joint defect cushion block and manufacturing method thereof | |
| Żmudzki et al. | Influence of model discretization density in FEM numerical analysis on the determined stress level in bone surrounding dental implants | |
| Khanoki et al. | Multiscale design and multiobjective optimization of orthopaedic cellular hip implants |
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 |