CN102262699A - Soft tissue deformation simulation method based on coupling of mesh-free Galerkin and mass spring - Google Patents

Soft tissue deformation simulation method based on coupling of mesh-free Galerkin and mass spring Download PDF

Info

Publication number
CN102262699A
CN102262699A CN 201110213387 CN201110213387A CN102262699A CN 102262699 A CN102262699 A CN 102262699A CN 201110213387 CN201110213387 CN 201110213387 CN 201110213387 A CN201110213387 A CN 201110213387A CN 102262699 A CN102262699 A CN 102262699A
Authority
CN
China
Prior art keywords
particle
subelement
grid
unit
particle spring
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN 201110213387
Other languages
Chinese (zh)
Other versions
CN102262699B (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.)
North China University of Water Resources and Electric Power
Original Assignee
North China University of Water Resources and Electric Power
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 North China University of Water Resources and Electric Power filed Critical North China University of Water Resources and Electric Power
Priority to CN201110213387A priority Critical patent/CN102262699B/en
Publication of CN102262699A publication Critical patent/CN102262699A/en
Application granted granted Critical
Publication of CN102262699B publication Critical patent/CN102262699B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention relates to an object deformation real-time simulation graphic processing technique, particularly a soft tissue deformation simulation method based on coupling of mesh-free Galerkin and mass spring, which comprises the following steps: in the pretreatment process, establishing a linear viscoelasticity biomechanical model for soft tissues; in the deformation computation process, dynamically partitioning a mesh-free region and a mass spring region according to the load carried by the soft tissues, establishing a transitional unit of the connection region between the mesh-free region and mass spring region, and constructing a transitional unit approximation displacement function, thereby implementing self-adapting coupling of a mesh-free Galerkin method and a mass spring method; and in the after-treatment process, outputting the state of the mass or node of each time step in the deformation process onto a screen, carrying out illumination rendering, and finally displaying the real-time deformation process of the soft tissue organ under stressed conditions on the screen, thereby implementing visualization effect of dynamic deformation. By utilizing the advantage of high efficiency in the mass spring method and the advantages of high precision and no need of mesh reconstruction in the mesh-free Galerkin method, the invention overcomes the defect that the Galerkin method is not suitable for solving a large-scale problem, thereby effectively lowering the complexity of computation in the soft tissue deformation simulation and enhancing the operation efficiency.

Description

