[go: up one dir, main page]

CN113284569B - A method, device, electronic device and storage medium for optimizing orthopedic implants - Google Patents

A method, device, electronic device and storage medium for optimizing orthopedic implants Download PDF

Info

Publication number
CN113284569B
CN113284569B CN202110406442.6A CN202110406442A CN113284569B CN 113284569 B CN113284569 B CN 113284569B CN 202110406442 A CN202110406442 A CN 202110406442A CN 113284569 B CN113284569 B CN 113284569B
Authority
CN
China
Prior art keywords
dimensional model
model
implant
femur
optimal solution
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110406442.6A
Other languages
Chinese (zh)
Other versions
CN113284569A (en
Inventor
高相胜
赵宇航
李紫茜
王辰露
郭子源
王�琦
刘成才
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN202110406442.6A priority Critical patent/CN113284569B/en
Publication of CN113284569A publication Critical patent/CN113284569A/en
Application granted granted Critical
Publication of CN113284569B publication Critical patent/CN113284569B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00Computational 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite 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)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Prostheses (AREA)

Abstract

本发明提供一种骨科内植物优化的方法、装置、电子设备及存储介质,该方法包括:获取人体股骨三维模型;根据所述人体股骨三维模型建立对应的内植物三维模型;将所述人体股骨三维模型和内植物三维模型进行嵌入处理,获取完整的股骨内植物三维模型;对所述股骨内植物三维模型进行参数设置及优化处理,获取所述股骨内植物三维模型的各个目标节段的最优解,并验证所述最优解;根据所述最优解对所述各个目标节段的多孔结构进行尺寸优化,获取符合所述最优解的多孔结构。本发明通过对内植物进行优化后得到的具有材料性能梯度的多孔结构,能够在保证强度的同时尽可能的削弱股骨应力遮蔽效应,提高了内植物植入的成功率。

The present invention provides a method, device, electronic device and storage medium for optimizing orthopedic implants, the method comprising: obtaining a three-dimensional model of a human femur; establishing a corresponding three-dimensional model of an implant based on the three-dimensional model of the human femur; embedding the three-dimensional model of the human femur and the three-dimensional model of the implant to obtain a complete three-dimensional model of the femoral implant; setting parameters and optimizing the three-dimensional model of the femoral implant to obtain the optimal solution of each target segment of the three-dimensional model of the femoral implant, and verifying the optimal solution; optimizing the size of the porous structure of each target segment according to the optimal solution to obtain a porous structure that meets the optimal solution. The present invention obtains a porous structure with a material property gradient after optimizing the implant, which can weaken the femoral stress shielding effect as much as possible while ensuring strength, thereby improving the success rate of implant implantation.

Description

