CN105550461B - A kind of union analogue system based on broken ends of fractured bone fine motion and blood supply - Google Patents

A kind of union analogue system based on broken ends of fractured bone fine motion and blood supply Download PDF

Info

Publication number
CN105550461B
CN105550461B CN201511027837.6A CN201511027837A CN105550461B CN 105550461 B CN105550461 B CN 105550461B CN 201511027837 A CN201511027837 A CN 201511027837A CN 105550461 B CN105550461 B CN 105550461B
Authority
CN
China
Prior art keywords
poroma
bone
union
tissue
module
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201511027837.6A
Other languages
Chinese (zh)
Other versions
CN105550461A (en
Inventor
王沫楠
杨宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin University of Science and Technology
Original Assignee
Harbin University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN201511027837.6A priority Critical patent/CN105550461B/en
Publication of CN105550461A publication Critical patent/CN105550461A/en
Application granted granted Critical
Publication of CN105550461B publication Critical patent/CN105550461B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/30Polynomial surface description

Abstract

A kind of union analogue system based on broken ends of fractured bone fine motion and blood supply, is related to biomedical engineering field.The present invention is used for predicting the complex process of union, seeks the mechanism that there is no method to determine in union.The system comprises union simulation subsystem, input module and display modules;The wherein described union simulation subsystem (1) includes:Geometric Modeling module, mesh generation module, context initialization module, biomethanics computing module and biological environment module.Biological environment module is used in fracture healing process, builds the biochemical environment at poroma, completes the tissue differentiation in poroma.The present invention is introduced into blood supply as dynamic variable in union analogue system, can more accurately simulate fracture healing process.Carry out experimental simulation can be infinitely repeated several times in the present invention, to provide beneficial to help finding the influence of best broken ends of fractured bone fine motion value and blood to union.

Description

A kind of union analogue system based on broken ends of fractured bone fine motion and blood supply
Technical field
The present invention relates to biomedical engineering field, more particularly to a kind of union based on broken ends of fractured bone fine motion and blood supply is imitative True system.
Background technology
Musculoskeletal system plays an important role in human body, and major function includes:Motor function, defencive function, support Function, hematopoiesis function, storage function.In the common diseases statistics that China endangers people's health, wound is by the early days of foundation The 9th, in recent years with industry and communication development, it is the 4th to have gone up.Fracture is a kind of common wound, fracture It is frequently-occurring so as to fracture mechanism and promote healing research it is particularly urgent, once fracture occur, repaiied with other tissue damages Unlike multiple, fracture does not lean on fibrous connective tissue to connect, the complete regeneration of bone tissue.Bone tissue regeneration is a pole Its complicated biology repair process, it is affected by various factors, these influence factors be broadly divided into biological factor and Biomechanical factor.The prerequisite of fracture healing is it with good biologically ability, it is therefore desirable in fractures Position forms a good biology and mechanical environment.
During union, the importance of mechanics influence repair of fractured bones has been taken seriously.The biomethanics of union Factor is based on wolff laws, i.e. increasing or decreasing to adapt to be applied to the stress on bone tissue by bone tissue, that is, It says in a certain range, can there is bon e formation when need to increase carrying, has bone information when not needing or reduce carrying.Meanwhile in stress Under the action of fracture end can fine motion, Goodship etc., Kenwright etc. respectively by zoopery and clinical research confirmation, Under external fixing rack is fixed, controllable fine motion and strong stationary phase ratio promote union.In addition to Biomechanical factor influences Outside union, biological factor also plays a crucial role union.The human hairs such as Carano, Filvaroff The growth factor in union part and surrounding tissue is showed to Chondrogenesis, the generation of bone and the generation tool of blood vessel It plays an important role.In addition, the generation of blood vessel and revascularization are also the key factor in fracture healing process.
Although the research about union is carrying out always, there are many shortcoming and deficiency:
1. fracture healing process is extremely complex and cannot directly observe;
2. the size of the best fine motion state such as fine motion of union can be promoted not yet clear;
3. blood is not yet clear to the Influencing Mechanism of union;
4. clinical test is mostly experiment made on the living, so experiment number and experiment species can be restricted.
Invention content
The union analogue system based on broken ends of fractured bone fine motion and blood supply that the object of the present invention is to provide a kind of, for predicting bone The complex process for rolling over healing seeks the mechanism that there is no method to determine in union.
The purpose of the invention is achieved by the following technical solution:It is a kind of to be emulated based on broken ends of fractured bone fine motion and the union of blood supply System, it is characterized in that union simulation subsystem software, input equipment and display equipment;
Wherein the union simulation subsystem software includes:
Geometric Modeling module:For the two-dimentional tomoscan image data according to importing, carried out after image preprocessing The foundation of the 3 d surface model of fracture poroma;
Mesh generation module:Surface grids for 3 d surface model divide;And for veil lattice model to be generated body net Lattice model realizes the discretization of continuous 3-D geometric model, guarantees to carry out subsequent FEM calculation;
Context initialization module:For carrying out parameter setting, including biology to established poroma three-dimensional finite element model Mechanical environment Initialize installation and biological environment Initialize installation.
Biomethanics computing module:Distance for calculating the broken ends of fractured bone fine motion in fracture healing process, and pass through finite element Analysis calculates the strain stimulation of unit at poroma.Additionally include different healing phases, the meter of poroma unit material attribute It calculates.
Biological environment module:For in fracture healing process, building the biochemical environment at poroma, complete in poroma Tissue differentiation.
Wherein the input equipment includes:
Image data input unit, the two-dimentional tomoscan image data for reading patient;
User operation unit is transmitted and is handled for the information between user and union simulation system software, information Personal essential information including patient and poroma model view;
Wherein, the display equipment realizes that the vision between user and system is handed over for observing union situation Mutually.
The beneficial effects of the invention are as follows:
1. union analogue system proposed by the present invention is to develop software based on windows development languages platform, lead to The dynamic analog that self-programming realizes fracture healing process is crossed, the form based on dialog box is easily operated, and the training period is short.
2. it is introduced into blood supply as dynamic variable in union analogue system, it being capable of more accurate analog bone Roll over agglutination.
3. carry out experimental simulation can be infinitely repeated several times, to finding best broken ends of fractured bone fine motion value and blood to fracture The influence of healing is provided beneficial to help.
4. by building union analogue system, optimal operation plan can be formulated doctor, guidance is provided, in turn It improves success rate of operation, improve the case where fracture is cured quality, reduces nonunion and delayed union.
5. by building union analogue system, multiplicating experimental study can be carried out to the simulation model of foundation, True Bioexperiment is reduced, the time is saved, efficiency is improved, saves money, avoid humanitarian dispute.
To sum up, the shortcomings that emulation platform of the invention overcomes the prior art and deficiency.
Description of the drawings
Fig. 1 is that (wherein dotted line frame 1 represents input equipment to flow diagram of the invention, and dotted line frame 2 represents union emulation Software, dotted line frame 3, which represents, shows equipment), Fig. 2 is that (dotted line frame 1 represents Geometric Modeling mould to union simulation software flow chart Block, dotted line frame 2 represent mesh generation module, and dotted line frame 3 represents context initialization module, and dotted line frame 4 represents biomethanics calculating Module, dotted line frame 5 represent biological environment module), Fig. 3 is union simulation software Geometric Modeling module, and Fig. 4 is that fracture is cured The mesh generation module of simulation software is closed, Fig. 5 is the biomethanics computing module of union simulation software, and Fig. 6 is that fracture is cured Close the biological environment module of simulation software.
Specific implementation mode
Specific implementation mode one:As shown in Figure 1, a kind of fracture based on broken ends of fractured bone fine motion and blood supply described in present embodiment Healing analogue system include:Union simulation system software 1, input equipment 2, display equipment 3.
Wherein the union simulation system software 1 includes:At the beginning of Geometric Modeling module, mesh generation module, environment Beginningization module, biomethanics computing module, biological environment module.
The Geometric Modeling module is used for the two-dimentional tomoscan image data (such as CT images) according to importing, by image The foundation of the 3 d surface model of fracture poroma is carried out after pretreatment;
Surface grids of the mesh generation module for 3 d surface model divide;And for veil lattice model to be generated body Grid model realizes the discretization of continuous 3-D geometric model, guarantees to carry out subsequent FEM calculation;
The context initialization module to established poroma three-dimensional finite element model for being disposed setting, including life Material resources context initialization is arranged and biological environment Initialize installation;
The biomethanics computing module is used to calculate the distance of the broken ends of fractured bone fine motion in fracture healing process, and by limited Meta analysis calculates the strain stimulation of unit at poroma.Unit material attribute also in different time step-length;
The biological environment module is used in fracture healing process, builds the biochemical environment at poroma, completes poroma Interior tissue differentiation.
Wherein, the input equipment 2 includes:Image data input unit and user operation unit.
Described image data input cell is for reading the formats such as DICOM, JPG, BMP of the image datas such as patient CT text Part is the basis of Geometric Modeling module and context initialization module in union simulation system software.
The user operation unit is transmitted for the information between user and union simulation system software, such as software circle The input of the contents such as patient basis, case and the input of healing operation plan, initial broken ends of fractured bone fine motion and blood supply state in face Input etc..
The display equipment is used for observing union situation, realizes the visual interactive between user and system.
Specific implementation mode two, as shown in Fig. 2, in present embodiment, the union simulation system software 1 is realized The detailed process of its function is:
1) Geometric Modeling module
As shown in figure 3, Geometric Modeling module includes threshold value screening, Interactive Segmentation, three-dimensional reconstruction:Pass through threshold value Screening, Interactive Segmentation and three-dimensional reconstruction and etc. realize the 3 d geometric modeling of poroma;It will by the process that threshold value is screened The tissues such as bone, muscle, cartilage, fat, the skin of patient carry out preliminary segmentation respectively;It is wiped manually by interactive segmentation process Single picture is removed and repaired to improve the quality and speed of image segmentation;It can will be in two-dimentional tomograph by the process of three-dimensional reconstruction The Area generation 3 d surface model split as in;3 can be seen that from figure, and the process of Geometric Modeling needs the aobvious of computer Show that equipment, input equipment (mouse, keyboard etc.) coordinate user to complete.
2) mesh generation module
It is divided with volume mesh stroke in two steps as shown in figure 4, mesh generation includes surface grids;Surface grids partition process is 3 d surface model is subjected to some optimizations, including:Surface model optimization, smoothing processing, repairing loophole etc.;Surface model it is excellent Change and realized by reducing the tri patch of surface model, which only need to be merged into a new summit by two connected vertex On, and continue original topological relation;During smoothing processing, denoising is carried out to three-dimensional veil lattice model;Repairing leakage system Process by the way that the cavity in model is extracted into space polygon, the method that trigonometric ratio is then carried out to empty polygon is real It is existing;The process that volume mesh divides be veil lattice model stretched, rotate and etc. realize;
3) context initialization module
Poroma is made of three kinds of organization types:Connective tissue, cartilaginous tissue and bone tissue.At union initial stage, poroma It is mainly made of connective tissue, so poroma material properties are the material properties of connective tissue.Initial period, at poroma Blood supply is made of three parts:The blood supply of cortex bone is 100%;Due to being organized at fracture site by major injury, So the blood supply at fracture site is 0;Other than the blood supply of fracture itself, the soft tissue of surrounding of fracturing also can be to Fracture carries out blood supply, sets it to 30%.In addition to this, further include setting and the cortex bone material of plus load The assignment of attribute.
4) biomethanics computing module
Main includes the calculating of fracture site fine motion distance, the calculating of poroma element stress and poroma unit material attribute It calculates.Fig. 5 provides the flow diagram for calculating poroma element stress.
A. the calculating of fracture site fine motion:
The calculating of fracture site fine motion includes the calculating of the calculating and broken ends of fractured bone fine motion of current poroma hardness.
To simplify poroma model, it is assumed that poroma has linear, isotropic property, so there is following relationship:
FN=kcalμcal
In formula, FNFor the load that poroma is subject to, kcalFor current poroma hardness, μcalFor poroma displacement distance;
Wherein FNIt can be acquired by following formula:
FN=σ A
In formula, FNFor the load that poroma is subject to, σ is poroma stress, and A is forced area;
Poroma Stress calculation formula is:
σ=E ε
In formula, σ is poroma stress, and E is poroma elasticity modulus, and ε strains for poroma;
WhereinFor poroma relative movement distance, L is fracture spacing;
Pass through the calculating of above-mentioned formula, you can acquire poroma hardness kcal
There is following relationship between plus load, external fixator and poroma:
F (μ)=Ffix(μ)+kcalμ
In formula, F is plus load, FfixThe load born by fixator, kcalFor current poroma hardness, μ is that the broken ends of fractured bone is micro- Dynamic distance;
Broken ends of fractured bone fine motion distance μ can be acquired by the inverse function μ (F) of F (μ);
B. the calculating of poroma element stress
The unit strain that poroma is indicated with bulk strain and distortion strain, is shown below:
In formula, ε0For bulk strain, γ0It is strained for distortion;
Utilize unit principal strain ε1, ε2, ε3Bulk strain ε is acquired respectively0γ is strained with distortion0, relational expression is as follows:
Wherein, ε0The change of volume is represented for bulk strain;γ0The change of shape is represented for distortion strain;
C. the calculating of poroma unit material attribute
With the progress of union, the concentration of connective tissue, cartilaginous tissue and bone tissue in poroma unit changes Become, to make the material properties of poroma unit change, calculation formula is as follows:
In formula, EeleFor unitary elasticity modulus;EconFor connective tissue elasticity modulus, cconFor connective tissue concentration;EcartFor Cartilaginous tissue elasticity modulus, ccartFor cartilaginous tissue concentration;EboneFor bone tissue elasticity modulus, cboneFor bone tissue concentration;
θeleconcconcartccartbonecbone
In formula, θeleFor unit Poisson's ratio;θconFor connective tissue Poisson's ratio, cconFor connective tissue concentration;θcartFor cartilage Organize Poisson's ratio, ccartFor cartilaginous tissue concentration;θboneFor bone tissue Poisson's ratio, cboneFor bone tissue concentration;
Since poroma is made of the mixture of connective tissue, cartilaginous tissue and bone tissue, so there is following formula:
ccon+ccart+cbone=1
In formula, cconFor connective tissue concentration, ccartFor cartilaginous tissue concentration, cboneFor bone tissue concentration;
Simultaneous above-mentioned formula can acquire poroma unit material attribute.
5) biological environment module
As shown in fig. 6, realizing the tissue differentiation in poroma by fuzzy control:Wherein fuzzy control is a pass In a kind of control process of seven input variables and three output variables.Seven input variables are respectively:Current blood supply in unit, Cartilage concentration and bone concentration, close on the blood supply of unit and the bulk strain of bone concentration and active cell and distortion strains;Three Output variable is respectively:The variation of active cell blood supply, the variation of cartilage concentration and the variation of bone concentration.The control process is specific For:It learns according to clinical tissue and obtains the membership function of input variable and defeated to the tissue distribution patterns in different phase poroma Input variable can be blurred by the membership function for going out variable according to the membership function of input variable;It is obscured by certain The membership function of rule and output variable, then updated blood supply, cartilaginous tissue and bone can be obtained by anti fuzzy method Structural state;Thus tissue differentiation can be realized.
Embodiment:
In order to better illustrate this method and meaning, lower mask body gives one example to illustrate the operation of the present invention Journey.
Simulate healing of tibial fractures process:
1) image data is obtained
In order to establish the 3-D geometric model of poroma at patients with fractures, need to carry out CT, MRI to tissue at the patients with fractures Scanning, obtains DICOM data formats, and be acquired to data.
2) input of fracture poroma 3-D geometric model and patient information is established
The CT of DICOM format and MRI image are imported into union simulation system software by USB interface of computer In, then by the information inputs such as the name of patient, gender, age, case, indagation result to operation interface.By to image into Row threshold value screening operation separates the destination organization in image such as bone, muscle, fat and skin histology.Pass through meter Calculation machine mouse interacts formula segmentation to image to be come after wiping and repair manually, and the three-dimensional table of poroma at patients with fractures is finally obtained Surface model.
3) mesh generation of fracture poroma
It is excellent that the three-dimensional surface grid model that previous step is generated carries out surface model optimization, smoothing processing and repairing loophole etc. Change step.Then the veil lattice model of generation is subjected to volume mesh division by stretching, rotating etc..
4) the biomethanics modeling of poroma threedimensional model
The biomethanics of poroma threedimensional model models:By the interfaces MFC to established poroma three-dimensional finite element mould Type carries out Initialize installation.Including:The load of external load, the poroma initial structure, that is, setting of connective tissue material properties, skin The poroma setting of blood supply state everywhere and setting for initial broken ends of fractured bone fine motion distance under the setting of matter bone material attribute, initial conditions It sets.Above-mentioned setting can be configured by the interfaces MFC in the form of parameter.
After above-mentioned parameter is provided with, the simulation of union can be carried out.
5) poroma biomethanics calculates
The biomethanics of poroma threedimensional model calculates:It is single in broken ends of fractured bone fine motion distance, poroma under the different time stage The bulk strain of member and the calculating of the material properties of distortion strain and the poroma unit under the different time stage.
6) the biology modeling of poroma threedimensional model
The biology of poroma unit models:Current poroma unit bone concentration, current poroma unit cartilage concentration, currently The blood supply state of poroma unit closes on the blood supply state of poroma unit, closes on the bone concentration of poroma unit, is asked in cooperation previous step The bulk strain and distortion of the current poroma unit obtained strain, you can complete seven input variables that biology models.In basis The fuzzy rule write, the membership function of above-mentioned input variable is blurred, and next state can be obtained in comprehensive function The knots modification of bone concentration, cartilage concentration and blood supply state, to be that organization type concentration is updated.
Above-mentioned steps 5), step 6) carries out in union simulation system software.
Step 5), step 6) iterative cycles carry out in union simulation system software, during cycle, by aobvious Show that equipment can see union situation, sees whether disunion phenomenon occur.When broken ends of fractured bone fine motion distance is reduced to 0, more It closes and completes, record healing time.Healing process in this way can simulate different healing schemes, pass through Change broken ends of fractured bone fine motion distance and blood supply parameter determines best fracture to obtain the union situation under different conditions Healing scheme and the broken ends of fractured bone fine motion distance and blood supply state that beneficial effect is generated to union.

