A Matlab toolbox for parametric identification of radiation-force models of ships and offshore structures

This article describes a Matlab toolbox for parametric identification of fluid-memory models associated with the radiation forces ships and offshore structures. Radiation forces are a key component of force-to-motion models used in simulators, motion control designs, and also for initial performance evaluation of wave-energy converters. The software described provides tools for preparing non-parmatric data and for identification with automatic model-order detection. The identification problem is considered in the frequency domain.

Time- vs. frequency-domain identification of parametric radiation force models for marine structures at zero speed

The dynamics describing the motion response of a marine structure in waves can be represented within a linear framework by the Cummins Equation. This equation contains a convolution term that represents the component of the radiation forces associated with fluid memory effects. Several methods have been proposed in the literature for the identification of parametric models to approximate and replace this convolution term. This replacement can facilitate the model implementation in simulators and the analysis of motion control designs. Some of the reported identification methods consider the problem in the time domain while other methods consider the problem in the frequency domain. This paper compares the application of these identification methods. The comparison is based not only on the quality of the estimated models, but also on the ease of implementation, ease of use, and the flexibility of the identification method to incorporate prior information related to the model being identified. To illustrate the main points arising from the comparison, a particular example based on the coupled vertical motion of a modern containership vessel is presented.

Air microjet system for non-contact force application and the actuation of micro-structures

We demonstrate a non-contact technique to apply calibrated and localized forces in the micro-Newton to milli-Newton range using an air microjet. An electromagnetically actuated diaphragm controlled by a signal generator is used to generate the air microjet. With a nozzle diameter of 150 mu m, the microjet diameter was maintained to a maximum of 1 mm at a distance of 5 mm from the nozzle. The force generated by the microjet was measured using a commercial force sensor to determine the velocity profile of the jet. Axial flow velocities of up to 25 m s(-1) were obtained at distances as long as 6 mm. The microjet exerted a force up to 1 mu N on a poly dimethyl siloxane (PDMS) micropillar (50 mu m in diameter, 157 mu m in height) and 415 mu N on a PDMS membrane (3 mm in diameter, 28 mu m thick). We also demonstrate that from a distance of 6 mm our microjet can exert a peak pressure of 187 Pa with a total force of about 84 mu N on a flat surface with 8 V operating voltage. Out of the cleanroom fabrication and robust design make this system cost effective and durable.

A new method for extending the range of conductive polymer sensors for contact force

This paper describes a technique for extending the force range of thin conductive polymer force sensors used for measuring contact force. These sensors are conventionally used for measuring force by changing electrical resistance when they are compressed. The new method involves measuring change in electrical resistance when the flexible sensor, which is sensitive to both compression and bending, is sandwiched between two layers of spring steel, and the structure is supported on a thin metal ring. When external force is applied, the stiffened sensor inside the spring steel is deformed within the annular center of the ring, causing the sensor to bend in proportion to the applied force. This method effectively increases the usable force range, while adding little in the way of thickness and weight. Average error for loads between 10 N and 100 N was 2.2 N (SD = 1.7) for a conventional conductive polymer sensor, and 0.9 N (SD = 0.4) using the new approach. Although this method permits measurement of greater loads with an error less than 1 N, it is limited since the modified sensor is insensitive to loads less than 5 N. These modified sensors are nevertheless useful for directly measuring normal force applied against handles and tools and other situations involving forceful manual work activities, such as grasp, push, pull, or press that could not otherwise be measured in actual work situations.

Endocardial versus epicardial ventricular radiofrequency ablation: utility of in vivo contact force assessment.

BACKGROUND Contact force (CF) is an important determinant of lesion formation for atrial endocardial radiofrequency ablation. There are minimal published data on CF and ventricular lesion formation. We studied the impact of CF on lesion formation using an ovine model both endocardially and epicardially. METHODS AND RESULTS Twenty sheep received 160 epicardial and 160 endocardial ventricular radiofrequency applications using either a 3.5-mm irrigated-tip catheter (Thermocool, Biosense-Webster, n=160) or a 3.5 irrigated-tip catheter with CF assessment (Tacticath, Endosense, n=160), via percutaneous access. Power was delivered at 30 watts for 60 seconds, when either catheter/tissue contact was felt to be good or when CF>10 g with Tacticath. After completion of all lesions, acute dimensions were taken at pathology. Identifiable lesion formation from radiofrequency application was improved with the aid of CF information, from 78% to 98% on the endocardium (P<0.001) and from 90% to 100% on the epicardium (P=0.02). The mean total force was greater on the endocardium (39±18 g versus 21±14 g for the epicardium; P<0.001) mainly because of axial force. Despite the force-time integral being greater endocardially, epicardial lesions were larger (231±182 mm(3) versus 209±131 mm(3); P=0.02) probably because of the absence of the heat sink effect of the circulating blood and covered a greater area (41±27 mm(2) versus 29±17 mm(2); P=0.03) because of catheter orientation. CONCLUSIONS In the absence of CF feedback, 22% of endocardial radiofrequency applications that are thought to have good contact did not result in lesion formation. Epicardial ablation is associated with larger lesions.