Soft tissue deformation emulating method based on no grid gal the Liao Dynasty's gold and the coupling of particle spring
Technical field
The present invention relates to a kind of object deformation real time modelling graph processing technique, particularly relate to a kind of soft tissue deformation emulating method based on no grid gal the Liao Dynasty's gold and the coupling of particle spring.
Background technology
Human body soft tissue deformation Study of calculation model can be traced back to the eighties in last century, early stage deformation model comes from Computer-aided Geometric Design (CAGD) field, be some non-physical models that use pure geometric techniques, this model is not considered the physics law of true deformation, after 1 year, the deformation model based on physical features that is proposed by Terzopoulos substitutes.At present, deformation computation model based on physical characteristics mainly is divided into particle spring model (Mass-Spring), finite element model (Finite-Element Model, be abbreviated as FEM) and boundary element model (Boundary Element Model is abbreviated as BEM) three major types.The particle spring model is separated into some spots with object, and using has the point of quality to come representation physique amount, influencing each other between expressing a little with the spring of no quality and damper.This modular concept is simple, calculated amount is little, and real-time is good, is widely used in the virtual operation emulation such as distortion, cutting, stitching of soft tissue by the scholar.But itself also has the shortcoming that himself can't overcome this model, owing to connect with spring between the supposition soft tissue, and the setting of the elastic coefficient does not have theoretical foundation, can only derive from laboratory technician's experience and debugging constantly, and when the soft tissue large deformation problem of dtmf distortion DTMF can appear, therefore, the finite element model that can overcome this shortcoming becomes the focus of research gradually.
Finite element method is the classical way of finding the solution Elasticity Problems, and its basic thought is to find the solution the assembly that discrete region is a group limited and the unit that is connected with each other by certain way with continuous.Come burst ground to represent to find the solution dependent variable undetermined in the territory with the approximate function of being supposed in each unit, utilize variational principle to set up the Algebraic Equation set of finding the solution dependent variable, thereby calculate deformation.Its advantage is to have higher precision, and adaptability is stronger, can be applicable to the emulation of various geometries and material behavior.Its shortcoming is: the pre-treatment for the 3 D complex structure is that 3D grid makes up difficulty, and is consuming time longer; The common displacement of its result of calculation is continuous, and stress, strain are discontinuous, need the aftertreatment smooth; " volume self-locking " phenomenon appears when structured material is incompressible; The mesh distortion problem can appear when large deformation.Because soft tissue organs has incompressibility, so when using finite element model, " volume self-locking " phenomenon can occur; In virtual operation, large deformation situation or cutting, sew application often take place, at this moment the topological structure of organ can change, and utilizes finite element method constantly to carry out grid reconstruction, not only takes time and effort, and occurs the mesh distortion problem easily; Because stress, the strain result of FEM (finite element) calculation are discontinuous, need the aftertreatment smooth, can influence the accuracy and the real-time of force feedback in the virtual operation.This shows that finite element method is not the preferred plan that soft tissue deformation is calculated, need to explore other computing method that are more suitable for.
Boundary element method is meant the elementary solution based on governing equation, and the boundary problem in zone is turned to equation on the boundary surface, and division unit on boundary surface is found the solution with Galerkin method or other numerical computation method more then.It can become two-dimensional problems to three-dimensional problem, and two-dimensional problems become one-dimensional problem, and the degree of freedom of finding the solution is descended, but because boundary element method also is based on the numerical method of grid, so also mesh distortion and distortion may occur.
Commonly used have the grid physical model need expend very big energy to construct grid model, and follow-up computation process mostly closely depends on this network structure; Mesh distortion might take place when the soft tissue large deformation; Cutting, during topology change such as stitching, need re-construct grid; Therefore there is grid method to be difficult to large deformation of accurate description soft tissue and topology change situation.For addressing these problems, produced the idea of in the numerical value processing procedure, breaking away from grid naturally, no gridding method is arisen at the historic moment.No gridding method is removed the network in the finite element method, replacing a series of node fully arranges, need not utilize predefined grid that field variable is carried out interpolation and approximate, and be based on a little approximate, adopt weight function to characterize nodal information, and the node on certain territory can influence the mechanical characteristic of any point on the research object.Like this, broken away from uncontinuity, guaranteed the precision found the solution to be particularly suitable for the topology change that cutting, sew application cause in the virtual operation the constraint of problem (as grid reconstruction etc.).In no grid method, Galerkin method only needs nodal information to describe with the border just can become one group of algebra system equation with formula formal transformation a little less than it, and adopt a kind of smooth mobile least squqre approximation function that can reach optimal approximation as the displacement approximate function with the value of each sampling spot, characterize nodal information with this, so computational accuracy is highly stable far above the strong formula method of no grid.
But, the golden shortcoming of no grid gal the Liao Dynasty is that calculated amount is bigger, the large-scale golden method of gal the Liao Dynasty of using can increase computing time greatly, therefore design the deformation emulating that a better model of performance is realized soft tissue, and the satisfied as much as possible requirement about real-time, robustness, accuracy, become the matter of utmost importance that the soft tissue deformation emulating faces.In fact, only the contact area at operating theater instruments and organ just can produce topology changes such as big distortion and cutting and stitching in virtual operation, and is less at other regional deformation, therefore there is no need the golden method in whole zone use no grid gal the Liao Dynasty.Therefore, we adopt the method for no grid gal the Liao Dynasty's gold and the coupling of particle spring, in the golden method of contact area employing no grid gal the Liao Dynasty of operating theater instruments and organ, in other zone employing particle spring method.
Summary of the invention
The present invention is directed to the prior art deficiency, a kind of soft tissue deformation emulating method based on no grid gal the Liao Dynasty's gold and the coupling of particle spring is proposed, to not have grid gal the Liao Dynasty gold and solve soft tissue deformation emulating problem in the virtual operation with the coupling of particle spring, purpose is to utilize the high and no grid Galerkin method precision height of particle SPRING METHOD OF STONE efficient, need not the advantage of reconstructed mesh, remedy the defective that Galerkin method is not suitable for finding the solution extensive problem, thereby effectively reduce the computation complexity in the soft tissue deformation emulating, improve operation efficiency, better solve the contradiction of accuracy and real-time in the soft tissue deformation emulating.
The technical solution adopted in the present invention:
A kind of soft tissue deformation emulating method based on no grid gal the Liao Dynasty's gold and the coupling of particle spring comprises the pre-treatment process, calculates deformation process, three steps of last handling process, in the pre-treatment process, for soft tissue is set up the linear viscoelasticity biomechanical model; In calculating deformation process, dynamically divide no net region and particle spring zone according to the load that soft tissue carried, and the join domain between no net region and particle spring zone is set up transition element, the approximate displacement function of structure transition element is realized the self-adaptation coupling of no grid gal golden method of the Liao Dynasty and particle spring method; In last handling process, according to deformation result of calculation, the particle of each time step of deformation process or the state of node are outputed on the screen, and carry out illumination and play up, finally, realize dynamic deformation effect of visualization showing the real-time deformation process of soft tissue organs under stressing conditions on the screen; The process of described calculating deformation specifically comprises following four steps:
1) function between design load and the distance, as the foundation of dividing no net region, and make no net region enough greatly with the assurance noncoherent boundary all in no net region, other zone is particle spring zone;
2) for no net region and particle spring zone, set up effective data structure respectively, the Classification Management data;
3) join domain between no net region and particle spring zone is set up transition element:
(1) the no net region that will divide and particle spring zone are as two entities, two entities contact part as transition boundary, on the basis of background grid, the background grid at segmentation transition boundary place and as subelement, making has a node or particle to exist at most in each subelement, with the unit that do not comprise node or particle as the gap unit;
(2) with the gap unit as object search, the left-right and front-back six direction is searched for subelement up and down respectively, through iteration, progressively will not have the gap unit that inside, net region meets switch condition and be converted into no grid cell, the gap unit that particle spring intra-zone is met switch condition is converted into the particle spring unit, no longer change until the gap unit, and with remaining gap unit as transition element;
4) in transition element, set up transition node, determine the approximate displacement function of transition node, realize the smooth transition between no net region and the particle spring zone; Two conditions are satisfied in the foundation of transition node:
(1) transition node is subjected to not having the acting force in net region and particle spring zone respectively, and it is made a concerted effort is zero;
(2) transition node is subjected to not having the acting force in net region and particle spring zone respectively, and its displacement equates;
According to above-mentioned condition, make up linear viscoelasticity dynamic motion equation respectively based on no grid gal the Liao Dynasty's gold and particle spring, find the solution the differential equation, draw the approximate displacement of two differences of same transition node, travel through all transition nodes, obtain the approximate displacement function that two kinds of model constructions go out, set up the funtcional relationship between the two, as the approximate displacement function of transition node.
Described soft tissue deformation emulating method based on no grid gal the Liao Dynasty's gold and the coupling of particle spring, the switch condition that inner gap unit, no net region is converted into no grid cell is:
1) top to bottom, left and right, front and rear three prescriptions are no grid cell to the subelement of middle existence direction more than two groups and two groups;
2) six subelements on the six direction of top to bottom, left and right, front and rear all are the gap unit;
3) for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situations that all exist, and every prescription is to getting a subelement, and three subelements on totally three directions are no grid cell;
4) for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situation that all exists, wherein the subelement that makes progress of a prescription is all got, in addition two prescriptions in every group of subelement of getting on one of them direction, four subelements that amount on the four direction are no grid cell;
5) for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situation that all exists, wherein the subelement that makes progress of two prescriptions is all got, another prescription is upwards got the subelement on one of them direction, and five subelements that amount on five directions are no grid cell.
Described soft tissue deformation emulating method based on no grid gal the Liao Dynasty's gold and the coupling of particle spring, the gap unit that particle spring intra-zone is met switch condition is converted into the particle spring unit, and its switch condition is:
1) top to bottom, left and right, front and rear three prescriptions are the particle spring unit to the subelement of middle existence direction more than two groups and two groups;
2) six subelements on the six direction of top to bottom, left and right, front and rear all are the gap unit;
3) for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situations that all exist, and every prescription is to getting a subelement, and three subelements on totally three directions are the particle spring unit;
4) for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situation that all exists, wherein the subelement that makes progress of a prescription is all got, in addition two prescriptions in every group of subelement of getting on one of them direction, four subelements that amount on the four direction are the particle spring unit;
5) for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situation that all exists, wherein the subelement that makes progress of two prescriptions is all got, another prescription is upwards got the subelement on one of them direction, and five subelements that amount on five directions are the particle spring unit.
Described soft tissue deformation emulating method based on no grid gal the Liao Dynasty's gold and the coupling of particle spring, remaining gap unit is meant as transition element:
1) do not meet no grid cell switch condition, do not meet the gap unit of particle spring unit switch condition yet;
2) six subelements on the six direction of top to bottom, left and right, front and rear comprise no grid cell and particle spring unit simultaneously, perhaps comprise no grid cell, particle spring unit and gap unit simultaneously.
Described soft tissue deformation emulating method based on no grid gal the Liao Dynasty's gold and the coupling of particle spring, in the pre-treatment process, suppose that soft tissue is even, isotropic material, have accurate incompressibility, linear viscoelasticity feature, and definite initialization data, definition starting condition, displacement boundary conditions and stress boundary condition; With the elasticity of spring simulation viscoelastic body, damper is simulated its viscosity then, and the distortion and the load of hypothesis spring are proportional, and the speed of the distortion of damper and distortion is proportional; Secondly the distortion of isotropic material can be separated into stereomutation and isopyknic shape distortion two parts, so the stress and strain tensor is resolved into its spherical tensor sum deviator, and then the constitutive equation of structure Kelvin model, simulation soft tissue stress relaxation and creep properties; Set up geometric equation, balance equation at last, and apply boundary condition and starting condition, for soft tissue is set up the linear viscoelasticity biomechanical model.
Described soft tissue deformation emulating method based on no grid gal the Liao Dynasty's gold and the coupling of particle spring, in calculating deformation process, at first generate initial background grid, divide no net region and particle spring zone according to external load, foundation from no grid to the data structure the particle spring data, to reduce the frequency that the data reusing phenomenon takes place in the adaptive process, secondly set up transition element at join domain, in transition element, set up the approximate displacement function of transition node and transition node, realize the smooth transition between no net region and the particle spring zone, apply boundary condition then, in conjunction with Kelvin normal linearity Viscous Elastic Mechanics model, adopt method of addition to find the solution displacement, ask stress according to displacement, strain.
Useful good effect of the present invention:
1, in large deformation and the topology change zone golden method of no grid gal the Liao Dynasty, the link information that does not need node, effectively avoid finite element and mesh distortion in the Element BEM and reconstruct commonly used at present, be fit to very much often take place in the operation emulation situation of large deformation and topology change; The golden method of no grid gal the Liao Dynasty has good continuity and form field function flexibly, has advantages of higher stability and computational accuracy; The ess-strain result of the golden method of no grid gal the Liao Dynasty has slickness preferably, can fine holding power feedback device.
2, adopt based on no grid gal the Liao Dynasty's gold and particle spring coupling process, adopt the golden method of high-precision no grid gal the Liao Dynasty in large deformation and topology change zone, adopt high efficiency particle spring method in other deformation than the zonule, not only can bring into play the plurality of advantages of the golden method of above-mentioned no grid gal the Liao Dynasty, and improve deformation emulating efficient.
3, use can effectively be simulated creep and the relaxation property that soft tissue is shown based on Kelvin normal linearity viscoelastic model in deformation process.This Elasticity model that generally uses at present based on viscoelastic model for coupling meets the biomechanics characteristic of human body soft tissue more.
Description of drawings
Fig. 1 is the soft tissue deformation emulating manner of execution process flow diagram that the present invention is based on no grid gal the Liao Dynasty's gold and the coupling of particle spring;
Fig. 2 dynamically divides the sectional view in no net region and particle spring zone for the present invention, among the figure, 1 is the operating theater instruments range of influence, the distant golden Model Calculation soft tissue deformation of no grid gal is used in this zone, 2 is transitional region, 3 is particle spring zone, and this zone uses the particle spring model to calculate soft tissue deformation;
Fig. 3 dynamically divides the stereographic map in no net region and particle spring zone for the present invention;
Fig. 4 is the sectional view of no grid gal the Liao Dynasty's gold of the present invention and particle spring coupling model;
Fig. 5 is that liver organ is respectively at the deformation design sketch that is subjected under pulling force and the pressure condition.
Embodiment
Embodiment one: referring to Fig. 1, the present invention is based on the soft tissue deformation emulating method of no grid gal the Liao Dynasty's gold and the coupling of particle spring, comprise the pre-treatment process, calculate deformation process, three parts of last handling process:
In the pre-treatment process, for soft tissue is set up the linear viscoelasticity biomechanical model.Suppose that soft tissue is even, isotropic material, have accurate incompressibility, linear viscoelasticity feature, structure constitutive equation, geometric equation and balance equation apply boundary condition and starting condition;
In calculating deformation process, dynamically divide no net region and particle spring zone according to the load that soft tissue carried, and the join domain between no net region and particle spring zone is set up transition element, the approximate displacement function of structure transition element is realized the self-adaptation coupling of no grid gal golden method of the Liao Dynasty and particle spring method;
In last handling process, according to deformation result of calculation, the particle of each time step of deformation process or the state of node are outputed on the screen, and carry out illumination and play up, finally, realize dynamic deformation effect of visualization showing the real-time deformation process of soft tissue organs under stressing conditions on the screen.
Embodiment two: present embodiment is based on the soft tissue deformation emulating method of no grid gal the Liao Dynasty gold with the coupling of particle spring, and different with embodiment one is: the process of calculating deformation specifically comprises following four steps:
1) function between design load and the distance, as the foundation of dividing no net region, and make no net region enough greatly with the assurance noncoherent boundary all in no net region, other zone is particle spring zone;
2) for no net region and particle spring zone, set up effective data structure respectively, the Classification Management data;
3) join domain between no net region and particle spring zone is set up transition element;
4) in transition element, set up the approximate displacement function of transition node and transition node, realize the smooth transition between no net region and the particle spring zone.
Embodiment three: present embodiment is based on the soft tissue deformation emulating method of no grid gal the Liao Dynasty's gold and the coupling of particle spring, different with embodiment two is: the join domain between no net region and particle spring zone is set up transition element, comprise following two steps: the no net region that 1) will divide and particle spring zone are as two entities, two entities contact part as transition boundary, on the basis of background grid, the background grid at segmentation transition boundary place and as subelement, making has a node or particle to exist at most in each subelement, with the unit that do not comprise node or particle as the gap unit;
2) with the gap unit as object search, the left-right and front-back six direction is searched for subelement up and down respectively, through iteration, progressively will not have the gap unit that inside, net region meets switch condition and be converted into no grid cell, the gap unit that particle spring intra-zone is met switch condition is converted into the particle spring unit, no longer change until the gap unit, and with remaining gap unit as transition element.
The present invention is based on the soft tissue deformation emulating method of no grid gal the Liao Dynasty's gold and the coupling of particle spring, in calculating deformation process, at first determine no net region and particle spring zone, node in two zones or particle information are stored respectively; Join domain segmentation background grid between no net region and particle spring zone guarantees at most only to comprise in each subelement a node or particle; Then, as object search, the application searches algorithm is no grid cell or particle spring unit with part gap cell translation with the gap unit, and with remaining gap unit as transition element.
Embodiment four: present embodiment is based on the soft tissue deformation emulating method of no grid gal the Liao Dynasty gold with the coupling of particle spring, and different with embodiment three is: set up transition node in transition element, determine the approximate displacement function of transition node as follows.Transition node satisfies two conditions:
1) transition node is subjected to not having the acting force in net region and particle spring zone respectively, and it is made a concerted effort is zero;
2) transition node is subjected to not having the acting force in net region and particle spring zone respectively, and its displacement equates;
According to above-mentioned condition, make up linear viscoelasticity dynamic motion equation respectively based on no grid gal the Liao Dynasty's gold and particle spring, find the solution the differential equation, draw the approximate displacement of two differences of same transition node, travel through all transition nodes, obtain the approximate displacement function that two kinds of model constructions go out, set up the funtcional relationship between the two, and as the approximate displacement function of transition node.
Embodiment five: referring to Fig. 1, present embodiment is based on the soft tissue deformation emulating method of no grid gal the Liao Dynasty gold with the coupling of particle spring, the pre-treatment process, calculate deformation process, last handling process is as follows respectively:
Pre-treatment process: at first set up the linear viscoelasticity biomechanical model for soft tissue.Suppose that soft tissue is even, isotropic material, have accurate incompressibility, linear viscoelasticity feature, with the elasticity of spring simulation viscoelastic body, damper is simulated its viscosity, and the distortion and the load of hypothesis spring are proportional, and the speed of the distortion of damper and distortion is proportional; Secondly the distortion of isotropic material can be separated into stereomutation and isopyknic shape distortion two parts, so the stress and strain tensor is resolved into its spherical tensor sum deviator, and then the constitutive equation of structure Kelvin model, simulate soft tissue stress relaxation and creep properties preferably; Set up geometric equation, balance equation at last, and apply boundary condition and starting condition.
Calculate deformation process: during operating theater instruments contact soft tissue, no matter push or lift or cut, all according to the suffered load of organ local surfaces of operating theater instruments contact, funtcional relationship between design load and the distance, thereby dynamically divide no net region and particle spring zone, make no net region and particle spring zone satisfy predetermined scope, and noncoherent boundary is finally still in no net region, continuum boundary is all in particle spring zone; Secondly establish the data structure of effect respectively in no net region and particle spring zone, and the Classification Management data; The no net region that to divide then and particle spring zone are as two entities, two entities contact part as transition boundary, the background grid at segmentation transition boundary place, and the background grid after will segmenting is as subelement, subelement satisfies the condition that has a node or particle at the most, obtains the border of two entities like this; Because the background grid after the segmentation exists and does not comprise the unit of node or particle, is defined as the gap unit so will not comprise the unit of node or particle; With the gap unit as object search, select for use and go up afterwards earlier, from left to right, vertical order is progressively searched for subelement, pass through iteration, the gap unit that progressively will meet switch condition is converted into no grid cell or particle spring unit, no longer change until the gap unit, and with remaining gap unit as the transition element between no grid and the particle spring; In transition element, set up transition node at last,, set up the approximate displacement function of transition node according to displacement and the dynamic balance condition that transition node satisfies.
The golden model of described no grid gal the Liao Dynasty is used for calculating the big zone of deformation that operating theater instruments contacts with soft tissue, adopts Pascal figure to set up three-dimensional basis function p (x), determines the number m of basis function; The background grid of the weak formula method of the structure overall situation is set point x QBe the calculation level x in the domain of influence, calculate the size d of the domain of influence s, definite field node number n that participates in the MFree shape function of structure calculation level guarantees m=n; Utilize selected node of the domain of influence, set up weight function W and derivative thereof; Adopt the node shape function φ (x) and the partial derivative thereof of all nodes in the mobile least square method structure influence territory
Figure BDA0000078959540000081
Energy functional by Viscous Elastic Mechanics governing equation correspondence obtains the golden weak variation form of gal the Liao Dynasty, and carry out totally dispersing, get and stay value, get the overall discrete governing equation of EFG method, apply essential boundary condition, utilize the Gauss of the two-dimentional equiband storage elimination method that circulates to obtain discrete nodal displacement.
Described particle spring model is used for calculating deformation than the zonule, according to Hooke's law, sets up the lagrange equation of motion of single particle movement:
m ∂ 2 u MS ∂ t 2 + γ ∂ u MS ∂ t + f int i = f ext i
In the formula, u MSThe position vector of expression particle i, m represents the quality of particle i, γ represents the stickiness density between adjacent particle,
Figure BDA0000078959540000083
The internal force that the expression spring is applied particle i,
Figure BDA0000078959540000084
The suffered external force of expression particle i according to the suffered external force of particle, is separated the differential equation and is got particle displacement.
Described no grid gal the Liao Dynasty gold and particle spring coupling model are used for realizing not having the natural transition between net region and the particle spring zone.The transition node of setting up in transition element satisfies two conditions: 1) transition node is subjected to not having the acting force in net region and particle spring zone respectively, and it is made a concerted effort is zero; 2) transition node is subjected to not having the acting force in net region and particle spring zone respectively, and its displacement equates.Now according to above-mentioned condition, make up the linear viscoelasticity dynamic simulation equation of motion respectively based on no grid and particle spring, find the solution the approximate displacement of two differences that the differential equation gets same transition node, travel through all transition nodes, obtain the approximate displacement function that uses two kinds of model constructions to go out, set up the funtcional relationship between the two, and as the approximate displacement function of transition node.
Suppose that t=0 displacement, speed and acceleration constantly is known, to find the solution the time territory 0～T be divided into n time interval Δ t (=T/n), utilize central difference method or Newmark method to calculate the useful load and the displacement of each time step, ask strain according to displacement, and by new bulk strain and the shape distortion of strain calculating, according to rheological equations, volume calculated stress, deviatoric stress and particle stress;
Last handling process: carry out playing up of illumination according to the deformation result calculated, finally real-time rendering goes out the particle of each time step or the state of node on screen, like this, just can realize dynamic deformation effect of visualization.
Embodiment six: present embodiment elaborates to technical solution of the present invention in conjunction with Fig. 1～Fig. 5:
As shown in Figure 1, when the deformation of simulation soft tissue, mainly be divided into three big modules and implement:
1, pre-treatment process
At first, set up the linear viscoelasticity biomechanical model for soft tissue.Suppose that soft tissue is even, isotropic material, have accurate incompressibility, linear viscoelasticity feature, this model combines by the Maxwell model is in parallel with Hookean spring, and the supposition Hookean spring produces the distortion that is directly proportional with load at once, damper can produce the speed that is directly proportional with load in arbitrary moment, so constitutive equation is:
σ + η E 2 σ · = E 1 ϵ + η ( 1 + E 1 E 2 ) ϵ ·
σ represents stress in the formula, and ε represents strain, and E represents the elasticity coefficient of Hookean spring, E 1The elasticity coefficient of representing the Hookean spring in parallel with the Maxwell model, E 2The elasticity coefficient of expression Maxwell model neutral line spring, η represents the coefficient of viscosity of damper, With
Figure BDA0000078959540000093
Represent the derivative of stress and strain respectively about the time.
The material parameter of Kelvin standard lines viscoelastic model is obtained more or less freely, and can realize somatometry, therefore can be by adjusting the lax and creep properties that is had when parameter is described soft tissue generation deformation, and wherein creep function is:
C ( t ) = 1 E 1 [ 1 - ( 1 - τ ϵ τ σ e - t / τ σ ) ] I ( t )
Relaxation function is:
K ( t ) = E 1 [ 1 - ( 1 - τ σ τ ϵ e - t / τ ϵ ) ] I ( t )
In the formula τ ϵ = η E 2 , τ σ = η E 1 ( 1 + E 1 E 2 ) ,
Figure BDA0000078959540000098
Generally speaking, the distortion of isotropic material can be separated into stereomutation and isopyknic shape distortion two parts, and promptly volumetric stress only changes volume, and the deviator of stress causes isopyknic shape distortion, so the stress and strain tensor is resolved into its spherical tensor sum deviator part, promptly
σ mn = s mn + 1 3 δ mn σ ii , ϵ mn = e mn + 1 3 δ mn ϵ ii , m，n＝1，2，3
Wherein, δ MnBe the Kronecker symbol; ε IiAnd σ IiBe respectively bulk strain and volumetric stress; e MnAnd s MnBe respectively the component of the strain deviation and the deviator of stress; e Mm=0 and s Mn=0.
Because viscous-elastic behaviour and effect under two kinds of situations can consider respectively, so according to Kelvin normal linearity Viscous Elastic Mechanics model, the three-dimensional viscoelastic constitutive equation between deviatoric stress tensor and the inclined to one side strain tensor, between volumetric stress and the bulk strain is:
s mn + η E 2 s · mn = E 1 e mn + η ( 1 + E 1 E 2 ) e · mn , m，n＝1，2，3
σ ii + η E 2 σ · ii = E 1 ϵ ii + η ( 1 + E 1 E 2 ) ϵ · ii
Wherein
Figure BDA0000078959540000102
Be respectively inclined to one side rate of strain and volumetric strain.
2, calculate deformation process
1) sets up the golden model of no grid gal the Liao Dynasty
The golden model of described no grid gal the Liao Dynasty is used for calculating the big zone of deformation that operating theater instruments contacts with soft tissue, adopts Pascal figure to set up three-dimensional basis function p (x), and the number m of definite basis function,
p T(x)=[1 x y z] TM=4, p=1 is linear basis function;
p T(x)=[1 x y z x 2Xy y 2Yz z 2Zx] TM=10, p=2 is the secondary base function;
p T(x)=[1 x y z x 2Xy y 2Yz z 2Zx x 3x 2Y xy 2y 3y 2Z yz 2z 3z 2X zx 2Xyz] TM=20, p=3 is three basis functions;
The background grid of the weak formula method of the structure overall situation is set point x QBe the calculation level x in the domain of influence, by setting the size d of the domain of influence s
d s = a s d c , d c = V s 3 n v s 3 - 1
Determine the field node n of the MFree shape function of participation structure calculation level, guarantee m=n.A in the formula sBe the dimensionless size of this domain of influence, d cFor being positioned at an x QNear nodal pitch, V sVolume for the pre-estimation domain of influence;
Figure BDA0000078959540000104
For being included in volume V sIn the pre-estimation domain of influence in the node number.
Utilize selected node of the domain of influence, set up weight function W and derivative thereof, existing is example with cubic spline function as weight function:
W i ( x ) = 2 / 3 - 4 r i 2 + 4 r i 3 r i &le; 0.5 4 / 3 - 4 r i + 4 r i 2 - 4 / 3 r i 3 0.5 < r i &le; 1 0 r i > 1
In the formula d i=| x-x i| be node x iAnd the distance between the calculation level x, r wBe weight function support region size.
In all directions differentiate, the first order derivative that gets weight function is to weight function:
W i ( x ) = - 8 r i + 12 r i 2 r i &le; 0.5 - 4 + 8 r i - 4 r i 2 0.5 < r i &le; 1 0 r i > 1
Adopt the node shape function φ (x) and the partial derivative thereof of all nodes in the mobile least square method structure influence territory
Figure BDA0000078959540000111
Energy functional by Viscous Elastic Mechanics governing equation correspondence obtains the golden weak variation form of gal the Liao Dynasty,
&Integral; &Omega; &delta; ( Lu ) T D ( Lu ) d&Omega; - &Integral; &Omega; &delta; u T bd&Omega; - &Integral; &Gamma; t &delta; u T t &OverBar; d&Gamma; - &delta; &Integral; &Gamma; u 1 2 ( u - u &OverBar; ) T &alpha; ( u - u &OverBar; ) d&Gamma; = 0
Above-mentioned equation is carried out totally dispersing, gets and stay value, get the overall discrete governing equation of EFG method,
[K+K α]U＝F+F α
In the formula, U is the nodal displacement parameter of all nodes in the whole Problem Areas, the global stiffness matrix of K for being assembled by the node stiffness matrix, the overall outer force vector of F for being assembled by the nodal force vector, K αThe added martix that the overall punishment stiffness matrix of serving as reasons causes, F αFor by the caused additional force vector of essential boundary condition, K, K α, F and F αDo three-dimensional Gauss integration or the Irons integration obtains by the point of the Gauss integration in the background grid; Apply essential boundary condition, utilize the Gauss of the two-dimentional equiband storage elimination method that circulates to obtain discrete nodal displacement.
2) set up the particle spring model
Described particle spring model is used for calculating deformation than the zonule, according to Hooke's law, sets up the lagrange equation of motion of single quality point motion
m &PartialD; 2 u MS &PartialD; t 2 + &gamma; &PartialD; u MS &PartialD; t + f int i = f ext i
Wherein: u MSThe position vector of expression particle i; M represents the quality of particle i, and γ represents the stickiness density between adjacent particle,
Figure BDA0000078959540000114
The internal force that the expression spring is applied particle i,
Figure BDA0000078959540000115
The suffered external force of expression particle i according to the suffered external force of particle, is separated partial differential equation and is got particle displacement.
3) set up no grid gal the Liao Dynasty's gold and particle spring coupling model
For a soft tissue, stressed when a certain position, then the deformation that takes place with the point of stress point close together is more violent, and more small away from the histogenetic deformation ratio of stress point.In order to guarantee computational accuracy, improve counting yield, the golden model of no grid gal the Liao Dynasty is adopted in the big zone of deformation, and deformation is adopted the particle spring model than the zonule.The existing load of being born according to the organ local surfaces that operating theater instruments contacted, function between design load and the distance, and with its foundation as no net region of division and particle spring zone, wherein there is not the net region and predetermined scope will be satisfied in particle spring zone, to guarantee noncoherent boundary in no net region, continuum boundary is in particle spring zone; Secondly in no net region and particle spring zone, set up effective data structure respectively, and the Classification Management data; The no net region that to divide then and particle spring zone are as two entities, two entities contact part as transition boundary, on the basis of background grid, the background grid at segmentation transition boundary place, and the background grid after will segmenting is as subelement, and the condition that subelement satisfies is that node or particle exist at the most in the unit; The unit that definition does not comprise node or particle is the gap unit, now with the gap unit as object search, according to the sequential search subelement before and after up and down, through iteration, progressively will not have the gap unit that inside, net region meets switch condition and be converted into no grid cell, the gap unit that particle spring intra-zone is met switch condition is converted into the particle spring unit, and with remaining gap unit as transition element, wherein
(1) the gap unit of inside, no net region is converted into the switch condition of no grid cell and is:
1. top to bottom, left and right, front and rear three prescriptions are no grid cell to the subelement of middle existence direction more than two groups and two groups;
2. six subelements on the six direction of top to bottom, left and right, front and rear all are the gap unit;
3. for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situations that all exist, and every prescription is to getting a subelement, and three subelements on totally three directions are no grid cell;
4. for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situation that all exists, wherein the subelement that makes progress of a prescription is all got, in addition two prescriptions in every group of subelement of getting on one of them direction, four subelements that amount on the four direction are no grid cell;
5. for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situation that all exists, wherein the subelement that makes progress of two prescriptions is all got, another prescription is upwards got the subelement on one of them direction, and five subelements that amount on five directions are no grid cell;
(2) the gap unit of particle spring intra-zone is converted into the switch condition of particle spring unit and is:
1. top to bottom, left and right, front and rear three prescriptions are the particle spring unit to the subelement of middle existence direction more than two groups and two groups;
2. six subelements on the six direction of top to bottom, left and right, front and rear all are the gap unit;
3. for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situations that all exist, and every prescription is to getting a subelement, and three subelements on totally three directions are the particle spring unit;
4. for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situation that all exists, wherein the subelement that makes progress of a prescription is all got, in addition two prescriptions in every group of subelement of getting on one of them direction, four subelements that amount on the four direction are the particle spring unit;
5. for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situation that all exists, wherein the subelement that makes progress of two prescriptions is all got, another prescription is upwards got the subelement on one of them direction, and five subelements that amount on five directions are the particle spring unit.
(3) being characterized as of remaining gap unit:
1. do not meet no grid cell switch condition, do not meet the gap unit of particle spring unit switch condition yet;
2. six subelements on the six direction of top to bottom, left and right, front and rear comprise no grid cell and particle spring unit simultaneously, perhaps comprise no grid cell, particle spring unit and gap unit simultaneously.
After determining transition element, in transition element, set up transition node, and transition node satisfies two equilibrium conditions:
(1) transition node is subjected to not having the acting force in net region and particle spring zone respectively, and it is made a concerted effort is zero;
(2) transition node is subjected to not having the acting force in net region and particle spring zone respectively, and its displacement equates;
According to above-mentioned condition, make up the linear viscoelasticity dynamic simulation equation of motion respectively based on no grid gal the Liao Dynasty's gold and particle spring, find the solution the differential equation, draw the approximate displacement of two differences of same transition node, travel through all transition nodes, obtain the approximate displacement function that two kinds of model constructions go out, set up the funtcional relationship between the two, and, realize the smooth transition between no net region and the particle spring zone as the approximate displacement function of transition node.
Suppose that at last t=0 displacement, speed and acceleration constantly is known, to find the solution the time territory 0～T be divided into n time interval Δ t (=T/n), utilize central difference method or Newmark method to calculate the useful load and the displacement of each time step, ask strain according to displacement, and by new bulk strain and the shape distortion of strain calculating, according to the constitutive equation of Kelvin standard lines viscoelasticity biomechanical model, volume calculated stress, deviatoric stress and particle stress.
3, last handling process
Carry out illumination according to the deformation result calculated and play up, finally real-time rendering goes out the particle of each time step or the state of node on screen, realizes dynamic deformation effect of visualization.
Shown in Fig. 2,3,4, the load of being born at first according to the organ local surfaces that operating theater instruments contacted, function between design load and the distance, dynamically divide no net region and particle spring zone, between no net region and particle spring zone, set up a series of transition element then, construct corresponding displacement approximate function at last, two kinds of methods are coupled to together.When the load that liver applied up or down, by based on the soft tissue deformation computation model of no grid gal the Liao Dynasty gold with the coupling of particle spring, simulation true to nature operative region soft tissue deformation effect, as shown in Figure 5.With respect to the finite element model that generally adopts at present, this model has reduced overall complexity, has accelerated computing velocity.In addition, this model is convenient to realize operation technique emulation such as cutting, stitching, can provide powerful support for force feedback equipment, can be used as soft tissue physical model in the virtual operation.