Claims (4)

1. a kind of union analogue system based on broken ends of fractured bone fine motion and blood supply, which is characterized in that the system comprises:
Union simulation subsystem (1), input module (2) and display module (3);
The wherein described union simulation subsystem (1) includes:
Geometric Modeling module:For the two-dimentional tomoscan image data according to importing, fracture after image preprocessing Locate the foundation of the 3 d surface model of poroma;
Mesh generation module:Surface grids for 3 d surface model divide, and for veil lattice model to be generated volume mesh mould Type realizes the discretization of continuous 3-D geometric model, guarantees to carry out subsequent FEM calculation, to obtain poroma three Tie up finite element model;
Context initialization module:For carrying out Initialize installation, including Biological Strength to established poroma three-dimensional finite element model Learn context initialization setting and biological environment Initialize installation;
Biomethanics computing module:Distance for calculating the broken ends of fractured bone fine motion in fracture healing process, and by finite element analysis, Calculate the strain stimulation of unit in poroma;Additionally include different healing phases, the calculating of poroma unit material attribute;
Biological environment module:For in fracture healing process, building the biochemical environment at poroma, the tissue in poroma is completed Atomization;
The union simulation subsystem (1) realizes that its function detailed process is:
1) Geometric Modeling module
The 3 d geometric modeling of poroma is realized by threshold value screening, Interactive Segmentation and three-dimensional reconstruction;It is screened by threshold value Process the bone of patient, muscle, cartilage, fat, skin histology are subjected to preliminary segmentation respectively;Pass through interactive segmentation process Single picture is wiped and repaired manually to improve the quality and speed of image segmentation;It will be in two-dimensional ct image by three-dimensional reconstruction In the Area generation 3 d surface model that splits;
2) mesh generation module
Mesh generation includes that surface grids are divided with volume mesh stroke in two steps;Surface grids partition process is used for three-dimensional table face mould Type optimizes, including:Surface model optimization, smoothing processing, repairing loophole;The optimization of surface model is by reducing surface model Tri patch realize that two connected vertex need to be only merged into a new summit by the process, and continue original open up Flutter relationship;During smoothing processing, denoising is carried out to three-dimensional veil lattice model;During repairing loophole, by by mould Cavity in type is extracted into space polygon, and the method for then carrying out trigonometric ratio to empty polygon is realized;Volume mesh divides Process be veil lattice model stretched, spin step is realized;
3) context initialization module
The poroma is made of three kinds of organization types:Connective tissue, cartilaginous tissue and bone tissue;At union initial stage, poroma It is mainly made of connective tissue, so poroma material properties are the material properties of connective tissue;
At union initial stage, the blood supply at poroma is made of three parts:Blood supply in cortex bone is 100%;By Tissue is by major injury at fracture site, so the blood supply at fracture site is 0;In addition to the blood of fracture itself Supply is outer, and the soft tissue for surrounding of fracturing also can carry out blood supply to fracture, set it to 30%;
Context initialization module is additionally operable to the setting of plus load and the assignment of cortex bone material properties;
4) biomethanics computing module:
The meter of the calculating and poroma unit material attribute of main calculating, the strain of poroma unit for including fracture site fine motion distance It calculates;
A. the calculating of fracture site fine motion:
The calculating of fracture site fine motion includes the calculating of the calculating and broken ends of fractured bone fine motion of current poroma hardness;
Assuming that poroma has linear, isotropic property;So there is following relationship:
FN=kcalμcal,
In formula, FNFor the load that poroma is subject to, kcalFor current poroma hardness, μcalFor poroma displacement distance;
Wherein FNIt can be acquired by following formula:
FN=σ A
In formula, FNFor the load that poroma is subject to, σ is poroma stress, and A is forced area;
Poroma Stress calculation formula is:
σ=E ε
In formula, σ is poroma stress, and E is poroma elasticity modulus, and ε strains for poroma;
Wherein△ L are poroma relative movement distance, and L is fracture spacing;
Pass through the calculating of above-mentioned formula, you can acquire poroma hardness kcal
There is following relationship between plus load, external fixator and poroma:
F (μ)=Ffix(μ)+kcalμ
In formula, F is plus load, FfixThe load born by fixator, kcalFor current poroma hardness, μ be broken ends of fractured bone fine motion away from From;
Broken ends of fractured bone fine motion distance μ can be acquired by the inverse function μ (F) of F (μ);
B. the calculating of poroma element stress
The element stress of poroma is indicated with bulk strain and distortion strain, is shown below:
In formula, ε0For bulk strain, γ0It is strained for distortion;
Utilize unit principal strain ε1, ε2, ε3Bulk strain ε is acquired respectively0γ is strained with distortion0, relational expression is as follows:
Wherein, ε0The change of volume is represented for bulk strain;γ0The change of shape is represented for distortion strain;
C. the calculating of poroma unit material attribute
With the progress of union, the concentration of connective tissue, cartilaginous tissue and bone tissue in poroma unit changes, from And the material properties of poroma unit is made to change, calculation formula is as follows:
In formula, EeleFor unitary elasticity modulus;EconFor connective tissue elasticity modulus, cconFor connective tissue concentration;EcartFor cartilage Tissue elasticity modulus, ccartFor cartilaginous tissue concentration;EboneFor bone tissue elasticity modulus, cboneFor bone tissue concentration;
θeleconcconcartccartbonecbone
In formula, θeleFor unit Poisson's ratio;θconFor connective tissue Poisson's ratio, cconFor connective tissue concentration;θcartFor cartilaginous tissue Poisson's ratio, ccartFor cartilaginous tissue concentration;θboneFor bone tissue Poisson's ratio, cboneFor bone tissue concentration;
Since poroma is made of the mixture of connective tissue, cartilaginous tissue and bone tissue, there is following formula:
ccon+ccart+cbone=1
In formula, cconFor connective tissue concentration, ccartFor cartilaginous tissue concentration, cboneFor bone tissue concentration;
Simultaneous above-mentioned formula can acquire poroma unit material attribute;
5) biological environment module
The tissue differentiation in poroma is realized by fuzzy control:Wherein fuzzy control is one about seven input variables With a kind of control process of three output variables;Seven input variables are respectively:Current blood supply, cartilage concentration and bone are dense in unit Degree, closes on the blood supply of unit and the bulk strain of bone concentration and active cell and distortion strains;Three output variables are respectively: The variation of active cell blood supply, the variation of cartilage concentration and the variation of bone concentration;The control process is specially:According to clinical tissue Obtain the degree of membership letter that the tissue distribution patterns in different phase poroma obtain the membership function and output variable of input variable Input variable, can be blurred by number according to the membership function of input variable;Pass through being subordinate to for fuzzy rule and output variable Function is spent, then updated blood supply, cartilaginous tissue and bone tissue state can be obtained by anti fuzzy method;
Thus tissue differentiation can be realized;
Wherein the input module (2) includes:
Image data input unit, the two-dimentional tomoscan image data for reading patient;
User operation unit, is transmitted for the information between user and union simulation system software and processing, information include The personal essential information and poroma model view of patient;
Wherein, the visual interactive between user and system is realized for observing union situation in the display module (3).
2. a kind of union analogue system based on broken ends of fractured bone fine motion and blood supply according to claim 1, it is characterised in that: The input module (2) realizes that the detailed process of its function is:
1) image data input unit
Two-dimentional tomoscan image is imported by computer input device in union simulation system software, is built for geometry Mould is prepared;
2) user operation unit
Operation relevant information is inputted by computer input device for operating unit, and is cured with fracture by operation input equipment The operation interface for closing simulation system software interacts.
3. a kind of union analogue system based on broken ends of fractured bone fine motion and blood supply according to claim 1, it is characterised in that: The display module (3) realizes that the detailed process of its function is:
The display equipment of computer is connected to computer, and for carrying out visual interaction.
4. a kind of union analogue system based on broken ends of fractured bone fine motion and blood supply according to claim 1, which is characterized in that Biomethanics computing module and biological environment functions of modules in union simulation software subsystem (1) are all by repeatedly Generation cycle is realized.
CN201511027837.6A 2015-12-30 2015-12-30 A kind of union analogue system based on broken ends of fractured bone fine motion and blood supply Expired - Fee Related CN105550461B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511027837.6A CN105550461B (en) 2015-12-30 2015-12-30 A kind of union analogue system based on broken ends of fractured bone fine motion and blood supply

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201511027837.6A CN105550461B (en) 2015-12-30 2015-12-30 A kind of union analogue system based on broken ends of fractured bone fine motion and blood supply