Method and device for optimizing orthopedic endophytes, electronic equipment and storage medium
Technical Field
The invention relates to the technical field of medical appliances, in particular to a method, a device, electronic equipment and a storage medium for optimizing orthopedic endophytes.
Background
With the gradual progress of China society to the aging direction, the incidence rate of the fracture disease of the femoral neck of the aged is increased year by year. According to relevant clinical reports, the artificial hip joint replacement is applied to the clinical treatment of senile femoral neck fracture diseases, so that satisfactory clinical curative effects can be obtained, adverse events can be reduced, and the safety of treatment is ensured.
At present, the biomechanical compatibility of the endophyte material is important in the application of the artificial hip joint field, the biomechanical compatibility has important influence on the effectiveness and reliability of the material after being implanted into a human body, the excellent biomechanical compatibility can ensure the friendly fusion of the endophyte and human bones, and the rejection reaction is avoided, wherein the elastic modulus is one of the most important physical properties for judging the biocompatibility of the biomedical metal material, particularly the bone substitute material. In the prior art, the used and mature endophyte 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 bone, and too high metal elastic modulus can cause loosening at the interface of endophyte and human bone, influence the function of an implanted device, or cause stress shielding effect, cause functional degradation or absorption of bone tissue and easily cause secondary injury such as postoperative osteoporosis; 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 the inner plant cannot be fixed and supported.
Disclosure of Invention
The invention provides a method, a device, electronic equipment and a storage medium for optimizing orthopedic endophytes, which are used for solving the technical problem that secondary injuries such as postoperative osteoporosis and the like are easily caused by stress shielding due to the fact that the elastic modulus of the endophytes is higher than that of human bones, so as to achieve the purposes of weakening the stress shielding effect and improving the success rate of endophytes implantation.
In a first aspect, the present invention provides a method of orthopedic endophyte optimization, comprising:
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 gradient conforming to the optimal solution.
In a second aspect, the present invention provides an orthopedic endophyte optimizing device comprising:
the acquisition module is used for acquiring the data of the electronic device, the method is used for obtaining a three-dimensional model of the femur of the human body;
the building module is used for building a corresponding endophyte three-dimensional model according to the human femur three-dimensional model;
the embedding module is used for carrying out embedding treatment on the human femur three-dimensional model and the endophyte three-dimensional model to obtain a complete femur endophyte 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 plant, obtaining an optimal solution of the equivalent elastic modulus of each target segment of the three-dimensional model of the femoral internal plant, and verifying the optimal solution;
And the optimization 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 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 complete communication 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 present invention provides a non-transitory computer readable storage medium storing computer instructions that cause a computer to perform a method as described in any of the above.
According to the method, the human femur three-dimensional model is obtained, then the corresponding endophyte three-dimensional model is built according to the human femur three-dimensional model, the human femur three-dimensional model and the endophyte three-dimensional model are subjected to embedding treatment, the complete endophyte three-dimensional model is obtained, parameter setting and optimization treatment are carried out on the complete endophyte three-dimensional model, the optimal solution of the equivalent elastic modulus of each target segment is obtained, and the initial porous structure of each target segment is subjected to size optimization treatment according to the obtained optimal solution, so that the porous structure meeting the requirements is obtained. According to the invention, through the mode of optimizing the porous structure of each target section of the three-dimensional model of the intra-femoral plant, the stress shielding effect of the femur is weakened while the strength of the intra-femoral plant is ensured, so that the success rate of the intra-femoral plant implantation is improved.
Drawings
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 an orthopedic implant optimization method provided by the invention;
fig. 2 is a schematic structural view of an orthopedic implant segmentation provided by the invention;
FIG. 3 is a schematic flow chart of the orthopedic implant optimization method provided by the invention;
Fig. 4 is a schematic structural diagram of a three-dimensional model of a human femur in the orthopedic implant optimization method provided by the invention;
FIG. 5 is a schematic structural view of an orthopedic implant optimizing device provided by the invention;
Fig. 6 is a schematic structural diagram of an electronic device provided by the present invention.
Reference numerals illustrate:
1-femoral cortical bone; 2-femur cancellous bone (bone cement);
3-titanium alloy inner plants.
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.

Claims (8)

