Interaction model for steel compliant riser on soft seabed

The use of catenary steel-compliant-riser (SCR) systems has increased as hydrocarbon production has moved progressively farther offshore and into deeper waters. The issue of fatigue damage caused by cyclic interaction of a riser with the seabed has gained prominence with the widespread use of SCRs and with the lengthening of the spans. The problem involves a number of complex factors, including trench configuration, nonlinear soil stiffness, breakaway of the riser from the seafloor, and degradation of soil resistance during cyclic loading. This paper presents a soilinteraction model capable of modeling these complexities, using input parameters that can be obtained with reasonable expenditure. Model simulations for typical offshore soft-soil conditions indicate that the model is capable of realistic predictions of cyclic bending moments. The degradation of soil resistance has a major effect on cyclic bending moments, particularly when uplift motions at the riser touchdown point (TDP) are large. © 2008 Society of Petroleum Engineers.

Friction-induced vibration: Model development and comparison with large-scale experimental tests

This paper presents a comparison between theoretical predictions and experimental results from a pin-on-disc test rig exploring friction-induced vibration. The model is based on a linear stability analysis of two systems coupled by sliding contact at a single point. Predictions are compared with a large volume of measured squeal initiations that have been post-processed to extract growth rates and frequencies at the onset of squeal. Initial tests reveal the importance of including both finite contact stiffness and a velocity-dependent dynamic model for friction, giving predictions that accounted for nearly all major clusters of squeal initiations from 0 to 5 kHz. However, a large number of initiations occurred at disc mode frequencies that were not predicted with the same parameters. These frequencies proved remarkably difficult to destabilise, requiring an implausibly high coefficient of friction. An attempt has been made to estimate the dynamic friction behaviour directly from the squeal initiation data, revealing complex-valued frequency-dependent parameters for a new model of linearised dynamic friction. These new parameters readily destabilised the disc modes and provided a consistent model that could account for virtually all initiations from 0 to 15 kHz. The results suggest that instability thresholds for a wide range of squeal-type behaviour can be predicted, but they highlight the central importance of a correct understanding and accurate description of dynamic friction at the sliding interface. © 2013 Elsevier Ltd. All rights reserved.

Extracting nanobelt mechanical properties from nanoindentation

A three-spring-in-series model is proposed for the nanobelt (NB) indentation test. Compared with the previous two-spring-in-series model, which considers the bending stiffness of atomic force microscope cantilever and the indenter/NB contact stiffness, this model adds a third spring of the NB/substrate contact stiffness. NB is highly flexural due to its large aspect ratio of length to thickness. The bending and lift-off of NB form a localized contact with substrate, which makes the Oliver-Pharr method [W. C. Oliver and G. M. Pharr, J. Mater. Res. 7, 1564 (1992)] and Sneddon method [I. N. Sneddon, Int. J. Eng. Sci. 3, 47 (1965)] inappropriate for NB indentation test. Because the NB/substrate deformation may have significant impact on the force-indentation depth data obtained in experiment, the two-spring-in-series model can lead to erroneous predictions on the NB mechanical properties. NB in indentation test can be susceptible to the adhesion influence because of its large surface area to volume ratio. NB/substrate contact and adhesion can have direct and significant impact on the interpretation of experimental data. Through the three-spring-in-series model, the influence of NB/substrate contact and adhesion is analyzed and methods of reducing such influence are also suggested. (C) 2010 American Institute of Physics. [doi:10.1063/1.3432748]

Delamination threshold load for dynamic impact on plates

A criterion is derived for delamination onset in transversely isotropic laminated plates under small mass, high velocity impact. The resulting delamination threshold load is about 21% higher than the corresponding quasi-static threshold load. A closed form approximation for the peak impact load is then used to predict the delamination threshold velocity. The theory is validated for a range of test cases by comparison with 3D finite element simulation using LS-DYNA and a newly developed interface element to model delamination onset and growth. The predicted delamination threshold loads and velocities are in very good agreement with the finite element simulations. Good agreement is also shown in a comparison with published experimental results. In contrast to quasi-static impacts, delamination growth occurs under a rapidly decreasing load. Inclusion of finite thickness effects and a proper description of the contact stiffness are found to be vital for accurate prediction of the delamination threshold velocity

