Simulated live wire treatment training method and system for high-risk orthopedic trauma

Abstract

The invention relates to the technical field of medical simulation training, in particular to a simulation live wire treatment training method and system for high-risk orthopedic wounds, and the method comprises the following steps: calling preset explosion shock wave parameters, calculating the influence of shock waves on bones of multiple parts in combination with the stress-strain characteristics of bone tissues, and analyzing the mechanical response of the bone tissues, and simulating fracture injury conditions to obtain injury simulation data. According to the invention, through refined simulation of fracture injury and fracture complications, injury condition data are close to a real battlefield environment, students can visually understand dynamic evolution of wounds, and by collecting joint movement data of the students in real time and combining mechanical response of fracture parts, direct influence of treatment actions on the wounds is analyzed, so that the treatment effect of the wounds is improved. In the training process, wound state mapping is adjusted in real time, students are helped to understand operation feedback in real time, in combination with assessment of treatment skills, the skill level of the students is helped to be understood, a data basis is provided for adjusting training courses, and the skill training effect and efficiency are improved.

Description

Technical Field

[0001] The present invention relates to the technical field of medical simulation training, and particularly to a high-risk orthopedic trauma simulation on-the-spot treatment training method and system. Background Art

[0002] The technical field of medical simulation training includes using computer simulation technology, virtual reality technology, and augmented reality technology to simulate different medical scenarios and treatment operations for medical operation training and skill improvement. The core content of this field is to construct a virtual environment or physical model so that trainees can perform practical operations in a realistic simulation environment and master medical treatment skills. Medical simulation training is applied to the training in multiple fields such as basic medical education, first aid, surgical operations, and emergency treatment. In this way, it helps trainees repeatedly perform operation training without actual risks and improve their emergency response ability and treatment ability in real situations.

[0003] Among them, the high-risk orthopedic trauma simulation on-the-spot treatment training method refers to using mixed reality technology, virtual reality technology, and motion capture technology to simulate high-risk orthopedic trauma scenarios to help trainees conduct trauma treatment skill training. The patent theme covers the virtual reconstruction of trauma scenarios, including simulating traumas such as limb fractures and massive bleeding, and guiding trainees to perform correct hemostasis, fracture fixation, etc. through establishing corresponding treatment processes. By using multiple means such as trauma scenario data input, human model simulation, motion capture, and target recognition technology, a highly restored training platform is provided to help trainees perform real-time operations in the simulation environment and adjust operation steps and techniques through system feedback.

[0004] Traditional orthopedic treatment training methods have deficiencies in the refined simulation of fracture injuries. The occurrence mechanism and injury degree of fracture traumas are often based on static modeling and cannot reflect the changes of bone tissues under different external force impacts in real time, resulting in limited authenticity of trauma scenarios. The evolution of trauma complications uses a preset model instead of real-time calculation based on hemodynamics, making it difficult to accurately predict complex complications such as ischemic necrosis and fracture dislocation, weakening the medical reference value of training. The evaluation of treatment actions relies on static comparison and cannot make full use of kinematic analysis for motion capture and deviation recognition, resulting in a lag in the operation feedback of trainees and making it difficult to immediately correct incorrect actions during the training process. The impact of treatment operations on the progress of trauma is not dynamically updated, and trainees cannot intuitively feel the immediate effect of operations on the injury condition during the training process, affecting the actual combat adaptability of battlefield emergency treatment skills. Summary of the Invention

[0005] In order to solve the technical problems existing in the prior art, embodiments of the present invention provide a high-risk orthopedic trauma simulation on-the-spot treatment training method and system. The technical solutions are as follows:

[0006] A training method for simulating on-site treatment of high-risk orthopedic trauma, comprising the following steps:

[0007] S1: Call the preset explosion shock wave parameters, combine with the stress-strain characteristics of bone tissue, calculate the impact of the shock wave on the bones of multiple parts, analyze the mechanical response of bone tissue, simulate fracture injuries, and obtain injury simulation data;

[0008] S2: Based on the injury simulation data, simulate the blood flow changes at the fracture trauma site, calculate the probability of local ischemia and blood interruption, and map the trauma complications to the dummy model to obtain injury mapping parameters;

[0009] S3: Call the injury mapping parameters, collect the joint movement data of the trainees in real time, obtain the angular change data of multiple joints, calculate the movement trajectory and acceleration of the joints, and obtain action deviation data by comparing with the standard treatment actions;

[0010] S4: Based on the action deviation data, combine with the mechanical response of the fracture site, analyze the impact of the treatment actions of multiple trainees on the trauma site, update the injury progress in real time and adjust the dummy mapping parameters to obtain injury update data;

[0011] S5: According to the injury update data, analyze the operation accuracy of each action of the trainees after the treatment is completed, combine with the treatment duration and treatment effect, calculate the trainee skill score, and obtain the treatment skill score.

[0012] As a further solution of the present invention, the injury simulation data specifically includes fracture type, bone tissue stress distribution, and fracture end displacement. The injury mapping parameters include local ischemia range, blood flow blockage probability, and complication type. The action deviation data specifically refers to joint angle error, movement trajectory deviation, and action stability score. The injury update data includes trauma deterioration degree, tissue repair rate, and dynamic injury status of the dummy model. The treatment skill score specifically refers to operation precision score, treatment efficiency score, and trauma recovery effect score.

[0013] As a further solution of the present invention, the step of calling the preset explosion shock wave parameters, combining with the stress-strain characteristics of bone tissue, calculating the impact of the shock wave on the bones of multiple parts, analyzing the mechanical response of bone tissue, simulating fracture injuries, and obtaining injury simulation data is specifically as follows:

[0014] S101: Obtain the preset explosion shock wave parameters, call the data of explosion wave pressure gradient, shock propagation speed, and shock duration, calculate the force conditions of the explosion shock wave at multiple bone sites, detect the stress distribution, deformation trend, and force area of bone tissue under the action of the shock wave, analyze the attenuation law of shock wave energy in bone tissue, and obtain the bone shock energy distribution data;

[0015] S102: Based on the bone impact energy distribution data, combined with the Young's modulus, yield stress, and trabecular bone structure parameters of bone tissue, calculate the stress concentration area, fracture occurrence probability, and fracture morphology during the force-bearing process of bone tissue, identify the fracture occurrence location, and obtain fracture stress characteristic parameters;

[0016] S103: According to the fracture stress characteristic parameters, combined with the mechanical response of bone tissue, simulate various degrees of fracture injuries, including fissure fractures, comminuted fractures, and complete fractures, and generate injury simulation data.

[0017] As a further solution of the present invention, the steps of simulating the blood flow change at the fracture trauma site based on the injury simulation data, calculating the probability of local ischemia and blood interruption, and mapping the trauma complications to the dummy model to obtain the injury mapping parameters are specifically as follows:

[0018] S201: Based on the injury simulation data, call the fracture morphology, tissue injury area, and bone stress distribution at the fracture site, calculate the degree of blood vessel damage, blood vessel diameter narrowing amount, and local blood flow velocity decrease rate in the fracture trauma area, and combine the blood flow resistance parameter and local venous return rate to calculate the blood flow supply change trend at the trauma site, and obtain the dynamic blood flow characteristics of the trauma site;

[0019] S202: According to the dynamic blood flow characteristics of the trauma site, calculate the local tissue perfusion reduction rate, microvascular occlusion probability, and oxygen supply level decrease amplitude, combine the tissue hypoxia threshold and blood perfusion recovery rate, judge the ischemic influence range of the fracture area, calculate the local blood flow interruption risk level, and obtain the local ischemia and blood flow blockage parameters;

[0020] S203: Call the local ischemia and blood flow blockage parameters, map the ischemic area, blood flow interruption point, and trauma complication category to the dummy model, and generate injury mapping parameters.

[0021] As a further solution of the present invention, the steps of calling the injury mapping parameters, collecting the joint movement data of the trainee in real time, obtaining the angle change data of multiple joints, calculating the movement trajectory and acceleration of the joints, and obtaining the action deviation data by comparing with the standard treatment actions are specifically as follows:

[0022] S301: Call the injury mapping parameters, collect the movement data of multiple joint parts of the trainee during the treatment process in real time, record the joint angle change, rotation rate, and acceleration characteristics during the treatment operation of the trainee, calculate the movement amplitude, angle change trend, and speed change rate of each joint, and obtain the trainee joint movement data;

[0023] S302: Calculate the motion trajectory, path curvature, and key time nodes of the treatment action based on the joint motion data of the trainee, evaluate the motion continuity, angular stability, and force application direction of each joint during the treatment operation, and obtain the treatment action trajectory features;

[0024] S303: Invoke the treatment action trajectory features, calculate the trajectory error of the trainee's operation by comparing with the standard treatment action, identify and record the abnormal postures existing during the treatment process, and generate action deviation data.

[0025] As a further solution of the present invention, the specific formula for calculating the trajectory error of the trainee's operation is:

[0026]

[0027] Calculate the trajectory error value, identify and record the abnormal postures existing during the treatment process, and generate action deviation data;

[0028] Among them, E traj is the trajectory error value, N' is the total number of trajectory sampling points, i' is the trajectory sampling point index, X i' is the X-axis coordinate value of the trainee at the i'-th trajectory sampling point, Y i' is the Y-axis coordinate value of the trainee at the i'-th trajectory sampling point, z i' is the Z-axis coordinate value of the trainee at the i'-th trajectory sampling point, is the X-axis coordinate value of the standard treatment action at the i'-th trajectory sampling point, is the Y-axis coordinate value of the standard treatment action at the i'-th trajectory sampling point, is the Z-axis coordinate value of the standard treatment action at the i'-th trajectory sampling point, M' is the total number of speed samplings during the treatment process, j' is the speed sampling point index, V j' is the hand movement speed of the trainee at the j'-th sampling, V std is the average speed of the standard treatment operation.

[0029] As a further solution of the present invention, based on the action deviation data, combined with the mechanical response of the fracture site, the steps of analyzing the influence of the treatment actions of multiple trainees on the trauma site, updating the injury progress in real time, and adjusting the dummy mapping parameters to obtain the injury update data are specifically as follows:

[0030] S401: Based on the action deviation data, combined with the stress state, displacement trend of the fracture end, and tissue stress distribution of the fracture site, calculate the influence of the treatment actions of multiple trainees on the trauma site, including local stress changes, misalignment amplitude of the fracture end, and tissue stretching deformation rate, and obtain the treatment operation influence parameters;

[0031] S402: Invoke the rescue operation impact parameters, analyze the tissue stress concentration area, the degree of blood flow supply obstruction, and the nerve injury spread trend of the fracture site, calculate the injury aggravation risk level, recovery delay rate, and local blood flow supply adjustment range of the trauma area, and obtain the injury progress status;

[0032] S403: Invoke the injury progress status, and adjust the mapping parameters of the dummy model in real time, including dynamic injury display, tissue necrosis process, and blood circulation recovery trend, to generate injury update data.

[0033] As a further solution of the present invention, according to the injury update data, after the rescue is completed, analyze the operation accuracy of each action of the trainee, and combine the rescue duration and rescue effect to calculate the trainee's skill score. The steps to obtain the rescue skill score are specifically as follows:

[0034] S501: Invoke the injury update data, calculate the operation precision and key action matching rate of the trainee in the rescue task by analyzing the joint angle adjustment range, force application direction deviation, and rescue posture stability of the trainee's rescue operation, and obtain the rescue operation accuracy data;

[0035] S502: Based on the rescue operation accuracy data, combine the trainee's rescue completion time, calculate the average rescue duration and emergency handling efficiency of the trainee in multiple rescue tasks, and obtain the rescue process efficiency parameter;

[0036] S503: Invoke the rescue process efficiency parameter, combine the trauma repair status, tissue recovery rate, and hemostasis stability of the injury progress during the rescue, obtain the trainee's rescue effect score, combine the operation accuracy and rescue efficiency, calculate the trainee's skill score, and generate the rescue skill score.

[0037] As a further solution of the present invention, the specific formula for calculating the trainee's skill score is:

[0038]

[0039] Calculate the rescue skill score;

[0040] Among them, S skill represents the trainee's skill score, E treat represents the rescue effect score, P acc represents the operation accuracy, T delay represents the delay time caused by stress or inexperience during the rescue process, W eff represents the rescue efficiency, R fatigue represents the correction value for the rescue performance affected by physiological fatigue.

[0041] On the other hand, a high-risk orthopedic trauma simulation and on-the-spot treatment training system is provided. This system is applied to the high-risk orthopedic trauma simulation and on-the-spot treatment training method, and the system includes:

[0042] The impact injury simulation module calculates the stress and strain effects of the shock wave on bone tissue under various pressure gradients, wave speed changes, and impact durations based on preset explosion shock wave parameters. Combining with the mechanical response characteristics of bone tissue, it analyzes the injury types of bones in multiple parts, including crack fractures and comminuted fractures, and obtains injury simulation data.

[0043] The complication mapping module simulates and analyzes the changes in the blood flow state at the trauma site based on the injury simulation data, predicts the occurrence probabilities of various complications, and maps them to the dummy model to display the pathological progress of the trauma site in real time, and obtains injury condition mapping parameters.

[0044] The treatment action capture module collects the joint angle change data of the trainee during the treatment operation in real time based on the injury condition mapping parameters, calculates the movement trajectories and accelerations of multiple joints, and compares them with the standard treatment action path to obtain action deviation data.

[0045] The treatment impact analysis module analyzes the impact of the trainee's treatment actions on the trauma area based on the action deviation data, combined with the mechanical response of the fracture site, adjusts the trauma mapping state of the dummy model, and displays the dynamic changes of the injury condition in real time to obtain injury condition update data.

[0046] The trainee skill assessment module analyzes the accuracy and standardization of each treatment action of the trainee after the treatment operation based on the injury condition update data. Combining with the treatment duration and the changes in the injury condition, it calculates the treatment ability score of the trainee to obtain the treatment skill score.

[0047] The beneficial effects brought by the technical solution provided by the embodiments of the present invention at least include:

[0048] In the embodiments of the present invention, through the refined simulation of fracture injuries and fracture complications, the injury condition data is close to the real battlefield environment, enabling trainees to intuitively understand the dynamic evolution of trauma. By collecting the joint movement data of trainees in real time and combining with the mechanical response of the fracture site, it analyzes the direct impact of treatment actions on trauma, and adjusts the trauma state mapping in real time during the training process to help trainees understand the operation feedback in real time. Combining with the assessment of treatment skills, it helps to understand the skill level of trainees, provides a data basis for adjusting the training curriculum, and improves the effect and efficiency of skill training. Description of the Drawings

[0049] To more clearly illustrate the technical solutions in the embodiments of the present invention, the following will briefly introduce the drawings required for the description of the embodiments. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0050] Figure 1 It is a schematic diagram of the working process of the present invention;

[0051] Figure 2 It is a system flowchart of the present invention. Specific embodiments

[0052] The following will describe the technical solutions in the present invention in conjunction with the drawings.

[0053] In the embodiments of the present invention, words such as "exemplarily" and "for example" are used to represent examples, illustrations or explanations. Any embodiment or design solution described as "example" in the present invention should not be construed as being more preferred or having more advantages than other embodiments or design solutions. Exactly speaking, the use of the word "example" is intended to present concepts in a specific way. In addition, in the embodiments of the present invention, the meaning expressed by "and / or" can be both, or either one of the two can be selected.

[0054] In the embodiments of the present invention, "image" and "picture" can sometimes be used interchangeably. It should be noted that when their differences are not emphasized, the meanings they express are the same. "(of)", "corresponding", and "corresponding" can sometimes be used interchangeably. It should be noted that when their differences are not emphasized, the meanings they express are the same.