1.一种骨科内植物优化的方法,其特征在于,包括:1. A method for optimizing an orthopedic implant, comprising: 获取人体股骨三维模型;Obtain a three-dimensional model of a human femur; 根据所述人体股骨三维模型建立对应的内植物三维模型;Establishing a corresponding three-dimensional model of the 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; 对所述股骨内植物三维模型进行参数设置及优化处理,获取所述股骨内植物三维模型的各个目标节段的等效弹性模量的最优解,并验证所述最优解;Performing parameter setting and optimization processing on the femoral implant three-dimensional model, obtaining an optimal solution of the equivalent elastic modulus of each target segment of the femoral implant three-dimensional model, and verifying the optimal solution; 根据所述最优解对所述各个目标节段的多孔结构进行尺寸优化,获取符合所述最优解的具有弹性模量梯度的多孔结构;其中,所述各个目标节段为所述股骨内植物三维模型需要处理的各个部分;According to the optimal solution, the size of the porous structure of each target segment is optimized to obtain a porous structure with an elastic modulus gradient that meets the optimal solution; wherein each target segment is a part of the femoral implant three-dimensional model that needs to be processed; 所述对所述股骨内植物三维模型进行优化处理,获取所述股骨内植物三维模型的各个目标节段的等效弹性模量的最优解,包括:The step of optimizing the femoral implant three-dimensional model to obtain an optimal solution of the equivalent elastic modulus of each target segment of the femoral implant three-dimensional model includes: 获取所述股骨内植物三维模型中股骨外表面七个格恩分区的初始应力和七个内植物目标节段对应的等效弹性模量;Obtaining the initial stress of seven Gern zones on the outer surface of the femur and the equivalent elastic modulus corresponding to seven implant target segments in the three-dimensional model of the femoral implant; 将获取的所述初始应力和所述等效弹性模量进行参数化建模,得到参数化模型;Performing parameterized modeling on the acquired initial stress and equivalent elastic modulus to obtain a parameterized model; 将所述参数化模型导入响应面优化模块中,并对所述参数化模型进行试验设计,获取响应面模型;Importing the parameterized model into a response surface optimization module, and performing experimental design on the parameterized model to obtain a response surface model; 基于多目标遗传优化算法对所述响应面模型进行优化设计,获取所述各个目标节段的等效弹性模量最优解;Optimizing 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; 其中,通过获取所述试验设计的样本点数和相关函数确定所述响应面模型;Wherein, the response surface model is determined by obtaining the number of sample points and the correlation function of the experimental design; 所述根据所述最优解对所述各个目标节段的初始多孔结构进行尺寸优化,包括:The step of optimizing the size of the initial porous structure of each target segment according to the optimal solution comprises: 获取所述股骨内植物三维模型中的各个目标节段内植物的初始多孔结构;Acquiring the initial porous structure of each target segment implant in the three-dimensional model of the femoral implant; 将各个目标节段内植物中需要调节的初始多孔结构的尺寸和等效弹性模量参数化;parameterizing the size and equivalent elastic modulus of the initial porous structure to be adjusted in the plant within each target segment; 基于序列二次规划法对参数化处理后的所述初始多孔结构进行尺寸优化,获取符合所述最优解的具有弹性模量梯度的多孔结构。The size of the initial porous structure after parameterization is optimized based on a sequential quadratic programming method to obtain a porous structure with an elastic modulus gradient that meets the optimal solution. 2.根据权利要求1所述的骨科内植物优化的方法,其特征在于,将所述人体股骨三维模型和内植物三维模型进行嵌入处理,获取完整的股骨内植物三维模型之前,包括:2. The method for optimizing orthopedic implants according to claim 1, characterized in that before embedding the human femur three-dimensional model and the implant three-dimensional model to obtain a complete femur implant three-dimensional model, the method comprises: 对所述内植物三维模型进行分段处理;Segmenting the three-dimensional model of the implant; 其中,将所述内植物切分为8个节段。Wherein, the implant is cut into 8 segments. 3.根据权利要求1所述的骨科内植物优化的方法,其特征在于,所述对所述股骨内植物三维模型进行参数设置,包括:3. The method for optimizing orthopedic implants according to claim 1, wherein the step of setting parameters of the femoral implant three-dimensional model comprises: 对所述股骨内植物三维模型进行材料设置、网格划分、接触设置、约束设置、载荷加载和求解设置,获取所述股骨内植物三维模型的初始应力分布。The material setting, mesh division, contact setting, constraint setting, load loading and solution setting are performed on the three-dimensional model of the femoral implant to obtain the initial stress distribution of the three-dimensional model of the femoral implant. 4.根据权利要求1所述的骨科内植物优化的方法,其特征在于,在所述对所述股骨内植物三维模型进行参数设置及优化处理之前,该方法还包括:4. The method for optimizing orthopedic implants according to claim 1, characterized in that before the parameters of the femoral implant three-dimensional model are set and optimized, the method further comprises: 对所述股骨内植物三维模型的股骨外表面进行曲面光滑处理,获取表面光滑的股骨内植物三维模型。The outer surface of the femur of the femoral implant three-dimensional model is subjected to curved surface smoothing processing to obtain a femoral implant three-dimensional model with a smooth surface. 5.根据权利要求1所述的骨科内植物优化的方法,其特征在于,所述获取人体股骨三维模型,包括:5. The method for optimizing orthopedic implants according to claim 1, wherein obtaining a three-dimensional model of a human femur comprises: 获取股骨模型,所述股骨模型通过扫描人体得到;Acquire a femur model, wherein the femur model is obtained by scanning a human body; 在所述股骨模型的股骨皮质骨内部建立与股骨松质骨属性相同的骨水泥模型,将所述股骨模型与所述骨水泥模型相互结合构成所述人体股骨三维模型。A bone cement model having the same properties as the femoral cancellous bone is established inside the femoral cortical bone of the femoral model, and the femoral model and the bone cement model are combined to form the human femur three-dimensional model. 6.一种骨科内植物优化的装置,其特征在于,包括:6. A device for optimizing orthopedic implants, comprising: 获取模块,用于获取人体股骨三维模型;An acquisition module, used for acquiring a three-dimensional model of a human femur; 建立模块,用于根据所述人体股骨三维模型建立对应的内植物三维模型;An establishing module, used for establishing a corresponding three-dimensional model of the implant according to the three-dimensional model of the human femur; 嵌入模块,用于将所述人体股骨三维模型和内植物三维模型进行嵌入处理,获取完整的股骨内植物三维模型;An embedding module, used for embedding the human femur three-dimensional model and the implant three-dimensional model to obtain a complete femur implant three-dimensional model; 处理模块,用于对所述股骨内植物三维模型进行参数设置及优化处理,获取所述股骨内植物三维模型的各个目标节段的等效弹性模量的最优解,并验证所述最优解;A processing module, used for setting and optimizing parameters of the femoral implant three-dimensional model, obtaining an optimal solution of the equivalent elastic modulus of each target segment of the femoral implant three-dimensional model, and verifying the optimal solution; 优化模块,用于根据所述最优解对所述各个目标节段的初始多孔结构进行尺寸优化,获取符合所述最优解的具有弹性模量梯度的多孔结构;其中,所述各个目标节段为所述股骨内植物三维模型需要处理的各个部分;An optimization module, used for optimizing the size of the initial porous structure of each target segment according to the optimal solution, and obtaining a porous structure with an elastic modulus gradient that meets the optimal solution; wherein each target segment is a part of the femoral implant three-dimensional model that needs to be processed; 所述处理模块,具体用于:The processing module is specifically used for: 获取所述股骨内植物三维模型中股骨外表面七个格恩分区的初始应力和七个内植物目标节段对应的等效弹性模量;Obtaining the initial stress of seven Gern zones on the outer surface of the femur and the equivalent elastic modulus corresponding to seven implant target segments in the three-dimensional model of the femoral implant; 将获取的所述初始应力和所述等效弹性模量进行参数化建模,得到参数化模型;Performing parameterized modeling on the acquired initial stress and equivalent elastic modulus to obtain a parameterized model; 将所述参数化模型导入响应面优化模块中,并对所述参数化模型进行试验设计,获取响应面模型;Importing the parameterized model into a response surface optimization module, and performing experimental design on the parameterized model to obtain a response surface model; 基于多目标遗传优化算法对所述响应面模型进行优化设计,获取所述各个目标节段的等效弹性模量最优解;Optimizing 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; 其中,通过获取所述试验设计的样本点数和相关函数确定所述响应面模型;Wherein, the response surface model is determined by obtaining the number of sample points and the correlation function of the experimental design; 所述优化模块,具体用于:The optimization module is specifically used for: 获取所述股骨内植物三维模型中的各个目标节段内植物的初始多孔结构;Acquiring the initial porous structure of each target segment implant in the three-dimensional model of the femoral implant; 将各个目标节段内植物中需要调节的初始多孔结构的尺寸和等效弹性模量参数化;parameterizing the size and equivalent elastic modulus of the initial porous structure to be adjusted in the plant within each target segment; 基于序列二次规划法对参数化处理后的所述初始多孔结构进行尺寸优化,获取符合所述最优解的具有弹性模量梯度的多孔结构。The size of the initial porous structure after parameterization is optimized based on a sequential quadratic programming method to obtain a porous structure with an elastic modulus gradient that meets the optimal solution. 7.一种电子设备,其特征在于,包括:处理器、存储器和总线,其中,7. An electronic device, comprising: a processor, a memory and a bus, wherein: 所述处理器和所述存储器通过所述总线完成相互间的通信;The processor and the memory communicate with each other via the bus; 所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行如权利要求1至5中任一所述的方法。The memory stores program instructions that can be executed by the processor, and the processor can execute any one of the methods according to claims 1 to 5 by calling the program instructions. 8.一种非暂态计算机可读存储介质,其特征在于,所述非暂态计算机可读存储介质存储计算机指令,所述计算机指令使所述计算机执行如权利要求1至5中任一所述的方法。8. A non-transitory computer-readable storage medium, characterized in that the non-transitory computer-readable storage medium stores computer instructions, and the computer instructions enable the computer to execute any one of the methods according to claims 1 to 5.
CN202110406442.6A 2021-04-15 2021-04-15 A method, device, electronic device and storage medium for optimizing orthopedic implants Active CN113284569B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110406442.6A CN113284569B (en) 2021-04-15 2021-04-15 A method, device, electronic device and storage medium for optimizing orthopedic implants

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110406442.6A CN113284569B (en) 2021-04-15 2021-04-15 A method, device, electronic device and storage medium for optimizing orthopedic implants

