Explicit solutions for depressed masonry arches loaded until collapse—Part II: a solution method for statically indeterminate systems

Meccanica - In this second part of the paper we use the results obtained in the first part in order to build the solution to two statically indeterminate cases of a depressed masonry arch. As in...     

A 2.5-D dynamic model for a saturated porous medium: Part I. Green’s function

Based on Biot’s theory, the dynamic 2.5-D Green’s function for a saturated porous medium is obtained using the Fourier transform and the potential decomposition methods. The 2.5-D Green’s function corresponds to the solutions for the following two problems: the point force applied to the solid skeleton, and the dilatation source applied within the pore fluid. By performing the Fourier transform on the governing equations for the 3-D Green’s function, the governing differential equations for the two parts of the 2.5-D Green’s function are established and then solved to obtain the dynamic 2.5-D Green’s function. The derived 2.5-D Green’s function for saturated porous media is verified through comparison with the existing solution for 2.5-D Green’s function for the elastodynamic case and the closed-form 3-D Green’s function for saturated porous media. It is further demonstrated that a simple form 2-D Green’s function for saturated porous media can be been obtained using the potential decomposition method.     

Analytical study of the nonlinear behavior of a shape memory oscillator: Part II—resonance secondary

In this work, the response of a single-degree-of-freedom shape memory oscillator subjected to the excitation harmonic has been investigated. Equation of motion is formulated assuming a polynomial constitutive model to describe the restitution force of the oscillator. Here the method of multiple scales is used to obtain an approximate solution to the equations of the motion describing the modulation equations of amplitude and phase, and to investigate theoretically its stability. This work is presented in two parts. In Part I of this study we showed the modeling of the problem where the free vibration of the oscillator at low temperature is analyzed, where martensitic phase is stable. Part I also presents the investigation dynamics of the primary resonance of the pseudoelastic oscillator. Part II of the work is focused on the study in the secondary resonance of a pseudoelastic oscillator using the model developed in Part I. The analysis of the system in Part I as well as in Part II is accomplished numerically by means of phase portraits, Lyapunov exponents, power spectrum and Poincare maps. Frequency-response curves are constructed for shape memory oscillators for various excitation levels and detuning parameter. A rich class of solutions and bifurcations, including jump phenomena and saddle-node bifurcations, is found.     

An enhanced beam-theory model of the mixed-mode bending (MMB) test—Part I: Literature review and mechanical model

The paper presents a mechanical model of the mixed-mode bending (MMB) test used to assess the mixed-mode interlaminar fracture toughness of composite laminates. The laminated specimen is considered as an assemblage of two sublaminates partly connected by an elastic–brittle interface. The problem is formulated through a set of 36 differential equations, accompanied by suitable boundary conditions. Solution of the problem is achieved by separately considering the two subproblems related to the symmetric and antisymmetric parts of the loads, which for symmetric specimens correspond to fracture modes I and II, respectively. Explicit expressions are determined for the interfacial stresses, internal forces, and displacements.     

A refined theory for laminated beams: Part I—A new high order approach

This paper presents a new displacement-based one-dimensional model for the analysis of multilayered composite beams. The kinematic restriction of cross sections rigid in their own planes is introduced. The axial displacements over the cross sections are represented in terms of explicitly defined piecewise polynomial warping functions with discontinuous derivatives at the interlaminae, whereas the amplitude of the displacements along the beam axis is established by means of a variational formulation. It is proved that the proposed representation of the axial displacements yields the exact solution of the interior domain problem for a beam subjected to a transverse load varying according to a polynomial law. It is shown that two or three coordinate functions are sufficient to yield continuous distributions of equilibrated stresses except for small neighborhoods of the constrained cross sections, where a higher number of warping functions could be used in order to obtain a better accuracy. The numerical results show excellent agreement with plane stress finite element and plane strain exact solutions.

---

Sommario In questo lavoro viene presentato un nuovo modello monodimensionale per l'analisi di travi composite multistrato. Viene introdotta l'ipotesi di indeformabilita delle sezioni nel proprio piano mentre gli spostamenti assiali nella sezione sono rappresentati facendo uso di funzioni ingobbamento definite sull'intera altezza e con derivata discontinua all'in erlamina. Infine, l'ampiezza degli spostamenti lungo l'asse della trave è determinata facendo uso di una formulazione variazionale. Si mostra come la rappresentazione degli spostamenti assiali proposta sia in grado di fornire la soluzione esatta, all'interno del dominio, per una trave soggetta ad un carico trasversale variabile con legge nolinomiale. Due o tre funzioni coordinate sono sufficienti a fornire distribuzioni di sforzi che verificano l'equilibrio anche all'interlamina, a meno di zone rislrette in vicinanza di sezioni vincolate. I risultati numerici mostrano un eccellente accordo con soluzioni agli elementi finiti in stato piano di tensione e con soluzioni esatte in stato piano di deformazione.

     