Multibody modelling of a TB31 and a TB32 crash test with vertical portable concrete barriers: model verification and sensitivity analysis.

In this article the multibody simulation software package MADYMO for analysing and optimizing occupant safety design was used to model crash tests for Normal Containment barriers in accordance with EN 1317. The verification process was carried out by simulating a TB31 and a TB32 crash test performed on vertical portable concrete barriers and by comparing the numerical results to those obtained experimentally. The same modelling approach was applied to both tests to evaluate the predictive capacity of the modelling at two different impact speeds. A sensitivity analysis of the vehicle stiffness was also carried out. The capacity to predict all of the principal EN1317 criteria was assessed for the first time: the acceleration severity index, the theoretical head impact velocity, the barrier working width and the vehicle exit box. Results showed a maximum error of 6% for the acceleration severity index and 21% for theoretical head impact velocity for the numerical simulation in comparison to the recorded data. The exit box position was predicted with a maximum error of 4°. For the working width, a large percentage difference was observed for test TB31 due to the small absolute value of the barrier deflection but the results were well within the limit value from the standard for both tests. The sensitivity analysis showed the robustness of the modelling with respect to contact stiffness increase of ±20% and ±40%. This is the first multibody model of portable concrete barriers that can reproduce not only the acceleration severity index but all the test criteria of EN 1317 and is therefore a valuable tool for new product development and for injury biomechanics research.

Análise estrutural dos efeitos da desadaptação em componentes de próteses odontológicas utilizando o método dos elementos finitos com elemento de contato e pré-carga

Pós-graduação em Engenharia Mecânica - FEB

Mechanical microcontacts of powder particles studied with the quartz crystal microbalance

Within this work, a particle-polymer surface system is studied with respect to the particle-surface interactions. The latter are governed by micromechanics and are an important aspect for a wide range of industrial applications. Here, a new methodology is developed for understanding the adhesion process and measure the relevant forces, based on the quartz crystal microbalance, QCM. rnThe potential of the QCM technique for studying particle-surface interactions and reflect the adhesion process is evaluated by carrying out experiments with a custom-made setup, consisting of the QCM with a 160 nm thick film of polystyrene (PS) spin-coated onto the quartz and of glass particles, of different diameters (5-20µm), deposited onto the polymer surface. Shifts in the QCM resonance frequency are monitored as a function of the oscillation amplitude. The induced frequency shifts of the 3rd overtone are found to decrease or increase, depending on the particle-surface coupling type and the applied oscillation (frequency and amplitude). For strong coupling the 3rd harmonic decreased, corresponding to an “added mass” on the quartz surface. However, positive frequency shifts are observed in some cases and are attributed to weak-coupling between particle and surface. Higher overtones, i.e. the 5th and 7th, were utilized in order to derive additional information about the interactions taking place. For small particles, the shift for specific overtones can increase after annealing, while for large particle diameters annealing causes a negative frequency shift. The lower overtones correspond to a generally strong-coupling regime with mainly negative frequency shifts observed, while the 7th appears to be sensitive to the contact break-down and the recorded shifts are positive.rnDuring oscillation, the motion of the particles and the induced frequency shift of the QCM are governed by a balance between inertial forces and contact forces. The adherence of the particles can be increased by annealing the PS film at 150°C, which led to the formation of a PS meniscus. For the interpretation, the Hertz, Johnson-Kendall-Roberts, Derjaguin-Müller-Toporov and the Mindlin theory of partial slip are considered. The Mindlin approach is utilized to describe partial slip. When partial slip takes place induced by an oscillating load, a part of the contact ruptures. This results in a decrease of the effective contact stiffness. Additionally, there are long-term memory effects due to the consolidation which along with the QCM vibrations induce a coupling increase. However, the latter can also break the contact, lead to detachment and even surface damage and deformation due to inertia. For strong coupling the particles appear to move with the vibrations and simply act as added effective mass leading to a decrease of the resonance frequency, in agreement with the Sauerbrey equation that is commonly used to calculate the added mass on a QCM). When the system enters the weak-coupling regime the particles are not able to follow the fast movement of the QCM surface. Hence, they effectively act as adding a “spring” with an additional coupling constant and increase the resonance frequency. The frequency shift, however, is not a unique function of the coupling constant. Furthermore, the critical oscillation amplitude is determined, above which particle detach. No movement is detected at much lower amplitudes, while for intermediate values, lateral particle displacement is observed. rnIn order to validate the QCM results and study the particle effects on the surface, atomic force microscopy, AFM, is additionally utilized, to image surfaces and measure surface forces. By studying the surface of the polymer film after excitation and particle removal, AFM imaging helped in detecting three different meniscus types for the contact area: the “full contact”, the “asymmetrical” and a third one including a “homocentric smaller meniscus”. The different meniscus forms result in varying bond intensity between particles and polymer film, which could explain the deviation between number of particles per surface area measured by imaging and the values provided by the QCM - frequency shift analysis. The asymmetric and the homocentric contact types are suggested to be responsible for the positive frequency shifts observed for all three measured overtones, i.e. for the weak-coupling regime, while the “full contact” type resulted in a negative frequency shift, by effectively contributing to the mass increase of the quartz..rnThe interplay between inertia and contact forces for the particle-surface system leads to strong- or weak-coupling, with the particle affecting in three mentioned ways the polymer surface. This is manifested in the frequency shifts of the QCM system harmonics which are used to differentiate between the two interaction types and reflect the overall state of adhesion for particles of different size.rn

