CN109766599A - A kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle - Google Patents
A kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle Download PDFInfo
- Publication number
- CN109766599A CN109766599A CN201811597044.1A CN201811597044A CN109766599A CN 109766599 A CN109766599 A CN 109766599A CN 201811597044 A CN201811597044 A CN 201811597044A CN 109766599 A CN109766599 A CN 109766599A
- Authority
- CN
- China
- Prior art keywords
- fusion device
- model
- centrum
- bone
- anterior intervertebral
- 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
- 230000004927 fusion Effects 0.000 title claims abstract description 230
- 210000000988 bone and bone Anatomy 0.000 title claims abstract description 101
- 238000013461 design Methods 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000009826 distribution Methods 0.000 claims abstract description 25
- 239000000463 material Substances 0.000 claims abstract description 22
- 238000002513 implantation Methods 0.000 claims abstract description 15
- 230000000694 effects Effects 0.000 claims abstract description 13
- 230000037182 bone density Effects 0.000 claims description 21
- 238000005457 optimization Methods 0.000 claims description 19
- 238000011156 evaluation Methods 0.000 claims description 10
- 238000002591 computed tomography Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 8
- 210000003041 ligament Anatomy 0.000 claims description 5
- 230000008676 import Effects 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 3
- 239000004744 fabric Substances 0.000 claims 1
- 239000000155 melt Substances 0.000 claims 1
- 238000004519 manufacturing process Methods 0.000 description 4
- 230000007774 longterm Effects 0.000 description 3
- 230000002980 postoperative effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 238000001356 surgical procedure Methods 0.000 description 2
- 239000004696 Poly ether ether ketone Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- JUPQTSLXMOCDHR-UHFFFAOYSA-N benzene-1,4-diol;bis(4-fluorophenyl)methanone Chemical compound OC1=CC=C(O)C=C1.C1=CC(F)=CC=C1C(=O)C1=CC=C(F)C=C1 JUPQTSLXMOCDHR-UHFFFAOYSA-N 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000001054 cortical effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- ZZUFCTLCJUWOSV-UHFFFAOYSA-N furosemide Chemical group C1=C(Cl)C(S(=O)(=O)N)=CC(C(O)=O)=C1NCC1=CC=CO1 ZZUFCTLCJUWOSV-UHFFFAOYSA-N 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 229920002530 polyetherether ketone Polymers 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
Landscapes
- Prostheses (AREA)
Abstract
The personalized anterior intervertebral fusion device design method based on bone remoulding principle that the invention discloses a kind of, comprising the following steps: A, centrum Image Acquisition and three-dimensional reconstruction obtain three-dimensional vertebral model；B, the initial model of anterior intervertebral fusion device is established；C, centrum-original fusion device model is established；D, centrum-original fusion device finite element model is established；E, optimize the Density Distribution of the original fusion device finite element model in centrum-original fusion device finite element model；The material properties of anterior intervertebral fusion device initial model, the Density Distribution cloud atlas of the original fusion device finite element model after being optimized are updated using bone remoulding method；F, it is formed, obtains anterior intervertebral fusion device model；G, the effect after estimating anterior intervertebral fusion device model implantation centrum.The present invention not only has the characteristics that can be improved fusion device intensity also have the characteristics that be convenient for and vertebral fusion.
Description
Technical field
The present invention relates to a kind of design technology of anterior intervertebral fusion device, especially a kind of individual character based on bone remoulding principle
Change anterior intervertebral fusion device design method.
Background technique
Waist interbody fusion is the common method for treating backbone retrogression pathological changes, for preventing the activity of ailing segment.?
In art, interverbebral disc is removed, instead implantation fusion device and bone graft.Invasive lumbar fusion device rises in waist interbody fusion
To recovery vertebra gap height, stablize the effect of backbone, and then provides suitable mechanical environment for the fusion of adjacent vertebral.It is existing to melt
Clutch is typically all that middle section is hole up and down, and to fill bone graft, periphery is entity part, to provide mechanics
Support.Currently, fusion device is mostly standardized fusion device on the market, due to the personalized difference of patient, standardized fusion device
It is applicable that not to be for all patients be all, often will appear some complication, serious that Fusion failure can occur.It is giving
After determining fusion device size, if the contact area of fusion device and terminal plate of vertebral body is too small, soleplate-fusion device interface will be caused to answer
Power is excessively high and causes the sedimentation of fusion device or the intensity of fusion device is caused not enough to be damaged, and this requires real in design
Body portion volume is as big as possible, to provide enough mechanical supports；On the contrary, if the contact area of fusion device and terminal plate of vertebral body
It is excessive, filling bone graft thus is left for without enough spaces, is also unfavorable for merging, and this requires the entity parts in design
Volume is as small as possible, to fill bone graft as much as possible.Therefore, being badly in need of one kind being capable of the above-mentioned contradictory fusion device of reasonable balance
Design method.
Summary of the invention
The object of the present invention is to provide a kind of personalized anterior intervertebral fusion device design side based on bone remoulding principle
Method.The present invention not only has the characteristics that can be improved fusion device intensity also have the characteristics that be convenient for and vertebral fusion.
Technical solution of the present invention: a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle, packet
Include following steps:
A, centrum Image Acquisition and three-dimensional reconstruction obtain three-dimensional vertebral model；
B, it establishes the initial model of anterior intervertebral fusion device: geometrical characteristic measurement being carried out to three-dimensional vertebral model, according to survey
The geometrical characteristic data measured are established to obtain anterior intervertebral fusion device initial model；
C, it establishes centrum-original fusion device model: anterior intervertebral fusion device initial model is implanted into three-dimensional vertebral model
Corresponding position obtains centrum-original fusion device model；
D, centrum-original fusion device finite element model is established: by finite element software in centrum-original fusion device model
Three-dimensional vertebral model and anterior intervertebral fusion device initial model carry out FEM meshing, be arranged centrum-original fusion device
The material properties and boundary condition of model, and corresponding Mechanical loading is applied to centrum-original fusion device model, generation contains vertebra
The centrum of body finite element model and original fusion device finite element model-original fusion device finite element model；
E, optimize the Density Distribution of the original fusion device finite element model in centrum-original fusion device finite element model: adopting
The material properties that anterior intervertebral fusion device initial model is updated with bone remoulding method, the original fusion device finite element after being optimized
The Density Distribution cloud atlas of model；
F, it is formed, obtains anterior intervertebral fusion device model: to the close of the original fusion device finite element model after adjusting and optimizing
The density value spent in cloud charts carries out binarization operation, and Density Distribution cloud atlas upper density value is close lower than specific threshold region
Angle value is set as 0, and the density value of remainder on Density Distribution cloud atlas is disposed as maximum value ρmax, it is limited to obtain fusion device
Meta-model；The position that density value is 0 in fusion device finite element model is emptied, the remainder of fusion device finite element model is retained,
Anterior intervertebral fusion device finite element mesh model is obtained, anterior intervertebral fusion device finite element mesh model is converted to preceding road vertebra
Between fusion device model；
G, anterior intervertebral fusion device model the effect after estimating anterior intervertebral fusion device model implantation centrum: is implanted to three
Between the upper and lower centrum for tieing up vertebral model, centrum-optimization fusion device finite element model is established, bone remoulding emulation is carried out, before estimating
Road Invasive lumbar fusion device model is implanted into the effect after centrum.
In a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle above-mentioned, the step A's
Specific steps are as follows: patient's lumbar vertebrae is subjected to CT scan, the vertebra initial data for then obtaining CT scan imports at medical image
It manages in software, the vertebra being partitioned into every CT film is reconstructed using volume drawing mode and needs to be implanted into anterior intervertebral fusion device
The three-dimensional structure of upper and lower two sections centrum, obtains three-dimensional vertebral model.
In a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle above-mentioned, in the step B
Geometrical characteristic data include centrum sagittal diameter, the horizontal shape diameter of centrum and disc height；The initial model of anterior intervertebral fusion device
For disc-shape, the horizontal shape diameter of centrum that the width of disc-shape is 2/3~3/4, the vertebra that the length of disc-shape is 2/3~3/4
Body sagittal diameter, the height of disc-shape are disc height.
It is limited in step D in a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle above-mentioned
First grid dividing method particularly includes: in finite element software to three-dimensional vertebral model and anterior intervertebral fusion device initial model into
Row FEM meshing, centrum grid cell use tetrahedron element or hexahedral element, anterior intervertebral fusion device introductory die
Type grid cell uses hexahedral element.
In a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle above-mentioned, bone weight in step E
The specific steps of construction method are as follows:
E1, the bone remoulding quantity of stimulus U/ ρ that anterior intervertebral fusion device initial model is calculated, the i.e. strain energy of unit mass
Density, and then judge bone remoulding quantity of stimulus U/ ρ and bone remoulding stimulus threshold k relationship；
E2, according to the relationship between the quantity of stimulus U/ ρ judged and bone remoulding stimulus threshold k, according to bone remoulding controlling party
The numerical value of Cheng Jinhang bone density calculates, and calculating is then terminated when bone density value reaches convergence state, is not restrained, preceding road vertebra is updated
Between fusion device initial model material properties continue calculate until bone density reach convergence；
The bone remoulding governing equation is expressed with piecewise function form,
Wherein, U is strain energy density；ρ is bone apparent density；B is bone remoulding rate；K is bone remoulding stimulus threshold；S is
The width of inertia area；T is incremental time.
In a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle above-mentioned, anterior intervertebral fusion
Formula in the update of the material properties of the initial model of device, between bone apparent density ρ and elastic modulus E are as follows: E=3790 ρ3；
Wherein, E unit is MPa, and the unit of ρ is g/cm3。
In a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle above-mentioned, the ρmaxIt indicates
The maximum value of bone apparent density ρ, the density value of selected material when being fusion device processing and fabricating.
In a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle above-mentioned, the centrum-is first
Beginning fusion device model includes three-dimensional vertebral model and anterior intervertebral fusion device initial model；The three-dimensional vertebral model includes upper
Centrum and the next centrum, are equipped with anterior intervertebral fusion device initial model between upper centrum and the next centrum；The upper centrum
Ligament is additionally provided between the next centrum.
In a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle above-mentioned, in the step G
Estimate the bone density and bone graft that evaluation index selected by effect after anterior intervertebral fusion device model implantation centrum is centrum
Bone density.
Compared with prior art, the present invention establishes three-dimensional vertebral model by the centrum image of acquisition patient, and as
According to the initial model for establishing anterior intervertebral fusion device, so as to carry out personal settings according to different patients, so as to
The effective syncretizing effect for improving anterior intervertebral fusion device；By establishing centrum-original fusion device finite element model, and to vertebra
The Density Distribution of original fusion device finite element model in body-original fusion device finite element model carries out reasonably optimizing, is based on bone
Method for reconstructing optimizes the Density Distribution of original fusion device finite element model, the anterior intervertebral fusion device enable with
Least material obtains maximum intensity, reaches optimal topological structure, stable structure, intensity is good, long service life.To sum up institute
It states, the present invention not only has the characteristics that can be improved fusion device intensity, and has the characteristics that be convenient for and vertebral fusion.
Detailed description of the invention
Fig. 1 is the top view of the next centrum；
Fig. 2 is the topology view of three-dimensional vertebral model；
Fig. 3 is the structural schematic diagram of anterior intervertebral fusion device initial model；
Fig. 4 is centrum-original fusion device finite element model structural schematic diagram；
Fig. 5 is the Density Distribution cloud atlas of the original fusion device finite element model after optimization；
Fig. 6 is anterior intervertebral fusion device model.
The label in accompanying drawing is: the upper centrum of 1-, the bottom 2- centrum, 3- anterior intervertebral fusion device initial model, 4- ligament.
Specific embodiment
The present invention is further illustrated with reference to the accompanying drawings and examples, but be not intended as to the present invention limit according to
According to.
Embodiment 1.A kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle is constituted such as Fig. 1 extremely
Shown in Fig. 6, comprising the following steps:
A, centrum Image Acquisition and three-dimensional reconstruction obtain three-dimensional vertebral model；
B, it establishes the initial model of anterior intervertebral fusion device: geometrical characteristic measurement being carried out to three-dimensional vertebral model, according to survey
The geometrical characteristic data measured are established to obtain anterior intervertebral fusion device initial model；
C, it establishes centrum-original fusion device model: anterior intervertebral fusion device initial model is implanted into three-dimensional vertebral model
Corresponding position obtains centrum-original fusion device model；
D, centrum-original fusion device finite element model is established: by finite element software in centrum-original fusion device model
Three-dimensional vertebral model and anterior intervertebral fusion device initial model carry out FEM meshing, be arranged centrum-original fusion device
The material properties and boundary condition of model, and corresponding Mechanical loading is applied to centrum-original fusion device model, generation contains vertebra
The centrum of body finite element model and original fusion device finite element model-original fusion device finite element model；
E, optimize the Density Distribution of the original fusion device finite element model in centrum-original fusion device finite element model: adopting
The material properties that anterior intervertebral fusion device initial model is updated with bone remoulding method, the original fusion device finite element after being optimized
The Density Distribution cloud atlas of model；
F, it is formed, obtains anterior intervertebral fusion device model: to the close of the original fusion device finite element model after adjusting and optimizing
The density value spent in cloud charts carries out binarization operation, and Density Distribution cloud atlas upper density value is close lower than specific threshold region
Angle value is set as 0, and the density value of remainder on Density Distribution cloud atlas is disposed as maximum value ρmax, it is limited to obtain fusion device
Meta-model；The position that density value is 0 in fusion device finite element model is emptied, the remainder of fusion device finite element model is retained,
Anterior intervertebral fusion device finite element mesh model is obtained, anterior intervertebral fusion device finite element mesh model is converted to preceding road vertebra
Between fusion device model；
Specific threshold is chosen between ρ min and ρ max.
G, anterior intervertebral fusion device model the effect after estimating anterior intervertebral fusion device model implantation centrum: is implanted to three
Between the upper and lower centrum for tieing up vertebral model, centrum-optimization fusion device finite element model is established, bone remoulding emulation is carried out, before estimating
Road Invasive lumbar fusion device model is implanted into the effect after centrum.
The specific steps of the step A are as follows: patient's lumbar vertebrae is subjected to CT scan, the vertebra for then obtaining CT scan is original
Data import in Medical Image Processing software, and the vertebra being partitioned into every CT film is reconstructed using volume drawing mode and needs to plant
The three-dimensional structure for entering the upper and lower two sections centrum of anterior intervertebral fusion device, obtains three-dimensional vertebral model.
Geometrical characteristic data in the step B include centrum sagittal diameter, the horizontal shape diameter of centrum and disc height；Preceding road vertebra
Between the initial model of fusion device be disc-shape, the horizontal shape diameter of centrum that the width of disc-shape is 2/3~3/4, disc-shape
The centrum sagittal diameter that length is 2/3~3/4, the height of disc-shape are disc height.
In computer aided design software, the lower soleplate of the upper centrum reconstructed is cut off, the next centrum it is upper
Soleplate excision, to simulate the operation of the excision soleplate in spinal fusion surgery.Specific way is to cut the lower surface of upper centrum
Flat, the upper surface of the next centrum is cut flat with, that is, thinks to complete soleplate excision.On the next centrum upper surface for completing excision operation
Find the anterior-posterior axis for dividing vertebral model equally into left and right two halves, measure anterior-posterior axis between two intersection point of vertebral model edge at a distance from,
The distance is centrum sagittal diameter.It measures at a distance between two intersection point of straight line and vertebral model edge perpendicular with anterior-posterior axis,
It is middle apart from it is maximum be the horizontal shape diameter of centrum.Measure lower surface and the interplanar of upper surface two of the next centrum of upper centrum
Distance, as disc height.
FEM meshing in step D method particularly includes: to three-dimensional vertebral model and preceding road in finite element software
Invasive lumbar fusion device initial model carries out FEM meshing, and centrum grid cell uses tetrahedron element or hexahedral element,
Anterior intervertebral fusion device initial model grid cell uses hexahedral element.
The specific steps of bone remoulding method in step E are as follows:
E1, the bone remoulding quantity of stimulus U/ ρ that anterior intervertebral fusion device initial model is calculated, the i.e. strain energy of unit mass
Density, and then judge bone remoulding quantity of stimulus U/ ρ and bone remoulding stimulus threshold k relationship；
E2, according to the relationship between the quantity of stimulus U/ ρ judged and bone remoulding stimulus threshold k, according to bone remoulding controlling party
The numerical value of Cheng Jinhang bone density calculates, and calculating is then terminated when bone density value reaches convergence state, is not restrained, preceding road vertebra is updated
Between fusion device initial model material properties continue calculate until bone density reach convergence；
When calculating, cause to calculate mistake as denominator in order to avoid occurring 0 value in division calculation, setting bone is apparently close
The minimum value of degree ρ is slightly over 0 ρ min.The maximum value that bone apparent density ρ is set selected material when being fusion device processing and fabricating
The density value ρ max of material.When bone remoulding is completed, the density of all units will be in one of following three situation: being located at the inertia phase
Between, reach minimum value ρ min, reach maximum value ρ max, thus be adjusted optimization after original fusion device finite element model it is close
Spend cloud charts.
The bone remoulding governing equation is expressed with piecewise function form,
Wherein, U is strain energy density；ρ is bone apparent density；B is bone remoulding rate；K is bone remoulding stimulus threshold；S is
The width of inertia area；T is incremental time.
In the update of the material properties of the initial model of anterior intervertebral fusion device, between bone apparent density ρ and elastic modulus E
Formula are as follows: E=3790 ρ3；Wherein, E unit is MPa, and the unit of ρ is g/cm3。
The ρmaxIndicate the maximum value of bone apparent density ρ, the density value of selected material when being fusion device processing and fabricating.
The centrum-original fusion device model includes three-dimensional vertebral model and anterior intervertebral fusion device initial model；It is described
Three-dimensional vertebral model includes upper centrum 1 and the next centrum 2, is equipped with anterior intervertebral fusion between upper centrum 1 and the next centrum 2
Device initial model 3；Ligament 4 is additionally provided between the upper centrum 1 and the next centrum 2.
Evaluation index selected by effect is centrum after estimating anterior intervertebral fusion device model implantation centrum in the step G
Bone density and bone graft bone density.
By adjusting specific threshold, different fusion device designs can be obtained.Threshold range is between ρ min and ρ max.
Between the fusion device model of obtained different designs is implanted to upper and lower centrum respectively, centrum-optimization is established respectively
Fusion device finite element model carries out implantation long-term postoperative evaluation.In evaluation, it is imitative that bone remoulding is carried out to centrum and bone graft
Very.Assuming that bone remoulding only occurs in vertebral cancellous bone, centrum trailing unit and bone graft region, without in vertebral body cortical bone
Region.After the completion of reconstruction, the bone density of centrum and the bone density of bone graft are chosen as evaluation index, analyzes melting for different designs
Bone density situation of change at a specified future date after clutch implantation, so that optimal fusion device be selected to design.The personalized intervertebral of final design
Fusion device is used for subsequent processing and manufacturing with the output of STL format.
Embodiment 2.A kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle is constituted such as Fig. 1 extremely
Shown in Fig. 6, comprising the following steps:
Step 101: centrum Image Acquisition and three-dimensional reconstruction
Patient's lumbar vertebrae is subjected to CT scan, setting thickness is 0.63mm, and resolution ratio is 0.39mm × 0.39mm, every tomographic image
Size be 512 × 512 pixels, the initial data for the vertebra for then obtaining CT scan imports the Medical Image Processings such as Mimics
In software, the vertebra in every CT film is gone out using manual automatic Segmentation, is finally reconstructed using volume drawing mode and needs to plant
Enter the three-dimensional structure of the upper and lower two sections centrum of fusion device.In this example, reconstruction is two section centrum of L3 and L4.
Step 102: the initial model of anterior intervertebral fusion device is established
In computer aided design software Solidworks, the lower soleplate of the upper centrum L3 reconstructed is cut off,
The upper soleplate of the next centrum L4 is cut off, to simulate the operation of the excision soleplate in spinal fusion surgery.Specific way is by L3 vertebra
The lower surface of body is cut flat with, and the upper surface of the next centrum L4 is cut flat with, that is, thinks to complete soleplate excision.In the case where completing excision operation
Position measures centrum sagittal diameter AB and the horizontal shape diameter CD value (Fig. 1) of centrum on the upper surface centrum L4.Measure the lower surface of upper centrum L3
With the distance EF (Fig. 2) of the interplanar of upper surface two of the next centrum L4, as disc height.
As shown in figure 3, shape is discoid, width for the initial model of the anterior intervertebral fusion device in the embodiment
Degree is in 2/3~3/4 range of the horizontal shape diameter CD of centrum, the length is in 2/3~3/4 range of centrum sagittal diameter AB, height
For disc height EF.Since Invasive lumbar fusion device initial model size is larger, can not be implanted into from posterior approach, it must be from preceding road hand
Art implantation, therefore, the fusion device to be designed are anterior intervertebral fusion device.
Step 103: the foundation of centrum-fusion device finite element model
Finite element net is carried out to upper and lower centrum and anterior intervertebral fusion device initial model in finite element software Abaqus
Lattice divide, and centrum grid cell uses tetrahedron element, and anterior intervertebral fusion device initial model grid cell uses hexahedron list
Member.After completing FEM meshing, corresponding material properties and boundary condition are set, corresponding Mechanical loading is applied, completes
The foundation (Fig. 4) of centrum-original fusion device finite element model.The finite element model by anterior intervertebral fusion device initial model, on
Position centrum, the next centrum and ligament composition.
Step 104: the anterior intervertebral fusion device Density Distribution optimization based on bone remoulding principle
Regard anterior intervertebral fusion device to be optimized as one piece of " bone " for deferring to bone remoulding rule, the mistake of topological optimization
Journey regards bone remoulding process as, and the attribute of " bone " material is updated using bone remoulding method, until balance is optimal close to obtain
Degree distribution.
Each parameter setting is as follows in bone remoulding governing equation: B=1.00 (g/cm3)2/ (MPa unit time), k=
0.004J/g, s=10%.The elasticity modulus that initial anterior intervertebral fusion device (" bone ") is arranged is 100MPa, passes through formula meter
Calculate initial apparent density is 0.2977g/cm3.The minimum value that apparent density is arranged is 0.01g/cm3, maximum value is set as springform
Amount corresponding density 0.9830g/cm when being 3600MPa (PEEK elasticity modulus of materials)3.When bone remoulding is completed, all units
Density will be one of in following three situation: positioned at during inertia, reach minimum value 0.01g/cm3, reach maximum value
0.9830g/cm3, to obtain the Density Distribution cloud atlas (Fig. 5) after fusion device optimization.
Step 105: anterior intervertebral fusion device design typification
By step 101-104, the optimization of anterior intervertebral fusion device Density Distribution, each portion of Invasive lumbar fusion device at this time are completed
The density divided is all inconsistent, and such fusion device processing and manufacturing is got up very difficult.Therefore, one is carried out to the density of fusion device
Binarization operation, on the Density Distribution cloud atlas after fusion device optimization, the density by density lower than the position of specific threshold is arranged
Being 0, emptying this position in design, the density of remainder is both configured to maximum value by implantable graft bone when for performing the operation
ρmax, retain the part in design.By adjusting the specific threshold, different fusion device designs can be obtained.Binarization operation
Fusion device mould shapes afterwards are possible and irregular, need to carry out appropriate adjustment, so that its shape is become more regular, to process
Manufacture.As shown in figure 5, for the Density Distribution cloud atlas after the fusion device optimization of the embodiment.As can be seen from Figure 5, density is past from centre
It is gradually increased outside, and Density Distribution is roughly divisible into high, medium and low three parts.Based on the distribution map, using the geometry of rule
Carry out centering, low density area is described, by low density area or in+low density area empties, to obtain two for filling bone graft
Different optimization fusion device models: the fusion device model S of optimization, optimization fusion device model L.
Step 106: the fusion device implantation long-term postoperative evaluation based on bone remoulding emulation
By the fusion device models of different designs obtained in the previous step (the fusion device model S of optimization, optimization fusion device model
L it) is implanted between upper and lower centrum respectively, establishes centrum-optimization fusion device S, centrum-optimization fusion device L finite element mould respectively
Type carries out implantation long-term postoperative evaluation.In evaluation, bone remoulding emulation is carried out to L3, L4 centrum and bone graft.It chooses
The bone density of upper and lower centrum and the bone density of bone graft are as evaluation index, after the fusion device implantation for analyzing different designs
Bone density situation of change at a specified future date, so that optimal fusion device be selected to design.The personalized Invasive lumbar fusion device of final design is with STL lattice
Formula output, is used for subsequent processing and manufacturing.
Claims (9)
1. a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle, which is characterized in that including following step
It is rapid:
A, centrum Image Acquisition and three-dimensional reconstruction obtain three-dimensional vertebral model；
B, it establishes the initial model of anterior intervertebral fusion device: geometrical characteristic measurement being carried out to three-dimensional vertebral model, according to measuring
To geometrical characteristic data establish to obtain anterior intervertebral fusion device initial model；
C, it establishes centrum-original fusion device model: anterior intervertebral fusion device initial model is implanted into the corresponding of three-dimensional vertebral model
Position obtains centrum-original fusion device model；
D, centrum-original fusion device finite element model is established: by finite element software to three in centrum-original fusion device model
It ties up vertebral model and anterior intervertebral fusion device initial model carries out FEM meshing, centrum-original fusion device model is set
Material properties and boundary condition, and corresponding Mechanical loading is applied to centrum-original fusion device model, generation has containing centrum
Limit centrum-original fusion device finite element model of meta-model and original fusion device finite element model；
E, optimize the Density Distribution of the original fusion device finite element model in centrum-original fusion device finite element model: using bone
Method for reconstructing updates the material properties of anterior intervertebral fusion device initial model, the original fusion device finite element model after being optimized
Density Distribution cloud atlas；
F, it is formed, obtains anterior intervertebral fusion device model: to the density point of the original fusion device finite element model after adjusting and optimizing
Density value on cloth cloud atlas carries out binarization operation, and Density Distribution cloud atlas upper density value is lower than to the density value in specific threshold region
It is set as 0, the density value of remainder on Density Distribution cloud atlas is disposed as maximum value ρmax, obtain fusion device finite element mould
Type；The position that density value is 0 in fusion device finite element model is emptied, retains the remainder of fusion device finite element model, obtains
Anterior intervertebral fusion device finite element mesh model is converted to preceding road intervertebral and melted by anterior intervertebral fusion device finite element mesh model
Clutch model；
G, anterior intervertebral fusion device model the effect after estimating anterior intervertebral fusion device model implantation centrum: is implanted to three-dimensional vertebra
Between the upper and lower centrum of body Model, centrum-optimization fusion device finite element model is established, bone remoulding emulation is carried out, estimates preceding road vertebra
Between fusion device model implantation centrum after effect.
2. a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle according to claim 1,
It is characterized in that, the specific steps of the step A are as follows: patient's lumbar vertebrae is subjected to CT scan, the vertebra for then obtaining CT scan is former
Beginning data import in Medical Image Processing software, and the vertebra being partitioned into every CT film reconstructs needs using volume drawing mode
It is implanted into the three-dimensional structure of the upper and lower two sections centrum of anterior intervertebral fusion device, obtains three-dimensional vertebral model.
3. a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle according to claim 1,
Be characterized in that: the geometrical characteristic data in the step B include centrum sagittal diameter, the horizontal shape diameter of centrum and disc height；Preceding road
The initial model of Invasive lumbar fusion device is disc-shape, the horizontal shape diameter of centrum that the width of disc-shape is 2/3~3/4, disc-shape
Length be 2/3~3/4 centrum sagittal diameter, the height of disc-shape is disc height.
4. a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle according to claim 1,
It is characterized in that, FEM meshing in step D method particularly includes: to three-dimensional vertebral model and preceding road in finite element software
Invasive lumbar fusion device initial model carries out FEM meshing, and centrum grid cell uses tetrahedron element or hexahedral element,
Anterior intervertebral fusion device initial model grid cell uses hexahedral element.
5. a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle according to claim 1,
It is characterized in that, the specific steps of bone remoulding method in step E are as follows:
E1, the bone remoulding quantity of stimulus U/ ρ that anterior intervertebral fusion device initial model is calculated, the i.e. strain energy of unit mass are close
Degree, and then judge bone remoulding quantity of stimulus U/ ρ and bone remoulding stimulus threshold k relationship；
E2, according to the relationship between the quantity of stimulus U/ ρ judged and bone remoulding stimulus threshold k, according to bone remoulding governing equation into
The numerical value of row bone density calculates, and calculating is then terminated when bone density value reaches convergence state, is not restrained, preceding road intervertebral is updated and melts
The material properties of clutch initial model continue to calculate until bone density reaches convergence；
The bone remoulding governing equation is expressed with piecewise function form,
Wherein, U is strain energy density；ρ is bone apparent density；B is bone remoulding rate；K is bone remoulding stimulus threshold；S is inertia
The width in region；T is incremental time.
6. a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle according to claim 1,
Be characterized in that, in the update of the material properties of the initial model of anterior intervertebral fusion device, bone apparent density ρ and elastic modulus E it
Between formula are as follows: E=3790 ρ3；Wherein, E unit is MPa, and the unit of ρ is g/cm3。
7. a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle according to claim 6,
It is characterized in that, the ρmaxIndicate the maximum value of bone apparent density ρ, the density value of selected material when being fusion device processing and fabricating.
8. a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle according to claim 1,
Be characterized in that: the centrum-original fusion device model includes three-dimensional vertebral model and anterior intervertebral fusion device initial model；It is described
Three-dimensional vertebral model includes upper centrum (1) and the next centrum (2), is equipped with preceding road between upper centrum (1) and the next centrum (2)
Invasive lumbar fusion device initial model (3)；Ligament (4) are additionally provided between the upper centrum (1) and the next centrum (2).
9. a kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle according to claim 1,
Be characterized in that: evaluation index selected by effect is vertebra after estimating anterior intervertebral fusion device model implantation centrum in the step G
The bone density of body and the bone density of bone graft.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811597044.1A CN109766599A (en) | 2018-12-26 | 2018-12-26 | A kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811597044.1A CN109766599A (en) | 2018-12-26 | 2018-12-26 | A kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109766599A true CN109766599A (en) | 2019-05-17 |
Family
ID=66450310
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811597044.1A Pending CN109766599A (en) | 2018-12-26 | 2018-12-26 | A kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109766599A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110368147A (en) * | 2019-08-05 | 2019-10-25 | 北京爱康宜诚医疗器材有限公司 | Invasive lumbar fusion device and preparation method thereof |
| CN110384572A (en) * | 2019-07-16 | 2019-10-29 | 天津市天津医院 | CAD personalization lumbar intervertebral fusion device and its designing and manufacturing method |
| CN110464514A (en) * | 2019-07-16 | 2019-11-19 | 太原理工大学 | The preparation method of application and implantation material of the honeycomb sandwich structure in vertebral implant |
| CN110495951A (en) * | 2019-07-31 | 2019-11-26 | 同济大学 | A kind of production method of the oblique lateral approach reset model of lumbar vertebra |
| CN110929379A (en) * | 2019-10-25 | 2020-03-27 | 西安交通大学 | Topology optimization artificial vertebral body and design method thereof |
| CN111317868A (en) * | 2020-02-19 | 2020-06-23 | 西安交通大学 | Directionally-enhanced 3D-printed polymer proximal femur substitute and design method thereof |
| CN113408172A (en) * | 2021-06-28 | 2021-09-17 | 大连理工大学 | Intervertebral fusion device manufacturing method and manufacturing system |
| CN117421940A (en) * | 2023-12-19 | 2024-01-19 | 山东交通学院 | Global mapping method and device between digital twin lightweight model and physical entity |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060276925A1 (en) * | 2003-04-23 | 2006-12-07 | The Regents Of The University Of Michigan | Integrated global layout and local microstructure topology optimization approach for spinal cage design and fabrication |
| CN104166775A (en) * | 2014-08-27 | 2014-11-26 | 上海理工大学 | Component distribution method for automobile seat framework |
| CN104318009A (en) * | 2014-10-20 | 2015-01-28 | 北京航空航天大学 | Method for designing personalized interbody fusion cage |
| CN104462723A (en) * | 2014-12-25 | 2015-03-25 | 北京航空航天大学 | Personalized interbody fusion cage design method based on topological optimization and bony reconstitution simulation |
| CN107252373A (en) * | 2017-06-14 | 2017-10-17 | 北京航空航天大学 | A kind of personalized concrete dynamic modulus Invasive lumbar fusion device and design method |
| CN107609277A (en) * | 2017-09-15 | 2018-01-19 | 上海凯利泰医疗科技股份有限公司 | The structural topological optimization method of fusion device |
-
2018
- 2018-12-26 CN CN201811597044.1A patent/CN109766599A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060276925A1 (en) * | 2003-04-23 | 2006-12-07 | The Regents Of The University Of Michigan | Integrated global layout and local microstructure topology optimization approach for spinal cage design and fabrication |
| CN104166775A (en) * | 2014-08-27 | 2014-11-26 | 上海理工大学 | Component distribution method for automobile seat framework |
| CN104318009A (en) * | 2014-10-20 | 2015-01-28 | 北京航空航天大学 | Method for designing personalized interbody fusion cage |
| CN104462723A (en) * | 2014-12-25 | 2015-03-25 | 北京航空航天大学 | Personalized interbody fusion cage design method based on topological optimization and bony reconstitution simulation |
| CN107252373A (en) * | 2017-06-14 | 2017-10-17 | 北京航空航天大学 | A kind of personalized concrete dynamic modulus Invasive lumbar fusion device and design method |
| CN107609277A (en) * | 2017-09-15 | 2018-01-19 | 上海凯利泰医疗科技股份有限公司 | The structural topological optimization method of fusion device |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110384572A (en) * | 2019-07-16 | 2019-10-29 | 天津市天津医院 | CAD personalization lumbar intervertebral fusion device and its designing and manufacturing method |
| CN110464514A (en) * | 2019-07-16 | 2019-11-19 | 太原理工大学 | The preparation method of application and implantation material of the honeycomb sandwich structure in vertebral implant |
| CN110464514B (en) * | 2019-07-16 | 2020-09-04 | 太原理工大学 | Application of honeycomb sandwich structure in vertebral body implant and preparation method of implant |
| CN110495951A (en) * | 2019-07-31 | 2019-11-26 | 同济大学 | A kind of production method of the oblique lateral approach reset model of lumbar vertebra |
| CN110368147A (en) * | 2019-08-05 | 2019-10-25 | 北京爱康宜诚医疗器材有限公司 | Invasive lumbar fusion device and preparation method thereof |
| CN110929379A (en) * | 2019-10-25 | 2020-03-27 | 西安交通大学 | Topology optimization artificial vertebral body and design method thereof |
| CN111317868A (en) * | 2020-02-19 | 2020-06-23 | 西安交通大学 | Directionally-enhanced 3D-printed polymer proximal femur substitute and design method thereof |
| CN113408172A (en) * | 2021-06-28 | 2021-09-17 | 大连理工大学 | Intervertebral fusion device manufacturing method and manufacturing system |
| CN113408172B (en) * | 2021-06-28 | 2024-05-14 | 大连理工大学 | Manufacturing method and manufacturing system of interbody fusion cage |
| CN117421940A (en) * | 2023-12-19 | 2024-01-19 | 山东交通学院 | Global mapping method and device between digital twin lightweight model and physical entity |
| CN117421940B (en) * | 2023-12-19 | 2024-03-19 | 山东交通学院 | Global mapping method and device between digital twin lightweight model and physical entity |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109766599A (en) | A kind of personalized anterior intervertebral fusion device design method based on bone remoulding principle | |
| CN104318009B (en) | A kind of personalized Invasive lumbar fusion device design method | |
| CN105877875A (en) | Personalized thyroid cartilage prosthesis and production method thereof | |
| US9250620B2 (en) | 3D design and fabrication system for implants | |
| CN105997306B (en) | A kind of bone, which implants, fills the design method of perforated grill structure | |
| CN109009572B (en) | Hip prosthesis stem and method of making same | |
| CN104462723A (en) | Personalized interbody fusion cage design method based on topological optimization and bony reconstitution simulation | |
| CN107411855A (en) | Atlas and axis fusion of intervertebral joints device | |
| US9138329B2 (en) | Vertebral disc prosthesis | |
| CN105581860A (en) | 3D print type artificial bone trabecula structure cervical fusion device and preparing method and application thereof | |
| CN110929379B (en) | Topology optimization artificial vertebral body and design method thereof | |
| WO2010140036A1 (en) | A method of designing a knee prosthesis | |
| JPWO2021141849A5 (en) | ||
| CN108601550B (en) | System and method for manufacturing personalized fibrous cartilage implants | |
| CN110236741B (en) | Personalized condylar prosthesis design method with topological optimization fixing unit and porous condylar-head unit and personalized condylar prosthesis | |
| CN109101686A (en) | A kind of porous structure design method based on trabecular bone structure form and mechanical property | |
| JP2023536518A (en) | Nasal prosthesis design method for patient-customized nasal prosthesis fabrication | |
| CN205322556U (en) | Backbone anterior cervical spine 3D prints bionical bone trabecula structure cervical vertebra and fuses device | |
| Wan et al. | Three‐Dimensional Biomechanical Finite Element Analysis of Lumbar Disc Herniation in Middle Aged and Elderly | |
| CN103340705B (en) | Personalized lumbar intervertebral disc artificial nucleus prosthesis manufacturing device and method | |
| CN115212014A (en) | Bionic scaffold design method for bone repair based on topological structure optimization | |
| Ni et al. | The bony labyrinth of the early platyrrhine primate Chilecebus | |
| CN113768668A (en) | Modeling method for designing personalized medical mandible model based on TPMS | |
| CN105662660A (en) | 3D printing bionic bone trabecular structure lumbar vertebra fusion device, manufacturing method and application | |
| CN107374785B (en) | A kind of pre- production method of orbital blow-out fracture titanium net |
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 | ||
| CB02 | Change of applicant information |
Address after: 100176 1 ronghua Middle Road, Daxing District economic and Technological Development Zone, Beijing Applicant after: NATIONAL RESEARCH CENTER FOR REHABILITATION TECHNICAL AIDS Address before: 100176 National Rehabilitation Aids Research Center, No. 1, Ronghua Middle Road, Daxing District Economic and Technological Development Zone, Beijing Applicant before: NATIONAL RESEARCH CENTER FOR REHABILITATION TECHNICAL AIDS |
|
| CB02 | Change of applicant information | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190517 |
|
| RJ01 | Rejection of invention patent application after publication |