Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite-element”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Thermal Buckling Analysis of Rectangular Panels Subjected to Humped Temperature Profile Heating

This research investigates thermal buckling characteristics of rectangular panels subjected to different types of humped temperature profile heating. Minimum potential energy and finite-element methods are used to calculate the panel buckling temperatures. The two methods give fairly close thermal buckling solutions. 'Buckling temperature magnification factor of the first kind, eta' is established for the fixed panel edges to scale up the buckling solution of uniform temperature loading case to give the buckling solution of the humped temperature profile loading cases. Also, 'buckling temperature magnification factor of the second kind, xi' is established for the free panel edges to scale up the buckling solution of humped temperature profile loading cases with unheated boundary heat sinks to give the buckling solutions when the boundary heat sinks are heated up.

Ko, William I.↗
NASA NTRS2001

Dynamic Analysis of a Building Under Rocket Engine Plume Acoustic Load

Studies have been performed to develop finite-element modeling and simulation techniques to predict the dynamic structural response of Building 4010 to the acoustic load from the plume of high-thrust rocket motors. The building is the Test Control Center and general office space for the E-complex at Stennis Space Center. It is a large single span; light-structured building located approximately 1,000 feet from the E-1 test stand. A three-dimensional shell/beam combined model of the building was built using Pro/Engineer platform and imported into Pro/Mechanica for analysis. An Equivalent Shell technique was developed to simplify the highly complex building structure so that the calculation is more efficient and accurate. A deterministic approach was used for the dynamic analysis. A pre-stressed modal analysis was performed to simulate the weight stiffening of the structure, through which about 200 modes ranging from 0 to 35 Hz were identified. In an initial dynamic frequency analysis, the maximum response over the model was found. Then the complete 3-D distributions of the displacement, as well as the stresses, were calculated through a final frequency analysis. The results were compared to a strain gage and accelerometer recordings from rocket engine tests and showed reasonable agreement.

Qian, Z.↗
NASA NTRS2002

Manufacturing Distortions of Curved Composite Panels

This papa briefly discusses the influences of through-thickness thermal expansion, a misaligned ply, and a resin-rich slightly thicker ply on the deformations of a curved composite laminate during cool down from tbc cure temperature. Both two-dimensional and three-dimensional level finite-element analyses are used. The deformations are categorized as to radial and tangential deformations and twist, and for each of the three influences, these deformations are quantified. An additional outcome of the study is an indication of the level of analysis needed to study each of these three influences.

Starnes, James H., Jr.↗
NASA NTRS2003

An Active Smart Material Control System for F/A-18 Buffet Alleviation

The vertical tail buffet problem of fighter aircraft occurs at high angles of attack when the vortical flow breaks down ahead of the vertical tails resulting in unsteady and unbalanced pressure loads on the vertical tails. The buffet loads imposed upon the vertical tails resulted in a premature fatigue failure of the tails, and consequently limits the performance and super maneuverability of twin-tail fighter aircraft. An active smart material control system using distributed piezoelectric actuators has been developed for buffet alleviation and is presented. The inboard and outboard surfaces of the vertical tail are equipped with piezoelectric actuators to control the buffet responses in the first bending and torsion modes. The electrodynamics of the piezoelectric actuators are expressed with a three-dimensional finite-element model. A single-input-single-output controller is designed to drive the active piezoelectric actuators. High-fidelity multidisciplinary analysis modules for the fluid dynamics, structure dynamics, electrodynamics of the piezoelectric actuators, control law, fluid structure interfacing, and grid motion are integrated into a multidisciplinary computing environment that controls the temporal synchronization of the analysis modules. At 30 degree angle of attack, RMS values of tip acceleration are reduced by as much as 12%. The peak values of the power spectral density of tail-tip acceleration are reduced by as much as 22% in the first bending mode and by as much as 82% in the first torsion mode. The actively controlled piezoelectric actuators were also effective in adding damping at wide range of angles of attack.

Sheta, Essam F.↗
NASA NTRS2000

Multidisciplinary High-Fidelity Analysis and Optimization of Aerospace Vehicles: Preliminary Results - Part 2