Publications (2)

Publication Number Publication Date
CN105550461A CN105550461A (en) 2016-05-04
CN105550461B true CN105550461B (en) 2018-10-02

Family

ID=55829650

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511027837.6A Expired - Fee Related CN105550461B (en) 2015-12-30 2015-12-30 A kind of union analogue system based on broken ends of fractured bone fine motion and blood supply

Country Status (1)

Country Link
CN (1) CN105550461B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106227993B (en) * 2016-07-14 2017-10-13 哈尔滨理工大学 A kind of union dynamic process simulation method based on Biological Mechanism
CN106557665A (en) * 2016-12-01 2017-04-05 哈尔滨理工大学 A kind of union emulation mode that algorithm is adjusted based on power biology
CN106777582B (en) * 2016-12-01 2018-01-19 哈尔滨理工大学 A kind of long bone fracture healing analogue system based on tissue differentiation
CN106777583B (en) * 2016-12-01 2018-01-19 哈尔滨理工大学 A kind of union emulation mode based on cellular activity
CN106777584B (en) * 2016-12-01 2017-12-08 哈尔滨理工大学 A kind of analogue system for simulating fracture healing process
CN109147048B (en) * 2018-07-23 2021-02-26 复旦大学 Three-dimensional mesh reconstruction method by utilizing single-sheet colorful image
CN111067532B (en) * 2019-12-16 2022-09-06 天津市威曼生物材料有限公司 Parameter acquisition device for fracture internal fixation system
CN111311740B (en) * 2020-03-23 2023-09-01 北京工业大学 Traction and expansion osteogenesis numerical simulation method based on tissue viscoelastic-plastic mechanical properties

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7642234B2 (en) * 2002-05-17 2010-01-05 Board Of Regents The University Of Texas Systems Beta-2-glycoprotein 1 is an inhibitor of angiogenesis
CN103530466A (en) * 2013-10-21 2014-01-22 哈尔滨理工大学 Method for optimally selecting thighbone prostheses based on material performance multi-objective optimization

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7642234B2 (en) * 2002-05-17 2010-01-05 Board Of Regents The University Of Texas Systems Beta-2-glycoprotein 1 is an inhibitor of angiogenesis
CN103530466A (en) * 2013-10-21 2014-01-22 哈尔滨理工大学 Method for optimally selecting thighbone prostheses based on material performance multi-objective optimization

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
实验性兔下颌骨骨折愈合过程的研究—生物力学仿真模型的建立;史俊;《万方数据知识服务平台》;20070807;第28-53页 *