Models for Lateral Dynamic Interaction of High-Speed Trains and Bridges

In this work a methodology for analysing the lateral coupled behavior of large viaducts and high-speed trains is proposed. The finite element method is used for the structure, multibody techniques are applied for vehicles and the interaction between them is established introducing wheel-rail nonlinear contact forces. This methodology is applied for the analysis of the railway viaduct of the R´ıo Barbantino, which is a very long and tall bridge in the north-west spanish high-speed line.

Efectos dinámicos debidos al tráfico de ferrocarril sobre la infraestructura de vía y las estructuras

La vía tradicional sobre balasto sigue siendo una selección para las líneas de alta velocidad a pesar de los problemas técnicos y la prestación del funcionamiento. El problema de la vía sobre balasto es el proceso continuo del deterioro de éste debido a las cargas asociadas al tráfico ferroviario. En consecuencia es imprescindible un mantenimiento continuado para mantener un alineamiento adecuado de la vía. Por eso se surge la necesidad de comprender mejor el mecanismo involucrado en el deterioro de la vía y los factores claves que rigen su progresión a lo largo de ciclos de carga con el fin de reducir los costos del mantenimiento de la vía y mejorar el diseño de las nuevas vías. La presente tesis intenta por un lado desarrollar los modelos más adecuados y eficientes del vehículo y de la vía para los cálculos de los efectos dinámicos debido al tráfico de ferrocarril sobre la infraestructura de la vía sobre balasto, y por otro evaluar estos efectos dinámicos sobre el deterioro de la vía sobre balasto a largo plazo, empleando un adecuado modelo de predicción del deterioro de la misma. Se incluye en el trabajo una recopilación del estado del arte en lo referente a la dinámica de la vía, a la modelización del vehículo, de la vía y de la interacción entre ambos. También se hace un repaso al deterioro de la vía y los factores que influyen en su proceso. Para la primera línea de investigación de esta tesis, se han desarrollado los diferentes modelos del vehículo y de la vía y la modelización de la interacción entre ambos para los cálculos dinámicos en dos y tres dimensiones. En la interacción vehículo-vía, se ha empleado la formulación de contacto nodo-superficie para establecer la identificación de las superficies en contacto y el método de los multiplicadores de Lagrange para imponer las restricciones de contacto. El modelo de interacción se ha contrastado con los casos reportados en la literatura. Teniendo en cuenta el contacto no lineal entre rueda-carril y los perfiles de irregularidades distribuidas de la vía, se han evaluado y comparado los efectos dinámicos sobre el sistema vehículo-vía en la interacción de ambos, para distintas velocidades de circulación del vehículo, en los aspectos como la vibración del vehículo, fuerza de contacto, fuerza transmitida en los railpads, la vibración del carril. También se hace un estudio de la influencia de las propiedades de los componentes de la vía en la respuesta dinámica del sistema vehículo-vía. Se ha desarrollado el modelo del asiento de la vía que consiste en la implementación del modelo de acumulación de Bochum y del modelo de hipoplasticidad en la subrutina del usuario \UMAT" del programa ABAQUS. La implementación numérica ha sido comprobado al comparar los resultados de las simulaciones numéricas con los reportados en la literatura. Se ha evaluado la calidad geométrica de la vía sobre balasto de los tramos de estudio con datos reales de la auscultación proporcionados por ADIF (2012). Se ha propuesto una metodología de simulación, empleando el modelo de asiento, para reproducir el deterioro de la geometría de la vía. Se usan los perfiles de la nivelación longitudinal de la auscultación como perfiles de irregularidades iniciales de la vía en las simulaciones numéricas. También se evalúa la influencia de la velocidad de circulación sobre el deterioro de la vía. The traditional ballast track structures are still being used in high speed railways lines with success, however technical problems or performance features have led to ballast track solution in some cases. The considerable maintenance work is needed for ballasted tracks due to the track deterioration. Therefore it is very important to understand the mechanism of track deterioration and to predict the track settlement or track irregularity growth rate in order to reduce track maintenance costs and enable new track structures to be designed. This thesis attempts to develop the most adequate and efficient models for calculation of dynamic track load effects on railways track infrastructure, and to evaluate these dynamic effects on the track settlement, using a track settlement prediction model, which consists of the vehicle/track dynamic model previously selected and a track settlement law. A revision of the state of the knowledge regarding the track dynamics, the modelling of the vehicle, the track and the interaction between them is included. An overview related to the track deterioration and the factors influencing the track settlement is also done. For the first research of this thesis, the different models of vehicle, track and the modelling of the interaction between both have been developed. In the vehicle-track interaction, the node-surface contact formulation to establish the identification of the surfaces in contact and the Lagrange multipliers method to enforce contact constraint are used. The interaction model has been verified by contrast with some benchmarks reported in the literature. Considering the nonlinear contact between wheel-rail and the track irregularities, the dynamic effects on the vehicle-track system have been evaluated and compared, for different speeds of the vehicle, in aspects as vehicle vibration, contact force, force transmitted in railpads, rail vibration. A study of the influence of the properties of the track components on the the dynamic response of the vehicle-track system has been done. The track settlement model is developed that consist of the Bochum accumulation model and the hipoplasticity model in the user subroutine \UMAT" of the program ABAQUS. The numerical implementation has been verified by comparing the numerical results with those reported in the literature. The geometric quality of the ballast track has been evaluated with real data of auscultation provided by ADIF (2012). The simulation methodology has been proposed, using the settlement model for the ballast material, to reproduce the deterioration of the track geometry. The profiles of the longitudinal level of the auscultation is used as initial profiles of the track irregularities in the numerical simulation. The influence of the running speed on the track deterioration is also investigated.

