Numerical study of geometric constraint and cohesive parameters in steady-state viscoelastic crack growth

This paper presents a finite element study of cohesive crack growth in a thin infinite viscoelastic strip to investigate the effects of viscoelastic properties, strip height, and cohesive model parameters on the crack growth resistance. The results of the study show that the dependence of the fracture energy on the viscoelastic properties for the strip problem is similar to that obtained for the infinite body problem even when the cohesive zone length is large compared to the height of the strip. The fracture energy also depends on the crack speed v through the dimensionless parameter v τ/L _∞ where L _∞ is the characteristic length of the cohesive zone and τ is the characteristic relaxation time of the bulk material. This relationship confirms that at least two properties of the fracture process must be prescribed accurately to model viscoelastic crack growth. In contrast, the fracture energy and crack speed are insensitive to the strip height even in situations where the growth of the dissipation zone is severely constrained by the strip boundaries. We observe that at high speeds, where the fracture energy asymptotically approaches the maximum value, the material surrounding the cohesive zone is in the rubbery (equilibrium) state and not the glassy state.     

An adaptive hybrid time‐stepping scheme for highly non‐linear strongly coupled problems

This paper deals with the design and implementation of an adaptive hybrid scheme for the solution of highly non‐linear, strongly coupled problems. The term ‘hybrid’ refers to a composite time stepping scheme where a controller decides whether a monolithic scheme or a fractional step (splitting) scheme is appropriate for a given time step. The criteria are based on accuracy and efficiency. The key contribution of this paper is the development of a framework for incorporating error criteria for stepsize selection and a mechanism for choosing from splitting or monolithic possibilities. The resulting framework is applied to silylation, a highly non‐linear, strongly coupled problem of solvent diffusion and reaction in deforming polymers. Numerical examples show the efficacy of our new hybrid scheme on both two‐ and three‐dimensional silylation simulations in the context of microlithography. Copyright © 2005 John Wiley & Sons, Ltd.     

Unilateral Buckling Restrained by Initial Force Supports

The common computation for understanding the buckling of reinforcement bars is associated with a unilateral rigid foundation with a set of uniformly spaced spring supports. In this technical note we examine the related situation of a set of uniformly spaced supports that require a finite force to initiate the displacement of the support.     

Inorganic anions induce state changes in spinach thylakoid membranes

The role of cations in excitation energy distribution between the two photosystems of photosynthesis is well established. This paper provides evidence, for the first time, for an important role of anions in the regulation of distribution of absorbed light energy between the two photosystems. Inorganic anions caused redistribution of energy more in favour of photosystem I, as judged from measurements of chlorophyll a fluorescence transients, rates of electron transport in low light and 77 K fluorescence emission spectra: the Fv/Fm ratio was decreased by inorganic anions even in the presence of DCMU, the PS II electron transport was decreased whereas PS I electron transport was increased and the F735 (77 K emission from PS I)/F685 (77 K emission from PS II) ratio was increased. Such changes were observed with inorganic anions having different valencies (Cl- , SO4(2-), PO4(3-)): the higher the valency of the inorganic anion, the more the energy transferred towards PS I. Change in the valency of the inorganic anions thus regulates distribution of absorbed light energy between the two photosystems. However, organic anions like acetate, succinate, and citrate caused no significant changes in the Fv/Fm ratio, and in rates of PS I and PS II electron transport, showing their ineffectiveness in regulating light energy distribution.     

Cyclic steady states of nonlinear electro‐mechanical devices excited at resonance

We present an efficient numerical method to solve for cyclic steady states of nonlinear electro‐mechanical devices excited at resonance. Many electro‐mechanical systems are designed to operate at resonance, where the ramp‐up simulation to steady state is computationally very expensive – especially when low damping is present. The proposed method relies on a Newton–Krylov shooting scheme for the direct calculation of the cyclic steady state, as opposed to a naïve transient time‐stepping from zero initial conditions. We use a recently developed high‐order Eulerian–Lagrangian finite element method in combination with an energy‐preserving dynamic contact algorithm in order to solve the coupled electro‐mechanical boundary value problem. The nonlinear coupled equations are evolved by means of an operator split of the mechanical and electrical problem with an explicit as well as implicit approach. The presented benchmark examples include the first three fundamental modes of a vibrating nanotube, as well as a micro‐electro‐mechanical disk resonator in dynamic steady contact. For the examples discussed, we observe power law computational speed‐ups of the form S  = 0.6·ξ  ^? 0.8, where ξ is the linear damping ratio of the corresponding resonance frequency. Copyright © 2016 John Wiley & Sons, Ltd.     

Inhibition of the reoxidation of the secondary electron acceptor of photosystem II by bicarbonate depletion