[0055] In the embodiments of the present invention, sometimes subscripts such as W1 may be written in a non-subscript form such as W1. When their differences are not emphasized, the meanings they express are the same.

[0056] To make the technical problems, technical solutions and advantages to be solved by the present invention clearer, the following will be described in detail in conjunction with the drawings and specific embodiments.

[0057] Please refer to Figure 1 , the present invention provides a technical solution, a high-risk orthopedic trauma simulation and on-the-spot treatment training method, including the following steps:

[0058] S1: Call the preset blast wave parameters, combine with the stress-strain characteristics of bone tissue, calculate the influence of the blast wave on the bones of multiple parts, analyze the mechanical response of bone tissue, simulate fracture injuries, and obtain injury simulation data;

[0059] S2: Based on the injury simulation data, simulate the blood flow changes at the fracture trauma site, calculate the probabilities of local ischemia and blood interruption, and map the trauma complications onto the dummy model to obtain the injury mapping parameters;

[0060] S3: Invoke the injury mapping parameters, collect the joint movement data of the trainees in real time, obtain the angular change data of multiple joints, calculate the movement trajectories and accelerations of the joints, and obtain the action deviation data by comparing with the standard treatment actions;

[0061] S4: Based on the action deviation data, combined with the mechanical response of the fracture site, analyze the impacts of the treatment actions of multiple trainees on the trauma site, update the injury progress in real time and adjust the dummy mapping parameters to obtain the injury update data;

[0062] S5: According to the injury update data, analyze the operation accuracy of each action of the trainees after the treatment is completed, and combined with the treatment duration and treatment effect, calculate the trainee skill scores to obtain the treatment skill scores.

[0063] The injury simulation data specifically includes fracture type, bone tissue stress distribution, and fracture end displacement. The injury mapping parameters include local ischemia range, blood flow interruption probability, and complication type. The action deviation data specifically refers to joint angle error, movement trajectory deviation, and action stability score. The injury update data includes the degree of trauma deterioration, tissue repair rate, and dynamic injury state of the dummy model. The treatment skill scores specifically refer to operation accuracy score, treatment efficiency score, and trauma recovery effect score.

[0064] The steps of invoking the preset explosion shock wave parameters, combined with the stress-strain characteristics of the bone tissue, calculating the impacts of the shock wave on the bones of multiple parts, analyzing the mechanical response of the bone tissue, and simulating the fracture injury conditions to obtain the injury simulation data are specifically as follows:

[0065] S101: Obtain the preset explosion shock wave parameters, invoke the data of explosion wave pressure gradient, shock propagation speed, and shock duration, calculate the forces on multiple bone parts by the explosion shock wave, detect the stress distribution, deformation trend, and force-bearing area of the bone tissue under the action of the shock wave, analyze the attenuation law of the shock wave energy in the bone tissue, and obtain the bone shock energy distribution data;

[0066] Obtain preset blast shock wave parameters, set the explosion center point, explosion equivalent, and shock radius based on the battlefield simulation environment, determine the maximum pressure value, peak rise time, and waveform decay rate of the blast shock wave, call the data of the pressure gradient, shock propagation speed, and shock duration of the blast shock wave, calculate the decay characteristics of the blast wave within different radius ranges, use the logarithmic decay model to calculate the curve of the pressure of the shock wave in the air changing with distance, take the explosion point as the center, and calculate the impact force magnitude at the bone-loaded position according to the shock wave velocity distribution equation. Set the normal pressure and shear force exerted by the shock wave on the bone, call the Young's modulus, density, and shear modulus of the bone tissue, calculate the stress transfer coefficient when the shock wave acts on the bone tissue, use the equilibrium mechanics equation to calculate the stress distribution inside the bone, detect the stress concentration area, deformation trend, and force-bearing area of the bone tissue under the action of the shock wave, combine the biomechanical characteristics, calculate the energy decay rate of the shock wave in the bone tissue, calculate the stress transfer situation at different parts inside the bone based on the wave propagation equation, calculate the absorption ratio of the shock wave energy in the bone tissue, determine the maximum force-bearing point of the bone after being impacted, and establish a simulation of the bone-loading scenario under the action of the blast shock in combination with the force distribution characteristics of the bone tissue.

[0067] The pressure P of the blast shock wave can be calculated by the shock wave decay formula:

[0068]

[0069] where P is the shock wave pressure at a certain point (Pa), P0 is the initial pressure at the explosion point (Pa), r0 is the reference distance from the explosion point to the calculation point (m), and r is the actual distance from the current calculation point to the explosion point (m).

[0070] If the initial pressure at the explosion point P0 = 500000 Pa, the reference distance r0 = 1 m, and the distance of the target calculation point r = 5 m, then the shock wave pressure at this point is calculated as follows:

[0071]

[0072] Finally, obtain the data of the bone impact energy distribution.

[0073] S102: Based on the bone impact energy distribution data, combine the Young's modulus, yield stress, and trabecular bone structure parameters of the bone tissue, calculate the stress concentration area, fracture occurrence probability, and fracture morphology during the force-bearing process of the bone tissue, identify the fracture occurrence location, and obtain the fracture stress characteristic parameters;

[0074] Based on the bone impact energy distribution data, call the Young's modulus, yield stress, and trabecular bone structure parameters of the bone tissue, calculate the stress concentration area, fracture probability, and fracture morphology during the force application process of the bone tissue, call the VonMises stress criterion, calculate the maximum principal stress distribution inside the bone tissue, and judge the possibility of fracture according to the bone yield limit, and set the yield stress threshold σ y , when the local stress σ max ≥σ y , it is judged that a fracture occurs at this part, calculate the fracture location, fracture angle, and displacement of the fracture ends, combine the spatial arrangement characteristics of the trabecular bone, calculate the distribution of trabecular bone fracture points after fracture, identify the fracture mode, analyze the fracture propagation trend according to the stress propagation path inside the bone tissue, call the fracture toughness coefficient to calculate the smoothness of the fracture surface, and use the stress intensity factor K IC to calculate the fracture propagation rate, judge the fracture propagation range, call the bone density data, analyze the area where the bone density is lower than the critical value, calibrate the fracture location, calculate the damage range of the bone tissue, and classify crack fractures, comminuted fractures, and complete fractures in combination with parameters such as stress gradient, local stress concentration, and trabecular bone arrangement.

[0075] The VonMises stress calculation formula is as follows:

[0076]

[0077] Among them, σ v is the VonMises stress (Pa), σ1 is the first principal stress (Pa), σ2 is the second principal stress (Pa), and σ3 is the third principal stress (Pa).

[0078] Assume that the principal stresses at a certain bone tissue are σ1 = 120 MPa, σ2 = 80 MPa, and σ3 = 40 MPa, then the VonMises stress at this point is calculated as follows:

[0079]

[0080] If the material yield stress σ y = 60 MPa, then since σ v >σ y , a fracture will occur here. Finally, obtain the fracture stress characteristic parameters.

[0081] S103: According to the fracture stress characteristic parameters, combined with the mechanical response of the bone tissue, simulate various degrees of fracture injuries, including crack fractures, comminuted fractures, and complete fractures, and generate damage simulation data;

[0082] Based on the fracture stress characteristic parameters, combined with the mechanical response characteristics of bone tissue, simulate the biomechanical changes after fracture, call the force redistribution model in the fracture area, calculate the mechanical adjustment of the stressed part after fracture, adopt the force calculation equation in the fracture healing stage, analyze the adjustment of stress distribution during the fracture healing process, calculate the relative displacement of the fracture ends, the contact area of the fracture surface, and the stress transfer efficiency of the fracture site, call the intramedullary pressure parameter, calculate the blood flow impact of the fracture on the surrounding tissues, analyze the microcirculation changes and the expansion trend of the ischemic area at the fracture site, combine the fracture morphology and stress state, simulate the trauma states of fissure fractures, comminuted fractures, and complete fractures, and use three-dimensional fracture model visualization to establish a fracture injury scenario simulation.

[0083] The formula for calculating the relative displacement of the fracture ends is as follows:

[0084]

[0085] Among them, ΔL is the relative displacement of the fracture ends (mm), F is the force at the fracture (N), L is the length of the fracture segment (mm), A is the cross-sectional area of the fracture (mm 2 ), and E is the Young's modulus of bone tissue (MPa).

[0086] Assume that the force at the fracture F = 500 N, the length of the fracture segment L = 50 mm, the cross-sectional area of the fracture A = 100 mm 2 , and the Young's modulus E = 15000 MPa. Then the relative displacement of the fracture ends is calculated as follows:

[0087]

[0088] Finally, generate injury simulation data.

[0089] Based on the injury simulation data, simulate the blood flow changes at the fracture trauma site, calculate the probabilities of local ischemia and blood interruption, and map the trauma complications to the dummy model. The specific steps to obtain the injury mapping parameters are as follows:

[0090] S201: Based on the injury simulation data, call the fracture morphology, tissue injury area, and bone stress distribution at the fracture site, calculate the degree of blood vessel damage, the amount of blood vessel diameter narrowing, and the local blood flow velocity decrease rate in the fracture trauma area, combine the blood flow resistance parameter and the local venous return rate, calculate the change trend of blood flow supply at the trauma site, and obtain the dynamic blood flow characteristics of the trauma site;

[0091] Based on the damage simulation data, call the fracture morphology, tissue damage area, and bone stress distribution of the fracture site, calculate the degree of blood vessel damage, the amount of blood vessel diameter narrowing, and the local blood flow velocity decrease rate in the fracture trauma area. Combine the blood flow resistance parameters and the local venous return rate to calculate the changing trend of blood flow supply at the trauma site. For the fracture morphology of the fracture site, detect the morphological characteristics of the fracture ends, including the fracture angle, crack width, and fracture surface roughness, calculate the shear force of the fracture ends on the surrounding tissues, and calculate the degree of microvascular rupture in this area according to the contact stress formula of the fracture ends. Call the blood vessel compliance parameter to calculate the deformation amount after blood vessel damage, and combine the blood vessel wall thickness and blood flow shear force to judge the possibility of blood vessel damage.

[0092] The calculation formula for the amount of blood vessel diameter narrowing is as follows:

[0093] ΔD = D0 - D f ;

[0094] where ΔD is the amount of blood vessel diameter narrowing (mm), D0 is the original blood vessel diameter (mm), and D f is the blood vessel diameter after damage (mm).

[0095] Suppose the original diameter of a certain blood vessel D0 = 2.5 mm, and the blood vessel diameter after damage D f = 1.8 mm, then the calculation of the amount of blood vessel diameter narrowing is as follows:

[0096] ΔD = 2.5 - 1.8 = 0.7 mm;

[0097] According to the calculated blood vessel narrowing situation, combine the local blood flow velocity decrease rate, call the Poiseuille equation to calculate the blood flow velocity change, analyze the local tissue blood supply situation, combine the blood vessel resistance coefficient and the venous return rate to calculate the changing trend of blood flow supply, and finally obtain the blood flow dynamic characteristics of the trauma site.

[0098] S202: According to the blood flow dynamic characteristics of the trauma site, calculate the local tissue perfusion reduction rate, the probability of microvascular occlusion, and the decrease amplitude of oxygen supply level. Combine the tissue hypoxia threshold and the blood perfusion recovery rate to judge the ischemic influence range of the fracture area, calculate the local blood flow interruption risk level, and obtain the local ischemia and blood flow interruption parameters;

[0099] According to the blood flow dynamic characteristics of the trauma site, calculate the perfusion reduction rate of local tissues, the probability of microvascular occlusion, and the decline amplitude of oxygen supply level. Combine the tissue hypoxia threshold and the blood perfusion recovery rate to judge the ischemic influence range of the fracture area, and calculate the risk level of local blood flow interruption. Based on the blood flow dynamic characteristics of the trauma site, call the local blood flow pressure distribution data, calculate the blood flow resistance values of each damaged area, use the local perfusion index calculation formula to calculate the blood perfusion level of tissues per unit time, and judge the degree of tissue ischemia by comparing with the normal tissue perfusion threshold.

[0100] The calculation formula for the tissue perfusion reduction rate is as follows:

[0101]

[0102] Among them, R p is the tissue perfusion reduction rate (%), Q0 is the blood flow of normal tissue (mL / min), and Q f is the blood flow of the tissue after damage (mL / min).

[0103] Assume that the blood flow of normal tissue Q0 = 50 mL / min and the blood flow after damage Q f = 30 mL / min, then the calculation is as follows:

[0104]

[0105] Combined with the calculation results, call the tissue hypoxia threshold to judge whether the area enters the hypoxic state, and combine the blood perfusion recovery rate to calculate the possible recovery time, and finally obtain the local ischemia and blood flow interruption parameters.

[0106] S203: Call the local ischemia and blood flow interruption parameters, map the ischemic area, blood flow interruption point, and trauma complication category to the dummy model to generate injury mapping parameters;

[0107] Call the local ischemia and blood flow interruption parameters, map the ischemic area, blood flow interruption point, and trauma complication category to the dummy model, and combine the degree of soft tissue injury and the blood vessel function attenuation rate to adjust the tissue color change, subcutaneous blood flow simulation, and trauma deterioration process of the dummy model. Based on the ischemic area data, set the blood flow dynamic display parameters of the dummy model, call the three-dimensional blood vessel network data, simulate the expansion of local ischemic points, and combine the blood flow interruption point to dynamically adjust the local color change of the dummy model to reflect the influence range of tissue ischemia.

[0108] The calculation formula for the ischemic area mapping ratio is as follows:

[0109]

[0110] Among them, R m is the ischemic area mapping ratio (%), Ai is the area of the ischemic region in the dummy model (cm 2 ), and A t is the total area of the overall simulation region (cm 2 ).

[0111] Assume that the total simulation region area A of the dummy model t = 200 cm 2 , and the ischemic region area A i = 50 cm 2 , then the calculation is as follows:

[0112]

[0113] Combined with the calculation results, simulate the local tissue color change trend, adjust the ischemic display layer of the dummy model, and finally generate the injury mapping parameters.

[0114] The steps of calling the injury mapping parameters, collecting the joint movement data of the trainees in real time, obtaining the angular change data of multiple joints, calculating the movement trajectory and acceleration of the joints, and obtaining the action deviation data by comparing with the standard treatment actions are as follows:

[0115] S301: Call the injury mapping parameters, collect the movement data of multiple joint parts of the trainees during the treatment process in real time, record the joint angle changes, rotation rates, and acceleration characteristics of the trainees during the treatment operation, calculate the movement amplitude, angle change trend, and speed change rate of each joint, and obtain the joint movement data of the trainees.

[0116] Call the injury mapping parameters, collect the movement data of multiple joint parts of the trainees during the treatment process in real time, record the joint angle changes, rotation rates, and acceleration characteristics of the trainees during the treatment operation, calculate the movement amplitude, angle change trend, and speed change rate of each joint. Collect the real-time coordinate data of key joint points such as the elbow joint, wrist joint, shoulder joint, and knee joint, call the inertial measurement unit (IMU) sensor to measure the three-axis angular velocity and acceleration data of each joint, calculate the instantaneous angular change amount of each joint, and based on the time step Δt, solve the rotation rate change trend of each joint, and call the angular velocity calculation formula to calculate the angular change situation within a unit time.

[0117] The angular velocity calculation formula is as follows:

[0118]

[0119] Among them, ω is the angular velocity (rad / s), θ f is the final angle (°), θ i is the initial angle (°), and Δt is the time interval (s).

[0120] Suppose the initial angle θ of a trainee's elbow joint i = 30°, and after 0.5 seconds, the angle becomes θ f = 90°. Then the angular velocity is calculated as follows:

[0121]

[0122] Combined with the calculation results, further analyze the movement amplitude, angle change trend, and speed change rate of each joint, and finally obtain the trainee's joint movement data.

