CN109992820A - The mechanical property evaluating method of porous bone implant based on dodecahedron rod structure - Google Patents

The mechanical property evaluating method of porous bone implant based on dodecahedron rod structure Download PDF

Info

Publication number
CN109992820A
CN109992820A CN201910081513.2A CN201910081513A CN109992820A CN 109992820 A CN109992820 A CN 109992820A CN 201910081513 A CN201910081513 A CN 201910081513A CN 109992820 A CN109992820 A CN 109992820A
Authority
CN
China
Prior art keywords
porous bone
bone implant
elastic modulus
equivalent elastic
porous
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910081513.2A
Other languages
Chinese (zh)
Other versions
CN109992820B (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.)
National Research Center for Rehabilitation Technical Aids
Original Assignee
National Research Center for Rehabilitation Technical Aids
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 National Research Center for Rehabilitation Technical Aids filed Critical National Research Center for Rehabilitation Technical Aids
Priority to CN201910081513.2A priority Critical patent/CN109992820B/en
Publication of CN109992820A publication Critical patent/CN109992820A/en
Application granted granted Critical
Publication of CN109992820B publication Critical patent/CN109992820B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P10/00Technologies related to metal processing
    • Y02P10/25Process efficiency

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Prostheses (AREA)

Abstract

The invention discloses a kind of mechanical property evaluating method of porous bone implant based on dodecahedron rod structure, porous bone implant is spliced to form by least one dodecahedron rod structure basic unit using modeling method;Dodecahedron rod structure basic unit includes the granatohedron being made of 24 girders, and the girder extended outward respectively by opposite two vertex in 12 faces of granatohedron, hole is formed between adjacent girder, dodecahedron rod structure basic unit is formed with multiple holes in three-dimensional spatial distribution;The the first quantitative analysis relationship established between the concrete moduli difference of porous bone implant and girder size predicts the mechanical property of porous bone implant according to the first quantitative analysis relationship;Wherein, concrete moduli difference is the difference of the equivalent elastic modulus and the equivalent elastic modulus for the porous bone implant tested using exemplar of the porous bone implant obtained using simulating analysis.

Description