Intravenous infusion of 600 ng/kg/min of 1-sarcosine, 8-isoleucine-angiotensin II, an angiotensin II antagonist, caused a marked blood pressure fall and a decrease in plasma aldosterone in 3 patients with Bartter's syndrome. These results indicate that proximal cause of Bartter's syndrome is an arteriolar hyporesponsiveness to angiotensin II and that this angiotensin II analogue has an antagonist activity on peripheral arterioles as well as adrenal cortex.     

Anisotropic modelling and numerical simulation of brittle damage in concrete

A framework for damage mechanics of brittle solids is presented and exploited in the design and numerical implementation of an anisotropic model for the tensile failure of concrete. The key feature exploited in the analysis is the hypothesis of maximum dissipation, which specifies a unique damage rule for the elastic moduli of the solid once a failure surface is specified. A complete algorithmic treatment of the resulting model is given which renders a useful tool for large-scale inelastic finite element calculations. A rather simple three-surface failure model for concrete, containing essentially no adjustable parameters, is shown to produce results in remarkably good agreement with sample experimental data.     

A high‐order immersed boundary discontinuous‐Galerkin method for Poisson's equation with discontinuous coefficients and singular sources

We adopt a numerical method to solve Poisson's equation on a fixed grid with embedded boundary conditions, where we put a special focus on the accurate representation of the normal gradient on the boundary. The lack of accuracy in the gradient evaluation on the boundary is a common issue with low‐order embedded boundary methods. Whereas a direct evaluation of the gradient is preferable, one typically uses post‐processing techniques to improve the quality of the gradient. Here, we adopt a new method based on the discontinuous‐Galerkin (DG) finite element method, inspired by the recent work of [A.J. Lew and G.C. Buscaglia. A discontinuous‐Galerkin‐based immersed boundary method. International Journal for Numerical Methods in Engineering, 76:427‐454, 2008]. The method has been enhanced in two aspects: firstly, we approximate the boundary shape locally by higher‐order geometric primitives. Secondly, we employ higher‐order shape functions within intersected elements. These are derived for the various geometric features of the boundary based on analytical solutions of the underlying partial differential equation. The development includes three basic geometric features in two dimensions for the solution of Poisson's equation: a straight boundary, a circular boundary, and a boundary with a discontinuity. We demonstrate the performance of the method via analytical benchmark examples with a smooth circular boundary as well as in the presence of a singularity due to a re‐entrant corner. Results are compared to a low‐order extended finite element method as well as the DG method of [1]. We report improved accuracy of the gradient on the boundary by one order of magnitude, as well as improved convergence rates in the presence of a singular source. In principle, the method can be extended to three dimensions, more complicated boundary shapes, and other partial differential equations. Copyright © 2014 John Wiley & Sons, Ltd.     

An efficient time‐domain perfectly matched layers formulation for elastodynamics on spherical domains

Many practical applications require the analysis of elastic wave propagation in a homogeneous isotropic media in an unbounded domain. One widely used approach for truncating the infinite domain is the so‐called method of perfectly matched layers (PMLs). Most existing PML formulations are developed for finite difference methods based on the first‐order velocity‐stress form of the elasticity equations, and they are not straight‐forward to implement using standard finite element methods (FEMs) on unstructured meshes. Some of the problems with these formulations include the application of boundary conditions in half‐space problems and in the treatment of edges and/or corners for time‐domain problems. Several PML formulations, which do work with FEMs have been proposed, although most of them still have some of these problems and/or they require a large number of auxiliary nodal history/memory variables. In this work, we develop a new PML formulation for time‐domain elastodynamics on a spherical domain, which reduces to a two‐dimensional formulation under the assumption of axisymmetry. Our formulation is well‐suited for implementation using FEMs, where it requires lower memory than existing formulations, and it allows for natural application of boundary conditions. We solve example problems on two‐dimensional and three‐dimensional domains using a high‐order discontinuous Galerkin (DG) discretization on unstructured meshes and explicit time‐stepping. We also study an approach for stabilization of the discrete equations, and we show several practical applications for quality factor predictions of micromechanical resonators along with verifying the accuracy and versatility of our formulation. Copyright © 2014 John Wiley & Sons, Ltd.     

Non-linear B-stability and symmetry preserving return mapping algorithms for plasticity and viscoplasticity

A class of second order accurate return mapping algorithms is presented which lead to symmetric algorithmic tangent moduli and contain the classical backward-Euler return maps as a particular case. More importantly, it is shown that this class of return maps is contractive relative to the natural norm defined by the complementary Helmholz free energy function (B-stability). Since the equations of classical plasticity and viscoplasticity are shown to be contractive relative to this natural norm, the requirement of B-stability furnishes the appropriate notion of unconditionally stable algorithms for plasticity and viscoplasticity. The analysis that follows depends critically on the assumption of convexity. In particular, the models of plasticity and viscoplasticity considered obey the principle of maximum plastic dissipation. The proposed algorithms obey the discrete counterpart of this classical principle.