[0123] S302: Based on the trainee's joint movement data, calculate the movement trajectory, path curvature, and key time nodes of the rescue action, evaluate the movement continuity, angle stability, and force application direction of each joint during the rescue operation, and obtain the trajectory characteristics of the rescue action;

[0124] Based on the trainee's joint movement data, calculate the movement trajectory, path curvature, and key time nodes of the rescue action, evaluate the movement continuity, angle stability, and force application direction of each joint during the rescue operation. Call the trajectory path calculation formula to calculate the movement trajectories of the arms and legs during the trainee's rescue process, identify the key points of each joint during the movement process, combine the three-dimensional coordinate conversion data, calculate the spatial curvature of the trajectory, call the joint angle stability judgment formula, calculate the angle change rate within different time steps, judge whether there are discontinuous changes in the action, and analyze the stability of the wrists and shoulders during the rescue process.

[0125] The path curvature calculation formula is as follows:

[0126]

[0127] where κ is the curvature, x', y' are the first derivatives of the trajectory, and x″, y″ are the second derivatives of the trajectory.

[0128] Suppose during the operation of a certain trainee, the derivatives of the arm movement trajectory are x' = 2, y' = 3, and the second derivatives are x″ = -1, y″ = 2. Then the curvature is calculated as follows:

[0129]

[0130] Combined with the calculation results, judge whether there are abrupt posture changes during the trainee's rescue process, and finally obtain the trajectory characteristics of the rescue action.

[0131] S303: Call the trajectory characteristics of the rescue action, calculate the trajectory error of the trainee's operation by comparing with the standard rescue action, identify and record the abnormal postures during the rescue process, and generate action deviation data;

[0132] The specific formula for calculating the trajectory error of the trainee's operation is:

[0133]

[0134] Calculate the trajectory error value, identify and record the abnormal postures during the treatment process, and generate action deviation data;

[0135] Among them, E traj is the trajectory error value, N' is the total number of trajectory sampling points, i' is the trajectory sampling point index, X i' is the X-axis coordinate value of the trainee at the i'-th trajectory sampling point, Y i' is the Y-axis coordinate value of the trainee at the i'-th trajectory sampling point, z i' is the Z-axis coordinate value of the trainee at the i'-th trajectory sampling point, is the X-axis coordinate value of the standard treatment action at the i'-th trajectory sampling point, is the Y-axis coordinate value of the standard treatment action at the i'-th trajectory sampling point, is the Z-axis coordinate value of the standard treatment action at the i'-th trajectory sampling point, M' is the total number of speed samplings during the treatment process, j' is the speed sampling point index, V j' is the hand movement speed of the trainee at the j'-th sampling, V std is the average speed of the standard treatment operation.

[0136] Formula:

[0137]

[0138] Detailed explanation of the formula and the derivation process of the formula calculation:

[0139] The formula is used to calculate the trajectory error value, and the obtained result is used to evaluate the stability of the trainee's treatment operation and the deviation degree from the standard action;

[0140] Parameter meanings and setting values:

[0141] N' is the total number of trajectory sampling points, and the setting value is 3, which reflects the number of coordinate points recorded during the treatment operation;

[0142] i' is the trajectory sampling point index, indicating the i'-th sampling point;

[0143] X i' is the X-axis coordinate value of the trainee at the i'-th trajectory sampling point, Y i' is the Y-axis coordinate value of the trainee at the i'-th trajectory sampling point, z i' is the Z-axis coordinate value of the trainee at the i'-th trajectory sampling point;

[0144] is the X-axis coordinate value of the standard treatment action at the i'-th trajectory sampling point, is the Y-axis coordinate value of the standard treatment action at the i'-th trajectory sampling point, is the Z-axis coordinate value of the standard treatment action at the i'-th trajectory sampling point;

[0145] M' is the total number of speed samplings during the treatment process, with a set value of 2, reflecting the number of speed measurement points recorded during the sampling process;

[0146] j' is the speed sampling point index, representing the j'-th speed measurement point;

[0147] V j' is the hand movement speed of the trainee during the j'-th sampling, in meters per second;

[0148] V std is the average speed of the standard treatment operation, with a set value of 0.8 m / s;

[0149] Assume that the data of three trajectory sampling points are (X1, Y1, z1) = (0.5, 1.2, 0.3), (X2, Y2, z2) = (0.6, 1.4, 0.5), (X3, Y3, z3) = (0.7, 1.6, 0.6), The speed data are V1 = 0.75 m / s, V2 = 0.85 m / s;

[0150] Substitute the parameters into the formula for calculation:

[0151] Calculate the error of the first sampling point:

[0152]

[0153] Calculate the error of the second sampling point:

[0154]

[0155] Calculate the error of the third sampling point:

[0156]

[0157] Calculate the mean value of the trajectory error:

[0158]

[0159] Calculate the mean value of the speed:

[0160]

[0161] Calculate the speed deviation:

[0162] |0.8 - 0.8| = 0;

[0163] Calculate the error value:

[0164] E traj = 0.1152 + 0 = 0.1152;

[0165] The result 0.1152 indicates that the average error between the trainee's operation trajectory and the standard trajectory is 0.1152. The higher the value, the more the trainee's operation trajectory deviates from the standard trajectory. The calculation process is used to evaluate the operation stability and accuracy of the trainee.

[0166] Based on the action deviation data, combined with the mechanical response of the fracture site, analyze the impact of multiple trainees' treatment actions on the trauma site, update the injury progress in real time, and adjust the dummy mapping parameters. The steps to obtain the injury update data are as follows:

[0167] S401: Based on the action deviation data, combined with the stress state, displacement trend of the fracture end, and tissue stress distribution of the fracture site, calculate the impact of multiple trainees' treatment actions on the trauma site, including local stress changes, fracture end dislocation amplitude, and tissue tensile deformation rate, and obtain the treatment operation impact parameters;

[0168] Based on the action deviation data, combined with the stress state, displacement trend of the fracture end, and tissue stress distribution of the fracture site, calculate the impact of multiple trainees' treatment actions on the trauma site, including local stress changes, fracture end dislocation amplitude, and tissue tensile deformation rate. Call the joint movement trajectory data to obtain the angle changes of key joints such as the arm, knee, and elbow during the trainee's rescue process, calculate the direction of the external force applied by the trainee, calculate the stress change of the damaged bone during the treatment process according to the fracture end force calculation formula, set the contact force distribution of the fracture end, combine the moment calculation in the fracture area, analyze the secondary force on the fractured bone when the trainee rescues, call the fracture displacement calculation formula to calculate the relative displacement of the fracture end, combine the elastic modulus of the tissue, analyze the tensile strain of the treatment action on the trauma area, and judge whether the tissue enters the stress overload state.

[0169] The fracture end dislocation calculation formula is as follows:

[0170]

[0171] Among them, ΔL is the fracture end dislocation amount (mm), F is the external force applied to the fracture (N), L is the length of the stress-bearing area (mm), A is the stress-bearing area (mm 2 ), and E is the Young's modulus of bone tissue (MPa).

[0172] Assume that during a trainee's rescue process, the external force F = 300N, the fracture area length L = 40mm, the stress-bearing area A = 80mm 2 , and the Young's modulus E = 15000MPa, then the calculation is as follows:

[0173]

[0174] Combined with the calculation results, analyze the external forces during the treatment of the trainee, and finally obtain the influencing parameters of the treatment operation.

[0175] S402: Invoke the influencing parameters of the treatment operation, analyze the tissue stress concentration area, the degree of blood flow supply obstruction, and the trend of nerve injury spread at the fracture site, calculate the risk level of injury aggravation, the recovery delay rate, and the local blood flow supply adjustment range in the trauma area, and obtain the injury progression status;