Also Published As

Publication number Publication date
CN105550461A (en) 2016-05-04

Similar Documents

Publication Publication Date Title
CN105550461B (en) A kind of union analogue system based on broken ends of fractured bone fine motion and blood supply
CN105303605B (en) A kind of bone surgery analogue system based on force feedback
CN105405167A (en) Finite element modeling method based on complete human head
CN106777584B (en) A kind of analogue system for simulating fracture healing process
CN108511076B (en) A kind of union analogue system based on mechanical stimulation and bio combined stimulation
CN106777583B (en) A kind of union emulation mode based on cellular activity
CN108133755A (en) A kind of atlas pivot surgery simulation system and its analogy method based on three-dimensional visualization
CN106227993A (en) A kind of union of fracture dynamic process simulation method based on Biological Mechanism
CN106777582A (en) A kind of long bone fracture healing analogue system based on tissue differentiation
Martin-Guerrero et al. Machine learning for modeling the biomechanical behavior of human soft tissue
Dao et al. Image-based skeletal muscle coordination: case study on a subject specific facial mimic simulation
CN107134010A (en) A kind of pattern effect prediction method of the elastic soft tissue based on finite element
Levine et al. Mesh processing in medical-image analysis—a tutorial
Dornelles et al. New virtual tool for accurate evaluation of facial volume
Shi et al. A new volumetric geometric model for cutting procedures in surgical simulation
Cheng et al. Facial morphology prediction after complete denture restoration based on principal component analysis
Wang et al. Real-time FEA-based breast deformation simulation using artificial neural network
Wu A computational framework for modelling the biomechanics of human facial expressions
Qiu et al. Biomechanical analysis of reduction malarplasty with L-shaped osteotomy
Yu et al. Creating and simulating a realistic physiological tongue model for speech production
Bian et al. A survey of tongue modeling methods in speech visualization
Dong et al. A Nonlinear Viscoelastic Meshless Model for Soft Tissue Deformation
Kumar et al. Modeling and analysis of human liver by using finite element analysis
Li et al. High fidelity 3D anatomical visualization of the fibre bundles of the muscles of facial expression as in situ
Inoue et al. Development of surgical simulator with high‐quality visualization based on finite‐element method and deformable volume rendering

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20181002

Termination date: 20191230