Explicit mixed strain–displacement finite elements for compressible and quasi-incompressible elasticity and plasticity

The final publication is available at Springer via http://dx.doi.org/ 10.1007/s00466-016-1305-zThis paper presents an explicit mixed finite element formulation to address compressible and quasi-incompressible problems in elasticity and plasticity. This implies that the numerical solution only involves diagonal systems of equations. The formulation uses independent and equal interpolation of displacements and strains, stabilized by variational subscales. A displacement sub-scale is introduced in order to stabilize the mean-stress field. Compared to the standard irreducible formulation, the proposed mixed formulation yields improved strain and stress fields. The paper investigates the effect of this enhancement on the accuracy in problems involving strain softening and localization leading to failure, using low order finite elements with linear continuous strain and displacement fields (P1P1 triangles in 2D and tetrahedra in 3D) in conjunction with associative frictional Mohr–Coulomb and Drucker–Prager plastic models. The performance of the strain/displacement formulation under compressible and nearly incompressible deformation patterns is assessed and compared to analytical solutions for plane stress and plane strain situations. Benchmark numerical examples show the capacity of the mixed formulation to predict correctly failure mechanisms with localized patterns of strain, virtually free from any dependence of the mesh directional bias. No auxiliary crack tracking technique is necessary.Peer ReviewedPostprint (author's final draft

A geometric approach to phase response curves and its numerical computation through the parameterization method

The final publication is available at link.springer.comThe phase response curve (PRC) is a tool used in neuroscience that measures the phase shift experienced by an oscillator due to a perturbation applied at different phases of the limit cycle. In this paper, we present a new approach to PRCs based on the parameterization method. The underlying idea relies on the construction of a periodic system whose corresponding stroboscopic map has an invariant curve. We demonstrate the relationship between the internal dynamics of this invariant curve and the PRC, which yields a method to numerically compute the PRCs. Moreover, we link the existence properties of this invariant curve as the amplitude of the perturbation is increased with changes in the PRC waveform and with the geometry of isochrons. The invariant curve and its dynamics will be computed by means of the parameterization method consisting of solving an invariance equation. We show that the method to compute the PRC can be extended beyond the breakdown of the curve by means of introducing a modified invariance equation. The method also computes the amplitude response functions (ARCs) which provide information on the displacement away from the oscillator due to the effects of the perturbation. Finally, we apply the method to several classical models in neuroscience to illustrate how the results herein extend the framework of computation and interpretation of the PRC and ARC for perturbations of large amplitude and not necessarily pulsatile.Peer ReviewedPostprint (author's final draft

An overlay J2 viscoelastic viscoplastic viscodamage model for stable shear localization problems

This work formulates a relatively simple isotropic local Overlay J2-Viscoelastic-Viscoplastic-Viscodamage constitutive model  (O-J2-VVV) which encompasses the merits of both the plastic and continuum damage formulations. The plastic component of the model account for inelastic permanent strains, while the damage component account for loss of stiffness. The plas- tic and damage softening moduli are regularized according to the material mode II fracture energy and the element size. The Orthogonal SubGrid Stabilization Method (OSGS ) is used to ensure existance and uniquess of the solution for strain shear strain localization processes, attaining global and local stability of the corresponding discrete finite element formulation. Consistent residual viscosity is used to enhance robustness and convergence of the formulation. Numerical examples show that the formulation derived is versatily, fully stable and remarkably robust, The solutions obtained are completely mesh independent, unlike those obtained with the ill-posed standard approaches. &nbsp

The uncoupling limit of identical Hopf bifurcations with an application to perceptual bistability

We study the dynamics arising when two identical oscillators are coupled near a Hopf bifurcation where we assume a parameter $\epsilon$ uncouples the system at $\epsilon=0$. Using a normal form for $N=2$ identical systems undergoing Hopf bifurcation, we explore the dynamical properties. Matching the normal form coefficients to a coupled Wilson-Cowan oscillator network gives an understanding of different types of behaviour that arise in a model of perceptual bistability. Notably, we find bistability between in-phase and anti-phase solutions that demonstrates the feasibility for synchronisation to act as the mechanism by which periodic inputs can be segregated (rather than via strong inhibitory coupling, as in existing models). Using numerical continuation we confirm our theoretical analysis for small coupling strength and explore the bifurcation diagrams for large coupling strength, where the normal form approximation breaks down