An objective of the High Performance Computing and Communication Program at the NASA Langley Research Center is to demonstrate multidisciplinary shape and sizing optimization of a complete aerospace vehicle configuration by using high-fidelity finite-element structural analysis and computational fluid dynamics aerodynamic analysis in a distributed, heterogeneous computing environment that includes high performance parallel computing. A software system has been designed and implemented to integrate a set of existing discipline analysis codes, some of them computationally intensive, into a distributed computational environment for the design of a high-speed civil transport configuration. The paper describes both the preliminary results from implementing and validating the multidisciplinary analysis and the results from an aerodynamic optimization. The discipline codes are integrated by using the Java programming language and a Common Object Request Broker Architecture compliant software product. A companion paper describes the formulation of the multidisciplinary analysis and optimization system.

Walsh, J. L.↗
NASA NTRS1998

Design-Oriented Analysis of Aircraft Fuselage Structures

A design-oriented analysis capability for aircraft fuselage structures that utilizes equivalent plate methodology is described. This new capability is implemented as an addition to the existing wing analysis procedure in the Equivalent Laminated Plate Solution (ELAPS) computer code. The wing and fuselage analyses are combined to model entire airframes. The paper focuses on the fuselage model definition, the associated analytical formulation and the approach used to couple the wing and fuselage analyses. The modeling approach used to minimize the amount of preparation of input data by the user and to facilitate the making of design changes is described. The fuselage analysis is based on ring and shell equations but the procedure is formulated to be analogous to that used for plates in order to take advantage of the existing code in ELAPS. Connector springs are used to couple the wing and fuselage models. Typical fuselage analysis results are presented for two analytical models. Results for a ring-stiffened cylinder model are compared with results from conventional finite-element analyses to assess the accuracy of this new analysis capability. The connection of plate and ring segments is demonstrated using a second model that is representative of the wing structure for a channel-wing aircraft configuration.

Giles, Gary L.↗
NASA NTRS1998

Comparison of Response Surface and Kriging Models for Multidisciplinary Design Optimization

In this paper, we compare and contrast the use of second-order response surface models and kriging models for approximating non-random, deterministic computer analyses. After reviewing the response surface method for constructing polynomial approximations, kriging is presented as an alternative approximation method for the design and analysis of computer experiments. Both methods are applied to the multidisciplinary design of an aerospike nozzle which consists of a computational fluid dynamics model and a finite-element model. Error analysis of the response surface and kriging models is performed along with a graphical comparison of the approximations, and four optimization problems m formulated and solved using both sets of approximation models. The second-order response surface models and kriging models-using a constant underlying global model and a Gaussian correlation function-yield comparable results.

Simpson, Timothy W.↗
NASA NTRS1998

Ducted-Fan Engine Acoustic Predictions using a Navier-Stokes Code

A Navier-Stokes computer code is used to predict one of the ducted-fan engine acoustic modes that results from rotor-wake/stator-blade interaction. A patched sliding-zone interface is employed to pass information between the moving rotor row and the stationary stator row. The code produces averaged aerodynamic results downstream of the rotor that agree well with a widely used average-passage code. The acoustic mode of interest is generated successfully by the code and is propagated well upstream of the rotor; temporal and spatial numerical resolution are fine enough such that attenuation of the signal is small. Two acoustic codes are used to find the far-field noise. Near-field propagation is computed by using Eversman's wave envelope code, which is based on a finite-element model. Propagation to the far field is accomplished by using the Kirchhoff formula for moving surfaces with the results of the wave envelope code as input data. Comparison of measured and computed far-field noise levels show fair agreement in the range of directivity angles where the peak radiation lobes from the inlet are observed. Although only a single acoustic mode is targeted in this study, the main conclusion is a proof-of-concept: Navier-Stokes codes can be used both to generate and propagate rotor/stator acoustic modes forward through an engine, where the results can be coupled to other far-field noise prediction codes.

Rumsey, C. L.↗
NASA NTRS1985

Subtle alternating electrocardiographic morphology as an indicator of decreased cardiac electrical stability

Observations from finite-element computer models, together with analytic developments based on percolation theory have suggested that subtle fluctuations of ECG morphology might serve as an indicator diminished cardiac electrical stability. With fixed-rate atrial pacing in canines, we have previously observed a pattern of alternation in T wave energy which correlated with cardiac electrical stability. We report here on a series of 20 canine experiments in which cardiac electrical stability (measured via Ventricular Fibrillation Threshold determination) was compared to a non-degenerate, multidimensional measurement of the degree of alternating activity present in the ECG complex morphology. The decrease in cardiac electrical stability brought on by both coronary artery occlusion and systemic hypothermia was consistently accompanied by subtle alternation in ECG morphology, with the absolute degree of alternating activity being significantly (negatively) correlated with cardiac electrical stability.

