CN108428256B - Soft tissue deformation simulation method based on self-adaptive grid refinement of softness - Google Patents

Soft tissue deformation simulation method based on self-adaptive grid refinement of softness Download PDF

Info

Publication number
CN108428256B
CN108428256B CN201810288388.8A CN201810288388A CN108428256B CN 108428256 B CN108428256 B CN 108428256B CN 201810288388 A CN201810288388 A CN 201810288388A CN 108428256 B CN108428256 B CN 108428256B
Authority
CN
China
Prior art keywords
soft tissue
softness
grid
area
deformation
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.)
Active
Application number
CN201810288388.8A
Other languages
Chinese (zh)
Other versions
CN108428256A (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.)
Fujian Polytechnic Normal University
Original Assignee
Fujian Normal University
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 Fujian Normal University filed Critical Fujian Normal University
Priority to CN201810288388.8A priority Critical patent/CN108428256B/en
Publication of CN108428256A publication Critical patent/CN108428256A/en
Application granted granted Critical
Publication of CN108428256B publication Critical patent/CN108428256B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Image Generation (AREA)

Abstract

The invention relates to a soft tissue deformation simulation method based on self-adaptive mesh refinement of softness, which comprises the steps of firstly determining the refinement range of meshes in a soft tissue model according to the conduction range of an external force, then improving on the basis of a common self-adaptive refinement rule, formulating a refinement rule based on softness, realizing that the most appropriate length is used as the side length threshold of a triangular mesh, and avoiding the occurrence of a larger triangular mesh in the refinement range. The triangular mesh is encrypted by using a mesh refining method in a soft tissue operation deformation area, so that the sense of reality of deformation simulation is improved, an original triangular mesh model with lower precision is reserved in a non-operation area, the calculation efficiency is improved, and the requirements of the sense of reality and the real-time property of the deformation simulation are met.

Description

