A simple and efficient BEM implementation of quasistatic linear visco-elasticity

A simple, yet efficient procedure to solve quasistatic problems of special linear visco-elastic solids at small strains with equal rheological response in all tensorial components, utilizing boundary element method (BEM), is introduced. This procedure is based on the implicit discretisation in time (the so-called Rothe method) combined with a simple "algebraic" transformation of variables, leading to a numerically stable procedure (proved explicitly by discrete energy estimates), which can be easily implemented in a BEM code to solve initial-boundary value visco-elastic problems by using the Kelvin elastostatic fundamental solution only. It is worth mentioning that no inverse Laplace transform is required here. The formulation is straightforward for both 2D and 3D problems involving unilateral frictionless contact. Although the focus is to the simplest Kelvin-Voigt rheology, a generalization to Maxwell, Boltzmann, Jeffreys, and Burgers rheologies is proposed, discussed, and implemented in the BEM code too. A few 2D and 3D initial-boundary value problems, one of them with unilateral frictionless contact, are solved numerically

Self-consistent model for the saturation mechanism of the response to harmonic forcing in the backward-facing step flow

Certain flows denominated as amplifiers arc characterized by their global linear stability while showing large linear amplifications to sustained perturbations. As the forcing amplitude increases, a strong saturation of the response appears when compared to the linear prediction. However, a predictive model that describes the saturation of the response to higher amplitudes of forcing in stable laminar flows is still missing. While an asymptotic analysis based on the weakly nonlinear theory shows qualitative agreement only for very small forcing amplitudes, the linear response to harmonic forcing around the mean flow computed by direct numerical simulations presents a good prediction of the saturation also at higher forcing amplitudes. These results suggest that the saturation process is governed by the Reynolds stress and thus motivate the introduction of a simple self-consistent model. The model consists of a decomposition of the full nonlinear Navier-Stokes equations in a wan flow equation together with a linear perturbation equation around the mean flow, which arc coupled through the Reynolds stress. The full fluctuating response and the resulting Reynolds stress are approximated by the first harmonic calculated from the linear response to the forcing around the aforementioned mean flow. This closed set of coupled equations is solved in an iterative manner as partial nonlinearity is still preserved in the mean flow equation despite the assumed simplifications. The results show an accurate prediction of the response energy when compared to direct numerical simulations. The approximated coupling is strong enough to retain the main nonlinear effects of the saturation process. Hence, a simple physical picture is formalized, wherein the response modifies the mean flow through the Reynolds stress in such a way that the correct response energy is attained

Crack arrest through branching at curved weak interfaces: an experimental and numerical study

The phenomenon of arrest of an unstably-growing crack due to a curved weak interface is investigated. The weak interface can produce the deviation of the crack path, trapping the crack at the interface, leading to stable crack growth for certain interface geometries. This idea could be used as a technical solution for a new type of crack arrester, with a negligible impact on the global stiffness, strength and weight of the structure. In order to exploit this concept, an experimental campaign based on photo-elasticity and digital image correlation is carried out, showing the capability of curved weak interfaces to arrest cracks. The experiment is repeated for several geometrical configurations through the modification of the interface curvature radii. The phenomenon of crack deviation and subsequent arrest at the interface is also investigated with the assistance of a computational model based on the finite element method. The computational predictions provide the rationale for the interpretation of the experimental observations, and distinguish between the different behaviour of concave and convex interfaces. Consequently, as is shown in the present study, the curved interface concept fosters new routes for the attainment of structures with enhanced fracture resistance capacities, which are of paramount importance for materials and components used in extreme conditions.Comment: 19 pages, 13 figure

Beer tapping: dynamics of bubbles after impact

Beer tapping is a well known prank where a bottle of beer is impacted from the top by a solid object, usually another bottle, leading to a sudden foam overflow. A description of the shock-driven bubble dynamics leading to foaming is presented based on an experimental and numerical study evoking the following physical picture. First, the solid impact produces a sudden downwards acceleration of the bottle creating a strong depression in the liquid bulk. The existing bubbles undergo a strong expansion and a sudden contraction ending in their collapse and fragmentation into a large amount of small bubbles. Second, the bubble clouds present a large surface area to volume ratio, enhancing the CO2 diffusion from the supersaturated liquid, hence growing rapidly and depleting the CO2. The clouds of bubbles migrate upwards in the form of plumes pulling the surrounding liquid with them and eventually resulting in the foam overflow. The sudden pressure drop that triggers the bubble dynamics with a collapse and oscillations is modelled by the Rayleigh-Plesset equation. The bubble dynamics from impact to collapse occurs over a time (t(b) similar or equal to 800 mu s) much larger than the acoustic time scale of the liquid bulk (t(ac) = 2H/c similar or equal to 80 mu s), for the experimental container of height H = 6 cm and a speed of sound around c similar or equal to 1500 m/s. This scale separation, together with the comparison of numerical and experimental results, suggests that the pressure drop is controlled by two parameters: the acceleration of the container and the distance from the bubble to the free surface

Multi-level fast multipole BEM for 3-D elastodynamics

To reduce computational complexity and memory requirement for 3-D elastodynamics using the boundary element method (BEM), a multi-level fast multipole BEM (FM-BEM) based on the diagonal form for the expansion of the elastodynamic fundamental solution is proposed and demonstrated on numerical examples involving single-region and multi-region configurations where the scattering of seismic waves by a topographical irregularity or a sediment-filled basin is examined

Crack onset in stretched open hole PMMA plates considering linear and non-linear elastic behaviours

Crack onset in PMMA holed plates subjected to tensile stresses is studied experimentally and by the coupled stress and energy criterion of the Finite Fracture Mechanics (CCFFM). The elastic, strength and fracture properties of PMMA are determined by the standard tests, a clearly nonlinear stress-strain relation being identified in the tensile tests. Thus, a novel numerical implementation of the CCFFM considering a non-linear elastic (NLE) material model, using the Ramberg-Osgood approximation, in addition to the usually used linear elastic (LE) model, is developed. Testing of plates with different hole sizes shows a hole size effect in the nominal failure load as expected. For a better fitting of the experimental results, higher strength values obtained by three point bending (TPB) flatwise and edgewise coupons (without any notch), for these material models, are used, apparently for the first time, in the CCFFM predictions. This approach reflects the observation that the strength values associated to smaller but highly stressed volumes, like those located at stress maxima in the holed plates and TPB specimens, are higher. For finite-width holed plates and both material behaviours, suitable FEM models are developed to implement the CCFFM for both LE and NLE models, considering plane stress state. Moreover, an inverse procedure is devised, using the experimental data for holed plates and predictions by CCFFM, to estimate the strength and fracture properties to be used in both material models, providing very good correlations of the CCFFM predictions with the experimental results.info:eu-repo/semantics/publishedVersio

Cohesive Crack Models and Finite Fracture Mechanics analytical solutions for FRP-concrete single-lap shear test: An overview

In the present paper we review and compare several analytical models describing the single-lap shear test, which is the most common test to determine the bonding behaviour between a strengthening FRP plates and the concrete substrate. The models are one-dimensional and formulated under the assumption that debonding occurs as a pure mode II cracking process throughout a zero-thickness interface between the FRP strip and the brittle substrate. As such, they are all amenable of an analytical treatment. The FRP-concrete interface is described by at most three parameters among the interfacial fracture energy, the tensile strength and the elastic stiffness. Particularly, we compare the effective bond length estimates provided by different models and compare them with the ones present in Design Codes. Finally, a comparison with experimental data sets available in the Scientific Literature is also given