Ventricular Fibrillation↗
NASA NTRS2003

Modeling of Biological Fluid Systems

This study presents a detailed methodology for modeling and validating three- dimensional fluid flow in realistic reconstructed vasculature geometry. Optical scanned data from a sectioned stenotic lumen having -65% area constriction within an arterial vessel of 9 mm diameter is processed and converted into a finite-element mesh. Mass and momentum transfer equations within the reconstructed stenotic vessel with a-flow rate of 200 ml/rnin is solved using the Galarkin Finite Element formulation. The calculated velocities at several axial locations along the length of the vessel are compared with the Particle Image Velocimetry (PIV) measurements. The computation results are in agreement with the initial steady state PIV measurements. In addition, initial results from pulsatile flow computations at different time points of a pulse cycle will be presented.

Banerjee, Rupak↗
NASA NTRS1998

Gel stretch method: a new method to measure constitutive properties of cardiac muscle cells

Diastolic dysfunction is an important cause of congestive heart failure; however, the basic mechanisms causing diastolic congestive heart failure are not fully understood, especially the role of the cardiac muscle cell, or cardiocyte, in this process. Before the role of the cardiocyte in this pathophysiology can be defined, methods for measuring cardiocyte constitutive properties must be developed and validated. Thus this study was designed to evaluate a new method to characterize cardiocyte constitutive properties, the gel stretch method. Cardiocytes were isolated enzymatically from normal feline hearts and embedded in a 2% agarose gel containing HEPES-Krebs buffer and laminin. This gel was cast in a shape that allowed it to be placed in a stretching device. The ends of the gel were held between a movable roller and fixed plates that acted as mandibles. Distance between the right and left mandibles was increased using a stepper motor system. The force applied to the gel was measured by a force transducer. The resultant cardiocyte strain was determined by imaging the cells with a microscope, capturing the images with a CCD camera, and measuring cardiocyte and sarcomere length changes. Cardiocyte stress was characterized with a finite-element method. These measurements of cardiocyte stress and strain were used to determine cardiocyte stiffness. Two variables affecting cardiocyte stiffness were measured, the passive elastic spring and viscous damping. The passive spring was assessed by increasing the force on the gel at 1 g/min, modeling the resultant stress vs. strain relationship as an exponential [sigma = A/k(ekepsilon - 1)]. In normal cardiocytes, A = 23.0 kN/m2 and k = 16. Viscous damping was assessed by examining the loop area between the stress vs. strain relationship during 1 g/min increases and decreases in force. Normal cardiocytes had a finite loop area = 1.39 kN/m2, indicating the presence of viscous damping. Thus the gel stretch method provided accurate measurements of cardiocyte constitutive properties. These measurements have allowed the first quantitative assessment of passive elastic spring properties and viscous damping in normal mammalian cardiocytes.

Non-NASA Center↗
NASA NTRS1992

Mathematical modelling of flow distribution in the human cardiovascular system

The paper presents a detailed model of the entire human cardiovascular system which aims to study the changes in flow distribution caused by external stimuli, changes in internal parameters, or other factors. The arterial-venous network is represented by 325 interconnected elastic segments. The mathematical description of each segment is based on equations of hydrodynamics and those of stress/strain relationships in elastic materials. Appropriate input functions provide for the pumping of blood by the heart through the system. The analysis employs the finite-element technique which can accommodate any prescribed boundary conditions. Values of model parameters are from available data on physical and rheological properties of blood and blood vessels. As a representative example, simulation results on changes in flow distribution with changes in the elastic properties of blood vessels are discussed. They indicate that the errors in the calculated overall flow rates are not significant even in the extreme case of arteries and veins behaving as rigid tubes.

NASA Center JSC↗
NASA NTRS2000

Propulsion Aeroelastic Analysis Developed for Flutter and Forced Response