Claims (6)

1. the soft tissue deformation emulating method based on no grid gal the Liao Dynasty's gold and the coupling of particle spring comprises the pre-treatment process, calculates deformation process, three steps of last handling process, in the pre-treatment process, for soft tissue is set up the linear viscoelasticity biomechanical model; In calculating deformation process, dynamically divide no net region and particle spring zone according to the load that soft tissue carried, and the join domain between no net region and particle spring zone is set up transition element, the approximate displacement function of structure transition element is realized the self-adaptation coupling of no grid gal golden method of the Liao Dynasty and particle spring method; In last handling process, according to deformation result of calculation, the particle of each time step of deformation process or the state of node are outputed on the screen, and carry out illumination and play up, finally, realize dynamic deformation effect of visualization showing the real-time deformation process of soft tissue organs under stressing conditions on the screen; It is characterized in that: the process of described calculating deformation specifically comprises following four steps:
1) function between design load and the distance, as the foundation of dividing no net region, and make no net region enough greatly with the assurance noncoherent boundary all in no net region, other zone is particle spring zone;
2) for no net region and particle spring zone, set up effective data structure respectively, the Classification Management data;
3) join domain between no net region and particle spring zone is set up transition element:
(1) the no net region that will divide and particle spring zone are as two entities, two entities contact part as transition boundary, on the basis of background grid, the background grid at segmentation transition boundary place and as subelement, making has a node or particle to exist at most in each subelement, with the unit that do not comprise node or particle as the gap unit;
(2) with the gap unit as object search, the left-right and front-back six direction is searched for subelement up and down respectively, through iteration, progressively will not have the gap unit that inside, net region meets switch condition and be converted into no grid cell, the gap unit that particle spring intra-zone is met switch condition is converted into the particle spring unit, no longer change until the gap unit, and with remaining gap unit as transition element;
4) in transition element, set up transition node, determine the approximate displacement function of transition node, realize the smooth transition between no net region and the particle spring zone; Two conditions are satisfied in the foundation of transition node:
(1) transition node is subjected to not having the acting force in net region and particle spring zone respectively, and it is made a concerted effort is zero;
(2) transition node is subjected to not having the acting force in net region and particle spring zone respectively, and its displacement equates;
According to above-mentioned condition, make up linear viscoelasticity dynamic motion equation respectively based on no grid gal the Liao Dynasty's gold and particle spring, find the solution the differential equation, draw the approximate displacement of two differences of same transition node, travel through all transition nodes, obtain the approximate displacement function that two kinds of model constructions go out, set up the funtcional relationship between the two, as the approximate displacement function of transition node.
2. the soft tissue deformation emulating method based on no grid gal the Liao Dynasty's gold and the coupling of particle spring according to claim 1, it is characterized in that: the switch condition that inner gap unit, no net region is converted into no grid cell is:
1) top to bottom, left and right, front and rear three prescriptions are no grid cell to the subelement of middle existence direction more than two groups and two groups;
2) six subelements on the six direction of top to bottom, left and right, front and rear all are the gap unit;
3) for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situations that all exist, and every prescription is to getting a subelement, and three subelements on totally three directions are no grid cell;
4) for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situation that all exists, wherein the subelement that makes progress of a prescription is all got, in addition two prescriptions in every group of subelement of getting on one of them direction, four subelements that amount on the four direction are no grid cell;
5) for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situation that all exists, wherein the subelement that makes progress of two prescriptions is all got, another prescription is upwards got the subelement on one of them direction, and five subelements that amount on five directions are no grid cell.
3. the soft tissue deformation emulating method based on no grid gal the Liao Dynasty's gold and the coupling of particle spring according to claim 1, it is characterized in that: the gap unit that particle spring intra-zone is met switch condition is converted into the particle spring unit, and its switch condition is:
1) top to bottom, left and right, front and rear three prescriptions are the particle spring unit to the subelement of middle existence direction more than two groups and two groups;
2) six subelements on the six direction of top to bottom, left and right, front and rear all are the gap unit;
3) for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situations that all exist, and every prescription is to getting a subelement, and three subelements on totally three directions are the particle spring unit;
4) for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situation that all exists, wherein the subelement that makes progress of a prescription is all got, in addition two prescriptions in every group of subelement of getting on one of them direction, four subelements that amount on the four direction are the particle spring unit;
5) for borderline unit, every group of subelement that top to bottom, left and right, front and rear three prescriptions make progress is not the situation that all exists, wherein the subelement that makes progress of two prescriptions is all got, another prescription is upwards got the subelement on one of them direction, and five subelements that amount on five directions are the particle spring unit.
4. the soft tissue deformation emulating method based on no grid gal the Liao Dynasty's gold and the coupling of particle spring according to claim 1, it is characterized in that: remaining gap unit is meant as transition element:
1) do not meet no grid cell switch condition, do not meet the gap unit of particle spring unit switch condition yet;
2) six subelements on the six direction of top to bottom, left and right, front and rear comprise no grid cell and particle spring unit simultaneously, perhaps comprise no grid cell, particle spring unit and gap unit simultaneously.
5. according to each described soft tissue deformation emulating method of claim 1～4 based on no grid gal the Liao Dynasty's gold and the coupling of particle spring, it is characterized in that: in the pre-treatment process, suppose that soft tissue is even, isotropic material, have accurate incompressibility, linear viscoelasticity feature, and definite initialization data, definition starting condition, displacement boundary conditions and stress boundary condition; With the elasticity of spring simulation viscoelastic body, damper is simulated its viscosity then, and the distortion and the load of hypothesis spring are proportional, and the speed of the distortion of damper and distortion is proportional; Secondly the distortion of isotropic material can be separated into stereomutation and isopyknic shape distortion two parts, so the stress and strain tensor is resolved into its spherical tensor sum deviator, and then the constitutive equation of structure Kelvin model, simulation soft tissue stress relaxation and creep properties; Set up geometric equation, balance equation at last, and apply boundary condition and starting condition, for soft tissue is set up the linear viscoelasticity biomechanical model.
6. the soft tissue deformation emulating method based on no grid gal the Liao Dynasty's gold and the coupling of particle spring according to claim 5, it is characterized in that: in calculating deformation process, at first generate initial background grid, divide no net region and particle spring zone according to external load, foundation from no grid to the data structure the particle spring data, to reduce the frequency that the data reusing phenomenon takes place in the adaptive process, secondly set up transition element at join domain, in transition element, set up the approximate displacement function of transition node and transition node, realize the smooth transition between no net region and the particle spring zone, apply boundary condition then, in conjunction with Kelvin normal linearity Viscous Elastic Mechanics model, adopt method of addition to find the solution displacement, ask stress according to displacement, strain.
CN201110213387A 2011-07-27 2011-07-27 Soft tissue deformation simulation method based on coupling of mesh-free Galerkin and mass spring Expired - Fee Related CN102262699B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110213387A CN102262699B (en) 2011-07-27 2011-07-27 Soft tissue deformation simulation method based on coupling of mesh-free Galerkin and mass spring

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110213387A CN102262699B (en) 2011-07-27 2011-07-27 Soft tissue deformation simulation method based on coupling of mesh-free Galerkin and mass spring

