Performance testing method of 3D printing titanium alloy bone defect repair structure

Abstract

The application discloses a kind of 3D printing titanium alloy bone defect repair structure performance test method, including the following steps: bone defect data establishes the three-dimensional model of bone defect repair structure, and titanium alloy material is prepared into titanium alloy structure entity;All stress points are obtained, and marked on bone defect repair structure entity;Compression fatigue life coefficient of each sub stress body under load T1 condition is calculated;The strength coefficient of tantalum coating is calculated;Compression fatigue life coefficient and the strength coefficient of tantalum coating are used to calculate the quality coefficient of bone defect repair structure entity, and the quality of bone defect repair structure entity is evaluated according to quality coefficient.The application carries out quality evaluation to the prepared bone defect repair structure body from multiple parameters, multiple performance angles, can qualitatively, quantitatively test the performance, quality of titanium alloy bone defect repair structure theme and tantalum coating, provide important data reference and basis for subsequent biological research.

Description

Technical Field

[0001] This invention relates to the field of 3D printing performance testing technology, specifically to a performance testing method for 3D printed titanium alloy bone defect repair structures. Background Technology

[0002] Bone diseases, bone trauma, and congenital malformations seriously affect human physical and mental health, and can even be life-threatening. Bone diseases, bone trauma, and congenital malformations are often accompanied by bone defects. With the development of bone repair or replacement surgery, artificial bone defect repair is frequently used to restore human bone function and improve patients' quality of life.

[0003] With the development and advancement of 3D printing technology, it has become possible to produce prosthetics for bone defect repair using titanium alloys, which offer advantages such as light weight and excellent performance. However, due to differences in the performance and production processes of 3D printing equipment, the quality of the produced prosthetic structures for bone defect repair varies significantly. In the initial stages of equipment debugging, the lack of reliable data parameters makes it impossible to assess the quality, leading to greater difficulty in subsequent repairs when applied to the patient's bone defect location, and increasing the patient's risks. Summary of the Invention

[0004] To address the aforementioned shortcomings of existing technologies, this invention provides a performance testing method for 3D-printed titanium alloy bone defect repair structures. This method can qualitatively and quantitatively test the performance and quality of titanium alloy bone defect repair structures, providing a basis for subsequent biological research.

[0005] To achieve the above-mentioned objectives, the technical solution adopted by this invention is as follows:

[0006] A performance testing method for 3D-printed titanium alloy bone defect repair structures is provided, which includes the following steps:

[0007] S1: Obtain bone defect data based on CT data, establish a three-dimensional model of the bone defect repair structure using the bone defect data, and correct the parameters of the three-dimensional model so that the three-dimensional model has an anatomical shape and internal porous structure that matches the location of the bone defect.

[0008] S2: Extract the corrected 3D model parameters and input them into a metal powder printer to prepare a titanium alloy structural entity using titanium alloy material;

[0009] S3: A tantalum coating is sprayed onto the surface of a titanium alloy structure using a high-vacuum magnetron sputtering coating system. The spraying process uses set sputtering power, sputtering angle, argon gas pressure and sputtering time to form a complete bone defect repair structure.

[0010] S4: Based on the stress characteristics and stress effectiveness of the bone defect location, perform biomechanical analysis on the bone defect repair structure to obtain all stress points and mark them on the bone defect repair structure.

[0011] S5: On the bone defect repair structure, cut out sub-force-bearing bodies for each force-bearing point location to obtain several sub-force-bearing bodies, so that each sub-force-bearing body has only one force-bearing point.

[0012] S6: Take each sub-body in sequence and conduct compression mechanical tests on the mechanical testing machine. During the experiment, the force point on the sub-body is directly facing the force application component of the mechanical testing machine. Based on the results of biomechanical analysis, set the actual load T1 on each force point, and apply the load T1 to each sub-body using the mechanical testing machine. Calculate the compression fatigue life coefficient of each sub-body under the load T1 condition.

[0013] S7: Obtain the peeling area and cracking area of ​​the tantalum coating for each sub-load body under load T1, test the indentation area at the contact point of the force-applying component, and calculate the strength coefficient of the tantalum coating.