The NASA Glenn Research Center at Lewis Field develops new technologies to increase the fuel efficiency of aircraft engines, improve the safety of engine operation, reduce emissions, and reduce engine noise. With the development of new designs for fans, compressors, and turbines to achieve these goals, the basic aeroelastic requirements are that there should be no flutter (self-excited vibrations) or high resonant blade stresses (due to forced response) in the operating regime. Therefore, an accurate prediction and analysis capability is required to verify the aeroelastic soundness of the designs. Such a three-dimensional viscous propulsion aeroelastic analysis capability has been developed at Glenn with support from the Advanced Subsonic Technology (AST) program. This newly developed aeroelastic analysis capability is based on TURBO, a threedimensional unsteady aerodynamic Reynolds-averaged Navier-Stokes turbomachinery code developed previously under a grant from Glenn. TURBO can model the viscous flow effects that play an important role in certain aeroelastic problems such as flutter with flow separation, flutter at high loading conditions near the stall line (stall flutter), flutter in the presence of shock and boundary-layer interaction, and forced response due to wakes and shock impingement. In aeroelastic analysis, the structural dynamics representation of the blades is based on normal modes. A finite-element analysis code is used to calculate these in-vacuum vibration modes and the associated natural frequencies. In an aeroelastic analysis using the TURBO code, flutter and forced response are modeled as being uncoupled. To calculate if a blade row will flutter, one prescribes the motion of the blade to be a harmonic vibration in a specified in-vacuum normal mode. An aeroelastic analysis preprocessor is used to generate the displacement field required for the analysis. The work done by aerodynamic forces on the vibrating blade during a cycle of vibration is calculated. If this work is positive, the blade is dynamically unstable, since it will extract energy from the flow, leading to an increase in the blade s oscillation amplitude. The forced-response excitations on a blade row are calculated by modeling the flow through two adjacent blade rows using the TURBO code. The blades are assumed to be rigid. As an option, a single blade row can be modeled with the upstream blade row influence represented by a time-varying disturbance (gust) at the inlet boundary. The unsteady forces on a blade row from such analyses are used in a structural analysis along with the blade structural dynamics characteristics and aerodynamic damping associated with blade vibration to calculate the resulting dynamic stresses on the blade.

Bakhle, Milind A.↗
NASA NTRS1998

Ice-shelf Dynamics Near the Front of Filchner-Ronne Ice Shelf, Antarctica, Revealed by SAR Interferometry: Model/Interferogram Comparison

We compare Earth Remote Sensing (ERS) satellite synthetic-aperture radar (SAR) interferograms with artificial interferograms constructed using output of a finite-element ice-shelf flow model to study the dynamics of Filchner-Ronne Ice Shelf (FRIS) near Hemmen Ice Rise (HIR) where the iceberg-calving front itersects Berkener Island (BI).

Ice-shelf Filchner-Ronne Antarctica SAR Interferom↗
NASA NTRS1995

NASA Tech Briefs, November 1995