[0176] Invoke the influencing parameters of the treatment operation, analyze the tissue stress concentration area, the degree of blood flow supply obstruction, and the trend of nerve injury spread at the fracture site, calculate the risk level of injury aggravation, the recovery delay rate, and the local blood flow supply adjustment range in the trauma area. Combine the mechanical data applied during the trainee's rescue process, calculate the stress concentration points in the fracture area, and based on the maximum principal stress criterion, analyze whether the tissue enters the ultimate stress state during the treatment. Invoke the hemodynamic parameters, calculate the degree of compression of local blood vessels, combine with the local perfusion pressure calculation formula, analyze the reduction in blood supply at the fracture site, invoke the nerve injury spread calculation formula, evaluate the injury progression in the damaged nerve area, set the recovery delay rate calculation formula, and analyze the degree of influence on the recovery process during the treatment.

[0177] The calculation formula for the local blood flow supply adjustment range is as follows:

[0178]

[0179] Where, R b is the local blood flow supply adjustment range (%), AQ is the blood flow before fracture (mL / min), and AQ f is the blood flow after rescue (mL / min).

[0180] Assume that the original blood flow in the fracture area is 45 mL / min and the blood flow after rescue is 30 mL / min, then the calculation is as follows:

[0181]

[0182] Combined with the calculation results, judge the influence on blood flow supply during the trainee's treatment, and finally obtain the injury progression status.

[0183] S403: Invoke the injury progression status, and adjust the mapping parameters of the dummy model in real time, including dynamic injury display, tissue necrosis process, and blood circulation recovery trend, to generate injury update data;

[0184] Call the injury progress status and adjust the mapping parameters of the dummy model in real time, including dynamic injury display, tissue necrosis process, and blood circulation recovery trend. According to the calculated injury aggravation risk level, set the trauma display parameters of the dummy model, call the dynamic tissue color mapping algorithm to simulate the color change of the ischemic area, combine the blood circulation recovery trend, adjust the blood perfusion rate of the dummy model, calculate the expansion rate of the ischemic area, call the necrosis tissue progress calculation formula, calculate the spatial distribution range of tissue necrosis, set the necrosis progress level of the dummy model, adjust the tissue damage diffusion rate in the simulation environment, and update the injury simulation data of the dummy model.

[0185] The tissue necrosis progress calculation formula is as follows:

[0186] A n = A0 × e -λt ;

[0187] where, A n is the area of necrotic tissue at a certain moment (cm 2 ), A0 is the initial area of necrotic tissue (cm 2 ), λ is the necrosis rate constant (h -1 ), and t is the time (h).

[0188] Assume that the initial area of necrotic tissue A0 = 10 cm 2 , the necrosis rate constant λ = 0.2 h -1 , after 5 hours, the calculation is as follows:

[0189] A n = 10 × e -0.2×5 = 10 × e -1 ≈ 10 × 0.368 = 3.68 cm 2 ;

[0190] Combined with the calculation results, adjust the simulation parameters of the necrotic tissue expansion of the dummy model, and finally generate the injury update data.

[0191] According to the injury update data, after the treatment is completed, analyze the operation accuracy of each action of the trainee, combine the treatment duration and treatment effect, and calculate the trainee's skill score. The specific steps to obtain the treatment skill score are as follows:

[0192] S501: Call the injury update data, calculate the operation accuracy and key action matching rate of the trainee in the treatment task by analyzing the adjustment range of joint angles, the deviation of force application direction, and the stability of the treatment posture during the trainee's treatment operation, and obtain the treatment operation accuracy data;

[0193] Call the injury update data. By analyzing the adjustment range of joint angles, the deviation of force application direction, and the stability of the treatment posture in the trainees' treatment operations, calculate the operation accuracy and the matching rate of key actions of the trainees in the treatment tasks. Obtain the angle change curve of each joint during the treatment process, calculate the angle adjustment range of the elbow joint, knee joint, and shoulder joint during the trainees' rescue process, call the standard action angle data, calculate the angle error of each joint during the trainees' treatment process, combine the angle deviation calculation formula, solve the deviation degree of each key action, set the allowable error range of the key action, compare the matching degree between the actual operation and the standard action, further calculate the operation accuracy of the trainees, and evaluate whether there are unstable postures during the trainees' treatment process.

[0194] The angle deviation calculation formula is as follows:

[0195]

[0196] Among them, D θ is the average angle deviation (°), θ i is the actual operation angle of the trainee (°), θ i * is the standard action angle (°), and n is the number of sampling points.

[0197] Assume that during the treatment process of the trainee, the elbow joint angles of 3 key points are collected. The trainee's angles are: 45°, 60°, 75°, and the standard angles are: 50°, 65°, 80°. Calculate the angle deviation:

[0198]

[0199] Combined with the calculation results, judge whether the trainee's actions exceed the error threshold, and finally obtain the treatment operation accuracy data.

[0200] S502: Based on the treatment operation accuracy data, combined with the treatment completion time of the trainee, calculate the average treatment duration and emergency handling efficiency of the trainee in multiple treatment tasks, and obtain the treatment process efficiency parameters;

[0201] Based on the treatment operation accuracy data, combined with the treatment completion time of the trainee, calculate the average treatment duration and emergency handling efficiency of the trainee in multiple treatment tasks. Call the task completion time of the trainee in different treatment tasks, calculate the execution duration of each action during the treatment process, set the standard time range of the treatment task, calculate the time deviation rate of each treatment, call the average treatment duration calculation formula, calculate the time mean of multiple tasks, and combine the time distribution of key nodes during the treatment process to calculate the trainee's emergency response ability. Set the emergency handling efficiency calculation formula to evaluate the trainee's reaction speed in case of emergencies.

[0202] The average treatment duration calculation formula is as follows:

[0203]

[0204] Among them, T avg is the average treatment duration (s), T i is the completion time (s) of the i-th treatment task, and m is the total number of tasks.

[0205] Suppose the completion times of the trainee in 3 treatment tasks are 120 s, 135 s, and 150 s respectively, then the calculation is as follows:

[0206]

[0207] Combined with the calculation results, evaluate whether the trainee's treatment time meets the standard, and finally obtain the efficiency parameter of the treatment process.

[0208] S503: Invoke the efficiency parameter of the treatment process, combine the wound repair status, tissue recovery rate, and hemostasis stability of the injury progress during the treatment process, obtain the treatment effect score of the trainee, combine the operation accuracy and treatment efficiency, calculate the trainee's skill score, and generate the treatment skill score;

[0209] The specific formula for calculating the trainee's skill score is:

[0210]

[0211] Calculate the treatment skill score;

[0212] Among them, S skill represents the trainee's skill score, E treat represents the treatment effect score, P acc represents the operation accuracy, T delay represents the delay time caused by stress or inexperience during the treatment process, W eff represents the treatment efficiency, R fatigue represents the correction value for the treatment performance affected by physiological fatigue.

[0213] Formula:

[0214]

[0215] Detailed explanation of the formula and the derivation process of the formula calculation:

[0216] The formula is used to calculate the trainee's skill score, and the obtained result is used to evaluate the trainee's overall ability in the treatment task, including treatment effect, operation accuracy, time management, and physiological tolerance.

[0217] Parameter meaning and setting value:

[0218] E treatRepresents the treatment effect score, which reflects the quality of the trainee's handling during the treatment process, including the correctness, stability, and comprehensive performance of operations such as wound dressing, hemostasis, and fixation. The set value is 85.

[0219] P acc Represents the operation accuracy, which reflects the precision of the trainee when performing the treatment task. The set value is 0.92.

[0220] T delay Represents the delay time caused by stress or inexperience during the treatment process, which reflects the additional time consumption of the trainee due to psychological pressure, slow thinking, or inexperience during the treatment process. The set value is 6 seconds.

[0221] W eff Represents the treatment efficiency, which reflects the trainee's ability to complete the treatment task per unit time. Set W eff =0.15

[0222] R fatigue Represents the correction value of the treatment performance affected by physiological fatigue, which reflects the impact of the trainee's physiological fatigue on the operation stability. The set value is 0.08.

[0223] Substitute the parameters into the formula for calculation:

[0224]

[0225] S skill =11.17 + 0.07;

[0226] S skill =11.24;