[0014] S8: Calculate the quality coefficient of the bone defect repair structure using the compressive fatigue life coefficient and the strength coefficient of the tantalum coating, and evaluate the quality of the bone defect repair structure based on the quality coefficient.

[0015] Further, step S4 includes:

[0016] S41: Use CT data to establish a complete model of the bone defect site, and implant the three-dimensional model of the bone defect repair structure into the bone defect site model so that the three-dimensional model of the bone defect repair structure matches the bone defect site model.

[0017] S42: Divide the outer surface of the three-dimensional model of the bone defect repair structure into several first analysis units and the outer surface of the bone defect site model into several second analysis units. Establish a three-dimensional coordinate system in the modeling space and obtain the coordinates (x, y, z) of the center of each first analysis unit in the three-dimensional coordinate system and the coordinates (x′, y′, z′) of the center of each second analysis unit in the three-dimensional coordinate system.

[0018] S43: Calculate the distance from the center of any first analysis unit to the center of each second analysis unit to obtain the distance data set.

[0019]

[0020] in, Let x be the distance between the center of the i-th first analysis unit and the center of the n-th second analysis unit, where i is the number of the first analysis unit and n is the number of the second analysis units.i ,y i ,z i Let (x) be the coordinates of the center of the i-th first analysis unit in the three-dimensional coordinate system. n ′,y n ′,z n (′) represents the coordinates of the center of the nth second analysis unit in the three-dimensional coordinate system;

[0021] S44: Compare distance data sets Calculate the value of each distance within the range and filter out the minimum value. Minimum value The corresponding second analysis unit is closest to the i-th first analysis unit;

[0022] S45: Repeat steps S43-S44 to calculate the distance between the center of the remaining first analysis unit and the center of each second analysis unit, and filter out the minimum distance to obtain the data set with the minimum distance. This represents the minimum distance corresponding to the m-th first analysis unit, where m is the number of first analysis units.

[0023] S45: Set a threshold d for determining the distance d required for contact between the bone defect site and the bone defect repair structure. yuzhi And each minimum distance value in the minimum distance data group is compared with the distance threshold d. yuzhi Comparison:

[0024] like Then it is determined that the i-th first analysis unit will be subjected to force, and the position corresponding to the i-th first analysis unit is taken as the force point;

[0025] like Then it is determined that the i-th first analysis unit will not be subjected to force, and the position corresponding to the i-th first analysis unit is not regarded as a force point;

[0026] S46: Based on the locations of all stress points marked on the three-dimensional model of the bone defect repair structure, mark the stress points on the actual defect repair structure.

[0027] Further, step S6 includes:

[0028] S61: Take each sub-body in sequence and perform compression mechanics test on the mechanical testing machine. During the test, the force point on the sub-body is directly facing the force-applying component of the mechanical testing machine.

[0029] S62: Based on biomechanical analysis, set the actual load T1 on each force point, and apply different loads T1 to each sub-force body using a mechanical testing machine;

[0030] S63: Set the duration v for applying load T1 to the sub-force-bearing body. After each application of load T1 to the sub-force-bearing body, record the distance l that the force-applying component moves and calculate the change in distance Δl = l - l0, where l0 is the initial length of the force-applying component.

[0031] S64: Based on the allowable deformation L of the sub-body, compare the deformation L with the change Δl. If Δl ≥ L, the sub-body has reached its deformation limit. Record the number of times load T1 has been applied to the sub-body at this point, and calculate the total duration for which the sub-body bears load T1.

[0032] S65: Repeat steps S61-S64 to calculate the total duration of the load T1 borne by all sub-load bodies, and obtain the total duration data set. A represents the number of sub-force-bearing bodies. This represents the total duration corresponding to the Ath sub-force-bearing body;

[0033] S66: Calculate the compressive fatigue life coefficient f1:

[0034]

[0035] in, Total duration data group The maximum value in, This represents the total duration corresponding to the a-th sub-force-bearing body.

[0036] Furthermore, the strength coefficient f2 of the tantalum coating is calculated as follows:

[0037]

[0038] Where s1, s2, and s3 are the average values ​​of the spalling area, cracking area, and indentation area of ​​all sub-force-bearing bodies under load T1, respectively, and k1, k2, and k3 are the influence coefficients of spalling, cracking, and indentation on the strength evaluation of tantalum coating, respectively.