Publications (2)

Publication Number Publication Date
CN102262699A true CN102262699A (en) 2011-11-30
CN102262699B CN102262699B (en) 2012-09-05

Family

ID=45009323

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110213387A Expired - Fee Related CN102262699B (en) 2011-07-27 2011-07-27 Soft tissue deformation simulation method based on coupling of mesh-free Galerkin and mass spring

Country Status (1)

Country Link
CN (1) CN102262699B (en)

Cited By (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102831280A (en) * 2012-09-10 2012-12-19 北京航空航天大学 Meshless physical deformation simulation method based on moving least squares
CN103310072A (en) * 2013-06-28 2013-09-18 哈尔滨理工大学 Thighbone biomechanics finite element analysis system based on force feedback
CN103400023A (en) * 2013-06-28 2013-11-20 华北水利水电大学 Soft tissue deformation simulation method
CN103699714A (en) * 2013-12-01 2014-04-02 北京航空航天大学 Flexible object real-time cutting simulation method based on finite element and meshless coupling
WO2014064066A1 (en) 2012-10-26 2014-05-01 Brainlab Ag Simulation of objects in an atlas and registration of patient data containing a specific structure to atlas data
CN103793552A (en) * 2013-12-18 2014-05-14 哈尔滨工程大学 Real-time dynamic generating method for local particle spring model with deformed soft tissues
CN103970960A (en) * 2014-05-23 2014-08-06 湘潭大学 Grid-free Galerkin method structural topology optimization method based on GPU parallel acceleration
CN105426339A (en) * 2015-11-06 2016-03-23 吉林大学 Line source time domain electromagnetic response numerical calculation method based on meshless method
CN105513130A (en) * 2016-02-01 2016-04-20 福建师范大学福清分校 Soft tissue deformation method based on mixing of gridding method and non-gridding method
CN105912859A (en) * 2016-04-11 2016-08-31 浙江工业大学义乌科学技术研究院有限公司 Tissue deformation method based on mass spring and hydromechanics
CN106156537A (en) * 2016-07-04 2016-11-23 南昌大学 The Force Feedback Model modeling method of radial direction base mesh free soft tissue data based on wheat quart algorithm
CN106570341A (en) * 2016-11-14 2017-04-19 南昌大学 Soft tissue stress deformable model modeling method based on mesh-free radial base data fitting
CN106650251A (en) * 2016-12-14 2017-05-10 南京信息工程大学 Modeling method of acupuncture force feedback deformable model
CN106777658A (en) * 2016-12-12 2017-05-31 南京信息工程大学 Based on the suture analogy method for loosening LRA bounding algorithms
CN106781941A (en) * 2016-11-24 2017-05-31 北京理工大学 A kind of method and its system for simulating microtrauma puncture operation
CN108446507A (en) * 2018-03-26 2018-08-24 中国人民解放军国防科技大学 Elastic body grid deformation method based on grid quality feedback optimization
CN108536936A (en) * 2018-03-27 2018-09-14 南京信息工程大学 A kind of mesh free soft tissue deformation analogy method of Multiple Optimization
CN108710735A (en) * 2018-05-08 2018-10-26 南京信息工程大学 A kind of mesh free soft tissue deformation analogy method of real-time, interactive
CN108877944A (en) * 2018-06-26 2018-11-23 南京信息工程大学 Virtual cropping method based on the grid model for being included in Kelvin's viscoelastic model
CN109344500A (en) * 2018-09-29 2019-02-15 南京信息工程大学 A kind of Soft Tissue Deformation analogy method of non-resilient mass-spring model and improved euler algorithm
CN109658796A (en) * 2018-11-14 2019-04-19 华中科技大学 A kind of soft tissue injury and abnormal connection simulation experiment method
US10388013B2 (en) 2012-10-26 2019-08-20 Brainlab Ag Matching patient images and images of an anatomical atlas
CN111317566A (en) * 2018-12-14 2020-06-23 西门子医疗有限公司 Planning support for interventional procedures
CN113343513A (en) * 2021-05-11 2021-09-03 南京信息工程大学 Method and device for simulating soft tissue deformation and path cutting
CN113470165A (en) * 2021-06-17 2021-10-01 南昌大学 Soft tissue modeling method based on radial basis point interpolation method and mass point spring method
CN116434885A (en) * 2023-03-20 2023-07-14 中国工程物理研究院研究生院 Viscoelastic fracture phase field calculation method considering tension and compression anisotropy of polymer composite material

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1975784A (en) * 2006-12-28 2007-06-06 上海交通大学 Point particle spring deformation simulating method based on skeleton linear net
CN102044086A (en) * 2010-11-30 2011-05-04 华北水利水电学院 Soft tissue deformation simulation method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1975784A (en) * 2006-12-28 2007-06-06 上海交通大学 Point particle spring deformation simulating method based on skeleton linear net
CN102044086A (en) * 2010-11-30 2011-05-04 华北水利水电学院 Soft tissue deformation simulation method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《2010 International Conference on Computer and Communication Technologies in Agriculture Engineering》 20100613 Xuemei Liu, Lei Mao. Visual Simulation of Soft Tissue Deformation 548-551 1-6 , *

Cited By (43)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102831280B (en) * 2012-09-10 2014-12-03 北京航空航天大学 Meshless physical deformation simulation method based on moving least squares
CN102831280A (en) * 2012-09-10 2012-12-19 北京航空航天大学 Meshless physical deformation simulation method based on moving least squares
US10388013B2 (en) 2012-10-26 2019-08-20 Brainlab Ag Matching patient images and images of an anatomical atlas
US10402971B2 (en) 2012-10-26 2019-09-03 Brainlab Ag Matching patient images and images of an anatomical atlas
US10417762B2 (en) 2012-10-26 2019-09-17 Brainlab Ag Matching patient images and images of an anatomical atlas
WO2014064066A1 (en) 2012-10-26 2014-05-01 Brainlab Ag Simulation of objects in an atlas and registration of patient data containing a specific structure to atlas data
CN103400023B (en) * 2013-06-28 2016-11-02 华北水利水电大学 Soft tissue deformation simulation method
CN103400023A (en) * 2013-06-28 2013-11-20 华北水利水电大学 Soft tissue deformation simulation method
CN103310072A (en) * 2013-06-28 2013-09-18 哈尔滨理工大学 Thighbone biomechanics finite element analysis system based on force feedback
CN103699714B (en) * 2013-12-01 2016-08-31 北京航空航天大学 A kind of flexible article real-time cutting simulation method coupled based on finite element and mesh free
CN103699714A (en) * 2013-12-01 2014-04-02 北京航空航天大学 Flexible object real-time cutting simulation method based on finite element and meshless coupling
CN103793552A (en) * 2013-12-18 2014-05-14 哈尔滨工程大学 Real-time dynamic generating method for local particle spring model with deformed soft tissues
CN103793552B (en) * 2013-12-18 2017-07-11 哈尔滨工程大学 A kind of real-time dynamic creation method of the local Mass-spring Model of soft tissue deformation
CN103970960A (en) * 2014-05-23 2014-08-06 湘潭大学 Grid-free Galerkin method structural topology optimization method based on GPU parallel acceleration
CN103970960B (en) * 2014-05-23 2016-11-23 湘潭大学 The element-free Galerkin structural topological optimization method accelerated parallel based on GPU
CN105426339A (en) * 2015-11-06 2016-03-23 吉林大学 Line source time domain electromagnetic response numerical calculation method based on meshless method
CN105426339B (en) * 2015-11-06 2018-05-29 吉林大学 A kind of line source time domain electromagnetic response numerical computation method based on gridless routing
CN105513130A (en) * 2016-02-01 2016-04-20 福建师范大学福清分校 Soft tissue deformation method based on mixing of gridding method and non-gridding method
CN105513130B (en) * 2016-02-01 2018-02-16 福建师范大学福清分校 A kind of soft tissue deformation method mixed based on grid with mesh free
CN105912859A (en) * 2016-04-11 2016-08-31 浙江工业大学义乌科学技术研究院有限公司 Tissue deformation method based on mass spring and hydromechanics
CN105912859B (en) * 2016-04-11 2018-07-17 浙江工业大学义乌科学技术研究院有限公司 One kind being based on mass spring and hydromechanical tissue deformation analogy method
CN106156537B (en) * 2016-07-04 2018-11-16 南昌大学 The Force Feedback Model modeling method of radial base mesh free soft tissue data based on wheat quart algorithm
CN106156537A (en) * 2016-07-04 2016-11-23 南昌大学 The Force Feedback Model modeling method of radial direction base mesh free soft tissue data based on wheat quart algorithm
CN106570341B (en) * 2016-11-14 2019-06-18 南昌大学 A kind of soft tissue stress and deformation model modelling approach based on the fitting of mesh free radial direction base data
CN106570341A (en) * 2016-11-14 2017-04-19 南昌大学 Soft tissue stress deformable model modeling method based on mesh-free radial base data fitting
CN106781941A (en) * 2016-11-24 2017-05-31 北京理工大学 A kind of method and its system for simulating microtrauma puncture operation
CN106777658A (en) * 2016-12-12 2017-05-31 南京信息工程大学 Based on the suture analogy method for loosening LRA bounding algorithms
CN106777658B (en) * 2016-12-12 2020-02-18 南京信息工程大学 Stitching simulation method based on relaxed LRA constraint algorithm
CN106650251B (en) * 2016-12-14 2018-08-14 南京信息工程大学 A kind of modeling method of acupuncture force feedback deformation model
CN106650251A (en) * 2016-12-14 2017-05-10 南京信息工程大学 Modeling method of acupuncture force feedback deformable model
CN108446507A (en) * 2018-03-26 2018-08-24 中国人民解放军国防科技大学 Elastic body grid deformation method based on grid quality feedback optimization
CN108536936A (en) * 2018-03-27 2018-09-14 南京信息工程大学 A kind of mesh free soft tissue deformation analogy method of Multiple Optimization
CN108710735A (en) * 2018-05-08 2018-10-26 南京信息工程大学 A kind of mesh free soft tissue deformation analogy method of real-time, interactive
CN108877944A (en) * 2018-06-26 2018-11-23 南京信息工程大学 Virtual cropping method based on the grid model for being included in Kelvin's viscoelastic model
CN109344500A (en) * 2018-09-29 2019-02-15 南京信息工程大学 A kind of Soft Tissue Deformation analogy method of non-resilient mass-spring model and improved euler algorithm
CN109658796A (en) * 2018-11-14 2019-04-19 华中科技大学 A kind of soft tissue injury and abnormal connection simulation experiment method
CN109658796B (en) * 2018-11-14 2021-04-20 华中科技大学 Soft tissue injury and abnormal connection simulation test method
CN111317566A (en) * 2018-12-14 2020-06-23 西门子医疗有限公司 Planning support for interventional procedures
CN111317566B (en) * 2018-12-14 2023-07-14 西门子医疗有限公司 Planning support for interventional procedures
CN113343513A (en) * 2021-05-11 2021-09-03 南京信息工程大学 Method and device for simulating soft tissue deformation and path cutting
CN113470165A (en) * 2021-06-17 2021-10-01 南昌大学 Soft tissue modeling method based on radial basis point interpolation method and mass point spring method
CN116434885A (en) * 2023-03-20 2023-07-14 中国工程物理研究院研究生院 Viscoelastic fracture phase field calculation method considering tension and compression anisotropy of polymer composite material
CN116434885B (en) * 2023-03-20 2023-09-12 中国工程物理研究院研究生院 Viscoelastic fracture phase field calculation method considering tension and compression anisotropy of polymer composite material

Also Published As

Publication number Publication date
CN102262699B (en) 2012-09-05

Similar Documents

Publication Publication Date Title
CN102262699B (en) Soft tissue deformation simulation method based on coupling of mesh-free Galerkin and mass spring
CN102044086B (en) Soft tissue deformation simulation method
Chen et al. Physically-based animation of volumetric objects
Krafczyk et al. Two-dimensional simulation of fluid–structure interaction using lattice-Boltzmann methods
Almeida et al. Finite element formulations for hyperelastic transversely isotropic biphasic soft tissues
Tsubota et al. Simulation study on effects of hematocrit on blood flow properties using particle method
CN104317985B (en) Fluid simulation method based on inter-belt finite element and Lagrange coordinate
CN103400023A (en) Soft tissue deformation simulation method
CN110298105A (en) The CCPDI-IMPM method of saturated porous media analysis on Large Deformation
CN110289104B (en) Soft tissue compression and deformation recovery simulation method
CN102156785A (en) Dynamic cloth simulation method
CN106504329A (en) A kind of gingiva deformity emulation mode of the Mass-spring Model based on tooth long axile
CN108984829A (en) The calculation method and system of rock-fill concrete rockfill banking process
CN112380764B (en) Gas scene end-to-end rapid reconstruction method under limited view
CN109101752A (en) A kind of complexity hydraulic structure partial structurtes natural frequency of vibration calculation method
Marinkovic et al. Towards real-time simulation of deformable structures by means of co-rotational finite element formulation
KR101350732B1 (en) Multi-Resolution Meshless Method for Real-Time Simulation of Deformable Objects
Arbain et al. Vibration analysis of Kenyir dam power station structure using a real scale 3D model
McDaniel et al. Efficient mesh deformation for computational stability and control analyses on unstructured viscous meshes
CN105912859B (en) One kind being based on mass spring and hydromechanical tissue deformation analogy method
Tan et al. An improved soft tissue deformation simulation model based on mass spring
Ge et al. Blending isogeometric and Lagrangian elements in three-dimensional analysis
Pederzani et al. A numerical method for the analysis of flexible bodies in unsteady viscous flows
CN106156537A (en) The Force Feedback Model modeling method of radial direction base mesh free soft tissue data based on wheat quart algorithm
Van Liedekerke et al. A particle based model to simulate plant cells dynamics

Legal Events

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

Granted publication date: 20120905

Termination date: 20130727