The mechanical property evaluating method of porous bone implant based on dodecahedron rod structure
Technical field
The present invention relates to bio-medical prosthese technical fields, particularly relate to a kind of porous bone based on dodecahedron rod structure The mechanical property evaluating method of implant.
Background technique
With deep and the traffic accident, the generation of natural calamity of aging society, Bone Defect Repari, the demand of bone displacement are more next It is more, how to provide safe and reliable, and there are excellent mechanical performances, (bone implant includes for the bone implant of regeneration performance Bone bracket and bone prosthese etc., bone prosthese include bone plate, bone nail etc.) have become one of clinical problem urgently to be resolved.Existing gold True body bone implant, is generally made of metal materials such as 316 stainless steels, cobalt-base alloys, titanium or titanium alloy, made of metal bone implant The rigidity of body is much larger than skeleton (3-30GPa of compact bone;0.02-2GPa of cancellous bone), easily cause bone plate and green bone power The problems such as performance mismatches, stress shielding occurs is learned, leads to adjacent bone osteoporosis, the problems such as prosthetic loosening, porous Bionics Bone Implant has rigidity low, and stress shielding is small, is easy to cell adherence growth and regeneration, the features such as personalized customization, has Huge clinical demand and wide application prospect.
Studies have shown that the form and Geometrical Parameter of porous structure are to influence mechanical property for porous bionical bone implant An important factor for energy.Currently, the porous structure of most of bionical bone implant is designed along one direction, structure is single, and performance is poor Not not less;Meanwhile the personalized designs of porous bionical bone implant depend on experience, still lack the scientific design of controllable quantization The problems such as method, there are mechanical property mismatches, and potential risk is big, and regeneration performance is bad greatly limits porous bionical The clinical application and popularization of bone implant;Person again is usually to prepare porous bionical bone implant using 3D printing technique at present, by Power of the size of the limitation of 3D printing technique, the process characteristic successively to add up and metal powder partial size for different porous structures It learns performance and regeneration performance is affected, the bone implant for causing the porous bionical bone implant of design and 3D printing to be prepared The mechanical property difference of body is larger, and bone implant is adapted to the mechanical property of host bone, and permanently effective military service is extremely not Benefit.
Summary of the invention
In view of this, it is an object of the invention to propose a kind of power of porous bone implant based on dodecahedron rod structure Learn performance evaluation methodology, porous bone implant by least one dodecahedron rod structure basic unit using modeling method splice and At several holes with three-dimensional spatial distribution can carry out quantitative analysis and prediction to the mechanical property of porous bone implant.
Based on above-mentioned purpose, the present invention provides a kind of mechanical property of porous bone implant based on dodecahedron rod structure Energy evaluating method, comprising:
The porous bone implant is spliced to form by least one dodecahedron rod structure basic unit using modeling method;
The dodecahedron rod structure basic unit includes: the granatohedron being made of 24 girders, and point The girder not extended outward by opposite two vertex in 12 faces of the granatohedron forms hole between adjacent girder Gap, the dodecahedron rod structure basic unit are formed with multiple holes in three-dimensional spatial distribution;
Establish the first quantitative analysis relationship between the concrete moduli difference and girder size of the porous bone implant, root According to the first quantitative analysis relationship, the mechanical property of the porous bone implant is predicted;
Wherein, the concrete moduli difference is the Equivalent Elasticity mould of the porous bone implant obtained using simulating analysis The difference of amount and the equivalent elastic modulus for the porous bone implant tested using exemplar.
Optionally, the constitutive relation between the concrete moduli difference of the porous bone implant and the girder size is y =25.091x3.1856, wherein x is the girder size, and y is the concrete moduli difference.
Optionally, the Geometrical Parameter of the porous bone implant and the equivalent elastic modulus of the porous bone implant are established Between the second quantitative analysis relationship determined according to the second quantitative analysis relationship described more according to preset Geometrical Parameter The equivalent elastic modulus of hole bone implant;Wherein, the Geometrical Parameter include: the volume of the porous bone implant, porosity, The aperture of surface area, specific surface area, maximum pore.
Optionally, it is obtained based on simulating analysis: the Equivalent Elasticity of the girder size and the porous bone implant Constitutive relation between modulus is y=30.321x2.589;Test to obtain based on exemplar: the girder size and the porous bone are planted Entering the constitutive relation between the equivalent elastic modulus of body is y=6.11x-07263, wherein y is the porous bone implant exemplar Equivalent elastic modulus, x be girder size.
Optionally, it is obtained based on simulating analysis: the porosity of the porous bone implant and the porous bone implant Constitutive relation between the equivalent elastic modulus of body is y=7E+12x-6.15;It tests to obtain based on exemplar: the porous bone implant Constitutive relation between the porosity of body and the equivalent elastic modulus of the porous bone implant is y=-0.0989x+9.5125, Wherein, y is the equivalent elastic modulus of the porous bone implant, and x is the porosity of the porous bone implant.
Optionally, it is obtained based on simulating analysis: the aperture of the maximum pore of the porous bone implant and described more Constitutive relation between the equivalent elastic modulus of hole bone implant is y=39.772x-4.005;It tests to obtain based on exemplar: described Constitutive relation between the aperture of the maximum pore of porous bone implant and the equivalent elastic modulus of the porous bone implant is y =-3.9826x+6.8216, wherein y is the equivalent elastic modulus of the porous bone implant, and x is the porous bone implant Maximum pore aperture.
Optionally, it is obtained based on simulating analysis: the volume of the porous bone implant and the porous bone implant Equivalent elastic modulus between constitutive relation be y=0.0024x1.4493;It tests to obtain based on exemplar: the porous bone implant Constitutive relation between the volume of body and the equivalent elastic modulus of the porous bone implant is y=0.0042x-0.3371, In, y is the equivalent elastic modulus of porous bone implant, and x is the volume of porous bone implant.
Optionally, it is obtained based on simulating analysis: the surface area of the porous bone implant and the porous bone implant Constitutive relation between the equivalent elastic modulus of body is y=5E-11x3.444;It tests to obtain based on exemplar: the porous bone implant Constitutive relation between the surface area of body and the equivalent elastic modulus of the porous bone implant is y=0.002x-2.8805, In, y is the equivalent elastic modulus of the porous bone implant, and x is the surface area of the porous bone implant.
Optionally, it is obtained based on simulating analysis: the specific surface area of the porous bone implant and porous bone implant Equivalent elastic modulus between constitutive relation be y=0.0096x2.4973, wherein y is the equivalent bullet of the porous bone implant Property modulus;It tests to obtain based on exemplar: the equivalent bullet of the specific surface area of the porous bone implant and the porous bone implant Property modulus between constitutive relation be y=0.1626x-1.5883, wherein y be the porous bone implant Equivalent Elasticity mould Amount, x are the specific surface area of the porous bone implant.
Optionally, the section of the girder is circle, and the material of the girder is titanium alloy.
The invention has the advantages that
1, the mechanical property evaluating method of the porous bone implant provided by the invention based on dodecahedron rod structure is established Quantitative analysis relationship between girder size and concrete moduli difference, can predict to produce more according to quantitative analysis relationship The mechanical property of hole bone implant exemplar assesses the mechanical property between the porous bone implant prepared and host bone in advance With degree, is conducive to the clinical adaptation of the porous bone implant of personalization of 3D printing preparation and promotes and applies;
2, the mechanical property evaluating method of the porous bone implant provided by the invention based on dodecahedron rod structure, also builds The Geometrical Parameter of the porous bone implant model under two kinds of measurement methods and the equivalent elastic modulus of porous bone implant are found Quantitative analysis relationship can provide data foundation for the controllable of porous bone implant, quantization, scientific design, be conducive to personalization The clinical application and popularization of bionical bone implant;
3, dodecahedron rod structure basic unit provided by the invention has several holes in three-dimensional spatial distribution, tool Have the mechanical characteristic of isotropic, may advantageously facilitate the conveying and exclusion of nutriment and metabolite, promote cell adherency, Differentiation, breeding and the regeneration of tissue;
4, dodecahedron rod structure basic unit provided by the invention, axisymmetricly structure, there are six Mosaic faces for setting, if Dry dodecahedron rod structure basic unit is by arrangement, the splicing porous bone implant of structure, convenient for setting according to actual needs Count out the diversified porous bone implant of personalization.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with It obtains other drawings based on these drawings.
Fig. 1 is the structural schematic diagram of the dodecahedron rod structure basic unit of the embodiment of the present invention;
Fig. 2 is the mechanotransduction schematic diagram of structure shown in Fig. 1;
Fig. 3 is the top view of structure shown in Fig. 1;
Fig. 4 A, 4B are the porous bone implant of two with different porosities kind that structure is arranged as shown in Figure 1, is spliced Body structural schematic diagram;
Fig. 5 A is the girder size of the embodiment of the present invention and the constitutive relation schematic diagram of porosity;
Fig. 5 B is the constitutive relation schematic diagram of the girder size of the embodiment of the present invention and the aperture of maximum pore;
Fig. 5 C is the girder size of the embodiment of the present invention and the constitutive relation schematic diagram of volume;
Fig. 5 D is the girder size of the embodiment of the present invention and the constitutive relation schematic diagram of surface area;
Fig. 5 E is the girder size of the embodiment of the present invention and the constitutive relation schematic diagram of specific surface area;
Fig. 6 A is the girder size of the embodiment of the present invention and the constitutive relation schematic diagram of equivalent elastic modulus;
Fig. 6 B is the porosity of the embodiment of the present invention and the constitutive relation schematic diagram of equivalent elastic modulus;
Fig. 6 C is the aperture of the embodiment of the present invention and the constitutive relation schematic diagram of equivalent elastic modulus;
Fig. 6 D is the volume of the embodiment of the present invention and the constitutive relation schematic diagram of equivalent elastic modulus;
Fig. 6 E is the surface area of the embodiment of the present invention and the constitutive relation schematic diagram of equivalent elastic modulus;
Fig. 6 F is the specific surface area of the embodiment of the present invention and the constitutive relation schematic diagram of equivalent elastic modulus;
Fig. 7 is the girder size of the embodiment of the present invention and is utilized respectively what finite element analysis was tested with mechanics machine The correspondence diagram of equivalent elastic modulus;
Fig. 8 is the girder size of the embodiment of the present invention and the constitutive relation schematic diagram of concrete moduli difference.
Specific embodiment
To make the objectives, technical solutions, and advantages of the present invention clearer, below in conjunction with specific embodiment, and reference Attached drawing, the present invention is described in more detail.
It should be noted that all statements for using " first " and " second " are for differentiation two in the embodiment of the present invention The non-equal entity of a same names or non-equal parameter, it is seen that " first " " second " only for the convenience of statement, does not answer It is interpreted as the restriction to the embodiment of the present invention, subsequent embodiment no longer illustrates this one by one.
Fig. 1 is the structural schematic diagram of the dodecahedron rod structure basic unit of the embodiment of the present invention, and Fig. 2 is knot shown in Fig. 1 The mechanotransduction schematic diagram of structure, Fig. 3 are the top view of structure shown in Fig. 1.As shown, 12 face provided in an embodiment of the present invention Body rod structural base member 1, including the granatohedron (centre of structure shown in Fig. 1 black being made of 24 girders 10 Line portion), and the girder 10 extended outward respectively by opposite two vertex in 12 faces of granatohedron, it is adjacent small Hole 12 is formed between beam 10, so that entire dodecahedron rod structure unit 1 is formed with several holes in three-dimensional spatial distribution 12。
By granatohedron, any one puts 13 centered on not setting the vertex for extending girder, adjacent with the central point 13 In four extension girders 101,102,103,104, with two opposite extension girders 101 of one pair of them, 103 place straight lines for S axis, Using two opposite extension girders 102 of another pair, 104 place straight lines as Z axis, the endpoint for extending girder 101,104 with adjacent two connects Straight line where the vertical line of line is Y-axis, using straight line where the vertical line of the endpoint line of adjacent two extensions girder 104,103 as X-axis, X, Y, Z, S axis intersect at central point 13, and the angle between adjacent two axis is 45 degree.Dodecahedron rod structure unit 1 respectively with X, Y, Z, S axis is symmetry axis, axisymmetricly structure, so that dodecahedron rod structure unit 1 has isotropic mechanical characteristic.
For the girder extended, its end is cut out to three smooth sections 11, so that entire dodecahedron rod structure The exterior contour of unit 1 is in square shape, and there are six Mosaic faces for the tool of dodecahedron rod structure unit 1 of square shape, respectively The Mosaic face of dodecahedron rod structure unit 1 forms the diversified porous bone implant of personalization by splicing, spread pattern.
Fig. 4 A, 4B are the porous bone implant of two with different porosities kind that structure is arranged as shown in Figure 1, is spliced The schematic diagram of body.As shown, dodecahedron rod structure basic unit 1 can be made when designing specific porous bone implant Entity bone implant is designed to there is difference by modeling methods such as linear array, mirror image, Boolean calculations for basic unit The bone implant of Geometrical Parameter, compared with entity bone implant, it is in three-dimensional space that the porous bone implant of the embodiment of the present invention, which has, Between several holes for being distributed, can both reduce the rigidity of bone implant, avoid stress shielding risk, and can promote cell adherence, promote Regeneration is grown into tissue, and convenient for designing the diversified porous bone implant of personalization according to actual needs.
In the embodiment of the present invention, the girder 10 of dodecahedron rod structure basic unit 1 is constituted, size and cross sectional shape can It is designed according to performance requirement.For example, the section of girder can be square, circle, equilateral triangle, ellipse etc..For The girder that section is square, the side length that girder size is square, section are the girder of equilateral triangle, and girder size is positive three Angular side length, section are circular girder, and for girder having a size of circular diameter, section is the girder of ellipse, girder size For the long axis and short axle of ellipse.
In the embodiment of the present invention, porous bone implant passes through modeling by least one dodecahedron rod structure basic unit 1 Method splicing is built-up, and a dodecahedron rod structure basic unit can be regarded as a small bone implant.Wherein, base In dodecahedron rod structure basic unit 1 constitute porous bone implant Geometrical Parameter include: porous bone implant volume, The aperture of the maximum pore of porous bone implant, the porosity of porous bone implant, the surface area of porous bone implant, porous bone The specific surface area etc. of implant.By adjusting the girder size of dodecahedron rod structure basic unit 1, according to corresponding quantitative point Analysis relationship, each Geometrical Parameter of adjustable dodecahedron rod structure basic unit 1, so as to adjust by dodecahedron rod structure base The Geometrical Parameter for the porous bone implant that this unit 1 designs.Specifically:
In a specific embodiment, the size of dodecahedron rod structure basic unit 1 are as follows: length × width × height be equal to 5 × 5 × 5mm, the material of girder are titanium alloy, and the section of girder is circle, and the size (i.e. circular diameter) that girder is respectively set is φ 0.5mm,φ0.7mm,φ0.9mm,φ1.1mm.Mirror image, array, Boolean calculation are carried out to dodecahedron rod structure basic unit 1 Equal modelling operabilities can be designed arbitrary shape, and porous bone scaffold or prosthese with different Geometrical Parameters.In a kind of implementation In example, it is square that design, which obtains outer profile, and length × width × height is equal to the first porous bone scaffold model of 10 × 10 × 20mm, should First porous bone scaffold model is used to obtain the relationship of girder size and Geometrical Parameter by experiment.In Solidworks or In Rhinoceros or other design softwares, to the volume of the first porous bone scaffold model, surface area, specific surface area, porosity, The Geometrical Parameters such as the aperture of maximum pore are measured and are calculated, and carry out data point using statistics software (such as SPSS software) Analysis, obtains according to data processed result, the Geometrical Parameter of girder size and porous bone scaffold there are apparent mathematics constitutive relation, Fitting degree value R2All close to 1.
It as shown in Figure 5A, is linear decrease relationship between girder size and the porosity of porous bone scaffold, between the two Constitutive relation is y=-42.6x+113.36, wherein x is girder size, and y is the porosity of porous bone scaffold.As shown in Figure 5 B, It is linear decrease relationship between girder size and the aperture of the maximum pore of porous bone scaffold, constitutive relation between the two is y =-1.065x+2.1445, wherein x is girder size, and y is the aperture of the maximum pore of porous bone scaffold.As shown in Figure 5 C, small It is linear increment relationship between beam size and the volume of porous bone scaffold, constitutive relation between the two is y=851.95x- 267.19, wherein x is girder size, and y is the volume of porous bone scaffold.As shown in Figure 5 D, girder size and porous bone scaffold It is linear increment relationship between surface area, constitutive relation between the two is y=1467.8x+676.28, wherein x is girder ruler Very little, y is the surface area of porous bone scaffold.It as shown in fig. 5e, is linear between girder size and the specific surface area of porous bone scaffold It is incremented by relationship, constitutive relation between the two is y=29.645x-2.4935, wherein x is girder size, and y is porous bone scaffold Specific surface area.
According to each Geometrical Parameter (porosity, the aperture of maximum pore, body of obtained girder size and porous bone scaffold Product, surface area, specific surface area) constitutive relation can go out the hole of porosity, maximum pore with careful design in the design process Diameter, volume, surface area, any one Geometrical Parameter in specific surface area, according to the specific Geometrical Parameter determined, according to small Constitutive relation between beam size and the Geometrical Parameter, the girder size designed needed for solving, after determining girder size, Allothimorph parameter can be determined according to the constitutive relation between girder size and corresponding Geometrical Parameter.It should be noted that due to Girder size and each Geometrical Parameter are once linear relationship, and are influenced each other between each Geometrical Parameter, can only be by preset Girder size meets a specific Geometrical Parameter.
The embodiment of the present invention provides a kind of mechanics of bone implant being made of above-mentioned dodecahedron rod structure basic unit 1 Performance evaluation methodology establishes the quantitative analysis relationship between girder size and equivalent elastic modulus, by adjusting dodecahedron bar The girder size of structural base member 1, its adjustable mechanical property, so as to adjust by dodecahedron rod structure basic unit 1 The mechanical property of the bone implant designed.
In a specific embodiment, it is cylindrical body, diameter × high in the second of 10 × 30mm of φ that design, which obtains outer profile, Porous bone scaffold model, the second porous bone scaffold model are used to carry out mechanical property assessment by experiment.By the second porous bone Stent model, using electron beam 3D printing equipment and titanium alloy (Ti6Al4V) metal powder, under the conditions of the thickness of 0.05mm, The porous bone scaffold exemplar of titanium alloy is prepared in printing, carries out mechanical property to porous bone scaffold exemplar using mechanics machine It tests (MPT, Mechanics Performance Testing), obtains the equivalent elastic modulus of porous bone scaffold.Meanwhile benefit With finite element analysis software (FEA, Finite Element Analysis), on the basis of the second porous bone scaffold model, into Row simulation analysis obtains the equivalent elastic modulus of porous bone scaffold.Further using SPSS statistics software to mechanical property The data obtained with finite element analysis can be tested and carry out statistical disposition, obtain the Geometrical Parameter and two methods of porous bone scaffold Mathematics constitutive relation between the equivalent elastic modulus for the porous bone scaffold that (finite element analysis and mechanics machine) measures, fitting Degree value R2All close to 1.
Fig. 6 A be the embodiment of the present invention girder size and equivalent elastic modulus constitutive relation schematic diagram, as shown, According to finite element simulation Measurement results (FEA data in such as figure), the equivalent elastic modulus of girder size and porous bone scaffold Between for power exponent be incremented by relationship, constitutive relation between the two be y=30.321x2.589, wherein x is girder size, and y is more The equivalent elastic modulus of hole bone bracket.According to the mechanics machine test result (MPT data in such as figure) of porous bone scaffold exemplar, It is linear increment relationship between girder size and the equivalent elastic modulus of porous bone scaffold, constitutive relation between the two is y= 6.11x-07263, wherein x is girder size, and y is the equivalent elastic modulus of porous bone scaffold.
Fig. 6 B is the porosity of the embodiment of the present invention and the constitutive relation schematic diagram of equivalent elastic modulus, as shown, root It is power between the porosity of porous bone scaffold and the equivalent elastic modulus of porous bone scaffold according to finite element simulation Measurement results Exponential decrease relationship, constitutive relation between the two are y=7E+12x-6.15, wherein x is the porosity of porous bone scaffold, and y is The equivalent elastic modulus of porous bone scaffold.According to the mechanics machine test result of porous bone scaffold exemplar, porous bone scaffold It is linear decrease relationship between porosity and the equivalent elastic modulus of porous bone scaffold, constitutive relation between the two is y=- 0.0989x+9.5125, wherein x is the porosity of porous bone scaffold, and y is the equivalent elastic modulus of porous bone scaffold.
Fig. 6 C be the embodiment of the present invention aperture and equivalent elastic modulus constitutive relation schematic diagram, as shown, according to The equivalent elastic modulus of finite element simulation Measurement results, the aperture of the maximum pore of porous bone scaffold and porous bone scaffold it Between successively decrease relationship for power exponent, constitutive relation between the two is y=39.772x-4.005, wherein x be porous bone scaffold most Macroporous aperture, y are the equivalent elastic modulus of porous bone scaffold.It is tested and is tied according to the mechanics machine of porous bone scaffold exemplar Fruit is linear decrease relationship between the aperture of the maximum pore of porous bone scaffold and the equivalent elastic modulus of porous bone scaffold, two Constitutive relation between person is y=-3.9826x+6.8216, wherein x is the aperture of the maximum pore of porous bone scaffold, and y is more The equivalent elastic modulus of hole bone bracket.
Fig. 6 D be the embodiment of the present invention volume and equivalent elastic modulus constitutive relation schematic diagram, as shown, according to Finite element simulation Measurement results are power exponent between the volume of porous bone scaffold and the equivalent elastic modulus of porous bone scaffold It is incremented by relationship, constitutive relation between the two is y=0.0024x1.4493, wherein x is the volume of porous bone scaffold, and y is porous The equivalent elastic modulus of bone bracket.According to the mechanics machine test result of porous bone scaffold exemplar, the volume of porous bone scaffold It is linear increment relationship between the equivalent elastic modulus of porous bone scaffold, constitutive relation between the two is y=0.0042x- 0.3371, wherein x is the volume of porous bone scaffold, and y is the equivalent elastic modulus of porous bone scaffold.
Fig. 6 E is the surface area of the embodiment of the present invention and the constitutive relation schematic diagram of equivalent elastic modulus, as shown, root It is power between the surface area of porous bone scaffold and the equivalent elastic modulus of porous bone scaffold according to finite element simulation Measurement results Exponential increasing relationship, constitutive relation between the two are y=5E-11x3.444, wherein x is the surface area of porous bone scaffold, and y is The equivalent elastic modulus of porous bone scaffold.According to the mechanics machine test result of porous bone scaffold exemplar, porous bone scaffold It is linear increment relationship between surface area and the equivalent elastic modulus of porous bone scaffold, constitutive relation between the two is y= 0.002x-2.8805, wherein x is the surface area of porous bone scaffold, and y is the equivalent elastic modulus of porous bone scaffold.
Fig. 6 F be the embodiment of the present invention specific surface area and equivalent elastic modulus constitutive relation schematic diagram, as shown, According to finite element simulation Measurement results, between the specific surface area of porous bone scaffold and the equivalent elastic modulus of porous bone scaffold It is incremented by relationship for power exponent, constitutive relation between the two is y=0.0096x2.4973, wherein x is the ratio table of porous bone scaffold Area, y are the equivalent elastic modulus of porous bone scaffold.It is porous according to the mechanics machine test result of porous bone scaffold exemplar It is linear increment relationship between the specific surface area of bone bracket and the equivalent elastic modulus of porous bone scaffold, this structure between the two closes System is y=0.1626x-1.5883, wherein x is the specific surface area of porous bone scaffold, and y is the Equivalent Elasticity mould of porous bone scaffold Amount.
According to Fig. 6 A to Fig. 6 F, for same second porous bone scaffold, the Equivalent Elasticity obtained using finite element analysis Constitutive relation between modulus and Geometrical Parameter is power exponent relationship;Second porous bone scaffold model is printed as porous bone scaffold Exemplar tests porous bone scaffold exemplar using mechanics machine, obtains the equivalent elastic modulus and Geometrical Parameter of porous bone scaffold Between constitutive relation be linear relationship.It is demonstrated experimentally that the second porous bone scaffold model designed is prepared with by 3D printing Porous bone scaffold exemplar compare, the mechanical property of the two is inconsistent, this is because the porous bone scaffold sample that 3D printing is prepared There are toughness deficiency, inclined brittleness etc. to influence the factor of bone implant mechanical property for part.
Fig. 7 is the girder size of the embodiment of the present invention and is utilized respectively what finite element analysis was tested with mechanics machine The correspondence diagram of equivalent elastic modulus.As shown, girder size under the same conditions, obtained using finite element analysis The equivalent elastic modulus of the porous bone scaffold arrived, it is equivalent greater than the porous bone scaffold exemplar tested using mechanics machine Elasticity modulus, this is because the problems such as there are powder to remain during 3D printing, sintering is incomplete, manufacturing defect, reduces system The mechanical property of standby bone implant out.
Fig. 8 is that the girder size of the embodiment of the present invention and this structure of concrete moduli difference (difference of equivalent elastic modulus) close It is schematic diagram.Mechanics is utilized as shown, the equivalent elastic modulus of the porous bone scaffold obtained using finite element analysis is subtracted The equivalent elastic modulus for the porous bone scaffold that testing machine is tested is obtained the difference of equivalent elastic modulus, is counted using SPSS It learns software and carries out data statistic analysis, obtain the constitutive relation between the concrete moduli difference of porous bone implant and girder size For y=25.091x3.1856, wherein x is girder size, and y is concrete moduli difference.It can be seen that, on the one hand, by manufacturing process, material The influence of the factors such as characteristic, the mechanical property for printing the porous bone implant exemplar produced are lower than the porous bone implant designed The mechanical property of body Model;On the other hand, the equivalent elastic modulus and power for the porous bone implant model that finite element analysis obtains The difference that the equivalent elastic modulus for the porous bone implant exemplar that testing machine is tested subtracts each other is learned, is existed with girder size Apparent constitutive relation can fast and accurately predict the porous bone implant sample produced according to this quantitative analysis relationship The mechanical property of part can also design adaptable porous bone according to the mechanical property of preset porous bone implant exemplar Implant cast.For example, as it is known that girder size, the Equivalent Elasticity mould of porous bone implant model can be obtained using finite element analysis Amount, utilizes above-mentioned constitutive relation y=25.091x3.1856, can fast and accurately predict the porous bone implant that printing produces Equivalent elastic modulus, assess the mechanical property matching degree between the porous bone implant prepared and host bone in advance, have Conducive to the clinical adaptation and popularization and application of the porous bone implant of personalization of 3D printing preparation;Furthermore it is intended to produce satisfaction centainly The porous bone implant exemplar of mechanical property condition is based on above-mentioned constitutive relation, by adjusting girder size, while utilizing limited The mechanical property of the porous bone implant model of meta analysis analog simulation can design the porous bone plant for meeting mechanical property condition Enter body exemplar.
The mechanical property evaluating method of the porous bone implant based on dodecahedron rod structure of the embodiment of the present invention, 12 Face body rod structural base member has several holes in three-dimensional spatial distribution, is conducive to promote cell adherence, tissue is promoted to grow into Regeneration;Dodecahedron rod structure basic unit axisymmetricly structure, and Mosaic face there are six settings, at least one dodecahedron bar Structural unit is by arrangement, the splicing porous bone implant of structure, convenient for designing diversified individual character according to actual needs Change porous bone implant;Porous bone implant that is arranged by several dodecahedron rod structure basic units, splicing and constitute, with Existing bone implant is compared, and can reduce the rigidity of bone implant, avoids stress shielding risk, has good mechanical property Energy;The mechanical property evaluating method of porous bone implant, establishes the quantitative analysis between girder size and concrete moduli difference Relationship can predict the mechanical property of porous bone implant exemplar produced according to quantitative analysis relationship, and assessment is prepared in advance Porous bone implant and host bone between mechanical property matching degree, the porous bone of personalization for being conducive to 3D printing preparation plants Enter the clinical adaptation of body and promotes and applies;Also set up porous bone implant model Geometrical Parameter and porous bone implant etc. The quantitative analysis relationship for imitating elasticity modulus can be the controllable of boniness implant, quantization, science by various quantitative analysis relationships Design provides data foundation, is conducive to the clinical application and popularization of personalized bionical bone implant.
It should be understood by those ordinary skilled in the art that: the discussion of any of the above embodiment is exemplary only, not It is intended to imply that the scope of the present disclosure (including claim) is limited to these examples;Under thinking of the invention, above embodiments Or can also be combined between the technical characteristic in different embodiments, step can be realized with random order, and be existed such as Many other variations of the upper different aspect of the invention, for simplicity, they are not provided in details.
The embodiment of the present invention be intended to cover fall into all such replacements within the broad range of appended claims, Modifications and variations.Therefore, all within the spirits and principles of the present invention, any omission, modification, equivalent replacement, the improvement made Deng should all be included in the protection scope of the present invention.

