Detailed Description
In order to more clearly illustrate the invention or the technical solutions of the prior art, the following description will briefly explain the drawings used in the embodiments or the description of the prior art, and it is obvious that the drawings in the following description are some embodiments of the invention, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.
Fig. 1 is a general flow diagram of the orthopedic implant optimization method provided by the invention. As shown in fig. 1, the orthopedic implant optimization method provided by the invention comprises the following steps:
Step 101: acquiring a three-dimensional model of the femur of a human body;
step 102: establishing a corresponding endophyte three-dimensional model according to the human femur three-dimensional model;
step 103: embedding the human femur three-dimensional model and the endophyte three-dimensional model to obtain a complete femur endophyte three-dimensional model;
Step 104: performing parameter setting and optimization treatment on the three-dimensional model of the intra-femoral plant, obtaining an optimal solution of the equivalent elastic modulus of each target segment of the three-dimensional model of the intra-femoral plant, and verifying the optimal solution;
Step 105: and performing size optimization on the initial porous structure of each target segment according to the optimal solution to obtain the porous structure with the elastic modulus gradient conforming to the optimal solution.
Specifically, the endophyte refers to a device for implantation into a femur of a human body, and the endophyte material may be stainless steel, cobalt-chromium alloy, titanium or titanium alloy, etc., and in this embodiment, the preferred endophyte is TC4 titanium alloy (Ti-6 Al-4V). Where each target segment refers to each part of the model that needs to be processed.
In the implementation, a three-dimensional model of a human femur is firstly obtained, an endophyte three-dimensional model is built according to the three-dimensional model of the human femur in SolidWorks software, the two models are embedded and combined to obtain a complete three-dimensional model of the endophyte of the femur, and the model is exported as an X_T file; then, the complete three-dimensional model of the intra-femur plant is imported into finite element analysis software WorkBench, parameter setting and optimization treatment are carried out on the three-dimensional model of the intra-femur plant, the optimal solution of each target segment of the three-dimensional model of the intra-femur plant is obtained, and the optimal solution of the equivalent elastic modulus of each target segment is verified; and finally, performing size optimization on the initial porous structure of each target segment according to the optimal solution to obtain the porous structure with the elastic modulus gradient, which accords with the optimal solution.
In the embodiment of the invention, a corresponding endophyte three-dimensional model is established by acquiring a human femur three-dimensional model, the endophyte three-dimensional model is embedded into the human femur three-dimensional model to acquire a complete femur endophyte three-dimensional model, parameter design and optimization processing are carried out on the femur endophyte three-dimensional model, the optimal solution of each target segment is acquired, and the porous structure of each target segment is optimized in size according to the optimal solution, so that the porous structure meeting the requirements is acquired. According to the invention, through optimizing the three-dimensional model and the porous structure of the intra-femur plant, the stress shielding effect of the femur can be weakened while the intensity of the intra-femur plant is ensured, and the success rate of the intra-femur plant implantation is improved.
In another embodiment of the present invention, the embedding process of the three-dimensional model of femur and the three-dimensional model of endophyte of human body is performed, before obtaining the complete three-dimensional model of endophyte of femur, comprising:
carrying out sectional treatment on the endophyte three-dimensional model;
Wherein the endophyte is cut into 8 segments.
Specifically, the endophyte three-dimensional model is a model which is obtained by three-dimensional software design according to a femur three-dimensional model which accords with an actual human body, and has more reality.
In this embodiment, as shown in fig. 2, the preferred endophyte is a titanium alloy endophyte, the titanium alloy endophyte model is segmented into eight sections a-H by adopting a segmentation design mode, the segmented endophyte 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 setting constraints.
In the embodiment of the invention, the three-dimensional model of the endophyte is subjected to sectional treatment, so that the individual treatment of the endophyte in each section is realized, the stress shielding effect is reduced, the endophyte still has higher strength in a low-sensitivity place by the characteristic of gradient, and the condition of insufficient strength caused by the integral regulation setting of the elastic modulus is avoided.
In one embodiment of the present invention, as shown in fig. 3, the parameter setting of the three-dimensional model of the intra-femur plant includes:
And performing material setting, grid division, contact setting, constraint setting, load loading and solving setting on the three-dimensional model of the intra-femur plant to obtain initial stress distribution of the three-dimensional model of the intra-femur plant.
Specifically, the material of the inner plant is titanium alloy; the parameter settings are operated in the finite element analysis software WorkBench, but may also be operated in other analysis software, and are not particularly limited herein.
In this embodiment, the material properties of each part of the three-dimensional model of the plant in the femur, such as the cortical bone of the outer layer of the femur, are set, the elastic modulus is set to 13.7GPa, and the poisson ratio is 0.33; the elastic modulus of the bone cement is set to 2.2GPa, and the Poisson ratio is 0.3; the plant material in the titanium alloy is Ti-6Al-4V, the elastic modulus is set to 114GPa, and the Poisson ratio is 0.36.
In the embodiment, the mesh division of the three-dimensional model of the plant in the femur is set, the preferred overall size in the embodiment is 2mm, wherein the plant in the titanium alloy is divided into meshes by adopting the methods MultiZone and Hex Dominant, and the mesh quality is controlled by using the surface size in a narrow part; dividing the bone cement into grids by Hex Dominant and Tetrahedrons methods, and adjusting the grid quality of a narrow area by using the surface size and the edge size; the femur cortical bone is meshed using Hex Dominant methods, others using global settings.
In the embodiment, contact between each part in the three-dimensional model of the plant in the femur is set, wherein the three-dimensional model of the plant in the titanium alloy is required to be set as an integral part, and all target sections are mutually influenced; the titanium alloy inner plant is in binding contact with bone cement, and the titanium alloy inner plant is in binding contact with femur cortical bone. The constraint condition of the distal end of the three-dimensional model of the plant in the femur is set, and the plane of distal end excision is set as fixed constraint.
In this embodiment, the load loading of the three-dimensional model of the intra-femoral plant is set, considering that the general weight of asian males is 70KG, the downward pressure on the femur is set to 1400N in the case of standing on one leg, and this pressure acts on the part of the titanium alloy intra-femoral plant tip that replaces the femoral mushroom head, in a direction down the femur.
In the embodiment, the stress solving of the three-dimensional model of the plant in the femur is set, and the set solving form is Von-Mises equivalent stress, wherein the equivalent stress of 7 Gruen areas on the outer surface of the femur is particularly required to be solved, so that the initial stress distribution condition of the plant in the titanium alloy before the optimization is obtained.
In the embodiment of the invention, the initial stress distribution of the three-dimensional model of the plant in the femur is obtained by setting parameters of the three-dimensional model of the plant in the femur, which mainly relate to 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 three-dimensional model of the intra-femoral plant can better optimize the three-dimensional model of the intra-femoral plant.
In another embodiment of the present invention, as shown in fig. 3, the optimizing the three-dimensional model of the intra-femur plant includes:
Acquiring initial stress of seven Gran partitions on the outer surface of a femur in the three-dimensional model of the femur inner plant and equivalent elastic modulus corresponding to seven target segments of the inner plant;
carrying out parameterized modeling on the initial stress and the equivalent elastic modulus to obtain a parameterized model;
the parameterized model is imported into a response surface optimization module, and experimental design is carried out on the parameterized model to obtain a response surface model;
and carrying out optimal 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 objective segment.
And determining the response surface model by acquiring the number of sample points and the correlation function of the test design.
Specifically, DOE (DESIGN OF EXPERIMENT, design of experiments) plays a very important role in the whole process of quality control, and is an important guarantee for improving product quality and improving process flow.
In the embodiment, initial stress and equivalent elastic modulus of each target segment of the three-dimensional model of the plant in the femur are obtained, then parameterized treatment is carried out on the obtained initial stress and equivalent elastic modulus to generate a parameterized model, the parameterized model is led into a response surface module, DOE test design is carried out on the parameterized model, and a response surface model is obtained; and optimizing the response surface model based on a multi-objective genetic optimization algorithm to obtain the optimal solution of each objective segment.
Wherein, firstly, the part for optimization is parameterized, including equivalent stress values of 7 Gruen areas and equivalent elastic modulus in material properties of the rear 7 sections in the titanium alloy endophyte, and the parameterized data is imported 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 response surface model optimization setting, firstly, sampling is needed by using DOE test design, 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, the coverage is uniform, the scale of the test can be effectively reduced, and the specific implementation steps are as follows:
step 1: dividing each dimension into m sections which are not overlapped with each other, so that each section has the same probability;
step 2: a point is randomly extracted in each interval in each dimension to form the total Latin hypercube.
Step3: the selected point component vectors in the population are randomly extracted in each dimension, and once each layer of samples is extracted, the layer is no longer 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. The Kriging model is a semi-parameterized interpolation technology based on a statistical theory, the accuracy of interpolation unknown information in known information of the model is very high, the model generally consists of a regression part and a random part, and based on the Kriging method, the output parameter is equal to the global design space plus local deviation, and the expression is as follows:
y(x)=F(β,x)+z(x)=fT(x)β+z(x)
Where β is the basis function regression coefficient, f (x) is a polynomial function of the variable x, representing the global model of the design space, and z (x) is a gaussian random function with a mean of 0 and a variance of σ 2.
Forming a correlation matrix by design points obtained based on Latin hypercube experimental design method:
The method comprises the steps of determining a final Kriging response surface model through sample points and related functions of test design, wherein n is the total number of data points, and the response surface model has good effect on solving nonlinear engineering optimization problems, and has the characteristics of high calculation efficiency, short time, good response surface fitting effect, high result accuracy and the like.
After the response surface model is established, optimizing design of the equivalent elastic modulus is carried out based on Optimization, a genetic algorithm used is a multi-objective genetic algorithm (MOGA), the maximum stress value of 7 Gruen areas on the outer surface of the femur is used as target stress, the equivalent elastic modulus of 7 sections of titanium alloy inner plants is used as an input variable, and the iteration times are set for final objective elastic modulus optimizing, so that the optimal reference point equivalent elastic modulus is obtained.
And the obtained optimal solution is guided into the three-dimensional model of the femoral internal plant after the parameter setting is completed for verification, the obtained optimal elastic modulus is guided into the initial titanium alloy internal plant model of stress analysis, the equivalent elastic modulus of the 7 sections of the lower end of the titanium alloy internal plant is reset, the rest conditions are unchanged, and the stress analysis verification and optimization model effect is carried out again.
In the embodiment of the invention, an initial stress and an equivalent elastic model of the model are obtained, the parameterized model is obtained, then the parameterized model is imported into a response surface module for DOE test design, 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 invention, the optimal solution of the equivalent elastic modulus of each target segment is obtained through optimizing the model, so that the elastic modulus of the inner plant can be reduced, and the femur stress shielding effect 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 includes:
acquiring an initial porous structure of plants in each target segment in the three-dimensional model of the intra-femoral plant;
parameterizing the size and equivalent elastic modulus of the initial porous structure to be regulated in the plants within each target segment;
And performing size optimization on the initial porous structure after parameterization 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, sequential quadratic programming algorithms are one of the most effective methods currently accepted to solve constrained nonlinear optimization problems, and specific implementations are not specifically stated herein.
In this embodiment, the size optimization process is performed on the initial porous structure of each target segment according to the obtained optimal solution value, and it should be noted that the size and the equivalent elastic modulus of the initial porous structure 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 cube cell, and the porous structure is characterized in that the whole porosity of the cube is influenced by changing one size, and the equivalent elastic modulus of plants in the titanium alloy is reduced by changing the porosity.
The size of the porous structure of the inner cell which is needed to be regulated by the plant in each section of titanium alloy and the equivalent elastic modulus are parameterized, and the size optimization treatment of the porous structure is also carried out by WorkBench so as to achieve the aim that the equivalent elastic modulus required by design is the optimal target value. The method comprises the steps of taking the internal size of a porous structure of each target segment of an endophyte as an input variable, using a direct optimization method (Direct Optimization), and using a sequence quadratic programming method (NLPQL) as an optimizing 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 conforming to the optimal solution is obtained based on a sequence quadratic programming method. According to the invention, the elastic modulus of the endophyte can be reduced, the femur stress shielding effect is weakened, and the success rate of endophyte implantation is improved by optimizing the size of the initial porous structure.
In another embodiment of the present invention, before the parameter setting and optimizing the three-dimensional model of the intra-femoral plant, the method further includes:
and carrying out curved surface smoothing treatment on the outer surface of the femur of the three-dimensional model of the intra-femur plant to obtain the three-dimensional model of the intra-femur plant with smooth surface.
Specifically, the curved surface smoothing is used for enabling the plant model in the femur to better perform finite element analysis and achieving the quality requirement of finite element mesh division.
In this embodiment, the outer surface of the femur of the three-dimensional model of the intra-femur plant also needs to be smoothed before the parameter setting and optimization process is performed. In this embodiment, the model is preferably imported into finite element analysis software WorkBench, and the three-dimensional model of the intra-femur plant is opened by using SPACECLAIM component, and the irregular curved surface and the narrow curved surface of the outer surface of the femur are combined by using the curved surface combining function, so as to obtain the three-dimensional model of the intra-femur plant with a smooth surface.
In the embodiment of the invention, before the parameter setting and optimizing treatment are carried out on the three-dimensional model of the plant in the femur, the model is subjected to the femur external surface smoothing treatment, 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 acquiring a three-dimensional model of a femur of a human body 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 attribute as that of the femur cancellous bone 2 in the femur cortical bone 1 of the femur model, and combining the femur model and the bone cement model to form the human femur three-dimensional model.
Specifically, the femur model refers to a real human model, obtained by CT scan.
In this embodiment, a femur model is obtained by scanning an actual human body, the inside of the femur model is hollowed, only the part of the outside corresponding to the femur cortical bone 1 is reserved, a three-dimensional model of the femur cancellous bone is built according to the hollowed shape, cancellous bone and bone cement are combined into bone cement to form a bone cement model, and the hollowed femur model and the bone cement model are combined to form the human body femur three-dimensional model. The material properties of the cancellous bone of the femur were similar to those of the cement of the adhesive in the replacement surgery.
The method comprises the steps of cutting out a femur mushroom head according to the operation requirement in SolidWorks software, and drawing out a cavity for implanting an endophytic three-dimensional model in the femur three-dimensional model according to the shape of the endophytic three-dimensional model along the position of the cut-out mushroom head.
In the embodiment of the invention, the human femur model is obtained through scanning, then the model is cut open, the model and the generated bone cement model are combined to form the femur three-dimensional model which accords with the actual condition of the human body, the cavity for implanting the endophyte three-dimensional model is cut open in the model, the construction of the femur three-dimensional model accords with the actual condition of the human body, and the model can be suitable for people with different physical characteristics through changing the size of parameters, so that the model has universality.
Fig. 5 is a schematic structural diagram of an orthopedic implant optimizing apparatus provided by the present invention, and as shown in fig. 5, the orthopedic implant optimizing apparatus provided by the present invention includes:
The acquisition module 501 is used for acquiring a three-dimensional model of the femur of the human body;
the building module 502 is configured to build a corresponding endophyte three-dimensional model according to the human femur three-dimensional model;
An embedding module 503, configured to perform an embedding process on the three-dimensional model of the femur of the human body and the three-dimensional model of the endophyte, so as to obtain a complete three-dimensional model of the endophyte of the femur;
The processing module 504 is configured to perform parameter setting and optimization processing on the three-dimensional model of the intra-femoral plant, obtain an optimal solution of an equivalent elastic modulus of each target segment of the three-dimensional model of the intra-femoral plant, and verify the optimal solution;
And the optimizing module 505 is configured to perform size optimization on the initial porous structure of each target segment according to the optimal solution, and obtain a porous structure with an elastic modulus grade according to the optimal solution.
Specifically, the endophyte three-dimensional model is subjected to segmentation processing.
According to the orthopedic endophyte optimizing device provided by the embodiment of the invention, the acquisition module is used for acquiring the human femur three-dimensional model, the corresponding endophyte three-dimensional model is established according to the model, the two models are embedded through the embedding module to obtain the complete femur endophyte three-dimensional model, the processing module is used for carrying out parameter setting and optimizing processing on the femur endophyte three-dimensional model to obtain the optimal solution of the equivalent elastic modulus of each target segment, and the optimizing module is used for carrying out size optimization on the initial porous structure of each target segment according to the optimal solution to obtain the porous structure with the elastic modulus grade, which accords with the optimal solution. The invention can weaken the stress shielding effect of femur while ensuring the strength, and improve the implantation success rate of endophytes.
Since the apparatus according to the embodiment of the present invention is the same as the method according to the above embodiment, the details of the explanation will not be repeated here.
Fig. 6 is a schematic diagram of an entity structure 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;
Wherein, the processor 601 and the memory 602 complete communication with each other through the bus 603;
The processor 601 is configured to invoke program instructions in the memory 602 to perform the methods provided by the method embodiments described above, including, for example: acquiring a three-dimensional model of the femur of a human body; establishing a corresponding endophyte three-dimensional model according to the human femur three-dimensional model; embedding the human femur three-dimensional model and the endophyte three-dimensional model to obtain a complete femur endophyte three-dimensional model; performing parameter setting and optimization treatment on the three-dimensional model of the intra-femoral plant, obtaining an optimal solution of the equivalent elastic modulus of each target segment of the three-dimensional model of the intra-femoral plant, and verifying the optimal solution; and performing size optimization on the initial porous structure of each target segment according to the optimal solution to obtain the porous structure with the elastic modulus grade, which accords with the optimal solution.
Embodiments of the present invention provide a non-transitory computer readable storage medium storing computer instructions that cause a computer to perform the methods provided by the above-described method embodiments, for example, including: acquiring a three-dimensional model of the femur of a human body; establishing a corresponding endophyte three-dimensional model according to the human femur three-dimensional model; embedding the human femur three-dimensional model and the endophyte three-dimensional model to obtain a complete femur endophyte three-dimensional model; performing parameter setting and optimization treatment on the three-dimensional model of the intra-femoral plant, obtaining an optimal solution of the equivalent elastic modulus of each target segment of the three-dimensional model of the intra-femoral plant, and verifying the optimal solution; and performing size optimization on the initial porous structure of each target segment according to the optimal solution to obtain the porous structure with the elastic modulus grade, which accords with the optimal solution.
Those of ordinary skill in the art will appreciate that: all or part of the steps for implementing the above method embodiments may be implemented by hardware associated with program instructions, where the foregoing program may be stored in a computer readable storage medium, and when executed, the program performs steps including the above method embodiments; and the aforementioned storage medium includes: various media that can store program code, such as ROM, RAM, magnetic or optical disks.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and are not limiting; although the 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 scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.