CN108461149A - A kind of blood analogy method based on PBF - Google Patents
A kind of blood analogy method based on PBF Download PDFInfo
- Publication number
- CN108461149A CN108461149A CN201810087356.1A CN201810087356A CN108461149A CN 108461149 A CN108461149 A CN 108461149A CN 201810087356 A CN201810087356 A CN 201810087356A CN 108461149 A CN108461149 A CN 108461149A
- Authority
- CN
- China
- Prior art keywords
- particle
- blood
- density
- pbf
- information
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 210000004369 blood Anatomy 0.000 title claims abstract description 41
- 239000008280 blood Substances 0.000 title claims abstract description 41
- 239000002245 particle Substances 0.000 claims abstract description 136
- 238000004088 simulation Methods 0.000 claims abstract description 25
- 230000008569 process Effects 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 4
- 238000006073 displacement reaction Methods 0.000 claims description 19
- 230000017531 blood circulation Effects 0.000 claims description 2
- 239000012530 fluid Substances 0.000 abstract description 35
- 230000004048 modification Effects 0.000 abstract description 4
- 238000012986 modification Methods 0.000 abstract description 4
- 230000008859 change Effects 0.000 description 3
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 210000002381 plasma Anatomy 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 235000019441 ethanol Nutrition 0.000 description 1
- 239000004519 grease Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000011499 joint compound Substances 0.000 description 1
- 239000008267 milk Substances 0.000 description 1
- 210000004080 milk Anatomy 0.000 description 1
- 235000013336 milk Nutrition 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 210000002966 serum Anatomy 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 229940034610 toothpaste Drugs 0.000 description 1
- 239000000606 toothpaste Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
Landscapes
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Public Health (AREA)
- Health & Medical Sciences (AREA)
- Pathology (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Biomedical Technology (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
The blood analogy method based on PBF that the embodiment of the invention discloses a kind of, wherein this method includes：The information data of particle is obtained, initialization process is carried out；According to the vessel position at each position of human body, the blood in blood vessel is subjected to sliding-model control, obtains particle discrete one by one；A variety of particles are subjected to calculation processing, obtain the multidate information of the field particle of each particle；The position of particle is corrected using PBD frames according to the data information of the field particle of each particle, continuous iteration is until terminate.By implementing the embodiment of the present invention, the limitation of conventional fluid analogue technique can be broken through so that the time step with bigger in simulation process has better real-time；In conjunction with a variety of particle hybrid analog-digital simulation blood, can modification appropriate be carried out according to the difference of application scenarios by using the quantity of particle, simultaneously the distribution situation of different densities particle is constrained using the density of blood entirety, the phenomenon that distribution situation that ensure that different type particle in unit volume is not in layering.
Description
Technical field
The present invention relates to computer graphics techniques field more particularly to a kind of blood analogy methods based on PBF.
Background technology
It is always the popular direction of research in the simulation of field of Computer Graphics, fluid, wherein two major classes are broadly divided into,
Lagrangian method and Euler's method.Euler's method lays particular emphasis on the research of " field ", the property of fluid, such as quality, density and temperature
Degree etc., the function being defined as between the added-time of spatial position.Opposite, Lagrangian method lays particular emphasis on the research of " particle ", fluid
Be abstracted as particle discrete one by one, each particle carries the physical attribute of oneself, by track the movement locus of particle come
Carry out the simulation of fluid.SPH methods are the Main Branches of Lagrangian method.
SPH methods are divided into two big directions, and one kind is micro- compressible SPH methods, such method is incompressible in simulation
When fluid, the density of fluid is allowed to have small fluctuation, so as to solve the pressure suffered by particle by density
Power；A kind of method is incompressible SPH methods, and when simulating incompressible fluid, fluid density will not be sent out such method
Raw fluctuation, the pressure of fluid is calculated by calculating the Poisson's equation of pressure.For micro- compressible SPH methods, it has
Intuitive feature is calculated, each individually particle can calculate stressing conditions using Newton's second law；It require that
The time step for having very little, to ensure situation that the density of fluid is not in out of control.For incompressible SPH methods, it
It can accurately ensure the constant of fluid density during calculating, but the calculating that it is often walked needs successive ignition, convergence
Speed is slow.And the implicit incompressible SPH methods proposed later are improved on the basis of incompressible SPH methods, are carried
High convergence rate when solving Poisson's equation.
Blood simulation is a big branch of fluid simulation, for this fluid of water, property that blood has its special
Matter.Fluid can be divided into two major class, and one kind is Newtonian fluid, and one kind is non-newtonian fluid.For Newtonian fluid, mainly
Feature, which is the viscosity of Newtonian fluid, will not cut the variation of variability with fluid and change, that is to say, that the coefficient of viscosity of Newtonian fluid
Be a definite value, such as water, ethyl alcohol, blood plasma, serum etc. all it is certain Newtonian fluid.And the viscosity of non-newtonian fluid then can be with stream
Body is cut the variation of variability and is changed, and typically winged Newtonian fluid includes grease, milk, toothpaste, blood and mud etc..So blood
As a kind of non-newtonian fluid, during simulation, we also need to the influence in view of its viscosity.
Among current analogue technique, there has been no a kind of methods can be good at simulating blood, wherein glutinous for blood
The simulation of degree feature is also seldom.
The basic thought of SPH algorithms is that continuous fluid is regarded as discrete one by one and interaction particle,
These particles influence each other, and together form complicated fluid.For each individual particle, they all follow most base
This Newton's second lawBy studying the stressing conditions of each particle, each particle can be calculated in we
Acceleration, and then speed and displacement is calculated, the movement of all particles together constitutes fluid.
But in simulation process, if time step is larger, the density fluctuation of fluid will be larger, has exceeded permission
Range, cause the at random of particle in simulation process, therefore this method requires time step very little, calculation amount in practical application
Greatly, the live effect of simulation is bad；Current techniques have used single particle to be simulated when simulating fluid, the party
When method is applied to simulation blood, it only will not be able to show the physiology characteristic of blood with a kind of particle, because of blood
Heterogeneity different effects is suffered from organism；If a variety of particles are added by force using the prior art,
The particle of different densities will be layered after the simulation of a period of time, and the larger particle of density can be deposited on bottom, and close
The smaller particle of degree will float.What present day analog technology was simulated is all Newtonian fluid, and viscosity will not be because of shearing rate etc.
Other external factor and change, and when applied to blood simulation, blood belongs to non-newtonian fluid, and viscosity can be because outer
The influence of portion's factor and change, while blood viscosity variation medically have diagnosis the state of an illness important function, very much
Disease can all be accompanied by the variation of blood viscosity, will not be able to accurately describe the viscosity of blood using current method.
Invention content
It is an object of the invention to overcome the deficiencies in the prior art, the blood simulation based on PBF that the present invention provides a kind of
Method can break through the limitation of conventional fluid analogue technique so that the time step with bigger in simulation process has more
Good real-time；In conjunction with a variety of particle hybrid analog-digital simulation blood, can by using particle quantity according to application scenarios it is different into
Row modification appropriate, while the distribution situation of different densities particle is constrained using the density of blood entirety, it ensure that unit bodies
The phenomenon that distribution situation of different type particle in product is not in layering.
In order to solve the above technical problem, the present invention provides a kind of blood analogy method based on PBF, the method packet
It includes：
The information data of particle is obtained, initialization process is carried out；
According to the vessel position at each position of human body, the blood in blood vessel is subjected to sliding-model control, obtain one by one from
Scattered particle；
A variety of particles are subjected to calculation processing, obtain the multidate information of the field particle of each particle；
The position of particle is corrected using PBD frames according to the data information of the field particle of each particle, continuous iteration is straight
To terminating.
Preferably, the information data after initialization process includes：Smooth kernel radius data, density, time step, density are missed
Poor threshold value, maximum iteration, Carson's year coefficient, Casson yield stress.
Preferably, described the step of a variety of particles are carried out calculation processing, includes：
To each particle, calculation processing is carried out, the field particle of particle is obtained；
The density of particle is calculated, the pressure that each particle externally generates is obtained；
According to the viscosity of particle, the pressure being subject to and viscous force, the displacement information of each particle is obtained.
Preferably, the data information of the field particle of each particle of the basis corrects the position of particle using PBD frames
The step of include：
According to the displacement information of particle, update generates the position of particle after displacement；
Updated particle position is obtained, calculation processing is carried out to its density, reacquires the density of the particle；
Integral particles density error is acquired according to all particle densities and preset global density, comparison particle generates displacement
Front and back integral particles density error is judged that it is close to be less than preset maximum for the density error of integral particles after generating displacement
Error threshold is spent, then is ignored；
The density error of integral particles is greater than or equal to preset maximal density error threshold after generating displacement, then carries out
It calculates, corrects the position of particle.
Preferably, by according to the position for constantly correcting particle, observing particle multidate information, entire simulation blood stream is obtained
Dynamic process.
In embodiments of the present invention, the limitation of conventional fluid analogue technique can be broken through so that have in simulation process
The time step of bigger has better real-time；It, can be by using the quantity of particle in conjunction with a variety of particle hybrid analog-digital simulation blood
Modification appropriate is carried out according to the difference of application scenarios, while constraining point of different densities particle using the density of blood entirety
The phenomenon that cloth situation ensure that the distribution situation of different type particle in unit volume, be not in layering.
Description of the drawings
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technology description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
Other attached drawings are obtained according to these attached drawings.
Fig. 1 is a kind of blood analogy method flow diagram based on PBF in the embodiment of the present invention.
Specific implementation mode
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation describes, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, those of ordinary skill in the art are obtained every other without creative efforts
Embodiment shall fall within the protection scope of the present invention.
Fig. 1 is a kind of blood analogy method flow diagram based on PBF of the embodiment of the present invention.As shown in Figure 1, described
Method includes：
S1 obtains the information data of particle, carries out initialization process；
Blood in blood vessel is carried out sliding-model control, obtained one by one by S2 according to the vessel position at each position of human body
Discrete particle；
A variety of particles are carried out calculation processing, obtain the multidate information of the field particle of each particle by S3；
S4 corrects the position of particle, continuous iteration according to the data information of the field particle of each particle using PBD frames
Until terminating.
Specifically, the data after the initialization obtained from S1 include：Smooth kernel radius data, density, time step, density
Error threshold, maximum iteration, Carson's year coefficient, Casson yield stress.
In a particular embodiment, we respectively represent blood plasma and haemocyte using two kinds of particles.
S3 is further described：
S31 carries out calculation processing, obtains the field particle of particle to each particle；
S32 calculates the density of particle, obtains the pressure that each particle externally generates；
S33 obtains the displacement information of each particle according to the viscosity of particle, the pressure being subject to and viscous force.
S4 is described further：
S41, according to the displacement information of particle, update generates the position of particle after displacement；
S42 obtains updated particle position, carries out calculation processing to its density, reacquires the density of the particle；
S43 acquires integral particles density error according to all particle densities and preset global density, and comparison particle generates
Integral particles density error before and after displacement judged, after generating displacement the density error of integral particles be less than it is preset most
Big density error threshold value, then ignore；
The density error of integral particles is greater than or equal to preset maximal density error threshold after generating displacement, then carries out
It calculates, corrects the position of particle.
Further, we obtain entire simulation by according to the position for constantly correcting particle, observing particle multidate information
The process of blood flow.
In embodiments of the present invention, the limitation of conventional fluid analogue technique can be broken through so that have in simulation process
The time step of bigger has better real-time；It, can be by using the quantity of particle in conjunction with a variety of particle hybrid analog-digital simulation blood
Modification appropriate is carried out according to the difference of application scenarios, while constraining point of different densities particle using the density of blood entirety
The phenomenon that cloth situation ensure that the distribution situation of different type particle in unit volume, be not in layering.
One of ordinary skill in the art will appreciate that all or part of step in the various methods of above-described embodiment is can
It is completed with instructing relevant hardware by program, which can be stored in a computer readable storage medium, storage
Medium may include：Read-only memory (ROM, Read Only Memory), random access memory (RAM, Random
Access Memory), disk or CD etc..
In addition, being provided for the embodiments of the invention a kind of blood analogy method based on PBF above has carried out detailed Jie
It continues, principle and implementation of the present invention are described for specific case used herein, and the explanation of above example is only
It is the method and its core concept for being used to help understand the present invention；Meanwhile for those of ordinary skill in the art, according to this hair
Bright thought, there will be changes in the specific implementation manner and application range, in conclusion the content of the present specification should not manage
Solution is limitation of the present invention.
Claims (5)
1. a kind of blood analogy method based on PBF, which is characterized in that the method includes：
The information data of particle is obtained, initialization process is carried out；
According to the vessel position at each position of human body, the blood in blood vessel is subjected to sliding-model control, is obtained discrete one by one
Particle；
A variety of particles are subjected to calculation processing, obtain the multidate information of the field particle of each particle；
The position of particle is corrected using PBD frames according to the data information of the field particle of each particle, continuous iteration is until knot
Beam.
2. a kind of blood analogy method based on PBF as described in claim 1, which is characterized in that the letter after initialization process
Ceasing data includes：Smooth kernel radius data, density, time step, density error threshold value, maximum iteration, Carson's year be
Number, Casson yield stress.
3. a kind of blood analogy method based on PBF as described in claim 1, which is characterized in that it is described by a variety of particles into
The step of row calculation processing includes：
To each particle, calculation processing is carried out, the field particle of particle is obtained；
The density of particle is calculated, the pressure that each particle externally generates is obtained；
According to the viscosity of particle, the pressure being subject to and viscous force, the displacement information of each particle is obtained.
4. a kind of blood analogy method based on PBF as described in claim 1, which is characterized in that each particle of basis
Field particle data information using PBD frames correct particle position the step of include：
According to the displacement information of particle, update generates the position of particle after displacement；
Updated particle position is obtained, calculation processing is carried out to its density, reacquires the density of the particle；
Integral particles density error is acquired according to all particle densities and preset global density, comparison particle generates before and after displacement
Integral particles density error judged, after generating displacement the density error of integral particles less than preset maximal density miss
Poor threshold value, then ignore；
The density error of integral particles is greater than or equal to preset maximal density error threshold after generating displacement, then is counted
It calculates, corrects the position of particle.
5. a kind of blood analogy method based on PBF as claimed in claim 4, which is characterized in that correct particle according to continuous
Position obtain entire simulation blood flow process by observing particle multidate information.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810087356.1A CN108461149A (en) | 2018-01-30 | 2018-01-30 | A kind of blood analogy method based on PBF |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810087356.1A CN108461149A (en) | 2018-01-30 | 2018-01-30 | A kind of blood analogy method based on PBF |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108461149A true CN108461149A (en) | 2018-08-28 |
Family
ID=63239342
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810087356.1A Pending CN108461149A (en) | 2018-01-30 | 2018-01-30 | A kind of blood analogy method based on PBF |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108461149A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109063375A (en) * | 2018-09-07 | 2018-12-21 | 中山大学 | The analogy method and system of incompressible fluid based on secrecy and without divergence |
| CN109271696A (en) * | 2018-09-07 | 2019-01-25 | 中山大学 | Blood clotting analogy method and system based on MPM |
| CN111047707A (en) * | 2019-11-11 | 2020-04-21 | 南昌大学 | Bleeding simulation method of mixed particle blood model based on SPH |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103093081A (en) * | 2012-12-15 | 2013-05-08 | 西南科技大学 | Simulation method of laser-induce biological tissue photothermal effect |
| CN103729555A (en) * | 2013-12-20 | 2014-04-16 | 深圳先进技术研究院 | Method and device for simulating blood flow and vascular wall effects |
| CN104462830A (en) * | 2014-12-12 | 2015-03-25 | 武汉大学 | GPU (graphics processing unit) acceleration based real-time hybrid particle blood flow-blood vessel coupling method |
-
2018
- 2018-01-30 CN CN201810087356.1A patent/CN108461149A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103093081A (en) * | 2012-12-15 | 2013-05-08 | 西南科技大学 | Simulation method of laser-induce biological tissue photothermal effect |
| CN103729555A (en) * | 2013-12-20 | 2014-04-16 | 深圳先进技术研究院 | Method and device for simulating blood flow and vascular wall effects |
| CN104462830A (en) * | 2014-12-12 | 2015-03-25 | 武汉大学 | GPU (graphics processing unit) acceleration based real-time hybrid particle blood flow-blood vessel coupling method |
Non-Patent Citations (2)
| Title |
|---|
| 徐梦佳等: "面向血管疾病诊断及预测分析的血流动力学模拟综述", 《中国图象图形学报》 * |
| 李忠态: "基于SPH和流体动力学的皮肤表面上三维流血模拟", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109063375A (en) * | 2018-09-07 | 2018-12-21 | 中山大学 | The analogy method and system of incompressible fluid based on secrecy and without divergence |
| CN109271696A (en) * | 2018-09-07 | 2019-01-25 | 中山大学 | Blood clotting analogy method and system based on MPM |
| CN109063375B (en) * | 2018-09-07 | 2019-05-31 | 中山大学 | The analogy method and system of incompressible fluid based on secrecy and without divergence |
| CN109271696B (en) * | 2018-09-07 | 2019-07-23 | 中山大学 | Blood clotting analogy method and system based on MPM |
| CN111047707A (en) * | 2019-11-11 | 2020-04-21 | 南昌大学 | Bleeding simulation method of mixed particle blood model based on SPH |
| CN111047707B (en) * | 2019-11-11 | 2021-09-28 | 南昌大学 | Bleeding simulation method of mixed particle blood model based on SPH |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Tozluoǧlu et al. | Planar differential growth rates initiate precise fold positions in complex epithelia | |
| US9251618B2 (en) | Skin and flesh simulation using finite elements, biphasic materials, and rest state retargeting | |
| ES2704080T3 (en) | Biometric simulation device, method to control the biometric simulation device, and program to control the biometric simulation device | |
| CN108461149A (en) | A kind of blood analogy method based on PBF | |
| Schumacher et al. | Efficient simulation of example-based materials | |
| Vella et al. | Effect of the alterations in contractility and morphology produced by atrial fibrillation on the thrombosis potential of the left atrial appendage | |
| Marchandise et al. | Quality open source mesh generation for cardiovascular flow simulations | |
| Walmsley et al. | Combining computer modelling and cardiac imaging to understand right ventricular pump function | |
| CN105022928A (en) | Digitized real-time determination method for center-of-gravity position of fuel system of aircraft | |
| Jagtap et al. | Simulation of buffered advection diffusion of calcium in a hepatocyte cell | |
| CN111047707B (en) | Bleeding simulation method of mixed particle blood model based on SPH | |
| CN103839280B (en) | A kind of human body attitude tracking of view-based access control model information | |
| WO2014205584A1 (en) | Methods and systems for computer-based skin animation | |
| CN106844994B (en) | Choroidal neovascularization growth prediction method combining constitutive model and finite element | |
| EP3179390A1 (en) | Method and apparatus for modeling movement of air bubble based on fluid particles | |
| CN106407932A (en) | Handwritten number recognition method based on fractional calculus and generalized inverse neural network | |
| CN111507304B (en) | Self-adaptive rigidity priori model training method, face tracking method and related devices | |
| CN110459262B (en) | Construction method of human atrial mathematical model | |
| KR101579447B1 (en) | Cellular signal transduction network analyzing method | |
| CN114424934A (en) | Apnea event screening model training method and device and computer equipment | |
| CN109036567B (en) | Soft tissue deformation simulation method based on subspace condensation algorithm | |
| Yan et al. | Soft tissue deformation simulation in virtual surgery using nonlinear finite element method | |
| CN108765432A (en) | Middle membrane boundary automatic division method and system in a kind of arteria carotis | |
| Ye et al. | Research on soft tissue deformation and cutting in the virtual surgery | |
| Hao et al. | MUE-CoT: multi-scale uncertainty entropy-aware co-training framework for left atrial segmentation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180828 |