Analytical study of the nonlinear behavior of a shape memory oscillator: Part I—primary resonance and free response at low temperatures

In this work, the response of a single-degree-of-freedom shape memory oscillator subjected to the excitation harmonic has been investigated. Equation of motion is formulated assuming a polynomial constitutive model to describe the restitution force of the oscillator. Here the method of multiple scales is used to obtain an approximate solution to the equations of the motion describing the modulation equations of amplitude and phase, and to investigate theoretically its stability. This work is presented in two parts. In Part I of this study we showed the modeling of the problem where the free vibration of the oscillator at low temperature is analyzed, where martensitic phase is stable. Part I also presents the investigation dynamics of the primary resonance of the pseudoelastic oscillator. Part II of the work is focused on the study in the secondary resonance of a pseudoelastic oscillator using the model developed in Part I. The analysis of the system in Part I as well as in Part II is accomplished numerically by means of phase portraits, Lyapunov exponents, power spectrum and Poincare maps. Frequency-response curves are constructed for shape memory oscillators for various excitation levels and detuning parameter. A rich class of solutions and bifurcations, including jump phenomena and saddle-node bifurcations, is found.     

Interaction between a submerged evacuated cylindrical shell and a shock wave—Part I: Diffraction–radiation analysis

The interaction between a submerged elastic circular cylindrical shell and an external shock wave is addressed. A linear, two-dimensional formulation of the problem is considered. A semi-analytical solution is obtained using a combination of the classical analytical approach based on the use of the Laplace transform and separation of variables, and finite difference methodology. The study consists of two parts. Part I focuses on the simulation and analysis of the acoustic fields induced during the interaction. Both the diffraction (absolutely rigid cylinder) and complete diffraction–radiation (elastic shell) are considered. Special attention is paid to the lower-magnitude shell-induced waves representing radiation by the elastic waves circumnavigating the shell. The focus of Part II is on the numerical analysis of the solution. The convergence of the series solution and finite-difference scheme is analysed. The computation of the response functions of the problem is discussed as well, as is the effect of the bending stiffness on the acoustic field. The membrane model of the shell is considered to analyse such an effect, which, in combination with the models addressed in Part I, allows for the analysis of the evolution of the acoustic field around the structure as its elastic properties change from an absolutely rigid cylinder to a membrane. The results of the numerical simulations are compared to available experimental data, and a good agreement is observed.     

A micromechanics-based strain gradient damage model for fracture prediction of brittle materials – Part I: Homogenization methodology and constitutive relations

In this paper, we first describe a homogenization methodology with the aim of establishing strain gradient constitutive relations for heterogeneous materials. The methodology presented in this work includes two main steps. The first one is the construction of the average strain-energy density for a well-chosen RVE by using a homogenization technique. The second one is the transformation of the obtained average strain-energy density to that for the continuum. An important characteristic of this method is its self-consistency with respect to the choice of the RVE: the strain gradient constitutive law built by using the present method is independent of the size and the form of the RVE. In the frame of this homogenization procedure, we have constructed a strain gradient constitutive relation for a two-dimensional elastic material with many microcracks by adopting the self-consistent scheme. It was shown that the effective behavior of cracked solids depends not only on the crack density but also on the average crack size with which the strain gradient is associated. The proposed constitutive relation provides a starting point for the development of an evolution law of damage including strain gradient effect, which will be presented in the second part of this work.     

Fluid velocity based simulation of hydraulic fracture: a penny shaped model—part I: the numerical algorithm

In the first part of this paper, a universal fluid velocity based algorithm for simulating hydraulic fracture with leak-off, previously demonstrated for the PKN and KGD models, is extended to obtain solutions for a penny-shaped crack. The numerical scheme is capable of dealing with both the viscosity and toughness dominated regimes, with the fracture being driven by a power-law fluid. The computational approach utilizes two dependent variables; the fracture aperture and the reduced fluid velocity. The latter allows for the application of a local condition of the Stefan type (the speed equation) to trace the fracture front. The obtained numerical solutions are carefully tested using various methods, and are shown to achieve a high level of accuracy.     