Claims (10)

1. a kind of mechanical property evaluating method of the porous bone implant based on dodecahedron rod structure, which is characterized in that
The porous bone implant is spliced to form by least one dodecahedron rod structure basic unit using modeling method;
The dodecahedron rod structure basic unit includes: the granatohedron being made of 24 girders, and respectively by The girder that opposite two vertex in 12 faces of the granatohedron extend outward forms hole between adjacent girder, The dodecahedron rod structure basic unit is formed with multiple holes in three-dimensional spatial distribution;
The first quantitative analysis relationship between the concrete moduli difference and girder size of the porous bone implant is established, according to institute The first quantitative analysis relationship is stated, predicts the mechanical property of the porous bone implant;
Wherein, the concrete moduli difference be the porous bone implant obtained using simulating analysis equivalent elastic modulus with The difference of the equivalent elastic modulus for the porous bone implant tested using exemplar.
2. the method according to claim 1, wherein the concrete moduli difference of the porous bone implant with it is described Constitutive relation between girder size is y=25.091x3.1856, wherein x is the girder size, and y is that the concrete moduli is poor Value.
3. the method according to claim 1, wherein establish the Geometrical Parameter of the porous bone implant with it is described The second quantitative analysis relationship between the equivalent elastic modulus of porous bone implant, according to preset Geometrical Parameter, according to described Second quantitative analysis relationship determines the equivalent elastic modulus of the porous bone implant;Wherein, the Geometrical Parameter includes: described The volume of porous bone implant, porosity, surface area, specific surface area, maximum pore aperture.
4. the method according to claim 1, wherein being obtained based on simulating analysis: the girder size with Constitutive relation between the equivalent elastic modulus of the porous bone implant is y=30.321x2.589;It tests to obtain based on exemplar: Constitutive relation between the girder size and the equivalent elastic modulus of the porous bone implant is y=6.11x-07263, In, y is the equivalent elastic modulus of the porous bone implant exemplar, and x is girder size.
5. according to the method described in claim 3, it is characterized in that, being obtained based on simulating analysis: the porous bone implant Constitutive relation between the porosity of body and the equivalent elastic modulus of the porous bone implant is y=7E+12x-6.15;Based on sample Part is tested to obtain: this structure between the porosity of the porous bone implant and the equivalent elastic modulus of the porous bone implant Relationship is y=-0.0989x+9.5125, wherein y is the equivalent elastic modulus of the porous bone implant, and x is the porous bone The porosity of implant.
6. according to the method described in claim 3, it is characterized in that, being obtained based on simulating analysis: the porous bone implant Constitutive relation between the aperture of the maximum pore of body and the equivalent elastic modulus of the porous bone implant is y=39.772x-4.005;It tests to obtain based on exemplar: the aperture of the maximum pore of the porous bone implant and the porous bone implant etc. Imitating the constitutive relation between elasticity modulus is y=-3.9826x+6.8216, wherein y is the equivalent bullet of the porous bone implant Property modulus, x be the porous bone implant maximum pore aperture.
7. according to the method described in claim 3, it is characterized in that, being obtained based on simulating analysis: the porous bone implant Constitutive relation between the volume of body and the equivalent elastic modulus of the porous bone implant is y=0.0024x1.4493;Based on sample Part is tested to obtain: this structure between the volume of the porous bone implant and the equivalent elastic modulus of the porous bone implant closes System is y=0.0042x-0.3371, wherein y is the equivalent elastic modulus of porous bone implant, and x is the body of porous bone implant Product.
8. according to the method described in claim 3, it is characterized in that, being obtained based on simulating analysis: the porous bone implant Constitutive relation between the surface area of body and the equivalent elastic modulus of the porous bone implant is y=5E-11x3.444;Based on sample Part is tested to obtain: this structure between the surface area of the porous bone implant and the equivalent elastic modulus of the porous bone implant Relationship is y=0.002x-2.8805, wherein y is the equivalent elastic modulus of the porous bone implant, and x is that the porous bone is planted Enter the surface area of body.
9. according to the method described in claim 3, it is characterized in that, being obtained based on simulating analysis: the porous bone implant Constitutive relation between the specific surface area of body and the equivalent elastic modulus of porous bone implant is y=0.0096x2.4973, wherein y For the equivalent elastic modulus of the porous bone implant;It tests to obtain based on exemplar: the specific surface area of the porous bone implant Constitutive relation between the equivalent elastic modulus of the porous bone implant is y=0.1626x-1.5883, wherein y is institute The equivalent elastic modulus of porous bone implant is stated, x is the specific surface area of the porous bone implant.
10. the method according to claim 1, wherein the section of the girder is circle, the material of the girder For titanium alloy.
CN201910081513.2A 2019-01-28 2019-01-28 Mechanical property evaluation method of porous bone implant based on dodecahedron rod structure Active CN109992820B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910081513.2A CN109992820B (en) 2019-01-28 2019-01-28 Mechanical property evaluation method of porous bone implant based on dodecahedron rod structure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910081513.2A CN109992820B (en) 2019-01-28 2019-01-28 Mechanical property evaluation method of porous bone implant based on dodecahedron rod structure