Soft tissue deformation simulation method based on self-adaptive grid refinement of softness
Technical Field
The invention relates to the technical field of computer graphics, in particular to a soft tissue deformation simulation method based on self-adaptive grid refinement of softness.
Background
The soft tissue deformation calculation model mainly comprises a spring mass point model, a finite element model and a non-grid model. Lim et al (Lim Y.J, De S. Real time relationship of non-linear tissue response using the point alignment-based method of properties [ J ]. Computer Methods in Applied Mechanics and engineering.2007,196 (31): 3011-3024) adopt a lattice-free method based on distribution points to establish a soft tissue deformation model of a deformation physical model of soft tissue, which has higher calculation efficiency than a finite element method but poorer stability. Li et al (Yan-Dong L I, Zhu L, Xiu-Fen YE, et al, Modeling and Simulation of Soft Tissue Deformation Based on Local Dynamic Model [ J ]. Computer Science, 2013, 40(10): 283-. The application number is CN201110213387.5, the name is soft tissue deformation simulation method based on non-grid Galerkin and mass point spring coupling, the method makes up the defect that the Galerkin method is not suitable for solving large-scale problems, but seamless coupling among multiple models is difficult to realize. Chen Weidong and the like (Chen Weidong, Chen Pan, Zhu Qiguang, research on deformation modeling and force feedback algorithm based on an improved spring-mass point model [ J ]. biomedicine engineering journal, 2015(5):989 + 996.) propose a new variable diamond topological structure model on the basis of the spring-mass point model, realize the deformation simulation of different organs by the model by changing the length, the spring coefficient and the initial included angle of a diamond-shaped edge spring, and compromise is performed on the two aspects of calculation precision and real-time property.
In summary, the simulation of soft tissue deformation plays an important role in virtual surgery research, and the simulation technology thereof is receiving more and more attention. In order to simulate the soft tissue deformation realistically in real time in the virtual surgery, many problems still need to be solved in the aspect of the soft tissue physical deformation model.
In the traditional refining method, a grid refining algorithm is used for carrying out global refining on an original three-dimensional grid model, so that the requirement of improving the simulation precision is met. However, a great disadvantage of the conventional refining method is that the data size is increased in a geometric progression during the grid refining, and the requirement of simulation instantaneity is difficult to meet.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a soft tissue deformation simulation method based on self-adaptive mesh refinement of softness, wherein a mesh refinement method is used for encrypting a triangular mesh in a soft tissue operation deformation area, so that the sense of reality of deformation simulation is improved, and an original triangular mesh model with lower precision is reserved in a non-operation area, so that the calculation efficiency is improved.
In order to achieve the purpose, the technical scheme of the invention is as follows: a soft tissue deformation simulation method based on self-adaptive mesh refinement of softness comprises the following steps:
step S1: detecting a collision point between the soft tissue model and an external force by using a collision detection algorithm, and addressing according to the stress point to use a mass point closest to the collision point as an external force action point;
step S2: calculating the conduction range of the external force, and determining the thinning area range of the deformation grid;
step S3: calculating the softness of the soft tissue model in the thinning area of the deformed grid, and calculating a grid thinning threshold according to the softness;
step S4: determining grids in a thinning area according to a grid thinning threshold, and further thinning a certain grid surface patch in the thinning area if the side length of the surface patch is greater than the grid thinning threshold;
step S5: marking a surface patch which does not need to be refined as a dead surface, marking a mesh which needs to be refined as a live surface, calculating normal vectors and three-dimensional coordinates of a new vertex and a related old vertex by adopting a loop refining method, and connecting the vertices according to a corresponding topological connection rule to finish mesh refinement;
step S6: the heterogeneous grids are transited by utilizing a crack elimination algorithm;
step S7: and repeating the steps S4-S6 until the grids all meet the requirements.
Further, the step S2 is specifically: and presetting a threshold, and when the stress of the external force action point is smaller than the threshold, assuming that the action force can not be conducted to the next mass point any more, and determining the boundary of the refining area by the force conduction range.
Further, the threshold is set to 0.01N.
Further, the step S3 specifically includes:
step S31: according to the acceleration difference between adjacent particles of the soft tissue model and the self curvature of the deformation area, the softness of the soft tissue in the deformation grid refining area is calculated:
adjacent particlesi、jThe acceleration difference between them is:
Figure DEST_PATH_IMAGE001
the curvature of the deformation area is calculated by the normal difference of the triangular patch where each vertex is located, and the normal difference of the triangular patch where adjacent particles i and j are located is as follows:
Figure 476757DEST_PATH_IMAGE002
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE003
,xiand xjRepresenting the coordinates of particle i and particle j respectively,a i anda j representing the acceleration of the particle i and the particle j respectively,arepresenting particle acceleration, n i And n j Respectively, the normal lines of the triangular patches where the particle i and the particle j are located are represented, n represents the particle normal line, and softness S is calculated as follows:
Figure 266334DEST_PATH_IMAGE004
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE005
Figure 921438DEST_PATH_IMAGE006
Figure DEST_PATH_IMAGE007
the maximum acceleration difference of the particle is represented,
Figure 910839DEST_PATH_IMAGE008
represents the maximum extremum of the normal difference;
step S32: calculating the maximum threshold of the side length of the triangular patch of the refined region
Figure DEST_PATH_IMAGE009
Figure 334999DEST_PATH_IMAGE010
Wherein S isiAnd SjRespectively representing the softness of the particle i and the softness of the particle j;
step S33: set the minimum side length of the triangular patch
Figure DEST_PATH_IMAGE011
Further, the air conditioner is provided with a fan,
Figure 767248DEST_PATH_IMAGE012
further, the step S6 is specifically: inserting a new point on the edge of the grid to be refined, and generating a crack on the edge adjacent to the transition grid;
and then connecting the new point on the edge with the opposite vertex of the adjacent transition mesh to divide the transition mesh into two, so that the refined area and the unrefined area are smoothly connected.
Compared with the prior art, the invention has the beneficial effects that:
(1) the mesh refining method is applied to the soft tissue organ model, when the soft tissue is deformed, the mesh is refined in the deformation area, not only can the geometric characteristics of the deformation area be shown, but also the deformation process can be simulated more truly;
(2) the triangular mesh is encrypted in the soft tissue operation deformation area by using a mesh refinement method, so that the reality of deformation simulation is improved. And the original triangular mesh model with lower precision is reserved in the non-operation area, so that the calculation efficiency is improved.
Drawings
FIG. 1 is a schematic flow chart of a soft tissue deformation simulation method based on softness adaptive mesh refinement according to the present invention;
FIG. 2 is a schematic diagram of a particle 1-neighborhood refinement region in an embodiment of the present invention;
FIG. 3(a) is a schematic diagram of an unrefined grid for fracture elimination in an embodiment of the invention;
FIG. 3(b) is a schematic diagram of crack elimination after one refinement in crack elimination in an embodiment of the present invention;
FIG. 3(c) is a schematic diagram of crack elimination after secondary refinement during crack elimination in the embodiment of the invention;
FIG. 4 is a schematic diagram of an experimental model for mesh adaptive refinement according to an embodiment of the present invention;
FIG. 5 is a schematic diagram of an experimental model for global grid refinement according to an embodiment of the present invention;
FIG. 6 is a drawing of the effect of the adaptive mesh primary refinement model on the tensile deformation in the embodiment of the present invention;
FIG. 7 is a drawing of the effect of global refinement model on tensile deformation in an embodiment of the present invention;
FIG. 8 is a drawing of the effect of the adaptive mesh quadratic refinement model on the tensile deformation in the embodiment of the present invention.
Detailed Description
The invention is further explained below with reference to the drawings and the embodiments.
As shown in fig. 1, a soft tissue deformation simulation method based on adaptive mesh refinement of softness includes the following steps:
step S1: detecting a collision point between the soft tissue model and an external force by using a collision detection algorithm, and addressing according to the stress point to use a mass point closest to the collision point as an external force action point;
step S2: calculating the conduction range of the external force, and determining the thinning area range of the deformation grid;
step S3: calculating the softness of the soft tissue model in the thinning area of the deformed grid, and calculating a grid thinning threshold according to the softness;
step S4: determining grids in a thinning area according to a grid thinning threshold, and further thinning a certain grid surface patch in the thinning area if the side length of the surface patch is greater than the grid thinning threshold;
step S5: marking a surface patch which does not need to be refined as a dead surface, marking a mesh which needs to be refined as a live surface, calculating normal vectors and three-dimensional coordinates of a new vertex and a related old vertex by adopting a loop refining method, and connecting the vertices according to a corresponding topological connection rule to finish mesh refinement;
step S6: the heterogeneous grids are transited by utilizing a crack elimination algorithm;
step S7: and repeating the steps S4-S6 until the grids all meet the requirements.
When the soft tissue is operated, no matter the soft tissue is pulled or pressed by an instrument, the contacted surface is a region with concentrated stress and a place with larger deformation, so the stress of mass points is taken as the criterion of local refinement. In this embodiment, the step S2 specifically includes: and presetting a threshold, and when the stress of the external force action point is smaller than the threshold, assuming that the action force can not be conducted to the next mass point any more, and determining the boundary of the refining area by the force conduction range.
The smaller the threshold, the larger the refined region. As shown in FIG. 2, assuming that an external force is applied to the position of mass point 0 in the figure, the displacement of mass point 0 will generate a pulling force on the mass point in the 1-neighborhood, and then may be transmitted to the 2-neighborhood, the 3-neighborhood, etc., and the force is less transmitted to the outside until the force is less than the threshold value, so that the calculation is not performed.
In the present embodiment, the threshold is set to 0.01N.
If the particle motion condition of the soft tissue is clear when the soft tissue is stressed, the most vivid deformation effect is represented by 'softness', and the softness is related to the fineness degree of the grid. Thereby establishing the relationship of the forced movement-softness-fine degree of the grid of the stomach inner wall of the soft tissue.
Second law of Newton
Figure 299861DEST_PATH_IMAGE014
It can be known that the object generates acceleration when being subjected to force movement, and the acceleration is in direct proportion to the force, so the acceleration of the mass point is used to represent the force-receiving movement state of the mass point in the grid. From the analysis, it is known that softness is related to the motion of the grid, and the difference in acceleration between adjacent particles in the grid requires softness.
In this embodiment, the step S3 specifically includes:
step S31: according to the acceleration difference between adjacent particles of the soft tissue model and the self curvature of the deformation area, the softness of the soft tissue in the deformation grid refining area is calculated:
adjacent particlesi、jThe acceleration difference between them is:
Figure 574460DEST_PATH_IMAGE001
except for acceleration difference of particle motion, the curvature of the deformation region needs certain softness, the curvature of the deformation region is calculated by the normal difference of a triangular patch where each vertex is located, and the normal difference of the triangular patches where adjacent particles i and j are located is as follows:
Figure 508918DEST_PATH_IMAGE002
wherein the content of the first and second substances,
Figure 33440DEST_PATH_IMAGE003
,xiand xjRepresenting the coordinates of particle i and particle j respectively,a i anda j representing the acceleration of the particle i and the particle j respectively,arepresenting particle acceleration, n i And n j Respectively, the normal lines of the triangular patches where the particle i and the particle j are located are represented, n represents the particle normal line, and softness S is calculated as follows:
Figure DEST_PATH_IMAGE015
wherein the content of the first and second substances,
Figure 397556DEST_PATH_IMAGE005
Figure 665726DEST_PATH_IMAGE006
Figure 286982DEST_PATH_IMAGE007
the maximum acceleration difference of the particle is represented,
Figure 264296DEST_PATH_IMAGE008
represents the maximum extremum of the normal difference; the maximum value is set to avoid one of themThe factors are infinitely amplified to influence the calculation of the softness of the area so as to better control the mesh refining process;
step S32: the triangular mesh needs to be finer as the softness is larger, and the side length of the triangle is smaller as the mesh is finer. Therefore, the side length threshold of the triangular mesh can be determined according to the softness value calculated above
Figure 568238DEST_PATH_IMAGE009
Calculating the maximum threshold of the side length of the triangular patch in the refined region
Figure 921990DEST_PATH_IMAGE009
Figure 831041DEST_PATH_IMAGE010
Wherein S isiAnd SjRespectively representing the softness of the particle i and the softness of the particle j;
step S33: set the minimum side length of the triangular patch
Figure 228524DEST_PATH_IMAGE016
Only the maximum side length for refining the triangular patch is specified to be insufficient, so that in order to avoid the phenomenon of excessive refinement of the refinement, the triangular patches with some edges close to the threshold and some edges short are generated, a plurality of triangles with small areas and poor quality in the region are caused, the simulation efficiency is influenced, the quality for refining the triangular patch is controlled, and therefore the minimum side length is set
Figure 36074DEST_PATH_IMAGE016
In the present embodiment, setting
Figure DEST_PATH_IMAGE017
Because the balanced adaptive refining strategy may cause excessive crack migration and refining, and the refining levels of adjacent surfaces change after the unbalanced adaptive refining strategy performs different refining for multiple times, the refining levels of the adjacent surfaces are different from each other finally, so that the refining quality is influenced, the invention provides an improved crack elimination algorithm on the basis of the two strategies, and the working schematic diagram of the improved crack elimination algorithm is shown in fig. 3(a), fig. 3(b) and fig. 3 (c). The unrefined meshes that are adjacent to the edge of the mesh being refined are referred to as transition meshes.
In this embodiment, the step S6 specifically includes: inserting a new point on the edge of the grid to be refined, and generating a crack on the edge adjacent to the transition grid; and then connecting the new point on the edge with the opposite vertex of the adjacent transition mesh to divide the transition mesh into two, so that the refined area and the unrefined area are smoothly connected. The same method is also used in the areas with different refinement levels to realize good transition of the grids between different refinement levels, so that the model grid is smoother, and a better visual effect is achieved.
The specific experiment is realized on a computer with Intel (R) core (TM) i7-5500U CPU 2.40GHz and 8G RAM, display card ATI Radeon R9M 375 and display memory 2G by taking OpenGL open source graphic library as a basis, VS2010 as a development platform and C + + as a development language. For better comparison, two models are set in the experiment, one is a soft tissue model with a mass point number of 715 and is used for the mesh adaptive refining experiment, as shown in fig. 4, and the other is a model after the model is globally refined, as shown in fig. 5.
In the experiment, the total mass of the model mass point is set to be 75, the elastic coefficient of the spring is 3, the damping coefficient is 0.5, the time step is 0.01, and the iteration times are 10 times. The mouse was used to control the surgical instrument to stretch the soft tissue model and the globally refined model of 715 particles, the results of which are shown in fig. 6-8.
FIG. 6 is a drawing of the effect of the soft tissue model to realize the once refinement of the adaptive mesh, the stressed particles on the soft tissue are obtained by collision detection and stress point addressing, the force transmission range is calculated by taking the particles as the center, the range is taken as the region for mesh refinement, and the mesh is refined according to the calculated mesh side length threshold, so that the deformation simulation of the soft tissue adaptive mesh refinement is realized, and the deformation effect is more real as can be seen from the drawing. FIG. 7 is a drawing of the effect of the tensile deformation of the global soft tissue refinement model under the action of a large force, and the deformation effect is relatively real. Fig. 8 is an effect diagram of the soft tissue model for realizing the secondary refinement of the tensile deformed grid under the action of a large force, and it can be seen from the diagram that the soft tissue generates large deformation, but the precision of the deformed area grid is improved due to the grid refinement, so that the deformation effect equivalent to the globally refined model is achieved.
The above are preferred embodiments of the present invention, and all changes made according to the technical scheme of the present invention that produce functional effects do not exceed the scope of the technical scheme of the present invention belong to the protection scope of the present invention.

Claims (5)

1. A soft tissue deformation simulation method based on self-adaptive mesh refinement of softness is characterized by comprising the following steps:
step S1: detecting a collision point between the soft tissue model and an external force by using a collision detection algorithm, and addressing according to the stress point to use a mass point closest to the collision point as an external force action point;
step S2: calculating the conduction range of the external force, and determining the thinning area range of the deformation grid;
step S3: calculating the softness of the soft tissue model in the thinning area of the deformed grid, and calculating a grid thinning threshold according to the softness;
step S4: determining grids in a thinning area according to a grid thinning threshold, and further thinning a certain grid surface patch in the thinning area if the side length of the surface patch is greater than the grid thinning threshold;
step S5: marking a surface patch which does not need to be refined as a dead surface, marking a mesh which needs to be refined as a live surface, calculating normal vectors and three-dimensional coordinates of a new vertex and a related old vertex by adopting a loop refining method, and connecting the vertices according to a corresponding topological connection rule to finish mesh refinement;
step S6: the heterogeneous grids are transited by utilizing a crack elimination algorithm;
step S7: repeating the steps S4-S6 until the grids all meet the requirements;
wherein, the step S3 specifically includes:
step S31: according to the acceleration difference between adjacent particles of the soft tissue model and the self curvature of the deformation area, the softness of the soft tissue in the deformation grid refining area is calculated:
the acceleration difference between adjacent particles i, j is:
Figure FDA0003255267640000011
the curvature of the deformation area is calculated by the normal difference of the triangular patch where each vertex is located, and the normal difference of the triangular patch where adjacent particles i and j are located is as follows:
Figure FDA0003255267640000012
wherein x isij=xi-xj,xiAnd xjCoordinates of a and j, respectively, are showniAnd ajThe acceleration of the mass point i and the mass point j are respectively shown, a is the acceleration of the mass point, and n isiAnd njRespectively, the normal lines of the triangular patches where the particle i and the particle j are located are represented, n represents the particle normal line, and softness S is calculated as follows:
Figure FDA0003255267640000013
wherein the content of the first and second substances,
Figure FDA0003255267640000021
Δamaxrepresents the maximum acceleration difference, Δ n, of the particlemaxRepresents the maximum extremum of the normal difference;
step S32: calculating the maximum side length threshold l of the triangular patch of the refined regionmax
Figure FDA0003255267640000022
Wherein S isiAnd SjRespectively representing the softness of the particle i and the softness of the particle j;
step S33: set the minimum side length l of the triangular patchmin
2. The soft tissue deformation simulation method according to claim 1, wherein the step S2 is specifically: and presetting a threshold, and when the stress of the external force action point is smaller than the threshold, assuming that the action force can not be conducted to the next mass point any more, and determining the boundary of the refining area by the force conduction range.
3. A soft tissue deformation simulation method according to claim 2, wherein the threshold value is set to 0.01N.
4. A soft tissue deformation simulation method according to claim 1, characterized in that/, ismin=0.3lmax
5. The soft tissue deformation simulation method according to claim 1, wherein the step S6 is specifically: inserting a new point on the edge of the grid to be refined, and generating a crack on the edge adjacent to the transition grid; and then connecting the new point on the edge with the opposite vertex of the adjacent transition mesh to divide the transition mesh into two, so that the refined area and the unrefined area are smoothly connected.
CN201810288388.8A 2018-04-03 2018-04-03 Soft tissue deformation simulation method based on self-adaptive grid refinement of softness Active CN108428256B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810288388.8A CN108428256B (en) 2018-04-03 2018-04-03 Soft tissue deformation simulation method based on self-adaptive grid refinement of softness

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810288388.8A CN108428256B (en) 2018-04-03 2018-04-03 Soft tissue deformation simulation method based on self-adaptive grid refinement of softness

Publications (2)

Publication Number Publication Date
CN108428256A CN108428256A (en) 2018-08-21
CN108428256B true CN108428256B (en) 2021-12-21

Family

ID=63160225

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810288388.8A Active CN108428256B (en) 2018-04-03 2018-04-03 Soft tissue deformation simulation method based on self-adaptive grid refinement of softness

Country Status (1)

Country Link
CN (1) CN108428256B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109992175B (en) * 2019-04-03 2021-10-26 腾讯科技(深圳)有限公司 Object display method, device and storage medium for simulating blind feeling
CN110096818B (en) * 2019-05-06 2023-07-25 南京信息工程大学 Soft tissue supervision deformation algorithm
CN113470165B (en) * 2021-06-17 2022-09-16 南昌大学 Soft tissue modeling method based on radial basis point interpolation method and mass point spring method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6760022B1 (en) * 2001-01-11 2004-07-06 Autodesk, Inc. Simple soft creases for subdivision surfaces
CN105513130A (en) * 2016-02-01 2016-04-20 福建师范大学福清分校 Soft tissue deformation method based on mixing of gridding method and non-gridding method
CN107256557A (en) * 2017-05-03 2017-10-17 华南理工大学 A kind of controllable subdivision curved surface image vector method of error

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6760022B1 (en) * 2001-01-11 2004-07-06 Autodesk, Inc. Simple soft creases for subdivision surfaces
CN105513130A (en) * 2016-02-01 2016-04-20 福建师范大学福清分校 Soft tissue deformation method based on mixing of gridding method and non-gridding method
CN107256557A (en) * 2017-05-03 2017-10-17 华南理工大学 A kind of controllable subdivision curved surface image vector method of error

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Cellular neural network modelling of soft tissue dynamics for surgical simulation;Jinao Zhang 等;《Technology and Health Care:official journal of the European Society for Engineering and Medicine》;20170531;第25卷(第S1期);第337-344页 *
基于改进的质点-弹簧模型的软组织形变仿真研究;王舒珍 等;《通化师范学院学报》;20170220(第1期);第64-69页 *
软组织形变模型的自适应网格细分算法;刘秀玲 等;《河北大学学报(自然科学版)》;20130525;第33卷(第3期);第312-316、323页 *

Also Published As

Publication number Publication date
CN108428256A (en) 2018-08-21

Similar Documents

Publication Publication Date Title
CN108428256B (en) Soft tissue deformation simulation method based on self-adaptive grid refinement of softness
Teran et al. Adaptive physics based tetrahedral mesh generation using level sets
CN106875462B (en) Real-time digital organ cutting method based on meta-sphere model and hybrid driving method
US20090135189A1 (en) Character animation system and method
CN106504329B (en) Gum deformation simulation method based on mass point spring model of tooth long axis
Guo et al. Meshless thin-shell simulation based on global conformal parameterization
CN107330972B (en) Real-time soft tissue deformation method and system for simulating biomechanical characteristics
Paloc et al. Online multiresolution volumetric mass spring model for real time soft tissue deformation
CN111488670B (en) Nonlinear mass point spring soft tissue deformation simulation method
Chen et al. Peridynamics‐Based Fracture Animation for Elastoplastic Solids
Zhong et al. A physically based method for triangulated surface flattening
Bender et al. Adaptive cloth simulation using corotational finite elements
CN109002571B (en) Cloth dynamic simulation method based on isogeometric spring particle model
CN109961514B (en) Cutting deformation simulation method and device, storage medium and terminal equipment
Mandal et al. Interactive physics-based virtual sculpting with haptic feedback
CN109829232B (en) Layered material distribution simulation method based on random forest algorithm
Wang et al. A novel virtual cutting method for deformable objects using high‐order elements combined with mesh optimisation
Zou et al. A new deformation simulation algorithm for elastic-plastic objects based on splat primitives
Shi et al. A mixed-depth visual rendering method for bleeding simulation
CN111047684A (en) Model simplification method based on three-dimensional model characteristics
CN113191018B (en) Cloth simulation method and system, electronic equipment and storage medium
Qi et al. Softness-based adaptive mesh refinement algorithm for soft tissue deformation
CN104156502B (en) A kind of location-based clothing fold geometry generation method
Völlinger et al. Evolutionary optimization of mass-spring models
Wang et al. An unfixed-elasticity mass spring model based simulation for soft tissue deformation

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
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20221128

Address after: No. 1, Campus New Village, Longrong Town, Fuqing City, Fuzhou City, Fujian Province 350300

Patentee after: Fujian Normal University of Technology

Address before: 350300 No. 1, campus new village, Longjiang street, Fuqing City, Fuzhou City, Fujian Province

Patentee before: FUQING BRANCH OF FUJIAN NORMAL University

TR01 Transfer of patent right
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20180821

Assignee: FUJIAN SHUBOXUN INFORMATION TECHNOLOGY CO.,LTD.

Assignor: Fujian Normal University of Technology

Contract record no.: X2023980036349

Denomination of invention: A Soft Tissue Deformation Simulation Method Based on Flexibility and Adaptive Mesh Refinement

Granted publication date: 20211221

License type: Common License

Record date: 20230613

EE01 Entry into force of recordation of patent licensing contract