[0039] Further, step S8 includes:

[0040] S81: Calculate the mass coefficient K of the bone defect repair structure using the compressive fatigue life coefficient and the strength coefficient of the tantalum coating.

[0041] K = e1 lnf1 - e2 lnf2 + b;

[0042] Where e1 and e2 are the influence coefficients of compressive fatigue life and tantalum coating strength on the quality of the bone defect repair structure, respectively, and b is the fluctuation coefficient.

[0043] S82: Quality coefficient of the bone defect repair structure under ideal conditions To assess the quality of bone defect repair structures produced under current production conditions;

[0044] like Therefore, the quality of the bone defect repair structure produced under the current production conditions meets the requirements;

[0045] like If the quality of the bone defect repair structure produced under the current production conditions does not meet the requirements, it is necessary to reset the parameters of the metal powder printer, as well as the sputtering power, sputtering angle, argon gas pressure, and sputtering time of the spraying process.

[0046] The beneficial effects of this invention are as follows: This invention evaluates the quality of the prepared bone defect repair structure from multiple parameter and performance perspectives. Mechanical performance tests are conducted at various stress points of the titanium alloy bone defect repair structure. Biomechanical evaluation and testing are performed from the perspectives of the mechanical safety, bone tissue safety, and stress-bearing effectiveness of the repair structure. Furthermore, based on the main body of the repair structure, magnetron sputtering tantalum coating technology is applied, and the physicochemical properties of the coated titanium alloy material are tested. The performance influence relationship between sputtering parameters (power, angle, gas pressure, and time) and coating characteristics (strength, thickness) is elucidated. Further research is conducted on whether the tantalum coating on the titanium alloy material is firm and does not peel off. This invention enables qualitative and quantitative testing of the performance and quality of the titanium alloy bone defect repair structure and the tantalum coating, providing important data references and basis for subsequent biological research. Attached Figure Description

[0047] Figure 1 A flowchart for performance testing methods of 3D-printed titanium alloy bone defect repair structures. Detailed Implementation

[0048] The specific embodiments of the present invention are described below to enable those skilled in the art to understand the present invention. However, it should be understood that the present invention is not limited to the scope of the specific embodiments. For those skilled in the art, various changes are obvious as long as they are within the spirit and scope of the present invention as defined and determined by the appended claims. All inventions utilizing the concept of the present invention are protected.

[0049] like Figure 1 As shown, the performance testing method for the 3D-printed titanium alloy bone defect repair structure of this scheme includes the following steps:

[0050] S1: Obtain bone defect data based on CT data, establish a three-dimensional model of the bone defect repair structure using the bone defect data, and correct the parameters of the three-dimensional model so that the three-dimensional model has an anatomical shape and internal porous structure that matches the location of the bone defect; in order to ensure the characteristics of light weight and high strength, a porous structure is generally used inside the titanium alloy bone defect repair structure.

[0051] S2: Extract the corrected 3D model parameters and input them into a metal powder printer to prepare a titanium alloy structural entity using titanium alloy material;

[0052] S3: A tantalum coating is sprayed onto the surface of a titanium alloy structure using a high-vacuum magnetron sputtering coating system. The spraying process uses set sputtering power, sputtering angle, argon gas pressure and sputtering time to form a complete bone defect repair structure.

[0053] S4: Based on the stress characteristics and stress effectiveness of the bone defect location, perform biomechanical analysis on the bone defect repair structure to obtain all stress points and mark them on the bone defect repair structure.

[0054] Step S4 includes:

[0055] S41: Use CT data to establish a complete model of the bone defect site, and implant the three-dimensional model of the bone defect repair structure into the bone defect site model so that the three-dimensional model of the bone defect repair structure matches the bone defect site model.

[0056] S42: Divide the outer surface of the 3D model of the bone defect repair structure into several first analysis units and the outer surface of the bone defect site model into several second analysis units. The areas of the first analysis units are the same. Establish a 3D coordinate system in the modeling space and obtain the coordinates (x, y, z) of the center of each first analysis unit in the 3D coordinate system and the coordinates (x′, y′, z′) of the center of each second analysis unit in the 3D coordinate system.

[0057] S43: Calculate the distance from the center of any first analysis unit to the center of each second analysis unit to obtain the distance data set.