The contents include: 1) Mission Accomplished; 2) Resource Report: Marshall Space Flight Center; 3) NASA 1995 Software of the Year Award; 4) Microbolometers Based on Epitaxial YBa2Cu3O(sub 7-x) Thin Films; 5) Garnet Random-Access Memory; 6) Fabrication of SNS Weak Links on SOS Substrates; 7) High-Voltage MOSFET Switching Circuit; 8) Asymmetric Switching for a PWM H-Bridge Power Circuit; 9) Better Ohmic Contacts for InP Semiconductor Devices; 10) Low-Bandgap Thermovoltaic Materials and Devices; 11) Digital Frequency-Differencing Circuit; 12) Imaging Magnetometer; 13) Computer-Assisted Monitoring of a Complex System; 14) Buffered Telemetry Demodulator; 15) Compact Multifunction Inspection Head; 16) Optical Detection of Fractures in Ceramic Diaphragms; 17) Eddy-Current Detection of Cracks in Reinforced Carbon/Carbon; 18) Apparent Thermal Conductivity of Multilayer Insulation; 19) Optimizing Misch-Metal Compositions in Metal Hydride Anodes; 20) Device for Sampling Surface Contamination; 21) Probabilistic Failure Assessment for Fatigue; 22) Probabilistic Fatigue and Flaw-Propagation Analysis; 23) Windows Program for Driving the TDU-850 Printer; 24) Subband/Transform MATLAB Functions for Processing Images; 25) Computing Equilibrium Chemical Compositions; 26) Program Processes Thermocouple Readings; 27) ICAN-Second-Generation Integrated Composite Analyzer; 28) Integrated Composite Analyzer with Damping Capabilities; 29) Computing Efficiency of Transfer of Microwave Power; 30) Program Calculates Power Demands of Electronic Designs; 31) Cost-Estimation Program; 32) Program Estimates Areas Required by Electronic Designs; 33) Program to Balance Mapped Turbopump Assemblies; 34) BiblioTech; 35) Controlling Mirror Tilt With a Bimorph Actuator; 36) Burst-Disk Device Simulates Effect of Pyrotechnic Device; 37) Bearing-Mounting Concept Accommodates Thermal Expansion; 38) Parallel-Plate Acoustic Absorbers for Hot Environments; 39) Adjustable-Length Strut Withstands Large Cyclic Loads; 40) Tool Indicates Contact Angles in Bearing Raceways; 41) Gravity Slides With Magnetic Braking; 42) High-Torque, Lightweight, Pneumatically Driven Wrench for Small Spaces; 43) Device for Testing Compatibility of an O-Ring; 44) Magnetic Heat Pump Containing Flow Diverters; 45) Variable-Tilt Helicopter Rotor Mast; 46) "Beach-Ball" Robotic Rovers; 47) Apparatus Would Measure Temperatures of Ball Bearings; 48) Flexible Borescope for Inspecting Ducts; 49) Texturing Copper To Reduce Secondary Emission of Electrons; 50) Automated Laser Cutting in Three Dimensions; 51) Algorithm Helps Monitor Engine Operation; 52) Flexible Revision of Data-Processing Communications; 53) Software for Managing the Use of Land; 54) Thermal Strap Increases Cryocooling Efficiency; 55) Reversible Nut With Engagement Indication; 56) Control Algorithms for Kinematically Redundant Manipulators; 57) Computed Hydrogen-Flow Splits in a Rocket Engine; 58) Pressure and Thermal Modeling of Rocket Launches; 59) Field of View of a Spacecraft Antenna: Analysis and Software; 60) Digital Controller for Laser-Beam-Steering Subsystem; 61) More About Beam-Steering Subsystem for Laser Communication; 62) Digital Controller for Laser-Beam-Steering Subsystem: Part 2; 63) Interface Circuit Board for Space-Shuttle Communications; 64) Automated Planning of Spacecraft Telecommunications; 65) Artifacts of Spectral Analysis of Instrument Readings; 66) Neural-Network Controller for Vibration Suppression; 67) Adaptive Finite-Element Computation in Fracture Mechanics; 68) Attitude Control for the Cassini Spacecraft; 69) Analytical Model for Fluid Dynamics in a Microgravity Environment; 70) Study of Rocket-Engine Joints Bonded by NVCU/NARloy-Z; 71) Improved Silicon Nitride for Advanced Heat Engines; 72) Parameters for Welding Aluminum/Lithium Alloys; 73) Lightweight Composite Intertank Structure; 74) Foil Patches Seal Small Vacuum Leaks; 75) Data Base on Cables and Connectors; 76) Effect of Clock Mode on Radiation Hardnessf an ADC; and 77) Fault-Tolerant Control for a Robotic Inspection System.

Source record↗
NASA NTRS2003

Explicit Pore Pressure Material Model in Carbon-Cloth Phenolic

An explicit material model that uses predicted pressure in the pores of a carbon-cloth phenolic (CCP) composite has been developed. This model is intended to be used within a finite-element model to predict phenomena specific to CCP components of solid-fuel-rocket nozzles subjected to high operating temperatures and to mechanical stresses that can be great enough to cause structural failures. Phenomena that can be predicted with the help of this model include failures of specimens in restrained-thermal-growth (RTG) tests, pocketing erosion, and ply lifting

Gutierrez-Lemini, Danton↗
NASA NTRS2003

NASA Tech Briefs, October 2003

