A whole-body-local multi-scale biomechanical model coupling construction method for human body impact simulation
By constructing a whole-body dynamic model of the human body and performing equivalent load conversion, the problem of low computational efficiency in existing technologies is solved. This enables the effective transfer of the overall dynamic effects of the human body to the local simulation model, reducing the scale of the simulation model and computational complexity. It is suitable for human injury assessment and protective equipment design.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- JINAN UNIVERSITY
- Filing Date
- 2025-12-30
- Publication Date
- 2026-05-29
AI Technical Summary
Existing technologies struggle to simultaneously achieve both the overall dynamic effects of the human body and the accuracy of local structural mechanical responses while meeting the computational efficiency requirements of engineering applications.
A whole-body dynamic model of the human body is constructed, the time history data of constraint forces and constraint torques are extracted, and equivalent transformations are performed to generate equivalent loads, which are then applied to local simulation models to achieve effective transfer of the overall dynamic effects of the human body to local simulation models.
It reduces the size and computational complexity of simulation models, improves computational efficiency, and is suitable for engineering applications, including human injury assessment and protective equipment design.
Smart Images
Figure CN122113352A_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of human biomechanical simulation and impact analysis technology, and more specifically, to a method for constructing a whole-body-local multi-scale biomechanical model coupling for human impact simulation. Background Technology
[0002] Under impact conditions such as falls and collisions, local anatomical structures (such as the hip, head, and spine) are prone to bearing large transient loads, which can lead to fractures or soft tissue injuries. Therefore, conducting mechanical response analysis during human impact is of great significance for the study of human injury mechanisms, the design of protective equipment, and the assessment of human safety performance.
[0003] Existing human impact simulation methods mainly include whole-body motion analysis methods based on multibody dynamics models and detailed simulation methods of local human structures based on finite element models. Multibody dynamics models are typically used to describe the overall motion state of the human body and have high computational efficiency, but their accuracy is limited in analyzing detailed mechanical responses such as local stress, strain, and contact forces. While whole-body human simulation methods based on finite element models can accurately describe the local mechanical response of the human body during impact, their large model size, high computational complexity, and low simulation efficiency make them difficult to meet the computational efficiency requirements of engineering applications.
[0004] To balance the overall dynamic effects of the human body with the accuracy of local structural responses, some whole-body-local simulation methods have been proposed in the prior art. These methods achieve fine analysis of local regions by transferring data between the whole-body model and local models. However, these methods usually rely on the transfer of kinematic information such as displacement or velocity, making it difficult to accurately reflect the dynamic effects of the overall mass distribution and inertial effects on local regions; or they require the construction of high-precision local structures within the whole-body model, thus still resulting in problems such as large model size and high computational cost.
[0005] Therefore, how to effectively transfer the whole-body dynamic effects of human impact to the local simulation model while preserving the overall dynamic effects of the human body, and at the same time reduce the size and computational complexity of the simulation model and improve the computational efficiency of human impact simulation, remains a technical problem that urgently needs to be solved in this field. Summary of the Invention
[0006] Given that the existing human impact simulation methods are difficult to simultaneously consider the overall dynamic effects of the human body and the accuracy of local structural mechanical response while meeting the computational efficiency requirements of engineering applications, the purpose of this invention is to provide a method for coupling and constructing a whole-body-local multi-scale biomechanical model for human impact simulation.
[0007] The method includes: S1: Construct a full-body dynamic model of the human body; S2: Set the impact condition parameters of the whole human body dynamics model, simulate the motion process of the human body under impact, and obtain the overall kinematic and dynamic response of the human body; S3: In the whole-body dynamics model of the human body, determine the connection constraint positions corresponding to the target local area; S4: Extract the time history data of the constraint forces and / or the time history data of the constraint moments at the connection constraint locations; S5: Perform an equivalent transformation on the constraint force and / or constraint moment to generate an equivalent load that characterizes the overall dynamic effect of the human body acting on the target local area; S6: Construct a local simulation model of the target local area, and apply the equivalent load to the local simulation model to obtain the mechanical response results of the target local area under impact. S7: Based on the mechanical response of the target local area under impact, determine whether it is within the preset error range; if yes, execute S8; if no, return to S3 and update the connection constraint position corresponding to the target local area. S8: Output a multi-scale biomechanical model of the whole body and local areas.
[0008] Preferably, the whole-body dynamics model is constructed using the finite element dynamics model.
[0009] Preferably, the construction of the whole-body dynamic model in S1 specifically includes: Human skeletal models were extracted using AnyBody biomechanical software; By combining reverse engineering with Geomagic Wrap's NURBS surface reconstruction technology, complex skeletal geometry features are optimized and mesh quality is improved. HyperMesh is used to refine the mesh of key bones, solving the problem of dividing irregularly shaped tetrahedral units; Human body point cloud data was acquired using Creaform Go!SCAN 3D white light scanner to construct a hip soft tissue model, and the nonlinear mechanical behavior of cartilage and soft tissue was simulated based on hyperelastic constitutive model. The material properties of the skeleton are defined using a homogeneous isotropic linear elastic model; The joint connector is used to simulate the dynamic constraints of the shoulder, knee, and lumbosacral joints; Finally, a whole-body finite element dynamic model of the human body was established.
[0010] Preferably, step S1 further includes applying preset initial velocity and gravity conditions after constructing the whole-body dynamics model of the human body to simulate the human body impact process, and verifying the effectiveness of the whole-body dynamics model of the human body by comparing it with the simulation results under verified working conditions.
[0011] Preferably, the equivalent load acting on the target local area in S5 is generated based on the dynamic d'Alembert principle and the static equivalence principle.
[0012] Preferably, S5 specifically comprises: Perform coordinate system transformation on the constraint forces and / or constraint moments; Based on the geometric features and connectivity of the target local region, the constraint forces and / or constraint moments are converted into nodal loads or distributed loads applied to the local simulation model. The assignment is completed according to the rule of "load value = amplitude × filled load", and the equivalent load used to characterize the overall dynamic effect of the human body on the target local area is obtained.
[0013] Preferably, the equivalent load is a time history load corresponding to the constraint force and constraint moment.
[0014] Preferably, the local simulation model is any one of the finite element model, discrete element model, or multibody dynamics model.
[0015] Preferably, the human body impact simulation includes transient impact simulation under the action of falling, collision, or external force.
[0016] Preferably, the target local area includes any one of the human hip, head, spine, knee joint, or ankle joint.
[0017] Compared with the prior art, the beneficial effects of the technical solution of the present invention are: 1) This invention achieves the effective transfer of the overall dynamic effects of the human body to the local simulation model by equivalently transforming the dynamic constraints in the whole-body dynamic model and applying them to the local simulation model. Without directly constructing a high-precision whole-body simulation model, it can still reflect the influence of the overall inertial effect of the human body on the target local area.
[0018] 2) This invention constructs a local simulation model only for the target local area, which helps to reduce the size and computational complexity of the simulation model, improve the computational efficiency of human body impact simulation, and is suitable for engineering application scenarios.
[0019] 3) The method of the present invention has good versatility and scalability, and can be applied to the analysis of human biomechanical response under different impact conditions and in different local areas of the human body, providing an effective simulation means for human injury assessment, protective equipment design and human safety performance research. Attached Figure Description
[0020] Figure 1 This is a schematic diagram of a method for constructing a whole-body-local multi-scale biomechanical model coupling for human body impact simulation as described in Example 1.
[0021] Figure 2 A logical diagram for constructing a whole-body dynamics model of the human body.
[0022] Figure 3 A logical diagram for extracting complete time history data of the target hip region.
[0023] Figure 4 A logical diagram for constructing a local hip simulation model.
[0024] Figure 5 A schematic diagram of the logic for extracting equivalent loads for the target region.
[0025] Figure 6 This is a logical diagram illustrating the comparison of the stress time history trends between the local model and the whole-body model at the key section of the femoral neck. Detailed Implementation
[0026] To better understand the above-mentioned objectives, features, and advantages of the present invention, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments. It should be noted that, unless otherwise specified, the embodiments and features described in these embodiments can be combined with each other.
[0027] Many specific details are set forth in the following description in order to provide a full understanding of the invention. However, the invention may also be practiced in other ways different from those described herein, and therefore the scope of protection of the invention is not limited to the specific embodiments disclosed below. Example 1
[0028] like Figure 1 As shown in the figure, this embodiment discloses a method for constructing a whole-body-local multi-scale biomechanical model coupling for human impact simulation, the method comprising: S1: Construct a full-body dynamic model of the human body; S2: Set the impact condition parameters of the whole human body dynamics model, simulate the motion process of the human body under impact, and obtain the overall kinematic and dynamic response of the human body; S3: In the whole-body dynamics model of the human body, determine the connection constraint positions corresponding to the target local area; S4: Extract the time history data of the constraint forces and / or the time history data of the constraint moments at the connection constraint locations; S5: Perform an equivalent transformation on the constraint force and / or constraint moment to generate an equivalent load that characterizes the overall dynamic effect of the human body acting on the target local area; S6: Construct a local simulation model of the target local area, and apply the equivalent load to the local simulation model to obtain the mechanical response results of the target local area under impact. S7: Based on the mechanical response of the target local area under impact, determine whether it is within the preset error range; if yes, execute S8; if no, return to S3 and update the connection constraint position corresponding to the target local area. S8: Output a multi-scale biomechanical model of the whole body and local areas.
[0029] In this embodiment, the whole-body dynamics model is constructed using the finite element dynamics model.
[0030] In this embodiment, the construction of the whole-body dynamic model in S1 specifically refers to: Human skeletal models were extracted using AnyBody biomechanical software; By combining reverse engineering with Geomagic Wrap's NURBS surface reconstruction technology, complex skeletal geometry features are optimized and mesh quality is improved. HyperMesh is used to refine the mesh of key bones, solving the problem of dividing irregularly shaped tetrahedral units; Human body point cloud data was acquired using Creaform Go!SCAN 3D white light scanner to construct a hip soft tissue model, and the nonlinear mechanical behavior of cartilage and soft tissue was simulated based on hyperelastic constitutive model. The material properties of the skeleton are defined using a homogeneous isotropic linear elastic model; The joint connector is used to simulate the dynamic constraints of the shoulder, knee, and lumbosacral joints; Finally, a whole-body finite element dynamic model of the human body was established.
[0031] In this embodiment, S1 further includes applying preset initial velocity and gravity conditions after constructing the whole-body dynamics model of the human body to simulate the human body impact process, and verifying the effectiveness of the whole-body dynamics model of the human body by comparing it with the simulation results under verified working conditions.
[0032] In this embodiment, the equivalent load acting on the target local area in S5 is generated based on the dynamic d'Alembert principle and the static equivalence principle.
[0033] In this embodiment, S5 specifically refers to: Perform coordinate system transformation on the constraint forces and / or constraint moments; Based on the geometric features and connectivity of the target local region, the constraint forces and / or constraint moments are converted into nodal loads or distributed loads applied to the local simulation model. The assignment is completed according to the rule of "load value = amplitude × filled load", and the equivalent load used to characterize the overall dynamic effect of the human body on the target local area is obtained.
[0034] In this embodiment, the equivalent load is the time history load corresponding to the constraint force and constraint moment.
[0035] It should be noted that the equivalent load described in this embodiment is a time history load corresponding to the constraint force and / or constraint moment, used to reflect the influence of the overall inertial effect of the human body on the target local area as time changes.
[0036] The equivalent load is used to characterize the effect of the inertial effect generated by the overall mass distribution and motion state of the human body on the target local area.
[0037] In this embodiment, the local simulation model is any one of the finite element model, discrete element model, or multibody dynamics model.
[0038] In this embodiment, the human body impact simulation includes transient impact simulation under the action of human body falling, collision impact, or external force.
[0039] In this embodiment, the target local area includes any one of the human hip, head, spine, knee joint, or ankle joint.
[0040] In summary, this embodiment discloses a method for constructing a coupled whole-body-local multi-scale biomechanical model for human impact simulation, comprising: constructing a whole-body dynamic model of the human body and obtaining the kinematic and dynamic responses of the human body under a preset impact condition; determining the connection constraint positions corresponding to the target local region in the whole-body dynamic model and extracting the constraint forces and / or constraint moments at the connection constraint positions; performing an equivalent transformation on the constraint forces and / or constraint moments to obtain an equivalent load characterizing the overall dynamic effects of the human body acting on the target local region; constructing a local simulation model of the target local region and applying the equivalent load to the local simulation model to achieve simulation analysis of the human body impact response. This invention, through the transformation from whole-body dynamic constraint forces to a local simulation model, reduces the computational complexity of human impact simulation while preserving the overall inertial effects of the human body, improving simulation efficiency and computational stability, and is applicable to application scenarios such as fall impact, protective equipment evaluation, and human safety analysis.
[0041] As a specific embodiment, the method described in this embodiment will be explained below with reference to specific examples: This embodiment uses the hip impact during a side fall in an elderly person as an example to illustrate a method for constructing a whole-body-local multi-scale biomechanical model coupling for human impact simulation.
[0042] 1. Construct and validate a whole-body dynamic model of the human body (S1) like Figure 2 As shown, a whole-body dynamic model of the human body is constructed based on the body size parameters of the target population to describe the overall mass distribution, joint connection relationships, and kinematic constraints of the human body.
[0043] In this embodiment, the whole-body dynamic model of the human body is constructed in the form of a finite element dynamic model. The model includes the human skeletal system and the soft tissue structure covering its surface. The bones and soft tissues related to the target hip are set as deformable bodies, while the remaining distal bones are set as rigid bodies to reduce computational complexity.
[0044] Connectors are installed at key joint locations such as the shoulder, lumbosacral, and knee joints to simulate the mechanical behavior of human joints in multiple degrees of freedom, and to assign corresponding stiffness and damping parameters to each degree of freedom.
[0045] By applying preset initial velocity and gravity conditions, the impact process of a human body falling sideways is simulated, and the effectiveness of the human body dynamics model is verified by comparing it with the simulation results under verified working conditions.
[0046] 2. Extract the time history of dynamic constraint forces (S2-S4) After completing the whole-body dynamics simulation, the dynamic output data at the connection constraint positions corresponding to the target hip region is extracted in the post-processing stage.
[0047] like Figure 3 As shown, in this embodiment, the lumbosacral joint and knee joint connectors adjacent to the target hip region are selected, and the complete time history data of their triaxial constraint reaction force and triaxial constraint torque as a function of time in the global coordinate system are extracted.
[0048] It should be noted that, in Figure 3 The continuous constraint positions correspond to the boundaries of the target local region and are used to extract dynamic boundary loads.
[0049] The time history data is used to characterize the dynamic effects of the overall inertial effect of the human body transmitted through the joints to the target hip region.
[0050] 3. Construct a local high-precision hip simulation model and implement dynamic constraint force conversion and equivalent load application (S5-S7).
[0051] like Figure 4 As shown, the target hip bone, femur and related soft tissues are isolated from the whole-body dynamic model of the human body to construct a local hip simulation model.
[0052] The bone and soft tissue meshes in the local model are locally refined to improve the calculation accuracy of stress and contact force during impact, and only the anatomical boundaries connected to the lumbosacral joint and knee joint are retained.
[0053] A reference point is set at the boundary of the anatomical structure, and the reference point is kept in line with the movement of the corresponding bone surface by coupling constraints.
[0054] like Figure 5 As shown, the extracted constraint reaction force and constraint moment time history data are transformed into coordinate systems and equivalently converted into concentrated force loads and concentrated moment loads that vary with time.
[0055] It should be noted that, Figure 5 The time history characteristics of the intermediate load are used to characterize the dynamic effect of the overall inertial effect of the human body on local areas.
[0056] The equivalent load is applied to a reference point set at the boundary of the anatomical structure and is transmitted to the corresponding anatomical structure boundary of the local simulation model through coupling constraints, so as to equivalently characterize the dynamic effect of the rest of the human body on the target hip region.
[0057] 4. Method Validation and Effect Analysis Run a local simulation model with equivalent load and compare the simulation results with those of a whole-body dynamics model of the human body under the same impact conditions.
[0058] like Figure 6 As shown, within the allowable error range in engineering, the stress time history variation trends of the local model and the whole-body model at the key section of the femoral neck are consistent, verifying the effectiveness of the proposed method in dynamic response transmission.
[0059] It should be noted that, Figure 6 The critical mechanical response can be the resultant force, contact force, or reaction force at the critical section of the target region; the curve in the figure is a schematic curve.
[0060] Meanwhile, compared to directly using a whole-body finite element model for impact simulation, the local model constructed using the method of this invention significantly reduces the computational scale and computation time, thereby improving simulation efficiency.
[0061] In the several embodiments provided in this application, it should be understood that the disclosed devices and methods can be implemented in other ways. The device embodiments described above are merely illustrative. For example, the division of units is only a logical functional division, and in actual implementation, there may be other division methods, such as: multiple units or components can be combined, or integrated into another system, or some features can be ignored or not executed. In addition, the coupling, direct coupling, or communication connection between the various components shown or discussed can be through some interfaces, and the indirect coupling or communication connection between devices or units can be electrical, mechanical, or other forms.
[0062] The units described above as separate components may or may not be physically separate. The components shown as units may or may not be physical units. They may be located in one place or distributed across multiple network units. Some or all of the units may be selected to achieve the purpose of this embodiment according to actual needs.
[0063] In addition, in the various embodiments of the present invention, each functional unit can be integrated into one processing unit, or each unit can be a separate unit, or two or more units can be integrated into one unit; the integrated unit can be implemented in hardware or in the form of hardware plus software functional units.
[0064] Those skilled in the art will understand that all or part of the steps of the above method embodiments can be implemented by hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium. When the program is executed, it performs the steps of the above method embodiments. The aforementioned storage medium includes various media capable of storing program code, such as mobile storage devices, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.
[0065] Alternatively, if the integrated units of this invention are implemented as software functional modules and sold or used as independent products, they can also be stored in a computer-readable storage medium. Based on this understanding, the technical solutions of the embodiments of this invention, or the parts that contribute to the prior art, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the methods described in the various embodiments of this invention. The aforementioned storage medium includes various media capable of storing program code, such as mobile storage devices, ROM, RAM, magnetic disks, or optical disks.
Claims
1. A method for coupling and constructing a whole-body-local multi-scale biomechanical model for human impact simulation, characterized in that, The method includes: S1: Construct a full-body dynamic model of the human body; S2: Set the impact condition parameters of the whole human body dynamics model, simulate the motion process of the human body under impact, and obtain the overall kinematic and dynamic response of the human body; S3: In the whole-body dynamics model of the human body, determine the connection constraint positions corresponding to the target local area; S4: Extract the time history data of the constraint forces and / or the time history data of the constraint moments at the connection constraint locations; S5: Perform an equivalent transformation on the constraint force and / or constraint moment to generate an equivalent load that characterizes the overall dynamic effect of the human body acting on the target local area; S6: Construct a local simulation model of the target local area, and apply the equivalent load to the local simulation model to obtain the mechanical response results of the target local area under impact. S7: Based on the mechanical response of the target local area under impact, determine whether it is within the preset error range; if yes, execute S8; if no, return to S3 and update the connection constraint position corresponding to the target local area. S8: Output a multi-scale biomechanical model of the whole body and local areas.
2. The method for constructing a whole-body-local multi-scale biomechanical model coupling for human impact simulation according to claim 1, characterized in that, The whole-body dynamics model is constructed using the finite element dynamics model.
3. The method for constructing a whole-body-local multi-scale biomechanical model coupling for human impact simulation according to claim 2, characterized in that, The construction of the whole-body dynamic model in S1 is specifically as follows: Human skeletal models were extracted using AnyBody biomechanical software; By combining reverse engineering with Geomagic Wrap's NURBS surface reconstruction technology, complex skeletal geometry features are optimized and mesh quality is improved. HyperMesh is used to refine the mesh of key bones, solving the problem of dividing irregularly shaped tetrahedral units; Human body point cloud data was acquired using Creaform Go!SCAN 3D white light scanner to construct a hip soft tissue model, and the nonlinear mechanical behavior of cartilage and soft tissue was simulated based on hyperelastic constitutive model. The material properties of the skeleton are defined using a homogeneous isotropic linear elastic model; The joint connector is used to simulate the dynamic constraints of the shoulder, knee, and lumbosacral joints; Finally, a whole-body finite element dynamic model of the human body was established.
4. The method for constructing a whole-body-local multi-scale biomechanical model coupling for human impact simulation according to claim 3, characterized in that, S1 further includes applying preset initial velocity and gravity conditions after constructing the whole-body dynamics model of the human body to simulate the human body impact process, and verifying the effectiveness of the whole-body dynamics model of the human body by comparing it with the simulation results under verified working conditions.
5. A method for constructing a whole-body-local multi-scale biomechanical model coupling for human impact simulation according to any one of claims 1 to 4, characterized in that, The equivalent load acting on the local area of the target in S5 is generated based on the dynamic d'Alembert principle and the static equivalence principle.
6. The method for constructing a whole-body-local multi-scale biomechanical model coupling for human impact simulation according to claim 5, characterized in that, Specifically, S5 is: Perform coordinate system transformation on the constraint forces and / or constraint moments; Based on the geometric features and connectivity of the target local region, the constraint forces and / or constraint moments are converted into nodal loads or distributed loads applied to the local simulation model. The assignment is completed according to the rule of "load value = amplitude × filled load", and the equivalent load used to characterize the overall dynamic effect of the human body on the target local area is obtained.
7. The method for constructing a whole-body-local multi-scale biomechanical model coupling for human impact simulation according to claim 6, characterized in that, The equivalent load is the time history load corresponding to the constraint force and constraint moment.
8. A method for constructing a whole-body-local multi-scale biomechanical model coupling for human impact simulation according to claim 1 or 7, characterized in that, The local simulation model can be any one of the following: finite element model, discrete element model, or multibody dynamics model.
9. The method for constructing a whole-body-local multi-scale biomechanical model coupling for human impact simulation according to claim 8, characterized in that, The human body impact simulation includes transient impact simulations of human body falling impact, collision impact, or external force action.
10. A method for constructing a whole-body-local multi-scale biomechanical model for human impact simulation according to claim 1 or 9, characterized in that, The target local area includes any one of the human hip, head, spine, knee, or ankle joints.