[0227] The calculation result of 11.24 indicates the comprehensive skill score of the trainee in this treatment training. Combining the treatment effect, operation accuracy, and efficiency, and considering the delay factor and the impact of physiological fatigue, a higher score indicates that the trainee has stronger treatment ability in a high-pressure environment. If the score is low, training optimization is required for treatment accuracy, time management, or physiological tolerance.

[0228] Please refer to Figure 2 , a high-risk orthopedic trauma simulation frontline treatment training system, including:

[0229] The impact injury simulation module calculates the stress and strain effects of the shock wave on bone tissue under various pressure gradients, wave speed changes, and shock durations based on preset explosion shock wave parameters. Combining the mechanical response characteristics of bone tissue, it analyzes the injury types of bones in multiple parts, including crack fractures and comminuted fractures, and obtains injury simulation data;

[0230] The complication mapping module simulates and analyzes the changes in the blood flow state at the trauma site based on the injury simulation data, predicts the occurrence probabilities of various complications, maps them to the dummy model, displays the pathological progress of the trauma site in real time, and obtains the injury mapping parameters;

[0231] The treatment action capture module collects the joint angle change data of the trainee during the treatment operation in real time based on the injury mapping parameters, calculates the movement trajectories and accelerations of multiple joints, and compares them with the standard treatment action path to obtain the action deviation data;

[0232] The treatment impact analysis module analyzes the impact of the trainee's treatment actions on the trauma area based on the action deviation data, combined with the mechanical response of the fracture site, adjusts the trauma mapping state of the dummy model, displays the dynamic changes of the injury in real time, and obtains the injury update data;

[0233] The trainee skill evaluation module analyzes the accuracy and standardization of each treatment action of the trainee after the treatment operation based on the injury update data, and combines the treatment duration and the injury change situation to calculate the treatment ability score of the trainee and obtain the treatment skill score.

[0234] The above embodiments can be implemented in whole or in part by software, hardware (such as circuits), firmware, or any other combination. When implemented using software, the above embodiments can be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions or computer programs. When the computer instructions or computer programs are loaded or executed on a computer, the processes or functions according to the embodiments of the present invention are generated in whole or in part. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices. The computer instructions can be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another. For example, the computer instructions can be transmitted from one website, computer, server, or data center to another website, computer, server, or data center by wired (such as infrared, wireless, microwave, etc.) means. The computer-readable storage medium can be any available medium that the computer can access or a data storage device such as a server or data center that contains one or more collections of available media. The available media can be magnetic media (such as floppy disks, hard disks, magnetic tapes), optical media (such as DVDs), or semiconductor media. The semiconductor media can be a solid-state drive.

[0235] It should be understood that the term "and / or" in this text is merely a description of the association relationship between associated objects, indicating that three relationships can exist. For example, A and / or B can represent three situations: A exists alone, A and B exist simultaneously, and B exists alone. Here, A and B can be singular or plural. Additionally, the character " / " in this text generally represents an "or" relationship between the associated objects before and after, but it may also represent an "and / or" relationship. The specific meaning can be understood by referring to the context before and after.

[0236] In the present invention, "at least one" means one or more, and "a plurality" means two or more. "At least one of the following" or its similar expressions refer to any combination of these items, including any combination of single items or plural items. For example, at least one of a, b, or c can represent: a, b, c, a - b, a - c, b - c, or a - b - c, where a, b, and c can be single or multiple.

[0237] It should be understood that in various embodiments of the present invention, the magnitudes of the serial numbers of the above processes do not imply the order of execution. The order of execution of each process should be determined by its function and internal logic, and should not constitute any limitation to the implementation process of the embodiments of the present invention.

[0238] Those of ordinary skill in the art can realize that the units and algorithm steps of each example described in combination with the embodiments disclosed in this text can be implemented by electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are executed in a hardware or software manner depends on the specific application and design constraints of the technical solution. Professional technicians can use different methods to implement the described functions for each specific application, but such implementation should not be considered to exceed the scope of the present invention.

[0239] Those skilled in the art can clearly understand that for the convenience and brevity of description, the specific working processes of the devices, apparatuses, and units described above can refer to the corresponding processes in the foregoing method embodiments, and will not be elaborated herein.

[0240] In several embodiments provided by the present invention, it should be understood that the disclosed devices, apparatuses, and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the division of the units is only a logical function division. In actual implementation, there can be other division methods. For example, multiple units or components can be combined or integrated into another device, or some features can be ignored or not executed. Another point is that the displayed or discussed coupling or direct coupling or communication connection between each other can be through some interfaces. The indirect coupling or communication connection of the devices or units can be in an electrical, mechanical, or other form.

[0241] The unit described as a separation component may or may not be physically separated. The component shown as a unit may or may not be a physical unit, that is, it may be located in one place or may be distributed to multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0242] In addition, each functional unit in various embodiments of the present invention may be integrated in a processing unit, may exist physically separately for each unit, or two or more units may be integrated in one unit.

[0243] If the above-mentioned function is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or a part of this technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical discs that can store program codes.

[0244] As described above, the above are only specific embodiments of the present invention, but the protection scope of the present invention is not limited thereto. Any person skilled in the art within the technical scope disclosed by the present invention can easily think of changes or substitutions, which should all be covered by the protection scope of the present invention. Therefore, the protection scope of the present invention should be subject to the protection scope of the claims.

Claims

1. A training method for simulated on-the-spot treatment of high-risk orthopedic trauma, characterized in that, The method includes: S1: Invoke the preset blast shock wave parameters, combine with the stress-strain characteristics of bone tissue, calculate the impact of the shock wave on the bones of multiple parts, analyze the mechanical response of bone tissue, simulate fracture injuries, and obtain injury simulation data; S2: Based on the injury simulation data, simulate the blood flow changes at the fracture trauma site, calculate the probabilities of local ischemia and blood interruption, and map the trauma complications onto the dummy model to obtain injury mapping parameters; S3: Invoke the injury mapping parameters, collect the joint movement data of the trainees in real time, obtain the angular change data of multiple joints, calculate the movement trajectories and accelerations of the joints, and obtain action deviation data by comparing with the standard treatment actions; S4: Based on the action deviation data, combine with the mechanical response of the fracture site, analyze the impact of the treatment actions of multiple trainees on the trauma site, update the injury progress in real time and adjust the dummy mapping parameters to obtain injury update data; S5: According to the injury update data, analyze the operation accuracy of each action of the trainees after the treatment is completed, combine with the treatment duration and treatment effect, calculate the trainee skill score, and obtain the treatment skill score.

2. The high-risk orthopedic trauma simulation on-site treatment training method according to claim 1, wherein The injury simulation data specifically includes fracture type, bone tissue stress distribution, and fracture end displacement. The injury mapping parameters include local ischemia range, blood flow blockage probability, and complication type. The action deviation data specifically refers to joint angle error, movement trajectory deviation, and action stability score. The injury update data includes the degree of trauma deterioration, tissue repair rate, and dynamic injury status of the dummy model. The treatment skill score specifically refers to operation precision score, treatment efficiency score, and trauma recovery effect score.

3. The high-risk orthopedic trauma simulation on-the-spot treatment training method according to claim 1, characterized in that The steps of invoking the preset blast shock wave parameters, combining with the stress-strain characteristics of bone tissue, calculating the impact of the shock wave on the bones of multiple parts, analyzing the mechanical response of bone tissue, simulating fracture injuries, and obtaining injury simulation data are specifically as follows: S101: Obtain the preset blast shock wave parameters, invoke the data of blast wave pressure gradient, shock propagation speed, and shock duration, calculate the force conditions of the blast shock wave at multiple bone sites, detect the stress distribution, deformation trend, and force area of bone tissue under the action of the shock wave, analyze the attenuation law of shock wave energy in bone tissue, and obtain the bone shock energy distribution data; S102: Based on the bone shock energy distribution data, combine with the Young's modulus, yield stress, and trabecular bone structure parameters of bone tissue, calculate the stress concentration area, fracture occurrence probability, and fracture morphology during the force application process of bone tissue, identify the fracture occurrence location, and obtain the fracture stress characteristic parameters; S103: According to the fracture stress characteristic parameters, combine with the mechanical response of bone tissue, simulate various degrees of fracture injuries, including fissure fractures, comminuted fractures, and complete fractures, and generate injury simulation data.