Linear latent force models using Gaussian Processes

Purely data-driven approaches for machine learning present difficulties when data are scarce relative to the complexity of the model or when the model is forced to extrapolate. On the other hand, purely mechanistic approaches need to identify and specify all the interactions in the problem at hand (which may not be feasible) and still leave the issue of how to parameterize the system. In this paper, we present a hybrid approach using Gaussian processes and differential equations to combine data-driven modeling with a physical model of the system. We show how different, physically inspired, kernel functions can be developed through sensible, simple, mechanistic assumptions about the underlying system. The versatility of our approach is illustrated with three case studies from motion capture, computational biology, and geostatistics.

Contact Force in Atrial Fibrillation: Role of Atrial Rhythm and Ventricular Contractions: Co-Force Atrial Fibrillation Study

In an experimental model, variable and intermittent contact force (CF) resulted in a significant decrease in lesion volume. In humans, variability of CF during pulmonary vein isolation has not been characterized. Methods and Results-In 20 consecutive patients undergoing CF-guided circumferential pulmonary vein isolation, 914 radiofrequency applications (530 in sinus rhythm and 384 in atrial fibrillation) were analyzed. The variability of the 60% CF range (CF60%) was 17 ± 9.6 g. Hundred seventy-one (19%) applications were delivered with constant, 717 (78%) with variable, and 26 (3%) with intermittent CF. The mean CF and force-time integral were significantly higher during applications with variable than with intermittent or constant CF. There was no significant difference in CF variability, CF60% variability, and force-time integral between applications delivered in sinus rhythm and atrial fibrillation. The main reasons for CF variability were systolo-diastolic heart movement (29%) and respiration (27%). In 10 additional patients, during adenosine-induced atrioventricular block, the minimum CF significantly increased at 19 sites (5.3 ± 4.4 versus 13.4 ± 5.9 g; P < 0.001) and at 16 sites intermittent or variable CF became constant. At only 1 site systolo-diastolic movement remained the main reason for variable CF. Conclusions-CF during pulmonary vein isolation remains highly variable despite efforts to optimize contact. CF and CF parameters were similar during sinus rhythm and atrial fibrillation. The main reasons for CF variability are systolodiastolic heart movement and respiration. The systolo-diastolic peaks and nadirs of CF are because of ventricular contractions at the large majority of pulmonary vein isolation sites.

Environmental and track perturbations on multiple pantograph interaction with catenaries in high-speed trains

The top velocity of high-speed trains is generally limited by the ability to supply the proper amount of energy through the pantograph-catenary interface. The deterioration of this interaction can lead to the loss of contact, which interrupts the energy supply and originates arcing between the pantograph and the catenary, or to excessive contact forces that promote wear between the contacting elements. Another important issue is assessing on how the front pantograph influences the dynamic performance of the rear one in trainsets with two pantographs. In this work, the track and environmental conditions influence on the pantograph-catenary is addressed, with particular emphasis in the multiple pantograph operations. These studies are performed for high speed trains running at 300 km/h with relation to the separation between pantographs. Such studies contribute to identify the service conditions and the external factors influencing the contact quality on the overhead system. (C) 2013 Elsevier Ltd. All rights reserved.

Adhesion Elastic Contact and Hysteresis Effect