Fluid velocity based simulation of hydraulic fracture—a penny shaped model. Part II: new,accurate semi-analytical benchmarks for an impermeable solid

In the first part of this paper a universal fluid velocity based algorithm for simulating hydraulic fracture with leak-off was created for a penny-shaped crack. The power-law rheological model of fluid was assumed and the final scheme was capable of tackling both the viscosity and toughness dominated regimes of crack propagation. The obtained solutions were shown to achieve a high level of accuracy. In this paper simple, accurate, semi-analytical approximations of the solution are provided for the zero leak-off case, for a wide range of values of the material toughness and parameters defining the fluid rheology. A comparison with other results available in the literature is undertaken.     

Tension–torsion fracture experiments—Part I: Experiments and a procedure to evaluate the equivalent plastic strain

Ductile failure experiments on a double notched tube (DNT) specimen subjected to a combination of tensile load and torque that was applied at a fixed ratio is presented. The experimental results extend those in Barsoum and Faleskog (2007a) down to zero stress triaxiality. A new and robust evaluation procedure for such tests is proposed, and a simple relation for the equivalent plastic strain at failure for combined normal and shear deformation, respectively, is developed. Tests were carried out on the medium strength medium hardening steel Weldox 420, and the high strength low hardening steel Weldox 960. The experimental results unanimously show that ductile failure not only depends on stress triaxiality, but is also strongly affected by the type of deviatoric stress state that prevails, which can be quantified by a stress invariant that discriminates between axisymmetric stressing and shear dominated stressing, e.g., the Lode parameter. Additional experiments on round notch bar (RNB) specimens are recapitulated in order to give a comprehensive account on how ductile failure depends on stress triaxiality, ranging from zero to more than 1.6, and the type of stress state for the two materials tested. This provides an extensive experimental data base that will be used to explore an extension of the Gurson model that incorporates damage development in shear presented in Xue et al. (2013) (Part II).     

A two dimensional Euler–Lagrangian model of wood gasification in a charcoal bed—Part II: Parameter influence and comparison

A Euler–Lagrangian simulation was employed for a comprehensive parameter study of wood gasification in a fluidized charcoal bed. The parameters that were varied include the initial bed temperature, fuel mass flow rate, inert tar fraction, and kinetic energy losses caused by particle–particle and particle–wall collisions. The results of each parameter variation are compared with a base scenario, previously described in detail in Part I of this study (Gerber & Oevermann, 2014). The results are interpreted by comparing the reactor outlet temperature, averaged particle temperature, overall wood mass, overall charcoal mass, concentrations of several gaseous species, and axial barycenter data for particles obtained with different sets of parameters. The inert tar fraction and fuel mass flow rate are the most sensitive parameter, while the particle–particle and particle–wall contact parameters have only a small impact on the results. Increasing the reactive tar components by 19% almost doubled the amount of reactive tars at the reactor outlet, while decreasing the restitution coefficients of the particle collisions by 0.2 results in higher overall gas production but almost no change in bed height. Herein, our numerical results are discussed in detail while assessing the model restrictions.     

The PELskin project—part I: fluid–structure interaction for a row of flexible flaps: a reference study in oscillating channel flow

Previous studies of flexible flaps attached to the aft part of a cylinder have demonstrated a favourable effect on the drag and lift force fluctuation. This observation is thought to be linked to the excitation of travelling waves along the flaps and as a consequence of that, periodic shedding of the von Kármán vortices is altered in phase. A more general case of such interaction is studied herein for a limited row of flaps in an oscillating flow; representative of the cylinder case since the transversal flow in the wake-region shows oscillating character. This reference case is chosen to qualify recently developed numerical methods for the simulation of fluid–structure interaction in the context of the EU funded ‘PELskin’ project. The simulation of the two-way coupled dynamics of the flexible elements is achieved via a structure model for the flap motion, which was implemented and coupled to two different fluid solvers via the immersed boundary method. The results show the waving behaviour observed at the tips of the flexible elements in interaction with the fluid flow and the formation of vortices in the gaps between the flaps. In addition, formation of vortices upstream of the leading and downstream of the trailing flap is seen, which interact with the formation of the shear-layer on top of the row. This leads to a phase shift in the wave-type motion along the row that resembles the observation in the cylinder case.