CN105760681B - The modeling method of more planar cone needle contacts - Google Patents

The modeling method of more planar cone needle contacts Download PDF

Info

Publication number
CN105760681B
CN105760681B CN201610110014.8A CN201610110014A CN105760681B CN 105760681 B CN105760681 B CN 105760681B CN 201610110014 A CN201610110014 A CN 201610110014A CN 105760681 B CN105760681 B CN 105760681B
Authority
CN
China
Prior art keywords
needle
contact
planar cone
planar
tissue
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201610110014.8A
Other languages
Chinese (zh)
Other versions
CN105760681A (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201610110014.8A priority Critical patent/CN105760681B/en
Publication of CN105760681A publication Critical patent/CN105760681A/en
Application granted granted Critical
Publication of CN105760681B publication Critical patent/CN105760681B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT 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

Abstract

The present invention discloses a kind of modeling method of general more planar cone needle contacts, belongs to biomethanics field.This method is derived from the contact of more planar cone needles according to the flexible deformation feature of soft tissue.Main feature of the invention are as follows: choosing more planar cone needles first is research object, and the interaction feature of more planar cone needles and elastic fibrous tissue is analyzed according to Elastic Contact Theory, establishes the contact model of more planar cone needles;And the contact problem of more planar cone needles is analyzed, by parameter transformation, study the Variation Features of different needle point outline plan taper needle contacts.Therefore analysis method proposed by the invention can be the profile optimization of planar cone needle needle, and foundation is instructed in the improvement offer of puncture resistance.

Description

The modeling method of more planar cone needle contacts
Technical field
The present invention relates to biomethanics fields, elaborate a kind of modeling method of more planar cone needle contacts.
Background technique
As clinical treatment technology develops to high, quasi-, thin direction, microtrauma puncture operation plays in treatment clinical course More and more important role.Current microtrauma puncture operation is a kind of emerging clinical means of medical domain, represents clinical medicine The new direction of development.The quality of the common tool that puncture needle is performed the operation as microtrauma puncture, puncture resistance determines Minimally Invasive Surgery Success or failure.As a kind of common puncture needle, plane puncture needle is widely used in clinical practice, in biopsy sampling, part fiber crops Important function is all played during liquor-saturated, blood sample extraction and drug injection etc..In the puncture process of needle, due to viscous by tissue Elastic and nonlinear influence, force analysis of the needle in contact process are restricted by correlation.
Currently, denaturation feature and relevant border condition of the domestic and foreign scholars according to elastic fibrous tissue, to cone, cylindrical body, The contact of the dependency rules body such as sphere has carried out detailed analysis, establishes the Elastic Contact model of dependency rule body.Due to rule Then certain regularity is presented to deformation caused by tissue contact action in body, this is mentioned for the Contact modeling analysis of dependency rule body Foundation is supplied.But the object of exposed earlier mechanics study is mostly confined to the dependency rules body such as cone, cylindrical body, sphere, but With the continuous expansion of practical application range, the contact problem between many irregular bodies is studied as current contact mechanics Problem.For more planar cone needles common during clinical diagnosis and treatment, since its geometric shape is not the axisymmetric shape of rule Body cannot directly adopt the contact analysis model of above-mentioned dependency rule body when analyzing the contact of puncture needle.
Therefore the technical solution for needing one kind new is to solve the above problems.
Summary of the invention
The technical problems to be solved by the invention are directed to the stress Variation Features of soft tissue contact phase, propose a kind of logical The modeling method of more planar cone needle contacts.
The present invention uses following technical scheme to solve above-mentioned technical problem:
A kind of modeling method of more planar cone needle contacts, when reflecting more planar cones for tissue penetration Stress, comprising the following steps:
(1) tissue penetration process is divided into non-penetration phase, stage after juxtaposition metamorphose stage and puncture;Wherein, it is contacting Deformation stage, due to the contact action by more planar cone needles, elastic fibrous tissue deforms;
(2) it is guidance with the puncture direction of more planar cone needles, establishes rectangular coordinate system in space;More planar cone needles with Elastic fibrous tissue starts point when contacting as coordinate origin, using the direction of more planar cone pin puncture elastic fibrous tissues as Z axis losing side To, and the plane perpendicular to Z axis is made as X/Y plane using coordinate origin;
(3) juxtaposition metamorphoses of more planar cone needles and elastic fibrous tissue is analyzed, according to the sShape features of more planar cone needles, When needle and elastic fibrous tissue are in contact effect, cross-sectional shape profile when elastosis is by a plurality of identical hyperbola It surrounds;Hyp item number in cross-sectional shape profile when wherein the molding flat number of taper needle is equal to metaplasia;
(4) parameter of more planar cone needle Elastic Contact models, the puncture including more planar cone needle points and tissue are set Depth D;Contact depth h for needle point with tissue;The molding flat of plane needle and angle, that is, semiapex angle θ of needle axis;More planes The normal force F that taper needle is subject to;Origin is to hyp minimum range a in metaplasia section;Origin is to double in deformed sections The maximum distance b of curve;
(5) it according to the flat pressure head of symmetric body and rib body pressure head contact problems, and is established using the equivalent law principle of function The contact model of more planar cone needles is obtained, and is derived from the expression formula of more planar cone needle normal force F;
(6) parameter transformation is carried out, it is hyp curved in metaplasia section when planar cone needle is in contact with tissue Qu Chengdu changes with the variation of paracentesis depth, since eccentricity of a hyperbola is unable to visual representation in contact model, thus it is fixed The ratio of minimum range a and maximum distance b is that parameter κ divides more planar cone needle Elastic Contact models in adopted deformed sections Analysis;The normal force of more planar cone needles is finally obtained with the Variation Features of parameter κ by parameter transformation.
Compared with the existing technology, the modeling method of more planar cone needle contacts provided by the invention, can propose one The universal modeling method to more planar cone needle contacts is planted, and considers eccentricity of a hyperbola simultaneously in the method The problem of visual representation is unable in contact model and be further carried out parameter transformation and obtained more accurately modeling.Energy Enough foundation is provided for the needle point profile optimization of more planar cone needles under different use environments and puncture object.
Detailed description of the invention
Fig. 1 is flow diagram of the invention;
Fig. 2 is the front view and top view when four planar cone needles and tissue of the invention contact;
Section deformation profile when Fig. 3 is three, four, six planar cone needles contact tissue of the invention;
Fig. 4 is the schematic diagram of planar cone needle contact resilient tissue of the invention.
Specific embodiment
The technical solution invented is described in detail with reference to the accompanying drawing:
Incorporated by reference to shown in Fig. 1.
1. the interaction feature of the needle of planar cone more than pair and elastic fibrous tissue is analyzed:
(1) tissue penetration process is analyzed.Tissue penetration process is divided into non-penetration phase, the juxtaposition metamorphose stage and Three processes of stage after puncture.Wherein, in the juxtaposition metamorphose stage, due to the contact action by needle, elastic fibrous tissue deforms.
(2) it is guidance with the puncture direction of more planar cone needles, establishes rectangular coordinate system in space.Needle and elastic fibrous tissue are opened Point when beginning to contact is as coordinate origin, using the direction of pin puncture elastic fibrous tissue as Z axis negative direction, and is made vertically with coordinate origin In Z axis plane be X/Y plane.
(3) it combines shown in Fig. 2 to Fig. 4, the juxtaposition metamorphose of more planar cone needles and elastic fibrous tissue is analyzed.According to more planes The sShape features of taper needle, when needle and elastic fibrous tissue are in contact effect, cross-sectional shape profile when elastosis is It is surrounded by a plurality of identical hyperbola.When wherein the molding flat number of taper needle is equal to metaplasia in cross-sectional shape profile Hyp item number.
2. being pushed away according to the interaction feature of needle and tissue to the establishment process of more planar cone needle contact models It leads, and its contact is analyzed as follows:
(1) bibliography " Elastic and viscoelastic indentation of flat surfaces by Pyramid indentors " in prismatic body contact model method for building up, wherein the Elastic Contact model of more planar cone needles Contain following parameter: D is the paracentesis depth of more planar cone needle points and tissue;H is contact depth of the needle point with tissue;Half pushes up Angle θ is the molding flat of plane needle and the angle of needle axis;F is the normal force that more planar cone needles are subject to;A cuts for metaplasia Origin is to hyp minimum range in face;B is the origin in deformed sections to hyp maximum distance.
(2) according to the characteristics of the geometrical relationship and needle and tissue juxtaposition metamorphose in the Elastic Contact model of planar cone needle pairs The establishment process of contact model makees following analysis.
Wherein there are following relationships by the minimum range a in contact depth h and deformed sections:
A=htan θ (1)
The general normal force F of needle are as follows:
Wherein, A is the contact area of planar cone needle and tissue;D is general puncture stress.
According to the contact of the flat pressure head of symmetric body and rib body pressure head Analysis of Contact Problem: in flat contact problems, D*, p*, F* are respectively the paracentesis depth generated in flat puncture process, puncture stress and act on the normal force of axisymmetric body; In rib body contact problems, D1,p1,F1The paracentesis depth respectively generated in rib body puncture process punctures stress and acts on rib The normal force of body.
It can be obtained by the equivalent law of function:
D*·F1=D1·F* (3)
It can be obtained by the contact problems of more planar cone bodies:
D1F*=∫ ∫Ap*·w·r·dr·dφ (4)
Wherein: r is distance of any point to origin on hyperbola in metaplasia section;W is permeability function, expression Formula are as follows:
The normal force F that can must act on more planar cone needles is solved by expression formula (1), (2), (3), (4) and (5) are as follows:
For the ease of expression, I is defined1,I2,I3It respectively indicates as follows:
Finally formula (7) Integration Solving can be obtained:
Formula (8) are substituted into formula (6) and obtain the normal force F of needle are as follows:
Wherein,
According to Barber variation conditionIt can obtain:
(3) parameter transformation is carried out.It is hyp curved in metaplasia section when planar cone needle and tissue are in contact Qu Chengdu changes with the variation of paracentesis depth.Since eccentricity of a hyperbola is unable to visual representation in contact model, thus it is fixed Adopted parameter κ analyzes model.Wherein, parameter κ is the ratio of minimum range a and maximum distance b in deformed sections, then:
Since the deformable contour of elastic fibrous tissue is that hyperbola is closed the region surrounded, then the value range of parameter κ can be obtained Are as follows: 0 < κ < cos Ф.Parameter κ is suitable for the above contact model in the case where being less than cos Ф.
Formula (13) are substituted into and obtain paracentesis depth D in formula (12):
Paracentesis depth D can be obtained by formula (14) to be codetermined by parameter κ and region angle Φ.
Formula (13), (14) are substituted into formula (9) and obtain the normal force F of more planar cone needles are as follows:
Wherein,
3. the contact of pair different needle point outline plan taper needles is analyzed.In the identical situation of paracentesis depth, The number of planes for studying aciculiform semiapex angle and needle influences the variation of the contact of taper needle.By reduce needle point semiapex angle θ or The method for reducing aciculiform number of planes can make the contact of planar cone needle become smaller.
Above-described specific embodiment has carried out further the purpose of the present invention, technical scheme and beneficial effects It is described in detail, it should be understood that being not limited to this hair the foregoing is merely a specific embodiment of the invention Bright, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should be included in the present invention Protection scope within.What the present invention was not disclosed in detail partly belongs to techniques known.

Claims (4)

1. a kind of modeling method of more planar cone needle contacts, when reflecting more planar cones for tissue penetration by Power, more planar cone needles are axisymmetric body, which comprises the following steps:
(1) tissue penetration process is divided into non-penetration phase, stage after juxtaposition metamorphose stage and puncture;Wherein, in juxtaposition metamorphose In the stage, due to the contact action by more planar cone needles, elastic fibrous tissue deforms;
(2) it is guidance with the puncture direction of more planar cone needles, establishes rectangular coordinate system in space;More planar cone needles and elasticity Organize point when starting contact as coordinate origin, using the direction of more planar cone pin puncture elastic fibrous tissues as Z axis negative direction, and Make the plane perpendicular to Z axis as X/Y plane using coordinate origin;
(3) juxtaposition metamorphose of more planar cone needles and elastic fibrous tissue is analyzed, according to the sShape features of more planar cone needles, works as needle When being in contact effect with elastic fibrous tissue, cross-sectional shape profile when elastosis is surrounded by a plurality of identical hyperbola It forms;Hyp item number in cross-sectional shape profile when wherein the molding flat number of taper needle is equal to metaplasia;
(4) parameter of more planar cone needle Elastic Contact models, the paracentesis depth including more planar cone needle points and tissue are set D;Contact depth h for needle point with tissue;The molding flat of plane needle and angle, that is, semiapex angle θ of needle axis;More planar cones The normal force F that needle is subject to;Origin is to hyp minimum range a in metaplasia section;Origin is to double in metaplasia section The maximum distance b of curve;
(5) it according to the flat pressure head of symmetric body and rib body pressure head contact problems, and establishes to obtain using the equivalent law principle of function The contact model of more planar cone needles, and it is derived from the expression formula of more planar cone needle normal force F;
(6) parameter transformation is carried out, when planar cone needle is in contact with tissue, hyp bending journey in metaplasia section Degree changes with the variation of paracentesis depth, since eccentricity of a hyperbola is unable to visual representation in contact model, therefore define a with The ratio of b is that parameter κ analyzes more planar cone needle Elastic Contact models;More plane cones are finally obtained by parameter transformation The normal force of shape needle with parameter κ Variation Features.
2. the modeling method of more planar cone needle contacts according to claim 1, which is characterized in that in step (5) Obtain the expression formula of normal force F are as follows:
Wherein,
G is the modulus of shearing of tissue;
ν is the Poisson's ratio of tissue;
A is origin in metaplasia section to hyp minimum range;
θ is the molding flat of plane needle and the angle of needle axis;
E is eccentricity of a hyperbola in metaplasia section;
Φ is region angle, i.e., the radial direction of two hyperbola crosspoints and the angle of X-axis positive direction.
3. the modeling method of more planar cone needle contacts according to claim 2, which is characterized in that in step (6), Parameter κ is the ratio of a and b, it may be assumed that
4. the modeling method of more planar cone needle contacts according to claim 1 or 2 or 3, which is characterized in that mostly flat The modeling method of face cone shape needle contact is applied to cuspidated axisymmetric body, by reducing needle point semiapex angle θ or reduction The method of aciculiform number of planes makes the contact of symmetric body become smaller.
CN201610110014.8A 2016-02-26 2016-02-26 The modeling method of more planar cone needle contacts Expired - Fee Related CN105760681B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610110014.8A CN105760681B (en) 2016-02-26 2016-02-26 The modeling method of more planar cone needle contacts

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610110014.8A CN105760681B (en) 2016-02-26 2016-02-26 The modeling method of more planar cone needle contacts

Publications (2)

Publication Number Publication Date
CN105760681A CN105760681A (en) 2016-07-13
CN105760681B true CN105760681B (en) 2018-12-14

Family

ID=56331436

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610110014.8A Expired - Fee Related CN105760681B (en) 2016-02-26 2016-02-26 The modeling method of more planar cone needle contacts

Country Status (1)

Country Link
CN (1) CN105760681B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107808698B (en) * 2017-11-08 2021-10-15 哈尔滨理工大学 Mechanical and bending deformation modeling method for soft tissue puncture by flexible needle

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105069301A (en) * 2015-08-14 2015-11-18 南通大学 Lumbar puncture virtual simulation and training system supporting haptic interaction

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105069301A (en) * 2015-08-14 2015-11-18 南通大学 Lumbar puncture virtual simulation and training system supporting haptic interaction

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Elastic and viscoelastic indentation of flat surfaces by pyramid indentors;A.E. Giannakopoulos;《Journal of the mechanics and physics of solids》;20060131(第54期);第1305-1332页 *
多平面锥形针接触建模与影响因素研究;李成刚 等;《中国机械工程》;20160630;第27卷(第12期);第1591-1597页 *
柔性针穿刺软组织的力学建模及变形测量;宿志亮;《中国优秀硕士论文数据库 医药卫生科技辑》;20150515;第2015年卷(第5期);E060-71 *
活检针几何建模与穿刺力研究;林家庆;《中国优秀硕士论文数据库 医药卫生科技辑》;20170315;第2017年卷(第3期);E060-237 *
针穿刺软组织实验与穿刺力力学建模;郭菁莘 等;《青海大学学报(自然科学版)》;20140630;第32卷(第3期);第37-42页 *

Also Published As

Publication number Publication date
CN105760681A (en) 2016-07-13

Similar Documents

Publication Publication Date Title
Enayati et al. Haptics in robot-assisted surgery: Challenges and benefits
JP6779610B2 (en) Beauty equipment
Vargas-Irwin et al. Linking objects to actions: encoding of target object and grasping strategy in primate ventral premotor cortex
JP2004024878A (en) Lancet needle with polygonal cross section
EP2184016A1 (en) Medical suturing needle
CN109273091A (en) A kind of percutaneous nephrolithy based on data in art takes stone system of virtual operation
CN108937973A (en) A kind of robotic diagnostic human body indignation mood method and device
CN105760681B (en) The modeling method of more planar cone needle contacts
Buganza Tepole et al. Multi-view stereo analysis reveals anisotropy of prestrain, deformation, and growth in living skin
CN105232120A (en) Ultrasonic guided in-plane puncture guide apparatus
Barua et al. Experimental analysis the tissue deformation of needle insertion process in tissue engineering
Hu et al. Real-time haptic feedback in laparoscopic tools for use in gastro-intestinal surgery
Barua et al. Analysis of surgical needle insertion modeling and viscoelastic tissue material interaction for minimally invasive surgery (MIS)
Lebert et al. Synchronization-based reconstruction of electromechanical wave dynamics in elastic excitable media
Lin et al. The investigation of laparoscopic instrument movement control and learning effect
US11291818B2 (en) Microneedle structure for intradermal injection
CN205251643U (en) Puncture guiding device in supersound guide plane
Li et al. Recovering stress distribution on deformable tissue for a magnetic actuated insertable laparoscopic surgical camera
CN104224141A (en) Pulse-taking signal sensor probe
WO2005065776A3 (en) Template grid needle/instrument guide and method for minimally invasive procedures
CN201094670Y (en) Puncture locating mould plate
CN214434347U (en) Non-standard double-pointed needle line movable fixing device
Fruhstorfer et al. Capillary blood sampling: how much pain is necessary? Part 4: comparison of lancets for automatic lancing devices
CN201139643Y (en) Type-B ultrasonic guidance puncture position finder
CN216933379U (en) Auxiliary forceps for laparoscope

Legal Events

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

Granted publication date: 20181214