On the contact stiffness and nonlinear vibration of an elastic body with a rough surface in contact with a rigid flat surface

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

On the undamped free vibration of a mass interacting with a Hertzian contact stiffness

Three methods are used to determine the natural frequency of undamped free vibration of a mass interacting with a Hertzian contact stiffness. The exact value is determined using the first integral of motion. The harmonic balance method is used on a transformed equation for an approximate solution, and the multiple scales method is used on an approximate equation. The maximum initial displacement avoiding contact loss is also determined, and the corresponding exact natural frequency is also obtained analytically. The methods are evaluated by studying the free vibration of an elastic sphere on a flat rigid surface. © 2011 Elsevier Ltd © 2011 Elsevier Ltd. All rights reserved.

The non-ideal problem behavior using a dynamic vibration absorber with nonlinear essential stiffness and timedependent damping properties

The purpose of this study is to develop a dynamic vibration absorber using viscoelastic material with nonlinear essential stiffness and time-dependent damping properties for a non-ideal vibrating system with Sommerfeld effect, resonance capture, and jump phenomenon. The absorber is a mass-bar subsystem that consists of a viscoelastic bar with memory attached to mass, in which the internal dissipative forces depend on current, deformations, and its operational frequency varies with limited temperature. The non-ideal vibrating system consists of a linear (nonlinear) oscillator (plane frame structure) under excitation, via spring connector, of a DC-motor with limited power supply. A viscoelastic dynamic absorber modeled with elastic stiffness essentially nonlinearities was developed to further reduce the Sommerfeld effect and the response of the structure. The numerical results show the performance of the absorber on the non-ideal system response through the resonance curves, time histories, and Poincarésections. Furthermore, the structure responses using the viscoelastic damper with and without memory were studied. © IMechE 2012.

Influence of nonlinear stiffness on the dynamics of a slender elastic beam under torsional oscillations

This study focuses on analysing the effects of nonlinear torsional stiffness on the dynam-ics of a slender elastic beam under torsional oscillations, which can be subject to helical buckling.The helical buckling of an elastic beam confined in a cylinder is relevant to many applications. Someexamples include oil drilling, medical cateters and even the conformation and functioning of DNAmolecules. A recent study showed that the formation of the helical configuration is a result of onlythe torsional load, confirming that there is a different path to helical buckling which is not related tothe sinusoidal buckling, stressing the importance of the geometrical behaviour of the beam. A lowdimensional model of an elastic beam under torsional oscillations is used to analyse its dynamical be-haviour with different stiffness characteristics, which are present before and after the helical buckling.Hardening and softening characteristics are present, as the effects of torsion and bending are coupled.With the use of numerical algorithms applied to nonlinear dynamics, such as bifurcation diagramsand basins of attraction, it is shown that the nonlinear stiffness can shift the bifurcations and inducechanges in the stability of the desirable and undesirable solutions. Therefore, the proper modellingof these stiffness nonlinearities seems to be important for a better understanding of the dynamicalbehaviour of such beams.

Dissipative contact force models

"Series: Solid mechanics and its applications, vol. 226"

Experimental Vibration Characteristics of a Beam with Nonlinear Support Using Receptance Coupling Analysis

This paper applies the Multi-Harmonic Nonlinear Receptance Coupling Approach (MUHANORCA) (Ferreira 1998) to evaluate the frequency response characteristics of a beam which is clamped at one end and supported at the other end by a nonlinear cubic stiffness joint. In order to apply the substructure coupling technique, the problem was characterised by coupling a clamped linear beam with a nonlinear cubic stiffness joint. The experimental results were obtained by a sinusoidal excitation with a special force control algorithm where the level of the fundamental force is kept constant and the level of the harmonics is kept zero for all the frequencies measured.

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.

Force and displacement transmissibility of a nonlinear isolator with high-static-low-dynamic-stiffness

Engineers often face the challenge of reducing the level of vibrations experienced by a given payload or those transmitted to the support structure to which a vibrating source is attached. In order to increase the range over which vibrations are isolated, soft mounts are often used in practice. The drawback of this approach is the static displacement may be too large for reasons of available space for example. Ideally, a vibration isolator should have a high-static stiffness, to withstand static loads without too large a displacement, and at the same time, a low dynamic stiffness so that the natural frequency of the system is as low as possible which will result in an increased isolation region. These two effects are mutually exclusive in linear isolators but can be overcome if properly configured nonlinear isolators are used. This paper is concerned with the characterisation of such a nonlinear isolator comprising three springs, two of which are configured to reduce the dynamic stiffness of the isolator. The dynamic behaviour of the isolator supporting a lumped mass is investigated using force and displacement transmissibility, which are derived by modelling the dynamic system as a single-degree-of-freedom system. This results in the system dynamics being approximately described by the Duffing equation. For a linear isolator, the dynamics of the system are the same regardless if the source of the excitation is a harmonic force acting on the payload (force transmissibility) or a harmonic motion of the base (displacement transmissibility) on which the payload is mounted. In this paper these two expressions are compared for the nonlinear isolator and it is shown that they differ. A particular feature of the displacement transmissibility is that the response is unbounded at the nonlinear resonance frequency unless the damping in the isolator is greater than some threshold value, which is not the case for force transmissibility. An explanation for this is offered in the paper. (C) 2011 Elsevier Ltd. All rights reserved.

On the shock performance of a nonlinear vibration isolator with high-static-low-dynamic-stiffness

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)