In this paper, we study the relationship between the pull-off force and the transition parameter (or Tabor number) as well as the variation of the pull-off radius with the transition parameter in the adhesion elastic contact. Hysteresis models are presented to describe the contact radius as a function of external loads in loading and unloading processes. Among these models, we verified the hysteresis model from Johnson{Kendall{Roberts theory, based on which the calculated results are in good agreement with experimental ones.

Contact Sensing: A Sequential Decision Approach to Sensing Manipulation Contact

This paper describes a new statistical, model-based approach to building a contact state observer. The observer uses measurements of the contact force and position, and prior information about the task encoded in a graph, to determine the current location of the robot in the task configuration space. Each node represents what the measurements will look like in a small region of configuration space by storing a predictive, statistical, measurement model. This approach assumes that the measurements are statistically block independent conditioned on knowledge of the model, which is a fairly good model of the actual process. Arcs in the graph represent possible transitions between models. Beam Viterbi search is used to match measurement history against possible paths through the model graph in order to estimate the most likely path for the robot. The resulting approach provides a new decision process that can be use as an observer for event driven manipulation programming. The decision procedure is significantly more robust than simple threshold decisions because the measurement history is used to make decisions. The approach can be used to enhance the capabilities of autonomous assembly machines and in quality control applications.

Semi-analytic method of contact modelling

An algorithm to improve the accuracy and stability of rigid-body contact force calculation is presented. The algorithm uses a combination of analytic solutions and numerical methods to solve a spring-damper differential equation typical of a contact model. The solution method employs the recently proposed patch method, which especially suits the spring-damper differential equations. The resulting semi-analytic solution reduces the stiffness of the differential equations, while performing faster than conventional alternatives.

Field measurement of wheel-rail impact force at insulated rail joint

When wheels pass over insulated rail joints (IRJs) a vertical impact force is generated. The ability to measure the impact force is valuable as the force signature helps understand the behaviour of the IRJs, in particular their potential for failure. The impact forces are thought to be one of the main factors that cause damage to the IRJ and track components. Study of the deterioration mechanism helps finding new methods to improve the service life of IRJs in track. In this research, the strain-gage-based wheel load detector, for the first time, is employed to measure the wheel–rail contact-impact force at an IRJ in a heavy haul rail line. In this technique, the strain gages are installed within the IRJ assembly without disturbing the structural integrity of IRJ and arranged in a full wheatstone bridge to form a wheel load detector. The instrumented IRJ is first tested and calibrated in the lab and then installed in the field. For comparison purposes, a reference rail section is also instrumented with the same strain gage pattern as the IRJ. In this paper the measurement technique, the process of instrumentation, and tests as well as some typical data obtained from the field and the inferences are presented.

Practical aspects of frequency-domain identification of dynamic models of marine structures from hydrodynamic data

The motion response of marine structures in waves can be studied using finite-dimensional linear-time-invariant approximating models. These models, obtained using system identification with data computed by hydrodynamic codes, find application in offshore training simulators, hardware-in-the-loop simulators for positioning control testing, and also in initial designs of wave-energy conversion devices. Different proposals have appeared in the literature to address the identification problem in both time and frequency domains, and recent work has highlighted the superiority of the frequency-domain methods. This paper summarises practical frequency-domain estimation algorithms that use constraints on model structure and parameters to refine the search of approximating parametric models. Practical issues associated with the identification are discussed, including the influence of radiation model accuracy in force-to-motion models, which are usually the ultimate modelling objective. The illustration examples in the paper are obtained using a freely available MATLAB toolbox developed by the authors, which implements the estimation algorithms described.

On the dynamic contact angle in simulation of impinging droplets with sharp interface methods

Effects of dynamic contact angle models on the flow dynamics of an impinging droplet in sharp interface simulations are presented in this article. In the considered finite element scheme, the free surface is tracked using the arbitrary Lagrangian-Eulerian approach. The contact angle is incorporated into the model by replacing the curvature with the Laplace-Beltrami operator and integration by parts. Further, the Navier-slip with friction boundary condition is used to avoid stress singularities at the contact line. Our study demonstrates that the contact angle models have almost no influence on the flow dynamics of the non-wetting droplets. In computations of the wetting and partially wetting droplets, different contact angle models induce different flow dynamics, especially during recoiling. It is shown that a large value for the slip number has to be used in computations of the wetting and partially wetting droplets in order to reduce the effects of the contact angle models. Among all models, the equilibrium model is simple and easy to implement. Further, the equilibrium model also incorporates the contact angle hysteresis. Thus, the equilibrium contact angle model is preferred in sharp interface numerical schemes.

Moderate energy impact analysis combining phenomenological contact law with localised damage and integral equation method

A computational impact analysis methodology has been developed, based on modal analysis and a local contact force-deflection model. The contact law is based on Hertz contact theory while contact stresses are elastic, defines a modified contact theory to take account of local permanent indentation, and considers elastic recovery during unloading. The model was validated experimentally through impact testing of glass-carbon hybrid braided composite panels. Specimens were mounted in a support frame and the contact force was inferred from the deceleration of the impactor, measured by high-speed photography. A Finite Element analysis of the panel and support frame assembly was performed to compute the modal responses. The new contact model performed well in predicting the peak forces and impact durations for moderate energy impacts (15 J), where contact stresses locally exceed the linear elastic limit and damage may be deemed to have occurred. C-scan measurements revealed substantial damage for impact energies in the range of 30-50 J. For this regime the new model predictions might be improved by characterisation of the contact law hysteresis during the unloading phase, and a modification of the elastic vibration response in line with damage levels acquired during the impact. © 2011 Elsevier Ltd. All rights reserved.