This invention relates to catheters and in particular a haptically controlled active catheter instrumented with shape memory alloy actuators and a three degrees of freedom force sensor.
Angioplasty is a minimally invasive procedure that involves the insertion of a catheter into the blood vessel for removal of blockages in blood flow. In the conventional approach, the catheter is inserted into the body through the femoral artery and is guided through the lumen of the blood vessel till it reaches a blockage. A stent (superelastic Shape Memory Alloy) is then deployed to open the blood vessel at the blockage and let normal blood flow resume. The surgeon is provided with images that are obtained either by X-ray imaging or by Magnetic Resonance Imaging. These images enable the healthcare provider to track the end point of the catheter (position in absolute coordinates) in real time and determine the future course of insertion.
There are a number of problems associated with the conventional way of performing angioplasty. Specifically, the catheter insertion completely depends on the expertise and dexterity of the surgeon. In the case of intravascular neurosurgery, the catheter is pushed through extremely delicate and complex cranial blood vessels to treat aneurisms in the brain. The repeated insertion of the catheter through several trials could tear the blood vessel at the junction and cause heavy bleeding. This could also result in prolonged operating times and fatigue to clinicians and patients.
The surgeon has no method of estimating the amount of force that is being applied by the tip of the catheter on the walls of the blood vessel. Excessive pressure could rupture the blood vessel with dire consequences. Plaque could also be dislodged which may block blood vessels in the brain or heart and cause a stroke or a myocardial infarction.
The healthcare provider could have prolonged exposure to radiation or be subjected to a high-level of noise caused by the machinery generating magnetic fields for MRI. These pose danger or discomfort to the healthcare providers who perform the procedure over a prolonged period of time.
Another problem with the present procedure of Angioplasty is restenosis. The deployment of a stent at the site of a blockage only provides a temporary solution for resuming blood flow. In 40% of the procedures already done, the plaque begins to build up after a couple of years, a process called restenosis. In addition, the stent, made of superelastic Shape Memory Alloy, has a life of only 10 years. For these reasons, there may arise a need to perform repeated angioplasties for a single patient.
Accordingly it would be advantageous to provide a catheter wherein the movement of the tip can be controlled. Further it would be advantageous to provide a three degree of freedom force sensor that can be attached to a catheter.
The present invention describes a system for controlling a thin flexible thermoplastic catheter. The system includes a plurality of shape memory alloy filaments attached to the distal end of the catheter, each filament having a phase and a temperature; a means for receiving a strain value for at least one of the filaments; a means for determining the phase change that will result in the strain value, whereby the phase change is dependent on the temperature and the temperature is dependent on a voltage; and a means for setting the voltage in each filament thereby resulting in the bending of the catheter.
In another aspect of the invention there is provided a three dimensional force sensor adapted to be attached to the distal end of a catheter. The sensor includes a two-dimensional optical position sensing detector, a spring, an optical fibre embedded within the catheter and a means for determining the position of the beam spot on the detector. The two-dimensional optical position sensing detector is spaced from the distal end of the catheter. The spring attaches the detector to the catheter whereby a force acting on the detector moves the detector relative to the distal end of the catheter. A light beam is emitted from the end of an optical fibre that is embedded within the lumen of the catheter.
Further features of the invention will be described or will become apparent in the course of the following detailed description.
The invention will now be described by way of example only, with reference to the accompanying drawings, in which:
Researchers in the past have tried to instrument the catheter with Shape Memory Alloy actuators to actively control the tip of the catheter. This way the catheter can be smoothly guided into a branch in the blood vessel. The problem, however, is that the active catheters developed so far consist of actuators without any position feedback. Since most actuators used for the active catheter have a non-linear behavior, the control of the catheter by means of a joy-stick is extremely difficult and inaccurate and could cause damage to the blood vessels.
In one embodiment, this invention describes a completely autonomous control of a catheter instrumented with Shape Memory Alloy (SMA) actuators 12 and a novel three Degree Of Freedom (DOE) force sensor 14 on the tip of the catheter 16, as shown in
In addition, the catheter 16 instrumented with SMA actuators 12 and the 3-DOF force sensor 14 can be used as a tool to perform minimally invasive surgery and therapy at remote locations. The instrumented catheter 12 can be interfaced to a haptic device 20 in a master-slave configuration, as shown in
The sections below describes the design for the 3-DOF force sensor, the modeling of the Shape Memory Alloy Actuators, the LQR and H∞ loop shaping controller design, the autonomous control of the catheter, and the force reflection to a haptic device.
Referring to
A force acting on the sensor 14 due to contact with the walls of the blood vessel causes the spring 30 to bend. Since the optical fibre 32 and one end of the spring 30 are fixed to the catheter 16, the spring bends 30 relative to the catheter 16 and as a result the position of the beam spot from the optical fibre 32 shifts along the surface of the 2d PSD 34. The deviation of the beam spot along the surface of the 2d PSD 34 can be used to calculate the amount of bending in the spring which is proportional to the force acting on the tip of the catheter, by Hooke's Law. Therefore, the magnitude and direction of the force acting on the sensor can be measured accurately. The total voltage from all four outputs is an indication of the total intensity of light falling laterally on the PSD.
The force sensor 14 is mounted on the tip of the catheter 16. The sensor 14 measures the force exerted by the artery walls on the tip of the catheter. This force reading can be used to autonomously guide the catheter to prevent damage to the epithelial cells of the artery or it can be reflected to a haptic device.
Referring to
The SMA actuators 12 on the instrumented catheter 16 could be either in the form of a wire, tube or sheet. Since the machinability of SMA is low, the wires or filaments 40 should be microwelded to the stainless steel pads 42 which can be glued to the catheter, as illustrated in
In order to control the bending angle of the catheter, the amount of strain generated in the SMA actuators 12 should be controlled. For this, the physical behavior of the SMA should be described in dynamic equations. The model of the SMA set out below is based on the laws of physics. The modeling equations are as follows:
Since an SMA exists in two states, it can be modelled as a two-state system. There is a similarity between this system and an electron, which can exist only in two states—the positive spin or the negative spin. The Fermi-Dirac statistics are used to describe the number of electrons in the two states depending on the energy of the electron. The state of an SMA in the Martensite and Austenite forms using the same statistics has been modeled.
Two modelling equations were used based on whether the alloy is being heated or cooled due to hysteresis with two different transition temperatures. Since the SMA is in the Martensite form at lower temperatures, the phase transformation equation during heating is described by analogy with the Fermi-Dirac statistics in the form:
where ξ is the fraction of the Austenite phase, is the fraction of the Martensite phase prior to the present transformation from Martensite to Austenite, T is the temperature, Tfa is the transition temperature from Martensite to Austenite, νa is an indication of the range of temperature around the transition temperature Tfa during which the phase change occurs, σ is the stress and Ka is the stress curve-fitting parameter which is obtained from the loading plateau of the stress-strain characteristic with no change in temperature.
On cooling, the Austenite phase gets converted to the Martensite phase and the modelling equation during cooling is described by analogy with the Fermi-Dirac statistics in the form:
where ξa is the fraction of the Austenite phase prior to the present transformation from Austenite to Martensite, T is the temperature, Tfm is the transition temperature from Austenite to Martensite, σm is an indication of the range of temperature around the transition temperature Tfm during which the phase change occurs, σ is the stress and Km is the stress curve-fitting parameter which is obtained from the unloading part of the stress-strain characteristic.
Since the SMA is modelled as a two-component system, at any given time, the sum of the mole fractions of the Austenite and Martensite phase is 1, i.e.,
ξa+ξm=1 (3)
The time derivatives of Eqns. (1) and (2) are as follows:
For heating:
For cooling:
An SMA actuator is heated by the process of Joules heating by applying a voltage across the SMA. The loss of heat from the SMA is through natural convection. Mathematically the dynamics of the temperature are given by equation (6).
where m is the mass per unit length, cp is the specific heat capacity, V is the voltage applied across the SMA, R is the resistance per unit length, h is the coefficient of convection, A=πd is the circumferential area of cooling, d is the diameter of the wire, T is the temperature and Ta is the ambient temperature. The coefficient h is assumed to have the characteristics of a second-order polynomial to enhance the rate of convection at higher temperatures as observed in open-loop results:
h=h
0
+h
2
T
2 (7)
The constitutive equation relating changes in stress, strain, temperature and mole fraction is given by the following equation:
{dot over (σ)}=D{dot over (ε)}+θt{dot over (T)}+Ω{dot over (ξ)} (8)
where σ is the stress in the SMA, D is the average of the Young's Moduli for the Martensite and Austenite phases, e is the strain, θt is the thermal expansion factor and Ω is the phase transformation contribution factor. The model explains the hysteresis as well as the minor loops in the hysteresis, as illustrated in
The dynamic characteristics of the SMA are completely defined by equation (4) or (5) (heating or cooling), together with equations (6) and (8). The σe is defined as the integral of the error, i.e.,
{dot over (σ)}e=ε−εref (9)
where ε is the strain of the SMA actuator and εref is the reference trajectory. The dynamic equations of the SMA along with equation (9) can be represented in the state-space form:
and u is the input voltage to the SMA wire. The nonlinear equations are linearized about a set of operating points (ε0, T0, ξ0, u0) on the reference trajectory.
To obtain the operating points, ε0 is chosen as the value of the reference strain at that instant of time. T0 and ξ0 are obtained by integrating equations (4) or (5)
and (8), depending on whether the SMA is being heated or cooled, for a given value of ε0. The value of u0 is obtained from equation (6) for a given value of T0, assuming steady-state conditions.
Equation (10) is linearized about the calculated operating points to obtain linear models in the form:
For the no-load case, which is the case considered here, σ and {dot over (σ)} are equal to zero. In this case, the model given by (11) is not controllable since the number of independently controllable states is only 2. On removing the uncontrollable modes, the following state space model is obtained:
The force measured by the sensor at the tip of the catheter is provided as input to the control strategy developed to provide an output to the SMA actuators to minimize the force at the tip of the catheter. The actuators could also be activated by the surgeon to guide the catheter into the branches of the blood vessel. The control scheme should provide fast and accurate control of the strain in the SMA. Two control strategies have been developed—a Gain Scheduled PI controller and a robust H∞ loop-shaping controller using normalized coprime stabilization.
For a PI controller, the feedback is of the form,
u=−K
P(ε−εref)−KIσe+u0 (14)
where KP is the proportional gain and K, is the integral gain. Writing K=[KP KI], the gains are computed such that the resulting controller minimizes the quadratic cost function,
J(u)=∫0∞(xTQx+uTRu)dt (15)
The choice of gains for minimizing J(u) are obtained by first solving for S in the algebraic Riccati equation:
A′
T
S+SA′−SB′R
−1
B′
T
S+Q=0 (16)
The matrix K is given by:
K=R−1B′TS (17)
The simulation result for a step reference input is shown in
The experimental verification of the controller was done on a 700 MHz Windows based PC at a sampling rate of about 65 Hz. Since the form of matrices A′ and B′ are fairly simple, the solution for the Ricatti equation was obtained in closed form using MAPLE. This also greatly reduced the computation time by avoiding matrix decompositions like the Schur decomposition. The response of the SMA to a step reference input is shown in
From the experimental results, it can be observed that the response time for heating is approximately 1.0 second and for cooling is 2.1 seconds for a 0.012″ diameter SMA wire. The rate of heating and cooling is much higher for a thinner wire since the ratio of surface area to volume increases as the wire diameter is reduced, thereby increasing the rate of convectional cooling.
The experimental result for a sinusoidal reference input is shown in
The LQR based PI controller shows excellent response to a step and a sinusoidal reference input. The controller also has a good stability margin, thereby making the control law fairly insensitive to uncertainties in the parameters of the model. However, the controller is sensitive to unmodelled dynamics at high frequencies and thereby could lead to instability. For this reason, an H∞ loop-shaping controller using normalized coprime stabilization is designed such that the gains are high when the model describes the SMA accurately and low at higher frequencies when the model is inaccurate. The loop-shaping controller, however, obtains a performance/robust stability tradeoff.
The gain matrix K=[KP KI] of the PI controller, as given by equation (14), is computed by the loop-shaping procedure using normalized coprime stabilization.
The nominal plant (SMA) is first shaped by using a precompensator W1 so that the singular values of the nominal plant G are shaped to a desired open-loop shape. The compensator also ensures that the shaped plant Gs=GW1 is square and Gs has no hidden modes.
For ensuring robust stabilization, εmax is calculated as,
where {tilde over (M)}s and Ñs define the normalized coprime factors of Gs. An ε=εmax is selected to form a stabilizing controller K∞, which satisfies:
The final feedback controller K is obtained as the product of the precompensator W1 and the H∞ controller K∞ i.e.
K=W1K∞ (21)
The simulation of the robust controller and the SMA was done on MATLAB™. The desired loop shape was chosen so that the gains are high at low frequencies where the model describes the SMA accurately and the gains roll off at −20 dB beyond the corner frequency, thereby ensuring a low gain at high frequencies.
The value of the corner frequency was chosen to be 1000.424 rad·sec−1. The simulation result for input consisting of step changes is shown in
The experimental verification of the controller was done on a 700 MHz Windows based PC at a sampling rate of about 100 Hz. The response of the SMA to step changes is shown in
In the previous section, modeling and robust control of SMA actuators have been described. The simulation and experimental results show excellent tracking response for the SMA, thereby validating both the model and the control scheme. In autonomous control of an active catheter 16 fitted with SMA actuators 12 and a 3-DOF force sensor 14, the force acting on the tip of the catheter 16 is regulated so that there is no damage to the epithelial cells of the artery. A force control algorithm is implemented to ensure that the tip interaction forces are minimized.
The miniature force sensor 14 at the tip of the catheter provides the magnitude and direction of the force acting on the tip. This force is provided as a negative feedback signal to the robust control system, as shown in
One of the major problems with the conventional way of performing angioplasty is restenosis. In nearly 40% of procedures, the plaque begins to rebuild at the site after a couple of years, causing an obstruction to blood flow. In addition, the superelastic stent used for reopening the blood vessel has a lifetime of only 10 years. Therefore, the patient would have to undergo multiple angioplasties during his/her life time. A solution to this problem is to perform microsurgery at the site of the blockage.
The active catheter 16 fitted with SMA actuators 12 and a 3-DOF force sensor 14 can be controlled by a master haptic device 20, as shown in
In addition to the application in Angioplasty, the instrumented catheter 16 can also be used to perform minimally invasive fetal surgery. An abnormality can arise in some fetuses where the diaphragm separating the chest cavity from the abdomen does not completely develop. As a result the intestines begin to grow in the chest cavity, preventing the development of the lungs. This abnormality is called Congenital Diaphragmatic Hernia (CDH). In such cases, the surgeon can perform “fetal endoscopic surgery” for CDH. In this procedure, a catheter along with surgical instruments is inserted through tiny ports in the mother's womb. The surgeon operates on the fetus entirely within the womb. The trachea of the fetus is clipped to block the windpipe. This causes the trapped liquid in the lungs to aid in the growth and expansion of the lungs. This is a high risk operation and requires high dexterity and maneuverability on the part of the surgeon to position the fetus and the mother accurately. In addition, it is extremely difficult to operate on a fetus floating in the amniotic fluid without any force or visual feedback. However, the procedure is safer for the mother and the fetus in comparison to open fetal surgery, where a large incision must be made in the uterus and in the amniotic sac, since it reduces the chances of pre-term labor. A premature birth would reduce the amount of time for the fetus to develop in the womb after prenatal surgery, and reduce its chances of survival.
An instrumented catheter 16 as described in the invention would be useful in guiding the catheter into the windpipe of the fetus under image and haptic guidance. This would reduce the time taken and the chances of injury to the fetus in inserting the instrumented catheter 16 into the fetus, thereby making the procedure more reliable and efficient and reducing the amount of stress to the surgeon.
The primary application of this invention is to insert a catheter instrumented 16 with a 3-DOF force sensor 14 and Shape Memory Alloy actuators 12 smoothly into an artery and guide it through the lumen of the blood vessel till the catheter 16 reaches a blockage. The stent is then deployed to resume normal blood flow. The 3-Degrees of Freedom (DOF) force sensor 14 measures the magnitude and direction of the force being applied. It could be used for a number of other applications. Specifically it could be used at the tip of a catheter to measure the forces exerted by the walls of the blood vessel on the tip of the catheter. It could also be used as a minimally invasive tool to measure the force exerted by the tool during surgery. As well it could be used to study soft tissue properties by measuring accurately the forces exerted by the soft tissue during deformation of the tissue.
The Shape Memory alloy actuators 12 have a number of advantages. Specifically, SMA actuators have large plastic deformations which can be recovered; can generate large load forces; have high power to weight ratio; have low driving voltages; have no moving parts; are noiseless; can be used in compact workspace; and are biocompatible and corrosion resistant.
The modeling and accurate control of the SMA actuators 12 enables its use for a number of applications. Specifically it could be used as an active catheter, biomedical forceps, tissue stapler, stent, and minimally invasive, tools.
As used herein, the terms “comprises” and “comprising” are to construed as being inclusive and opened rather than exclusive. Specifically, when used in this specification including the claims, the terms “comprises” and “comprising” and variations thereof mean that the specified features, steps or components are included. The terms are not to be interpreted to exclude the presence of other features, steps or components.
It will be appreciated that the above description related to the invention by way of example only. Many variations on the invention will be obvious to those skilled in the art and such obvious variations are within the scope of the invention as described herein whether or not expressly described.
Number | Date | Country | |
---|---|---|---|
Parent | 11510846 | Aug 2006 | US |
Child | 12986001 | US |