Determining the elastic-plastic properties of metallic materials through instrumented indentation

En los últimos años, y asociado al desarrollo de la tecnología MEMS, la técnica de indentación instrumentada se ha convertido en un método de ensayo no destructivo ampliamente utilizado para hallar las características elástico-plásticas de recubrimientos y capas delgadas, desde la escala macroscópica a la microscópica. Sin embargo, debido al complejo mecanismo de contacto debajo de la indentación, es urgente proponer un método más simple y conveniente para obtener unos resultados comparables con otras mediciones tradicionales. En este estudio, el objetivo es mejorar el procedimiento analítico para extraer las propiedades elástico-plásticas del material mediante la técnica de indentación instrumentada. La primera parte se centra en la metodología llevada a cabo para medir las propiedades elásticas de los materiales elásticos, presentándose una nueva metodología de indentación, basada en la evolución de la rigidez de contacto y en la curva fuerza-desplazamiento del ensayo de indentación. El método propuesto permite discriminar los valores de indentación experimental que pudieran estar afectados por el redondeo de la punta del indentador. Además, esta técnica parece ser robusta y permite obtener valores fiables del modulo elástico. La segunda parte se centra en el proceso analítico para determinar la curva tensión-deformación a partir del ensayo de indentación, empleando un indentador esférico. Para poder asemejar la curva tension-deformación de indentación con la que se obtendría de un ensayo de tracción, Tabor determinó empíricamente un factor de constricción de la tensión () y un factor de constricción de la deformación (). Sin embargo, la elección del valor de y  necesitan una derivación analítica. Se describió analíticamente una nueva visión de la relación entre los factores de constricción de tensión y la deformación basado en la deducción de la ecuación de Tabor. Un modelo de elementos finitos y un diseño experimental se realizan para evaluar estos factores de constricción. A partir de los resultados obtenidos, las curvas tension-deformación extraidas de los ensayos de indentación esférica, afectadas por los correspondientes factores de constricción de tension y deformación, se ajustaron a la curva nominal tensión-deformación obtenida de ensayos de tracción convencionales. En la última parte, se estudian las propiedades del revestimiento de cermet Inconel 625-Cr3C2 que es depositado en el medio de una aleación de acero mediante un láser. Las propiedades mecánicas de la matriz de cermet son estudiadas mediante la técnica de indentación instrumentada, haciendo uso de las metodologías propuestas en el presente trabajo. In recent years, along with the development of MEMS technology, instrumented indentation, as one type of a non-destructive measurement technique, is widely used to characterize the elastic and plastic properties of metallic materials from the macro to the micro scale. However, due to the complex contact mechanisms under the indentation tip, it is necessary to propose a more convenient and simple method of instrumented indention to obtain comparable results from other conventional measurements. In this study, the aim is to improve the analytical procedure for extracting the elastic plastic properties of metallic materials by instrumented indentation. The first part focuses on the methodology for measuring the elastic properties of metallic materials. An alternative instrumented indentation methodology is presented. Based on the evolution of the contact stiffness and indentation load versus the depth of penetration, the possibility of obtaining the actual elastic modulus of an elastic-plastic bulk material through instrumented sharp indentation tests has been explored. The proposed methodology allows correcting the effect of the rounding of the indenter tip on the experimental indentation data. Additionally, this technique does not seem too sensitive to the pile-up phenomenon and allows obtaining convincing values of the elastic modulus. In the second part, an analytical procedure is proposed to determine the representative stress-strain curve from the spherical indentation. Tabor has determined the stress constraint factor (stress CF), and strain constraint factor (strain CF), empirically but the choice of a value for and is debatable and lacks analytical derivation. A new insight into the relationship between stress and strain constraint factors is analytically described based on the formulation of Tabor’s equation. Finite element model and experimental tests have been carried out to evaluate these constraint factors. From the results, representative stress-strain curves using the proposed strain constraint factor fit better with the nominal stress-strain curve than those using Tabor’s constraint factors. In the last part, the mechanical properties of an Inconel 625-Cr3C2 cermet coating which is deposited onto a medium alloy steel by laser cladding has been studied. The elastic and plastic mechanical properties of the cermet matrix are studied using depth-sensing indentation (DSI) on the micro scale.