[0058]

[0059] in, Let x be the distance between the center of the i-th first analysis unit and the center of the n-th second analysis unit, where i is the number of the first analysis unit and n is the number of the second analysis units. i ,y i ,z i Let (x) be the coordinates of the center of the i-th first analysis unit in the three-dimensional coordinate system. n ′,y n ′,z n (′) represents the coordinates of the center of the nth second analysis unit in the three-dimensional coordinate system;

[0060] S44: Compare distance data sets Calculate the value of each distance within the range and filter out the minimum value. Minimum value The corresponding second analysis unit is closest to the i-th first analysis unit;

[0061] S45: Repeat steps S43-S44 to calculate the distance between the center of the remaining first analysis unit and the center of each second analysis unit, and filter out the minimum distance to obtain the data set with the minimum distance. This represents the minimum distance corresponding to the m-th first analysis unit, where m is the number of first analysis units.

[0062] S45: Set a threshold d for determining the distance d required for contact between the bone defect site and the bone defect repair structure. yuzhi And each minimum distance value in the minimum distance data group is compared with the distance threshold d. yuzhi Comparison:

[0063] like Then it is determined that the i-th first analysis unit will be subjected to force, and the position corresponding to the i-th first analysis unit is taken as the force point;

[0064] like Then it is determined that the i-th first analysis unit will not be subjected to force, and the position corresponding to the i-th first analysis unit is not regarded as a force point;

[0065] S46: Based on the locations of all stress points marked on the three-dimensional model of the bone defect repair structure, mark the stress points on the actual defect repair structure.

[0066] S5: On the bone defect repair structure, for each stress point, cut out sub-stress bodies to obtain several sub-stress bodies, so that each sub-stress body has only one stress point; the cut sub-stress bodies are more conducive to the mechanical analysis of the bone defect repair structure, and the different sub-stress bodies have different positions and different stress magnitudes.

[0067] S6: Each sub-body is sequentially subjected to compression mechanics tests on a mechanical testing machine. During the experiment, the force-bearing point on the sub-body is directly opposite the force-applying component of the mechanical testing machine. Based on the results of biomechanical analysis, the actual load T1 on each force-bearing point is set, and the load T1 is applied to each sub-body using the mechanical testing machine. The compressive fatigue life coefficient of each sub-body under load T1 is calculated. Through biomechanical analysis of different locations of bone defects, the magnitude of the force at different locations can be obtained.

[0068] Step S6 includes:

[0069] S61: Take each sub-body in sequence and perform compression mechanics test on the mechanical testing machine. During the test, the force point on the sub-body is directly facing the force-applying component of the mechanical testing machine.

[0070] S62: Based on biomechanical analysis, set the actual load T1 on each force point, and apply different loads T1 to each sub-force body using a mechanical testing machine;

[0071] S63: Set the duration v for applying load T1 to the sub-force-bearing body. After each application of load T1 to the sub-force-bearing body, record the distance l that the force-applying component moves and calculate the change in distance Δl = l - l0, where l0 is the initial length of the force-applying component.

[0072] S64: Based on the allowable deformation L of the sub-body, compare the deformation L with the change Δl. If Δl ≥ L, the sub-body has reached its deformation limit. Record the number of times load T1 has been applied to the sub-body at this point, and calculate the total duration for which the sub-body bears load T1.

[0073] S65: Repeat steps S61-S64 to calculate the total duration of the load T1 borne by all sub-load bodies, and obtain the total duration data set. A represents the number of sub-force-bearing bodies. This represents the total duration corresponding to the Ath sub-force-bearing body;

[0074] S66: Calculate the compressive fatigue life coefficient f1:

[0075]

[0076] in, Total duration data group The maximum value in, This represents the total duration corresponding to the a-th sub-force-bearing body.

[0077] S7: Obtain the peeling area and cracking area of ​​the tantalum coating for each sub-loaded body under load T1, test the indentation area at the contact point of the force-applying component, and calculate the strength coefficient of the tantalum coating.

[0078] The strength coefficient f2 of the tantalum coating is calculated as follows:

[0079]