Topics covered include: Cryogenic Temperature-Gradient Foam/Substrate Tensile Tester; Flight Test of an Intelligent Flight-Control System; Slat Heater Boxes for Thermal Vacuum Testing; System for Testing Thermal Insulation of Pipes; Electrical-Impedance-Based Ice-Thickness Gauges; Simulation System for Training in Laparoscopic Surgery; Flasher Powered by Photovoltaic Cells and Ultracapacitors; Improved Autoassociative Neural Networks; Toroidal-Core Microinductors Biased by Permanent Magnets; Using Correlated Photons to Suppress Background Noise; Atmospheric-Fade-Tolerant Tracking and Pointing in Wireless Optical Communication; Curved Focal-Plane Arrays Using Back-Illuminated High-Purity Photodetectors; Software for Displaying Data from Planetary Rovers; Software for Refining or Coarsening Computational Grids; Software for Diagnosis of Multiple Coordinated Spacecraft; Software Helps Retrieve Information Relevant to the User; Software for Simulating a Complex Robot; Software for Planning Scientific Activities on Mars; Software for Training in Pre-College Mathematics; Switching and Rectification in Carbon-Nanotube Junctions; Scandia-and-Yttria-Stabilized Zirconia for Thermal Barriers; Environmentally Safer, Less Toxic Fire-Extinguishing Agents; Multiaxial Temperature- and Time-Dependent Failure Model; Cloverleaf Vibratory Microgyroscope with Integrated Post; Single-Vector Calibration of Wind-Tunnel Force Balances; Microgyroscope with Vibrating Post as Rotation Transducer; Continuous Tuning and Calibration of Vibratory Gyroscopes; Compact, Pneumatically Actuated Filter Shuttle; Improved Bearingless Switched-Reluctance Motor; Fluorescent Quantum Dots for Biological Labeling; Growing Three-Dimensional Corneal Tissue in a Bioreactor; Scanning Tunneling Optical Resonance Microscopy; The Micro-Arcsecond Metrology Testbed; Detecting Moving Targets by Use of Soliton Resonances; and Finite-Element Methods for Real-Time Simulation of Surgery.

Source record↗
NASA NTRS2003

Simulation System for Training in Laparoscopic Surgery

A computer-based simulation system creates a visual and haptic virtual environment for training a medical practitioner in laparoscopic surgery. Heretofore, it has been common practice to perform training in partial laparoscopic surgical procedures by use of a laparoscopic training box that encloses a pair of laparoscopic tools, objects to be manipulated by the tools, and an endoscopic video camera. However, the surgical procedures simulated by use of a training box are usually poor imitations of the actual ones. The present computer-based system improves training by presenting a more realistic simulated environment to the trainee. The system includes a computer monitor that displays a real-time image of the affected interior region of the patient, showing laparoscopic instruments interacting with organs and tissues, as would be viewed by use of an endoscopic video camera and displayed to a surgeon during a laparoscopic operation. The system also includes laparoscopic tools that the trainee manipulates while observing the image on the computer monitor (see figure). The instrumentation on the tools consists of (1) position and orientation sensors that provide input data for the simulation and (2) actuators that provide force feedback to simulate the contact forces between the tools and tissues. The simulation software includes components that model the geometries of surgical tools, components that model the geometries and physical behaviors of soft tissues, and components that detect collisions between them. Using the measured positions and orientations of the tools, the software detects whether they are in contact with tissues. In the event of contact, the deformations of the tissues and contact forces are computed by use of the geometric and physical models. The image on the computer screen shows tissues deformed accordingly, while the actuators apply the corresponding forces to the distal ends of the tools. For the purpose of demonstration, the system has been set up to simulate the insertion of a flexible catheter in a bile duct. [As thus configured, the system can also be used to simulate other endoscopic procedures (e.g., bronchoscopy and colonoscopy) that include the insertion of flexible tubes into flexible ducts.] A hybrid approach has been followed in developing the software for real-time simulation of the visual and haptic interactions (1) between forceps and the catheter, (2) between the forceps and the duct, and (3) between the catheter and the duct. The deformations of the duct are simulated by finite-element and modalanalysis procedures, using only the most significant vibration modes of the duct for computing deformations and interaction forces. The catheter is modeled as a set of virtual particles uniformly distributed along the center line of the catheter and connected to each other via linear and torsional springs and damping elements. The interactions between the forceps and the duct as well as the catheter are simulated by use of a ray-based haptic-interaction- simulating technique in which the forceps are modeled as connected line segments.

Basdogan, Cagatay↗