On the possibility of elastic strain localisation in a fault

The phenomenon of strain localisation is often observed in shear deformation of particulate materials, e.g., fault gouge. This phenomenon is usually attributed to special types of plastic behaviour of the material (e.g., strain softening or mismatch between dilatancy and pressure sensitivity or both). Observations of strain localisation in situ or in experiments are usually based on displacement measurements and subsequent computation of the displacement gradient. While in conventional continua the symmetric part of the displacement gradient is equal to the strain, it is no longer the case in the more realistic descriptions within the framework of generalised continua. In such models the rotations of the gouge particles are considered as independent degrees of freedom the values of which usually differ from the rotation of an infinitesimal volume element of the continuum, the latter being described for infinitesimal deformations by the non-symmetric part of the displacement gradient. As a model for gouge material we propose a continuum description for an assembly of spherical particles of equal radius in which the particle rotation is treated as an independent degree of freedom. Based on this model we consider simple shear deformations of the fault gouge. We show that there exist values of the model parameters for which the displacement gradient exhibits a pronounced localisation at the mid-layers of the fault, even in the absence of inelasticity. Inelastic effects are neglected in order to highlight the role of the independent rotations and the associated additional parameters. The localisation-like behaviour occurs if (a) the particle rotations on the boundary of the shear layer are constrained (this type of boundary condition does not exist in a standard continuum) and (b) the contact moment-or bending stiffness is much smaller than the product of the effective shear modulus of the granulate and the square of the width of the gouge layer. It should be noted however that the virtual work functional is positive definite over the range of physically meaningful parameters (here: contact stiffnesses, solid volume fraction and coordination number) so that strictly speaking we are not dealing with a material instability.