Publications (2)

Publication Number Publication Date
CN109992820A true CN109992820A (en) 2019-07-09
CN109992820B CN109992820B (en) 2023-06-09

Family

ID=67129407

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910081513.2A Active CN109992820B (en) 2019-01-28 2019-01-28 Mechanical property evaluation method of porous bone implant based on dodecahedron rod structure

Country Status (1)

Country Link
CN (1) CN109992820B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111494060A (en) * 2020-04-28 2020-08-07 河北科技大学 Method for constructing bone implant and bone implant structural unit
CN112743088A (en) * 2020-12-28 2021-05-04 北京航星机器制造有限公司 Rhombic dodecahedron titanium alloy lattice structure, interlayer structure and manufacturing method
CN114398799A (en) * 2022-01-26 2022-04-26 华侨大学 Optimization design method of variable gradient implant structure of molar prosthesis

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN204971711U (en) * 2015-05-31 2016-01-20 西安赛隆金属材料有限责任公司 High strength aperture footpath metal bone trabecula
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
CN106650141A (en) * 2016-12-29 2017-05-10 北京航空航天大学 Uncertainty analysis method for predicting performance of periodic material

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN204971711U (en) * 2015-05-31 2016-01-20 西安赛隆金属材料有限责任公司 High strength aperture footpath metal bone trabecula
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
CN106650141A (en) * 2016-12-29 2017-05-10 北京航空航天大学 Uncertainty analysis method for predicting performance of periodic material

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王颖: "面向骨组织工程的多孔结构设计优化及其SLM制备实验研究", 《中国优秀硕士学位论文全文数据库医药卫生科技辑》 *
赵剑 等: "蜂窝芯体材料面外等效弹性模量预测与分析", 《航空材料学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111494060A (en) * 2020-04-28 2020-08-07 河北科技大学 Method for constructing bone implant and bone implant structural unit
CN112743088A (en) * 2020-12-28 2021-05-04 北京航星机器制造有限公司 Rhombic dodecahedron titanium alloy lattice structure, interlayer structure and manufacturing method
CN114398799A (en) * 2022-01-26 2022-04-26 华侨大学 Optimization design method of variable gradient implant structure of molar prosthesis

