Elliptical and circular FRP-confined concrete sections: A Mohr–Coulomb analytical model

An analytical stress–strain model is developed for predicting the compressive behavior of elliptical and circular fiber reinforced polymer (FRP)-confined concrete members. The model is based on a diagonal Poisson’s ratio formulation expressed as a function of the mechanical properties of the unconfined concrete and confining FRP jacket, the geometry of the concrete section, and the extent of internal damage in the confined concrete core. A Mohr–Coulomb yield criterion is introduced for analysis of the compressive behavior of confined concrete. Equilibrium and strain compatibility are used to obtain the ultimate compressive strength and strain of elliptical and circular FRP-confined concrete sections as a function of the effective confining stiffness of the FRP jacket. A simplified expression is derived for the FRP reinforcement ratio which precludes strain softening in elliptical and circular FRP-confined concrete sections.     

A combined experimental–numerical investigation of ductile fracture in bending of a class of ferritic–martensitic steel

We present a combined experimental–numerical study on fracture initiation at the convex surface and its propagation during bending of a class of ferritic–martensitic steel. On the experimental side, so-called free bending experiments are conducted on DP1000 steel sheets until fracture, realizing optical and scanning electron microscopy analyses on the post mortem specimens for fracture characterization. A blended Mode I – Mode II fracture pattern, which is driven by cavitation at non-metallic inclusions as well as martensitic islands and resultant softening-based intense strain localization, is observed. Phenomena like crack zig-zagging and crack alternation at the bend apex along the bending axis are introduced and discussed. On the numerical side, based on this physical motivation, the process is simulated in 2D plane strain and 3D, using Gurson’s dilatant plasticity model with a recent shear modification, strain-based void nucleation, and coalescence effects. The effect of certain material parameters (initial porosity, damage at coalescence and failure, shear modification term, etc.), plane strain constraint and mesh size on the localization and the fracture behavior are investigated in detail.     

Limit analysis and homogenization of porous materials with Mohr–Coulomb matrix. Part II: Numerical bounds and assessment of the theoretical model

This second part of the two-part study is devoted to the numerical Limit Analysis of a hollow sphere model with a Mohr–Coulomb matrix and its use for the assessment of theoretical results. Brief background and fundamental of the static and kinematic approaches in the context of numerical limit analysis are first recalled. We then present the hollow sphere model, together with its axisymmetric FEM discretization and its mechanical position. A conic programming adaptation of a previous iterative static approach, based on a piecewise linearization (PWL) of the plasticity criterion, was first realized. Unfortunately, the resulting code, no more than the PWL one, did not allow sufficiently refined meshes for loss of convergence of the conic optimizer. This problem was solved by using the projection algorithm of Ben Tal and Nemriovski (BTN) and the (interior point) linear programming code XA. For the kinematic approach, a first conic adaptation appeared also inefficient. Then, an original mixed (but fully kinematic) approach dedicated to the general Mohr–Coulomb axisymmetric problem was elaborated. The final conic mixed code appears much more robust than the classic one when using the conic code MOSEK, allowing us to take into account refined numerical meshes. After a fine validation in the case of spherical cavities and isotropic loadings (for which the exact solution is known) and comparison to previous (partial) results, numerical lower and upper bounds (a posteriori verified) of the macroscopic strength are provided. These bounds are used to assess and validate the theoretical results of the companion (part I) paper. Effects of the friction angle as well as that of the porosity are illustrated.     

Tension–torsion fracture experiments – Part II: Simulations with the extended Gurson model and a ductile fracture criterion based on plastic strain

An extension of the Gurson model that incorporates damage development in shear is used to simulate the tension–torsion test fracture data presented in Faleskog and Barsoum (2013) (Part I) for two steels, Weldox 420 and 960. Two parameters characterize damage in the constitutive model: the effective void volume fraction and a shear damage coefficient. For each of the steels, the initial effective void volume fraction is calibrated against data for fracture of notched round tensile bars and the shear damage coefficient is calibrated against fracture in shear. The calibrated constitutive model reproduces the full range of data in the tension–torsion tests thereby providing a convincing demonstration of the effectiveness of the extended Gurson model. The model reinforces the experiments by highlighting that for ductile alloys the effective plastic strain at fracture cannot be based solely on stress triaxiality. For nominally isotropic alloys, a ductile fracture criterion is proposed for engineering purposes that depends on stress triaxiality and a second stress invariant that discriminates between axisymmetric stressing and shear dominated stressing.     