Publications (2)

Publication Number Publication Date
CN113284569A CN113284569A (en) 2021-08-20
CN113284569B true CN113284569B (en) 2024-11-26

Family

ID=77276749

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110406442.6A Active CN113284569B (en) 2021-04-15 2021-04-15 A method, device, electronic device and storage medium for optimizing orthopedic implants

Country Status (1)

Country Link
CN (1) CN113284569B (en)

Citations (2)

* Cited by examiner, † Cited by third party
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

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107391784B (en) * 2017-06-03 2020-08-14 大连理工大学 Cancellous bone porous structure modeling method based on topology 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 吉林大学 A design method of a femoral stem prosthesis
CN110866354B (en) * 2019-11-08 2021-08-20 大连理工大学 Structure optimization design method of polymer vascular stent considering scale effect
CN111967172B (en) * 2020-07-06 2024-05-17 华南理工大学 Internal high-pressure forming process optimization design method based on kriging model

Patent Citations (2)

* Cited by examiner, † Cited by third party
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

Also Published As

Publication number Publication date
CN113284569A (en) 2021-08-20

Similar Documents

Publication Publication Date Title
San Cheong et al. Novel adaptive finite element algorithms to predict bone ingrowth in additive manufactured porous implants
CN105997306B (en) A kind of bone, which implants, fills the design method of perforated grill structure
US20190240029A1 (en) Method for adjusting mechanical properties of implant and patient specific surgical implants
US20160256279A1 (en) Patient-Specific Implant for Bone Defects and Methods for Designing and Fabricating Such Implants
Helgason et al. Risk of failure during gait for direct skeletal attachment of a femoral prosthesis: a finite element study
Jafari et al. A comparative study of bone remodeling around hydroxyapatite-coated and novel radial functionally graded dental implants using finite element simulation
CN112233801B (en) Mathematical model construction method for topology optimization of endoprosthesis and topology optimization design method
CN108577958A (en) A Design Method of Personalized 3D Printed Calcaneal Plate Based on Topology Optimization
CN110226991A (en) Method for producing the orthopaedic implant of customization
CN104933263A (en) Serial bone plate design method based on average thighbone model
Van Cleynenbreugel et al. Trabecular bone scaffolding using a biomimetic approach
CN113284569B (en) A method, device, electronic device and storage medium for optimizing orthopedic implants
KR20190100837A (en) An elastic modulus adjustment method for implants having lattice scaffolds structure and patient specific surgical implants which use the method
Virulsri et al. Femoral hip prosthesis design for Thais using multi-objective shape optimization
Karuppudaiyan et al. Finite element analysis of scaffold for large defect in femur bone
CN112220588B (en) A method and system for generating a controllable gradient bone tissue engineering scaffold
González-Carbonell et al. Bone remodeling simulation of subject-specific model of tibia under torque
US20220151789A1 (en) Method for adjusting mechanical properties of implant and patient specific surgical implants
Sadollah et al. Optimum functionally gradient materials for dental implant using simulated annealing
Gkaintes et al. Computational modelling of human femur after total hip arthroplasty
Żmudzki et al. Influence of model discretization density in FEM numerical analysis on the determined stress level in bone surrounding dental implants
Miodowska et al. Callus remodelling model
Ghias et al. Reliability-based design algorithm for artificially replaced hip prosthesis considering material property uncertainty
Slowinski Procedure of generating the individually matched bone scaffolds.
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
GR01 Patent grant
GR01 Patent grant