4. The high-risk orthopedic trauma simulation live-fire treatment training method according to claim 1, wherein, The steps of simulating the blood flow changes at the fracture trauma site based on the injury simulation data, calculating the probabilities of local ischemia and blood interruption, and mapping the trauma complications onto the dummy model to obtain injury mapping parameters are specifically as follows: S201: Based on the injury simulation data, call the fracture morphology, tissue injury area, and bone stress distribution of the fracture site, calculate the degree of blood vessel damage, the amount of blood vessel diameter narrowing, and the local blood flow velocity decrease rate in the fracture trauma area, combine the blood flow resistance parameter and the local venous return rate, calculate the changing trend of blood flow supply at the trauma site, and obtain the dynamic blood flow characteristics of the trauma site; S202: According to the dynamic blood flow characteristics of the trauma site, calculate the local tissue perfusion reduction rate, the probability of microvascular occlusion, and the decrease amplitude of oxygen supply level, combine the tissue hypoxia threshold and the blood perfusion recovery rate, judge the ischemic influence range of the fracture area, calculate the local blood flow interruption risk level, and obtain the local ischemia and blood flow blockage parameters; S203: Call the local ischemia and blood flow blockage parameters, map the ischemic area, blood flow interruption point, and trauma complication category to the dummy model to generate injury mapping parameters.

5. The high-risk orthopedic trauma simulation first-aid training method according to claim 1, wherein The steps of calling the injury mapping parameters, collecting the joint movement data of the trainee in real time, obtaining the angular change data of multiple joints, calculating the movement trajectory and acceleration of the joints, and obtaining the action deviation data by comparing with the standard treatment actions are as follows: S301: Call the injury mapping parameters, collect the movement data of multiple joint parts of the trainee during the treatment process in real time, record the joint angle changes, rotation rates, and acceleration characteristics of the trainee during the treatment operation, calculate the movement amplitude, angular change trend, and speed change rate of each joint, and obtain the trainee joint movement data; S302: Based on the trainee joint movement data, calculate the movement trajectory, path curvature, and key time nodes of the treatment action, evaluate the movement continuity, angular stability, and force application direction of each joint during the treatment operation, and obtain the treatment action trajectory characteristics; S303: Call the treatment action trajectory characteristics, calculate the trajectory error of the trainee's operation by comparing with the standard treatment action, identify and record the abnormal postures existing during the treatment process, and generate action deviation data.

6. The high-risk orthopedic trauma simulation live-fire treatment training method according to claim 5, wherein The specific formula for calculating the trajectory error of the trainee's operation is: Calculate the trajectory error value, identify and record the abnormal postures existing during the treatment process, and generate action deviation data; Among them, E traj is the trajectory error value, N' is the total number of trajectory sampling points, i' is the trajectory sampling point index, X i' is the X-axis coordinate value of the trainee at the i'-th trajectory sampling point, Y i' is the Y-axis coordinate value of the trainee at the i'-th trajectory sampling point, z i' is the Z-axis coordinate value of the trainee at the i'-th trajectory sampling point, is the X-axis coordinate value of the standard rescue action at the i'-th trajectory sampling point, is the Y-axis coordinate value of the standard rescue action at the i'-th trajectory sampling point, is the Z-axis coordinate value of the standard rescue action at the i'-th trajectory sampling point, M' is the total number of speed samplings during the rescue process, j' is the speed sampling point index, V j' is the hand movement speed of the trainee at the j'-th sampling, V std is the average speed of the standard rescue operation.

7. The high-risk orthopedic trauma simulation live-fire treatment training method according to claim 1, characterized in that The steps of analyzing the influence of the treatment actions of multiple trainees on the trauma site based on the action deviation data, combining the mechanical response of the fracture site, updating the injury progress in real time, and adjusting the dummy mapping parameters to obtain the injury update data are as follows: S401: Based on the action deviation data, combine the stress state, displacement trend of the fracture end, and tissue stress distribution of the fracture site, calculate the influence of the treatment actions of multiple trainees on the trauma site, including local stress changes, the dislocation amplitude of the fracture end, and the tissue stretching deformation rate, and obtain the treatment operation influence parameters; S402: Call the treatment operation influence parameters, analyze the tissue stress concentration area, blood flow supply obstruction degree, and nerve injury diffusion trend of the fracture site, calculate the injury aggravation risk level, recovery delay rate, and local blood flow supply adjustment range of the trauma area, and obtain the injury progress status; S403: Call the injury progress status, and adjust the mapping parameters of the dummy model in real time, including dynamic injury display, tissue necrosis process, and blood circulation recovery trend, to generate injury update data.

8. The high-risk orthopedic trauma simulation live-fire treatment training method according to claim 1, wherein, According to the injury update data, after the treatment is completed, analyze the operation accuracy of each action of the trainee, and calculate the trainee's skill score by combining the treatment duration and treatment effect. The steps for obtaining the treatment skill score are specifically as follows: S501: Call the injury update data, and calculate the operation accuracy and key action matching rate of the trainee in the treatment task by analyzing the adjustment range of joint angles, the deviation of force application direction, and the stability of treatment posture during the trainee's treatment operation, so as to obtain treatment operation accuracy data; S502: Based on the treatment operation accuracy data, combine the trainee's treatment completion time, and calculate the average treatment duration and emergency handling efficiency of the trainee in multiple treatment tasks to obtain treatment process efficiency parameters; S503: Call the treatment process efficiency parameters, combine the trauma repair status, tissue recovery rate, and hemostasis stability of the injury progress during the treatment process, obtain the trainee's treatment effect score, combine the operation accuracy and treatment efficiency, calculate the trainee's skill score, and generate a treatment skill score.

9. The high-risk orthopedic trauma simulation on-site treatment training method according to claim 8, characterized in that The specific formula for calculating the trainee's skill score is: Calculate the treatment skill score; Among them, S skill represents the trainee skill score, E treat represents the treatment effect score, P acc represents the operation accuracy, T delay represents the delay time caused by stress or inexperience during the treatment process, W eff represents the treatment efficiency, R fatigue represents the correction value of the treatment performance affected by physical fatigue.

10. A high-risk orthopedic trauma simulation and on-the-spot treatment training system, characterized in that, According to the high-risk orthopedic trauma simulation frontline treatment training method described in any one of claims 1-9, the system includes: The impact injury simulation module calculates the stress and strain effects of the shock wave on bone tissue under various pressure gradients, wave speed changes, and shock duration based on the preset explosion shock wave parameters, and analyzes the injury types of bones in multiple parts, including crack fractures and comminuted fractures, in combination with the mechanical response characteristics of bone tissue, to obtain injury simulation data; The complication mapping module simulates and analyzes the changes in the blood flow state of the trauma site based on the injury simulation data, predicts the occurrence probability of various complications, and maps them to the dummy model to display the pathological progress of the trauma site in real time, so as to obtain injury mapping parameters; The treatment action capture module collects the joint angle change data of the trainee during the treatment operation in real time based on the injury mapping parameters, calculates the movement trajectories and accelerations of multiple joints, and compares them with the standard treatment action path to obtain action deviation data; The treatment impact analysis module analyzes the impact of the trainee's treatment action on the trauma area based on the action deviation data, in combination with the mechanical response of the fracture site, adjusts the trauma mapping state of the dummy model, and displays the dynamic changes of the injury in real time to obtain injury update data; The trainee skill evaluation module analyzes the accuracy and standardization of each treatment action of the trainee after the treatment operation based on the injury update data, and calculates the trainee's treatment ability score by combining the treatment duration and the injury change situation to obtain a treatment skill score.

Images