Limit analysis and homogenization of porous materials with Mohr–Coulomb matrix. Part I: Theoretical formulation

The present two-part study aims at investigating the specific effects of Mohr–Coulomb matrix on the strength of ductile porous materials by using a kinematic limit analysis approach. While in the Part II, static and kinematic bounds are numerically derived and used for validation purpose, the present Part I focuses on the theoretical formulation of a macroscopic strength criterion for porous Mohr–Coulomb materials. To this end, we consider a hollow sphere model with a rigid perfectly plastic Mohr–Coulomb matrix, subjected to axisymmetric uniform strain rate boundary conditions. Taking advantage of an appropriate family of three-parameter trial velocity fields accounting for the specific plastic deformation mechanisms of the Mohr–Coulomb matrix, we then provide a solution of the constrained minimization problem required for the determination of the macroscopic dissipation function. The macroscopic strength criterion is then obtained by means of the Lagrangian method combined with Karush–Kuhn–Tucker conditions. After a careful analysis and discussion of the plastic admissibility condition associated to the Mohr–Coulomb criterion, the above procedure leads to a parametric closed-form expression of the macroscopic strength criterion. The latter explicitly shows a dependence on the three stress invariants. In the special case of a friction angle equal to zero, the established criterion reduced to recently available results for porous Tresca materials. Finally, both effects of matrix friction angle and porosity are briefly illustrated and, for completeness, the macroscopic plastic flow rule and the voids evolution law are fully furnished.     

Dynamic flight stability of hovering model insects: theory versus simulation using equations of motion coupled with Navier–Stokes equations

In the present paper, the longitudinal dynamic flight stability properties of two model insects are predicted by an approximate theory and computed by numerical sim- ulation. The theory is based on the averaged model （which assumes that the frequency of wingbeat is sufficiently higher than that of the body motion, so that the flapping wings＇ degrees of freedom relative to the body can be dropped and the wings can be replaced by wingbeat-cycle-average forces and moments）; the simulation solves the complete equations of motion coupled with the Navier-Stokes equations. Comparison between the theory and the simulation provides a test to the validity of the assumptions in the theory. One of the insects is a model dronefly which has relatively high wingbeat frequency （164 Hz） and the other is a model hawkmoth which has relatively low wingbeat frequency （26 Hz）. The results show that the averaged model is valid for the hawkmoth as well as for the dronefly. Since the wingbeat frequency of the hawkmoth is relatively low （the characteristic times of the natural modes of motion of the body divided by wingbeat period are relatively large） compared with many other insects, that the theory based on the averaged model is valid for the hawkmoth means that it could be valid for many insects.     

The onset of gas entrainment from a flowing stratified gas–liquid regime in dual discharging branches: Part II: Critical conditions at low to moderate branch Froude numbers

This is the second of a two part investigation. Experiments were performed in a 50.8 mm diameter horizontal pipe with three 6.35 mm diameter branches located at the test section mid-span. The inlet length was 1.8 m, and three branch orientations were tested at 0° (side), 45° (inclined), and 90° (bottom) from horizontal. Water and air, operating at 206 kPa, were used and both fluids flowed co-currently within the inlet in the smooth-stratified regime. The inlet superficial velocity of the liquid phase ranged between 0.04 and 0.15 m/s while in the gas phase values of 0.3, 0.4, and 1 m/s were tested. Three different dual discharge combinations were tested and included side-inclined, side-bottom, and inclined-bottom. The tested branch flow Froude numbers were limited between low to moderate values which ranged between 1 and 23. Extensive experimental data are reported for the critical conditions at the onset of gas entrainment during dual discharge. A novel map was developed for the inclined-bottom branch configuration showing the relationship between the inlet superficial liquid velocity and branch Froude numbers. This map was used to quantify the three observed modes of gas entrainment during dual discharge. These modes were classified as onset of gas entrainment in the inclined branch only, in the bottom branch only, or both the inclined and bottom branches simultaneously. The critical height at the onset of gas entrainment results are compared to published models and data sets and poor agreement was found with studies conducted in stratified gas–liquid reservoirs.