Investigating the mesh dependency and upscaling of 3D grain-based models for the simulation of brittle fracture processes in low-porosity crystalline rock

The application of 3D grain-based modelling techniques is investigated in both small and large scale 3DEC models, in order to simulate brittle fracture processes in low-porosity crystalline rock. Mesh dependency in 3D grain-based models (GBMs) is examined through a number of cases to compare Voronoi and tetrahedral grain assemblages. Various methods are used in the generation of tessellations, each with a number of issues and advantages. A number of comparative UCS test simulations capture the distinct failure mechanisms, strength profiles, and progressive damage development using various Voronoi and tetrahedral GBMs. Relative calibration requirements are outlined to generate similar macro-strength and damage profiles for all the models. The results confirmed a number of inherent model behaviors that arise due to mesh dependency. In Voronoi models, inherent tensile failure mechanisms are produced by internal wedging and rotation of Voronoi grains. This results in a combined dependence on frictional and cohesive strength. In tetrahedral models, increased kinematic freedom of grains and an abundance of straight, connected failure pathways causes a preference for shear failure. This results in an inability to develop significant normal stresses causing cohesional strength dependence. In general, Voronoi models require high relative contact tensile strength values, with lower contact stiffness and contact cohesional strength compared to tetrahedral tessellations. Upscaling of 3D GBMs is investigated for both Voronoi and tetrahedral tessellations using a case study from the AECL’s Mine-by-Experiment at the Underground Research Laboratory. An upscaled tetrahedral model was able to reasonably simulate damage development in the roof forming a notch geometry by adjusting the cohesive strength. An upscaled Voronoi model underestimated the damage development in the roof and floor, and overestimated the damage in the side-walls. This was attributed to the discretization resolution limitations.

Nonlinear diffusion and exclusion processes with contact interactions

Exclusion processes on a regular lattice are used to model many biological and physical systems at a discrete level. The average properties of an exclusion process may be described by a continuum model given by a partial differential equation. We combine a general class of contact interactions with an exclusion process. We determine that many different types of contact interactions at the agent-level always give rise to a nonlinear diffusion equation, with a vast variety of diffusion functions D(C). We find that these functions may be dependent on the chosen lattice and the defined neighborhood of the contact interactions. Mild to moderate contact interaction strength generally results in good agreement between discrete and continuum models, while strong interactions often show discrepancies between the two, particularly when D(C) takes on negative values. We present a measure to predict the goodness of fit between the discrete and continuous model, and thus the validity of the continuum description of a motile, contact-interacting population of agents. This work has implications for modeling cell motility and interpreting cell motility assays, giving the ability to incorporate biologically realistic cell-cell interactions and develop global measures of discrete microscopic data.

Suitability of a point-contact SAW-based stiffness transducer for biomedical applications

Methods of diagnosis in Biomedical applications can be broadly divided into contact and non-contact based methods. So far, ultrasound based methods have been found to be most favorable for non-contact, non-invasive diagnosis, especially in the case of tissue stiffness analysis. We report here, the fabrication and characterization details of a new contact based transducer system for qualitative determination of the stiffnesses of non-piezoelectric substrates using the phenomenon of Surface Acoustic Waves (SAW). Preliminary trials to study the functionality of this system were carried out on various metallic and non-metallic substrates, and the results were found to be satisfactory. To confirm the suitability of this system for biomedical applications, similar trials have been conducted on tissue mimicking phantoms with varying degrees of stiffness.

Nonlinear dynamics of atomic force microscopy with intermittent contact

When the atomic force microscopy (AFM) in tapping mode is in intermittent contact with a soft substrate, the contact time can be a significant portion of a cycle, resulting in invalidity of the impact oscillator model, where the contact time is assumed to be infinitely small. Furthermore, we demonstrate that the AFM intermittent contact with soft substrate can induce the motion of higher modes in the AFM dynamic response. Traditional ways of modeling AFM (one degree of freedom (DOF) system or single mode analysis) are shown to have serious mistakes when applied to this kind of problem. A more reasonable displacement criterion on contact is proposed, where the contact time is a function of the mechanical properties of AFM and substrate, driving frequencies/amplitude, initial conditions, etc. Multi-modal analysis is presented and mode coupling is also shown. (c) 2006 Published by Elsevier Ltd.