CN114078355B - Method and apparatus for simulating anal sphincter complex stress - Google Patents
Method and apparatus for simulating anal sphincter complex stress Download PDFInfo
- Publication number
- CN114078355B CN114078355B CN202111388883.4A CN202111388883A CN114078355B CN 114078355 B CN114078355 B CN 114078355B CN 202111388883 A CN202111388883 A CN 202111388883A CN 114078355 B CN114078355 B CN 114078355B
- Authority
- CN
- China
- Prior art keywords
- anal sphincter
- sphincter complex
- complex
- rigid body
- finite element
- 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
Links
- 210000005070 sphincter Anatomy 0.000 title claims abstract description 194
- 238000000034 method Methods 0.000 title claims abstract description 24
- 210000003205 muscle Anatomy 0.000 claims abstract description 29
- 238000004088 simulation Methods 0.000 claims description 58
- 238000004458 analytical method Methods 0.000 claims description 33
- 238000012216 screening Methods 0.000 claims description 10
- 239000007787 solid Substances 0.000 claims description 8
- 210000004197 pelvis Anatomy 0.000 claims description 7
- 210000001519 tissue Anatomy 0.000 claims description 5
- 238000007920 subcutaneous administration Methods 0.000 claims description 4
- 230000000452 restraining effect Effects 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 2
- 230000004064 dysfunction Effects 0.000 abstract description 8
- 238000002560 therapeutic procedure Methods 0.000 abstract description 8
- 206010021639 Incontinence Diseases 0.000 abstract description 5
- 210000003903 pelvic floor Anatomy 0.000 abstract description 5
- 230000006378 damage Effects 0.000 abstract description 4
- 238000003745 diagnosis Methods 0.000 abstract description 4
- 201000010099 disease Diseases 0.000 abstract description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 abstract description 4
- 239000003814 drug Substances 0.000 abstract description 4
- 230000007246 mechanism Effects 0.000 abstract description 4
- 230000008506 pathogenesis Effects 0.000 abstract description 4
- 208000004680 Rectal Fistula Diseases 0.000 description 10
- 206010002156 anal fistula Diseases 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 239000000463 material Substances 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 210000004872 soft tissue Anatomy 0.000 description 4
- 230000003993 interaction Effects 0.000 description 3
- 230000002452 interceptive effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000001356 surgical procedure Methods 0.000 description 3
- 206010055023 Anal injury Diseases 0.000 description 2
- 210000000436 anus Anatomy 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000002203 pretreatment Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 125000000174 L-prolyl group Chemical group [H]N1C([H])([H])C([H])([H])C([H])([H])[C@@]1([H])C(*)=O 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 210000002640 perineum Anatomy 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000284 resting effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
- G09B—EDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
- G09B23/00—Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes
- G09B23/28—Models for scientific, medical, or mathematical purposes, e.g. full-sized devices for demonstration purposes for medicine
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Medicinal Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Algebra (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Medical Informatics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Educational Technology (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The disclosed embodiments disclose methods and devices for simulating anal sphincter complex stress. According to the embodiment of the disclosure, the stress of the anal sphincter complex is simulated by combining non-invasive examination with a finite element model, so that the mechanical characteristics of the anal sphincter complex and the pathogenesis of anal incontinence caused by damage of the anal sphincter complex can be more intuitively known from the mechanical angle, a biomechanical basis is provided for the occurrence mechanism of pelvic floor dysfunction diseases, particularly the bowel control dysfunction caused by perianal muscle tissue, and theoretical support is provided for the establishment of anal wilt diagnosis and treatment schemes such as traditional Chinese medicine string therapy and the implementation of quantitative operations.
Description
Technical Field
The present disclosure relates to the field of computer technology, and in particular, to a method and apparatus for simulating anal sphincter complex stress.
Background
Anal incontinence is the most serious sequelae after anal fistula surgery, and in particular, anal sphincter muscle is destroyed to various degrees by anal fistula surgery with high complexity, so that the treatment of the anal sphincter muscle is critical during the surgery. At present, for high anal fistula crossing the anal sphincter, a clinician mainly adopts a hanging wire therapy according to imaging data and self experience during operation. The anal sphincter injury degree is accurately evaluated clinically, the advantages and disadvantages of anal fistula incision and string hanging operation are weighed, string hanging therapy is properly applied, the string hanging position, the thickness, the tightness and the string falling time of the muscle string are further standardized, the string hanging is not only blindly, but also the function of protecting anus and the comfort of a patient are considered, the anal fistula treatment device has important significance for treating anal fistula, and intensive research is urgently needed by adopting modern technology.
Disclosure of Invention
To address the problems in the related art, embodiments of the present disclosure provide methods and apparatus for simulating anal sphincter complex stress.
In a first aspect, embodiments of the present disclosure provide a method of simulating anal sphincter complex stress comprising:
constructing an anal sphincter complex model based on pelvic CT data and anal sphincter complex muscle group MR data, wherein the anal sphincter complex model comprises an anal sphincter complex circular hole;
acquiring the anal sphincter complex model in a finite element simulation system;
creating an analytical rigid body in the finite element simulation system, the analytical rigid body comprising a cylindrical portion and a hemispherical portion connected to the cylindrical portion;
and enabling the analysis rigid body to pass through the circular hole of the anal sphincter complex in the finite element simulation system to obtain the stress of the anal sphincter complex.
With reference to the first aspect, in a first implementation manner of the first aspect, the diameter of the cylindrical portion and the hemispherical portion of the analytical rigid body is 13mm-15mm.
With reference to the first aspect, in a second implementation manner of the first aspect, the elastic modulus of the analytical rigid body is 0.95KPa, the poisson ratio is 0.45, and the density is 1.12e -9 Kg/mm 3 。
With reference to any one of the first aspect, the first implementation manner of the first aspect, and the second implementation manner of the first aspect, in a third implementation manner of the first aspect, before passing the analytical rigid body through the anal sphincter complex circular hole in the finite element simulation system, the method further includes:
converting the anal sphincter complex model and the analytical rigid body into an instance in the finite element simulation system;
in the finite element simulation system, the anal sphincter complex model and the analysis rigid body are positioned at adjacent positions, the analysis rigid body is a contact main surface, the anal sphincter complex is negative in contact, the contact attribute is friction contact, and the contact type is surface-to-surface contact;
boundary conditions are defined in the finite element simulation system, constraining six degrees of freedom on both sides of the anal sphincter complex.
With reference to any one of the first aspect, the first implementation manner of the first aspect, and the second implementation manner of the first aspect, in a fourth implementation manner of the first aspect, before passing the analytical rigid body through the anal sphincter complex circular hole in the finite element simulation system, the method further includes:
and meshing the anal sphincter complex, wherein the mesh attribute is a tetrahedral freeform mesh, and the cell type of the mesh is an explicit nonlinear cell type.
With reference to the first aspect, in a fifth implementation manner of the first aspect, the causing the analytical rigid body to pass through the anal sphincter complex circular hole in the finite element simulation system, and obtaining stress of the anal sphincter complex includes:
rotating the analysis rigid body to enable the axis direction of the analysis rigid body to be perpendicular to the circular hole of the anal sphincter complex;
the analytical rigid body is moved down through the anal sphincter complex circular aperture along the Z-axis.
With reference to the first aspect, in a sixth implementation manner of the first aspect, the constructing an anal sphincter complex model based on pelvis CT data and anal sphincter complex muscle group MR data includes:
the pelvis CT data and the anal sphincter complex muscle group MR data are subjected to data reading, image positioning, image screening and region growing treatment to obtain an anal sphincter complex global geometric model;
converting the global geometric model of the anal sphincter complex into a solid, and selecting muscle tissue from the solid from the subcutaneous part of the external sphincter at the anal margin to the lower margin of the anal sphincter complex to obtain the model of the anal sphincter complex.
With reference to the sixth implementation manner of the first aspect, in a seventh implementation manner of the first aspect, the image screening includes screening the outlined anal sphincter complex tissue from a complex image; the region growing serves to segment out the anal sphincter complex musculature.
With reference to the sixth implementation manner of the first aspect, in an eighth implementation manner of the first aspect, the disclosure further includes performing a three-dimensional finite element pre-processing on the global geometric model of the anal sphincter complex before the converting the global geometric model into the entity.
In a second aspect, embodiments of the present disclosure provide an apparatus for simulating anal sphincter complex stress, comprising:
a model building unit configured to build an anal sphincter complex model based on pelvic CT data and anal sphincter complex muscle group MR data, wherein the anal sphincter complex model comprises an anal sphincter complex circular hole;
a stress simulation unit configured to acquire the anal sphincter complex model data in a finite element simulation system; creating an analytical rigid body in the finite element simulation system, the analytical rigid body comprising a cylindrical portion and a hemispherical portion connected to the cylindrical portion; and enabling the analytic rigid body to pass through the circular hole of the anal sphincter complex in the finite element simulation system to obtain the stress of the anal sphincter complex.
With reference to the second aspect, in a first implementation manner of the second aspect, the diameter of the cylindrical portion and the hemispherical portion of the analytical rigid body is 13mm-15mm.
Combine with the firstIn a second implementation manner of the second aspect, the elastic modulus of the analytical rigid body is 0.95KPa, the poisson ratio is 0.45, and the density is 1.12e -9 Kg/mm 3 。
With reference to any one of the second aspect, the first implementation manner of the second aspect, and the second implementation manner of the second aspect, in a third implementation manner of the second aspect, before passing the analytical rigid body through the anal sphincter complex circular hole in the finite element simulation system, the method further includes:
converting the anal sphincter complex model and the analytical rigid body into an instance in the finite element simulation system;
in the finite element simulation system, the anal sphincter complex model and the analysis rigid body are positioned at adjacent positions, the analysis rigid body is a contact main surface, the anal sphincter complex is negative in contact, the contact attribute is friction contact, and the contact type is surface-to-surface contact;
boundary conditions are defined in the finite element simulation system, constraining six degrees of freedom on both sides of the anal sphincter complex.
With reference to any one of the second aspect, the first implementation manner of the second aspect, and the second implementation manner of the second aspect, in a fourth implementation manner of the second aspect, before passing the analytical rigid body through the anal sphincter complex circular hole in the finite element simulation system, the method further includes:
and meshing the anal sphincter complex, wherein the mesh attribute is a tetrahedral freeform mesh, and the cell type of the mesh is an explicit nonlinear cell type.
With reference to the second aspect, in a fifth implementation manner of the second aspect, the causing the analytical rigid body to pass through the anal sphincter complex circular hole in the finite element simulation system, obtaining stress of the anal sphincter complex includes:
rotating the analysis rigid body to enable the axis direction of the analysis rigid body to be perpendicular to the circular hole of the anal sphincter complex;
the analytical rigid body is moved down through the anal sphincter complex circular aperture along the Z-axis.
With reference to the second aspect, in a sixth implementation manner of the second aspect, the constructing an anal sphincter complex model based on pelvic CT data and anal sphincter complex muscle group MR data includes:
the pelvis CT data and the anal sphincter complex muscle group MR data are subjected to data reading, image positioning, image screening and region growing treatment to obtain an anal sphincter complex global geometric model;
converting the global geometric model of the anal sphincter complex into a solid, and selecting muscle tissue from the solid from the subcutaneous part of the external sphincter at the anal margin to the lower margin of the anal sphincter complex to obtain the model of the anal sphincter complex.
With reference to the sixth implementation manner of the second aspect, in a seventh implementation manner of the second aspect, the image screening includes screening the outlined anal sphincter complex tissue from a complex image; the region growing serves to segment out the anal sphincter complex musculature.
With reference to the sixth implementation manner of the second aspect, in an eighth implementation manner of the second aspect, the present disclosure further includes performing three-dimensional finite element pretreatment on the global geometric model of the anal sphincter complex before the converting the global geometric model into the entity.
The technical scheme provided by the embodiment of the disclosure can comprise the following beneficial effects:
according to the technical scheme provided by the embodiment of the disclosure, the stress of the anal sphincter complex is simulated by combining non-invasive examination with the finite element model, so that the mechanical characteristics of the anal sphincter complex and the pathogenesis of anal incontinence caused by damage of the anal sphincter complex can be more intuitively known from the mechanical angle, a biomechanical basis is provided for pelvic floor dysfunction diseases, particularly the occurrence mechanism of bowel control dysfunction caused by perianal muscle tissue, and theoretical support is provided for the establishment of anal wilt diagnosis and treatment schemes such as traditional Chinese medicine on-line therapy and the implementation of quantitative operations.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the disclosure.
Drawings
Other features, objects and advantages of the present disclosure will become more apparent from the following detailed description of non-limiting embodiments, taken in conjunction with the accompanying drawings. In the drawings:
fig. 1 illustrates a flow chart of simulating anal sphincter complex stress in accordance with an embodiment of the present disclosure;
FIG. 2A illustrates a scan parameter map of a CT scan in accordance with an embodiment of the present disclosure;
FIG. 2B illustrates a scan parameter map of an MR scan in accordance with an embodiment of the present disclosure;
FIG. 3 illustrates a flowchart of constructing an anal sphincter complex model based on pelvic CT data and anal sphincter complex muscle group MR data in accordance with an embodiment of the present disclosure;
fig. 4A shows a schematic diagram of an anal sphincter complex global geometric model after three-dimensional finite element pre-treatment in accordance with an embodiment of the present disclosure;
FIG. 4B illustrates a close-up view of the global geometric model shown in FIG. 4A;
fig. 5 shows a schematic representation of an anal sphincter complex model obtained according to an embodiment of the present disclosure;
FIG. 6 illustrates a schematic diagram of a resolved rigid body created in accordance with an embodiment of the present disclosure;
FIG. 7 illustrates a flow chart for configuring an anal sphincter complex model and analytical rigid bodies in a finite element simulation system according to an embodiment of the present disclosure;
FIG. 8 illustrates a schematic model diagram after setting up an anal sphincter complex model and analyzing contact properties of a rigid body in a finite element simulation system according to an embodiment of the present disclosure;
FIG. 9 illustrates a schematic diagram of a model with two sides constrained to an anal sphincter complex model in a finite element simulation system according to an embodiment of the present disclosure;
fig. 10 shows a schematic diagram of stress curves and stress values of an anal sphincter complex simulated according to an embodiment of the present disclosure;
fig. 11 shows a block diagram of a simulated anal sphincter complex stress device in accordance with an embodiment of the present disclosure.
Detailed Description
Hereinafter, exemplary embodiments of the present disclosure will be described in detail with reference to the accompanying drawings so that those skilled in the art can easily implement them. In addition, for the sake of clarity, portions irrelevant to description of the exemplary embodiments are omitted in the drawings.
In this disclosure, it should be understood that terms such as "comprises" or "comprising," etc., are intended to indicate the presence of features, numbers, steps, acts, components, portions, or combinations thereof disclosed in this specification, and are not intended to exclude the possibility that one or more other features, numbers, steps, acts, components, portions, or combinations thereof are present or added.
In addition, it should be noted that, without conflict, the embodiments of the present disclosure and features of the embodiments may be combined with each other. The present disclosure will be described in detail below with reference to the accompanying drawings in conjunction with embodiments.
As mentioned above, in the present clinic, for the high anal fistula crossing the anal sphincter, the clinician adopts the hanging wire therapy based on the imaging data and self experience during the operation. The anal sphincter injury degree is accurately evaluated clinically, the advantages and disadvantages of anal fistula incision and string hanging operation are weighed, string hanging therapy is properly applied, the string hanging position, the thickness, the tightness and the string falling time of the muscle string are further standardized, the string hanging is not only blindly, but also the function of protecting anus and the comfort of a patient are considered, the anal fistula treatment device has important significance for treating anal fistula, and intensive research is urgently needed by adopting modern technology.
In view of the above, the present disclosure presents a method of simulating anal sphincter complex stress.
Fig. 1 shows a flowchart of simulating anal sphincter complex stress according to an embodiment of the present disclosure, and as shown in fig. 1, the method of simulating anal sphincter complex stress includes steps S110-S140:
in step S110, constructing an anal sphincter complex model comprising an anal sphincter complex circular aperture based on the pelvic CT data and the anal sphincter complex muscle group MR data;
in step S120, the anal sphincter complex model is acquired in a finite element simulation system;
in step S130, creating an analytical rigid body in the finite element simulation system;
in step S140, the analytical rigid body is passed through the anal sphincter complex circular hole in the finite element simulation system to obtain the stress of the anal sphincter complex.
Through the steps of the embodiment of the disclosure, the pathogenesis of anal incontinence caused by the mechanical property of the anal sphincter complex and the damage thereof can be more intuitively known from the mechanical angle, a biomechanical basis is provided for the occurrence mechanism of pelvic floor dysfunction diseases, particularly the bowel control dysfunction caused by perianal musculature, and theoretical support is provided for the establishment of anal wilt diagnosis and treatment schemes such as traditional Chinese medicine thread hanging therapy and the implementation of quantitative operations.
In one embodiment of the present disclosure, CT is used to scan pelvic bony structures and MR scans anal sphincter complex muscle groups because CT has a higher resolution of bony tissue and MR has a higher resolution of soft tissue. The CT and MR scanning range can be up to the level of sacrum 2 and down to 1.0cm below perineum, DICOM-format original image data of soft tissues CT and MR of pelvis and anal sphincter complex muscle groups of a subject are respectively obtained, and a CD.R optical disc is used for recording the disc. The CT can be scanned in a conventional supine position using U.S. GE MEDICAL SYSTEMS/Discovery CT 750HD, with scanning parameters shown in FIG. 2A. MR can use GE MEDICAL SYSTEMS/Discovery MR 750, wherein a thin layer scan is used for three-dimensional reconstruction, a conventional scan is used as a reference for anatomical recognition, a conventional supine position is taken, an abdomen coil is used, a subject is in a resting state, and images of the axial position, the coronal position and the sagittal position T2WI of pelvic floor soft tissue are acquired, a 3DMR sequence is adopted for MR scan, and scanning parameters are shown in FIG. 2B.
In one embodiment of the present disclosure, as shown in fig. 3, the constructing the anal sphincter complex model based on the pelvic CT data and the anal sphincter complex muscle group MR data further comprises steps S1101-S1103:
in step S1101, the CT data and MR data are processed based on the interactive medical image control system, so as to obtain a global geometric model of the anal sphincter complex. Specifically, the global geometric model of the anal sphincter complex can be obtained by performing data reading, image positioning, image screening, region growing and other processes on the pelvis CT data and the anal sphincter complex muscle group MR data.
In one embodiment of the present disclosure, first, the obtained raw image data of the subject's pelvis, anal sphincter complex muscle group soft tissue CT and MR is marked, resulting in an imagewise tomogram of the marked anal sphincter complex contour, which may be saved in a "tif" format. And then reading and processing the marked screenshot through an interactive medical image control system to construct an anal sphincter complex model. The interactive medical image control system can be Mimics software, the Mimics software reads in the marked screenshot in sequence, ultrasonic original scanning data are automatically set during reading, the interlayer spacing can be set to be 0.8mm, the sketched anal sphincter complex tissues are screened out from the complex images through positioning images based on the threading function of the Mimics software, the anal sphincter complex musculature is segmented through region growing, and the overall geometric model of the anal sphincter complex, namely the anal sphincter complex, is obtained by combining with pelvic data.
In step S1102, a three-dimensional finite element pre-treatment is performed on the anal sphincter complex global geometric model.
In step S1103, the global geometric model of the anal sphincter complex is converted into an entity, and the anal sphincter complex model is obtained according to the entity.
In one embodiment of the present disclosure, the global geometric model of the anal sphincter complex is first subjected to a three-dimensional finite element pre-process to generate a finite element mesh model. The finite element grid generation is a process of dispersing an object in a working environment into simple units, wherein a tetrahedral unit is adopted for dividing a three-dimensional space, a quadrilateral unit is adopted for dividing a two-dimensional unit, and for units with curved boundaries and curved surfaces, finite element analysis requires a plurality of points on each side or surface. The global geometric model of the anal sphincter complex after three-dimensional finite element pretreatment is shown in fig. 4A, and fig. 4B is a partial enlarged view of fig. 4A.
In one embodiment of the present disclosure, the global geometric model of the anal sphincter complex may be converted by digital model software into a solid body, and then the anal sphincter complex model as shown in fig. 5 may be obtained by selecting from the solid body the muscle tissue from the subcutaneous part of the external sphincter at the anal vernal site to the lower rim of the anal sphincter complex. As shown in fig. 5, the model of the anal sphincter complex obtained by the previous step includes the anal sphincter complex 1 and the anal sphincter complex circular hole 2, and the diameter of the circular hole 2 can be measured by parameterized three-dimensional software. Wherein the digital model software may be Geomagic Studio and the parametric three-dimensional software may be Pro/ENGINEER.
In one embodiment of the present disclosure, the obtaining the anal sphincter complex model in the finite element simulation system may be converting the anal sphincter complex model obtained in the previous step into a suitable format and importing the same into the finite element simulation system. Wherein the suitable format may be an igs format and the finite element simulation system may be ABAQUS simulated finite element software.
In one embodiment of the present disclosure, the creating an analytical rigid body in the finite element simulation system may be creating an analytical rigid body part in the finite element simulation system, defining the size and material properties of the analytical rigid body. Taking the example of creating an analytical rigid body in ABAQUS simulated finite element software, an analytical rigid body Part may first be created by a "Create Part" module and the dimensions of the analytical rigid body defined. In one example, the parsing rigid body as shown in fig. 6 may include a cylindrical portion, and a hemispherical portion connected to the cylindrical portion, and the cylindrical portion and hemispherical portion of the parsing rigid body may have diameters of 12mm-15mm, preferably 12mm; then through the Create Material in the Property ModuleThe sub-module imparts material properties to the analytical rigid body, wherein the material properties can be elastic modulus of 0.95KPa, poisson's ratio of 0.45 and density of 1.12e -9 Kg/mm 3 Meanwhile, the stress-strain relation of the anal sphincter complex obtained through experiments can be input into the material mechanical property of the sub-module; finally, a cross-Section is created by the "Create Section" module and the cross-Section properties are imparted to the anal sphincter complex component by the "Assign Section" module.
In one embodiment of the present disclosure, as shown in fig. 7, before the step of passing the analytical rigid body through the anal sphincter complex circular hole in the finite element simulation system, the method further includes steps S210-S240:
in step S210, the anal sphincter complex model and the analytical rigid body are transformed into an instance in the finite element simulation system.
In one embodiment of the present disclosure, the anal sphincter complex model and the analytical rigid body may be converted into an Instance by the "Instance Part" sub-module in the "Assembly" module in ABAQUS simulated finite element software.
In step S220, the anal sphincter complex model and the analytical rigid body are brought into close proximity in the finite element simulation system.
In one embodiment of the present disclosure, a Dynamic Explicit analysis Step, i.e., a "Dynamic, explicit" analysis Step, may be created by first simulating the "Create Step" sub-module in the "Step" module in the finite element software by ABAQUS, where the parameter "Time period" may be set to 1; then creating a contact attribute as friction contact through a Create Interaction Property submodule of an Interaction module so as to enable the anal sphincter complex model and the analysis rigid body to be in a close position; finally, setting the analysis rigid body as a contact main Surface through a 'Create Interaction' submodule, wherein the anal sphincter complex is negative in contact, and the contact type is Surface-to-surfacecontact (Explicit) Surface contact, and the analysis step selects the dynamic display analysis step created in the previous step. The model after setting the contact property is shown in fig. 8.
In step S230, boundary conditions are defined in the finite element simulation system.
In one embodiment of the present disclosure, boundary conditions may be defined by a "Create Boundary Condition" sub-module in the "Load" module in ABAQUS simulated finite element software. In one example, the analysis steps in the ABAQUS simulation finite element software are first selected to "Initial" to obtain a model with two sides constrained as described in FIG. 9; and then selecting and analyzing a Reference Point (RP) on the rigid body, and restraining all six degrees of freedom in an analysis step (Initial), so as to realize all restraint of six degrees of freedom on two sides of the anal sphincter complex. The analysis step may then select "Initial".
In step S240, the anal sphincter complex is network-classified.
In one embodiment of the present disclosure, the anal sphincter complex may be network partitioned by a "Mesh" module in ABAQUS simulation finite element software, where the parameter "Approximate global size" in the "Seed Part" sub-module may be set to 5, with other properties unchanged; in the grid attribute definition module Assign Mesh Controls, the grid attribute can be set as a tetrahedron Tet free grid, and other attributes are unchanged; the unit type definition module "Assign Element Type" can selectively display nonlinear unit types, namely, "Element Library" is set as "Explicit", "geometry Order" is set as "numeric", and "Family" is set as "3D Stress".
In one embodiment of the present disclosure, the step S140 further includes: rotating the analysis rigid body to enable the axis direction of the analysis rigid body to be perpendicular to the circular hole of the anal sphincter complex; and enabling the analysis rigid body to move downwards along the Z axis to pass through the circular hole of the anal sphincter complex, so as to obtain the stress of the anal sphincter complex.
Taking the example of simulating the stress of the anal sphincter complex in ABAQUS simulation finite element software, firstly, in the step S210, after the anal sphincter complex model and the analysis rigid body are converted into examples through an "Instance Part" sub-module in an "Assembly" module, the analysis rigid body is rotated so that the axial direction of the analysis rigid body is perpendicular to the circular hole of the anal sphincter complex. Then, after all six degrees of freedom are constrained in the analysis Step "Initial" in the aforementioned Step S230, the parameter "U3" in the "Step-1" module is set to a preset value, so that the analytical rigid body moves downward along the Z axis by the preset value to pass through the circular hole. Wherein, the preset value may be 20mm. Finally, starting the stress simulation operation of the anal sphincter complex through a 'Create Job' submodule in a 'Job' module, and creating a node Path1 at the round hole part after the simulation operation is completed, so as to obtain a stress curve formed by connecting stress points on the output Path1 as shown in fig. 10.
According to the technical scheme provided by the embodiment of the disclosure, the stress of the anal sphincter complex is simulated by combining non-invasive examination with the finite element model, so that the mechanical characteristics of the anal sphincter complex and the pathogenesis of anal incontinence caused by damage of the anal sphincter complex can be more intuitively known from the mechanical angle, a biomechanical basis is provided for pelvic floor dysfunction diseases, particularly the occurrence mechanism of bowel control dysfunction caused by perianal muscle tissue, and theoretical support is provided for the establishment of anal wilt diagnosis and treatment schemes such as traditional Chinese medicine on-line therapy and the implementation of quantitative operations.
A block diagram of a simulated anal sphincter complex stress device according to an embodiment of the present disclosure is described below with reference to fig. 11.
Fig. 11 shows a block diagram of a device 1100 for modeling anal sphincter complex stress in accordance with an embodiment of the present disclosure. The apparatus may be implemented as part or all of an electronic device by software, hardware, or a combination of both. The device comprises:
a model construction unit 1110 configured to construct an anal sphincter complex model based on pelvic CT data and anal sphincter complex muscle group MR data, wherein the anal sphincter complex model comprises an anal sphincter complex circular hole;
a stress simulation unit 1120 configured to acquire the anal sphincter complex model data in a finite element simulation system; creating an analytical rigid body in the finite element simulation system, the analytical rigid body comprising a cylindrical portion and a hemispherical portion connected to the cylindrical portion; and enabling the analytic rigid body to pass through the circular hole of the anal sphincter complex in the finite element simulation system to obtain the stress of the anal sphincter complex.
It will be appreciated by those skilled in the art that the technical solution described with reference to fig. 11 may be combined with the embodiments described with reference to fig. 1 to 10, thereby achieving the technical effects achieved by the embodiments described with reference to fig. 1 to 10. The specific content may be referred to the description above with reference to fig. 1 to 10, and the specific content is not repeated here.
The foregoing description is only of the preferred embodiments of the present disclosure and description of the principles of the technology being employed. It will be appreciated by those skilled in the art that the scope of the invention referred to in this disclosure is not limited to the specific combination of features described above, but encompasses other embodiments in which any combination of features described above or their equivalents is contemplated without departing from the inventive concepts described. Such as those described above, are mutually substituted with the technical features having similar functions disclosed in the present disclosure (but not limited thereto).
Claims (9)
1. A method of simulating anal sphincter complex stress comprising:
constructing an anal sphincter complex model based on pelvic CT data and anal sphincter complex muscle group MR data, wherein the anal sphincter complex model comprises an anal sphincter complex circular hole;
acquiring the anal sphincter complex model in a finite element simulation system;
creating an analytical rigid body in the finite element simulation system, the analytical rigid body comprising a cylindrical portion and a hemispherical portion connected to the cylindrical portion;
converting the anal sphincter complex model and the analytical rigid body into an instance in the finite element simulation system;
in the finite element simulation system, the anal sphincter complex model and the analysis rigid body are positioned at adjacent positions, the analysis rigid body is a contact main surface, the anal sphincter complex is negative in contact, the contact attribute is friction contact, and the contact type is surface-to-surface contact;
defining boundary conditions in the finite element simulation system, and restraining six degrees of freedom of two sides of the anal sphincter complex;
and enabling the analysis rigid body to pass through the circular hole of the anal sphincter complex in the finite element simulation system to obtain the stress of the anal sphincter complex.
2. The method of claim 1, wherein the diameter of the cylindrical and hemispherical portions of the analytical rigid body is 13mm-15mm.
3. The method of claim 1, wherein the analytical rigid body has an elastic modulus of 0.95KPa, a poisson's ratio of 0.45, and a density of 1.12e -9 Kg/mm 3 。
4. The method of any of claims 1-3, further comprising, prior to passing the analytical rigid body through the anal sphincter complex orifice in the finite element simulation system:
and meshing the anal sphincter complex, wherein the mesh attribute is a tetrahedral freeform mesh, and the cell type of the mesh is an explicit nonlinear cell type.
5. The method of claim 1, the passing the analytical rigid body through the anal sphincter complex orifice in the finite element simulation system, obtaining the stress of the anal sphincter complex comprising:
rotating the analysis rigid body to enable the axis direction of the analysis rigid body to be perpendicular to the circular hole of the anal sphincter complex;
the analytical rigid body is moved down through the anal sphincter complex circular aperture along the Z-axis.
6. The method of claim 1, the constructing an anal sphincter complex model based on pelvic CT data and anal sphincter complex muscle group MR data comprising:
the pelvis CT data and the anal sphincter complex muscle group MR data are subjected to data reading, image positioning, image screening and region growing treatment to obtain an anal sphincter complex global geometric model;
converting the global geometric model of the anal sphincter complex into a solid, and selecting muscle tissue from the solid from the subcutaneous part of the external sphincter at the anal margin to the lower margin of the anal sphincter complex to obtain the model of the anal sphincter complex.
7. The method of claim 6, wherein the image screening comprises screening the delineated anal sphincter complex tissue from a complex image; the region growing serves to segment out the anal sphincter complex musculature.
8. The method of claim 6, further comprising three-dimensional finite element pre-processing the anal sphincter complex global geometric model prior to said converting the global geometric model into an entity.
9. A device for simulating anal sphincter complex stress, comprising:
a model building unit configured to build an anal sphincter complex model based on pelvic CT data and anal sphincter complex muscle group MR data, wherein the anal sphincter complex model comprises an anal sphincter complex circular hole;
a stress simulation unit configured to acquire the anal sphincter complex model data in a finite element simulation system; creating an analytical rigid body in the finite element simulation system, the analytical rigid body comprising a cylindrical portion and a hemispherical portion connected to the cylindrical portion; converting the anal sphincter complex model and the analytical rigid body into an instance in the finite element simulation system; in the finite element simulation system, the anal sphincter complex model and the analysis rigid body are positioned at adjacent positions, the analysis rigid body is a contact main surface, the anal sphincter complex is negative in contact, the contact attribute is friction contact, and the contact type is surface-to-surface contact; defining boundary conditions in the finite element simulation system, and restraining six degrees of freedom of two sides of the anal sphincter complex; and enabling the analytic rigid body to pass through the circular hole of the anal sphincter complex in the finite element simulation system to obtain the stress of the anal sphincter complex.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111388883.4A CN114078355B (en) | 2021-11-22 | 2021-11-22 | Method and apparatus for simulating anal sphincter complex stress |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111388883.4A CN114078355B (en) | 2021-11-22 | 2021-11-22 | Method and apparatus for simulating anal sphincter complex stress |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114078355A CN114078355A (en) | 2022-02-22 |
| CN114078355B true CN114078355B (en) | 2023-09-26 |
Family
ID=80284236
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111388883.4A Active CN114078355B (en) | 2021-11-22 | 2021-11-22 | Method and apparatus for simulating anal sphincter complex stress |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114078355B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107405098A (en) * | 2015-10-08 | 2017-11-28 | 汉斯·葛根森 | A device for testing terminal colon and anorectal function |
| CN112169171A (en) * | 2020-10-27 | 2021-01-05 | 上海交通大学 | In situ implantation of bionic anal sphincter system based on functional reconstruction of convenience |
| CN113570723A (en) * | 2021-07-29 | 2021-10-29 | 中国人民解放军陆军军医大学 | Method for establishing female pelvic floor structure finite element model with mechanical simulation function |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9292636B2 (en) * | 2012-12-18 | 2016-03-22 | The Procter & Gamble Plaza | Method for modeling an exudate on a substrate |
| US9092585B2 (en) * | 2013-01-22 | 2015-07-28 | The Procter & Gamble Company | Computer based models for absorbent articles |
| US9715760B2 (en) * | 2014-05-16 | 2017-07-25 | University Of South Florida | Analysis-suitable geometry from discrete point sets using a mesh-free method |
-
2021
- 2021-11-22 CN CN202111388883.4A patent/CN114078355B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107405098A (en) * | 2015-10-08 | 2017-11-28 | 汉斯·葛根森 | A device for testing terminal colon and anorectal function |
| CN112169171A (en) * | 2020-10-27 | 2021-01-05 | 上海交通大学 | In situ implantation of bionic anal sphincter system based on functional reconstruction of convenience |
| CN113570723A (en) * | 2021-07-29 | 2021-10-29 | 中国人民解放军陆军军医大学 | Method for establishing female pelvic floor structure finite element model with mechanical simulation function |
Non-Patent Citations (3)
| Title |
|---|
| 基于MRI正常人肛管直肠环三维有限元建模的研究;贝绍生;李华山;刘义辉;赵艳哲;王晓锋;孙杰;孙涛;;结直肠肛门外科;第22卷(第05期);第451-454页 * |
| 有限元方法分析直肠电极在直肠角的生物电学效应;孙源;周奇;王志刚;;重庆理工大学学报(自然科学)(第12期);全文 * |
| 肛门直肠有限元模型黏膜节点位移量评价消痔灵注射的效果;柯敏辉;叶玲;陈立武;牛素生;郑鸣霄;张之玺;;中国组织工程研究(第49期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114078355A (en) | 2022-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP7100048B2 (en) | Systems and methods for generating simulated computer tomography (CT) images | |
| Li et al. | Methodology for generating a 3D computerized breast phantom from empirical data | |
| CN110458817B (en) | Medical image quality prediction method, device, equipment and storage medium | |
| CN100380404C (en) | Method and device for image processing using three-dimensional region of interest | |
| US20090149977A1 (en) | Methods, systems, and computer program products for shaping medical implants directly from virtual reality models | |
| WO2018119766A1 (en) | Multi-modal image processing system and method | |
| CN1666710A (en) | Method for producing resultant images of inspected objects | |
| CN1655193A (en) | Method and image processing system for segmenting tomographic image data | |
| JP2007172604A (en) | Method and apparatus for selecting computer-assisted algorithm based on protocol and/or parameter of acquisition system | |
| US10013797B1 (en) | Analysis-suitable geometry from discrete point sets using a mesh-free method | |
| CN114078355B (en) | Method and apparatus for simulating anal sphincter complex stress | |
| Gerasimov et al. | Application of a digital prototype for CT-based bone strength analysis | |
| US11955228B2 (en) | Methods and system for simulated radiology studies based on prior imaging data | |
| Victoria et al. | Real-time 3D ultrasound reconstruction using octrees | |
| US11227399B2 (en) | Analysis apparatus, ultrasound diagnostic apparatus, and analysis method | |
| Koc et al. | 3D Mesh Model Generation from CT and MRI data | |
| CN113538395A (en) | Image processing method, device, equipment and readable storage medium | |
| Ollé et al. | MedEdit: a computer assisted image processing and navigation system for orthopedic trauma surgery | |
| US20230394717A1 (en) | System and Method for Improved Artefact Correction in Reconstructed 2D/3D Images | |
| CN120125700B (en) | Ultrasonic CT (computed tomography) rapid reconstruction method, device and equipment based on sparse view angle | |
| CN118587304B (en) | Blood vessel imaging method and device, storage medium and electronic equipment | |
| Danping | Development of Three Dimensional Volumetric Rendering in Medical Image Application | |
| WO2010060170A1 (en) | Method for forming physical tridimensional models of a foetus in the womb using fast prototyping technology from archives generated in ultrasonography, magnetic resonance imaging and/or computer tomography equipment | |
| JP2025149688A (en) | Image processing device, image processing method, and image processing program | |
| EP1776669A1 (en) | Imaging system |
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 |