Also Published As

Publication number Publication date
CN109992820B (en) 2023-06-09

Similar Documents

Publication Publication Date Title
CN109992820A (en) The mechanical property evaluating method of porous bone implant based on dodecahedron rod structure
CN204971711U (en) High strength aperture footpath metal bone trabecula
EP2001381B1 (en) Method and device for producing a planar implant for a human or animal body, which planar implant is preformed corresponding to a desired anatomical shape
CN109938882A (en) The mechanical property evaluating method of porous bone implant based on octahedra rod structure
CN102087676B (en) Pore network model (PNM)-based bionic bone scaffold designing method
CN107563056A (en) Variable density porous metals orthopaedics implant preparation method based on topological optimization technology
CN104887351A (en) High-strength small-bore metal bone trabecula and preparation method of same
DE112012004953B4 (en) Oscillator device and method for measuring physical properties
EP3462129A1 (en) Device and method for measuring a surface topography and calibration method
EP2130402B1 (en) Method for producing an otoplasty
CN109938883A (en) Dodecahedron rod structure unit, bone implant and its mass-transfer performance evaluating method
CN105608744B (en) The internal void design method of imitative period minimal surface based on tetrahedral grid
CN109966027A (en) A kind of Bone Defect Repari Gradient Unit, porous support and preparation method
CN109925097B (en) Octahedral rod structure unit, bone implant and mass transfer performance evaluation method thereof
EP2322925B1 (en) Device for determining the impedance of cell layers
CN103655014A (en) Artificial joint wear evaluation method and artificial joint wear evaluation system
CN110037832B (en) Composite unit structure porous scaffold for bone repair and processing method
EP2443433B1 (en) Method and device for determining the location of corrosion sites in reinforced concrete
CN210664206U (en) Novel plant sample measuring device
CN107194913A (en) A kind of most suitable Research scale detection method and device of crop groups
US11162064B2 (en) Electrochemical measurement device and transducer
CN112468948B (en) Device and method for testing wearing comfort of earphone
Bengough et al. Image analysis of maize root caps—estimating cell numbers from 2-D longitudinal sections
CN108665511A (en) A kind of simulation speckle pattern production method based on secondary enthesis
CN211576776U (en) Soft tissue fixing splint

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