[0080] Where s1, s2, and s3 are the average values ​​of the spalling area, cracking area, and indentation area of ​​all sub-force-bearing bodies under load T1, respectively, and k1, k2, and k3 are the influence coefficients of spalling, cracking, and indentation on the strength evaluation of tantalum coating, respectively.

[0081] S8: Calculate the quality coefficient of the bone defect repair structure using the compressive fatigue life coefficient and the strength coefficient of the tantalum coating, and evaluate the quality of the bone defect repair structure based on the quality coefficient.

[0082] Step S8 includes:

[0083] S81: Calculate the mass coefficient K of the bone defect repair structure using the compressive fatigue life coefficient and the strength coefficient of the tantalum coating.

[0084] K = e1 lnf1 - e2 lnf2 + b;

[0085] Where e1 and e2 are the influence coefficients of compressive fatigue life and tantalum coating strength on the quality of the bone defect repair structure, respectively, and b is the fluctuation coefficient.

[0086] S82: Quality coefficient of the bone defect repair structure under ideal conditions To assess the quality of bone defect repair structures produced under current production conditions;

[0087] like Therefore, the quality of the bone defect repair structure produced under the current production conditions meets the requirements;

[0088] like If the quality of the bone defect repair structure produced under the current production conditions does not meet the requirements, it is necessary to reset the parameters of the metal powder printer, as well as the sputtering power, sputtering angle, argon gas pressure, and sputtering time of the spraying process.

[0089] This invention evaluates the quality of the prepared bone defect repair structure from multiple parameter and performance perspectives. Mechanical performance tests are conducted at various stress points of the titanium alloy bone defect repair structure. Biomechanical assessments and tests are performed from the perspectives of the repair structure's mechanical safety, bone tissue safety, and stress-bearing effectiveness. Furthermore, based on the main body of the repair structure, magnetron sputtering tantalum coating technology is applied, and the physicochemical properties of the coated titanium alloy material are tested. The performance influence relationship between sputtering parameters (power, angle, gas pressure, and time) and coating characteristics (strength, thickness) is elucidated. Further research is conducted on whether the tantalum coating on the titanium alloy material is firm and does not peel off. This invention enables qualitative and quantitative testing of the performance and quality of the titanium alloy bone defect repair structure and the tantalum coating, providing important data references and basis for subsequent biological research.

Claims

1. A performance testing method for a 3D-printed titanium alloy bone defect repair structure, characterized in that, Includes the following steps: S1: Obtain bone defect data based on CT data, establish a three-dimensional model of the bone defect repair structure using the bone defect data, and correct the parameters of the three-dimensional model so that the three-dimensional model has an anatomical shape and internal porous structure that matches the location of the bone defect. S2: Extract the corrected 3D model parameters and input them into a metal powder printer to prepare a titanium alloy structural entity using titanium alloy material; S3: A tantalum coating is sprayed onto the surface of a titanium alloy structure using a high-vacuum magnetron sputtering coating system. The spraying process uses set sputtering power, sputtering angle, argon gas pressure and sputtering time to form a complete bone defect repair structure. S4: Based on the stress characteristics and stress effectiveness of the bone defect location, perform biomechanical analysis on the bone defect repair structure to obtain all stress points and mark them on the bone defect repair structure. S5: On the bone defect repair structure, cut out sub-force-bearing bodies for each force-bearing point location to obtain several sub-force-bearing bodies, so that each sub-force-bearing body has only one force-bearing point. S6: Each sub-body is subjected to compression mechanics tests on a mechanical testing machine in sequence. During the experiment, the force application point on the sub-body is directly opposite the force application component of the mechanical testing machine. Based on the results of biomechanical analysis, the actual load on each force application point is set. T 1. And apply loads to each sub-body using a mechanical testing machine. T 1. Calculate the load on each sub-body. T Compression fatigue life coefficient under condition 1; S7: Obtain the load on each sub-body T Under the condition of 1, the area of ​​detachment and cracking of tantalum coating, the area of ​​indentation on the contact point of the force-applying component, and the strength coefficient of tantalum coating are calculated. S8: Calculate the quality coefficient of the bone defect repair structure using the compressive fatigue life coefficient and the strength coefficient of the tantalum coating, and evaluate the quality of the bone defect repair structure based on the quality coefficient. Step S4 includes: S41: Use CT data to establish a complete model of the bone defect site, and implant the three-dimensional model of the bone defect repair structure into the bone defect site model so that the three-dimensional model of the bone defect repair structure matches the bone defect site model. S42: Divide the outer surface of the 3D model of the bone defect repair structure into several first analysis units, and the outer surface of the bone defect site model into several second analysis units. Establish a 3D coordinate system in the modeling space, and obtain the coordinates of the center of each first analysis unit in the 3D coordinate system. The coordinates of the center of each second analysis unit in the three-dimensional coordinate system ; S43: Calculate the distance from the center of any first analysis unit to the center of each second analysis unit to obtain the distance data set. : ； in, For the first i The center of the first analysis unit and the first n The distance from the center of each second analysis unit i This is the number of the first analysis unit. n The number of the second analysis units. For the first i The coordinates of the center of the first analysis unit in the three-dimensional coordinate system For the first n The coordinates of the center of the second analysis unit in the three-dimensional coordinate system; S44: Compare distance data sets Calculate the value of each distance within the range and filter out the minimum value. minimum value The corresponding second analysis unit and the first i The first analysis unit is closest; S45: Repeat steps S43-S44 to calculate the distance between the center of the remaining first analysis unit and the center of each second analysis unit, and filter out the minimum distance to obtain the data set with the minimum distance. , For the first m The minimum distance corresponding to the first analysis unit m This represents the number of the first analysis units; S45: Set a threshold distance for determining whether the bone defect site and the bone defect repair structure can make contact. And each minimum distance value in the minimum distance data group is compared with the distance threshold. Comparison: like ≤ Then determine the corresponding number i The first analysis unit will be subjected to force, the first i The position corresponding to the first analysis unit is taken as the force point; like > Then determine the corresponding number i The first analysis unit will not be subjected to force, the second i The location corresponding to the first analysis unit is not considered as a stress point; S46: Based on the locations of all stress points marked on the three-dimensional model of the bone defect repair structure, mark the stress points on the actual defect repair structure. Step S6 includes: S61: Take each sub-body in sequence and perform compression mechanics test on the mechanical testing machine. During the test, the force point on the sub-body is directly facing the force-applying component of the mechanical testing machine. S62: Based on biomechanical analysis, set the actual load on each stress point. T 1. Different loads were applied to each sub-body using a mechanical testing machine. T 1; S63: Set the load to be applied to the sub-load body T Duration of 1 v Each time a load is applied to the sub-body T After step 1, record the distance the force-applying component moves. l And calculate the change in the distance traveled. , The initial length of the force-applying component; S64: Based on the allowable deformation of the sub-body. L Compare the amount of deformation L With change The size, if ≥ L If the sub-body reaches its deformation limit, record the load applied to the sub-body at this point. T The number of times 1 t Calculate the load borne by the sub-body. T Total duration of 1 ; S65: Repeat steps S61-S64 to calculate the load on all sub-load-bearing bodies in sequence. T After calculating the total duration, the total duration data set is obtained. , A For the number of sub-force-bearing bodies, For the first A The total duration corresponding to each load-bearing body; S66: Calculate the compressive fatigue life coefficient : ； in, Total duration data group The maximum value in, For the first a The total duration corresponding to each load-bearing body; The strength coefficient of the tantalum coating The calculation method is as follows: ； in, Each of the sub-body under load T The average value of the spalling area, crack area, and indentation area under condition 1. These are the influence coefficients of peeling, cracking, and indentation on the strength assessment of tantalum coatings; Step S8 includes: S81: Calculate the mass coefficient of the bone defect repair structure using the compressive fatigue life coefficient and the strength coefficient of the tantalum coating. K : ； in, These are the coefficients representing the influence of compressive fatigue life and tantalum coating strength on the quality of the bone defect repair structure, respectively. b This is the volatility coefficient; S82: Quality coefficient of the bone defect repair structure under ideal conditions To assess the quality of bone defect repair structures produced under current production conditions; like K ≥ If so, the quality of the bone defect repair structure produced under the current production conditions meets the requirements; like K < If the quality of the bone defect repair structure produced under the current production conditions does not meet the requirements, it is necessary to reset the parameters of the metal powder printer, as well as the sputtering power, sputtering angle, argon gas pressure, and sputtering time of the spraying process.

Images