Analysis of nonlinear stability and post-buckling for Euler-type beam-column structure

Based on the assumption of finite deformation, the Hamilton variational principle is extended to a nonlinear elastic Euler-type beam-column structure located on a nonlinear elastic foundation. The corresponding three-dimensional (3D) mathematical model for anaiyzing the nonlinear mechanical behaviors of structures is established, in which the effects of the rotation inertia and the nonlinearity of material and geometry are considered. As an application, the nonlinear stability and the post-buckling for a linear elastic beam with the equal cross-section located on an elastic foundation are analyzed.One end of the beam is fully fixed, and the other end is partially fixed and subjected to an axial force. A new numerical technique is proposed to calculate the trivial solution,bifurcation points, and bifurcation solutions by the shooting method and the Newton-Raphson iterative method. The first and second bifurcation points and the corresponding bifurcation solutions are calculated successfully. The effects of the foundation resistances and the inertia moments on the bifurcation points are considered.     

Multiple bifurcations and local energy minimizers in thermoelastic martensitic transformations

Thermoelastic martensitic transformations in shape memory alloys can be modeled on the basis of nonlinear elastic theory.Microstructures of fine phase mixtures are local energy minimizers of the total energy.Using a one-dimensional effective model,we have shown that such microstructures are inhomogeneous solutions of the nonlinear Euler-Lagrange equation and can appear upon loading or unloading to certain critical conditions,the bifurcation conditions.A hybrid numerical method is utilized to calculate the inhomogeneous solutions with a large number of interfaces.The characteristics of the solutions are clarified by three parameters:the number of interfaces,the interface thickness,and the oscillating amplitude.Approximated analytical expressions are obtained for the interface and inhomogeneity energies through the numerical solutions.     

Effect of bulk deformation and surface tension on the growth of a coherent inclusion

The object of this paper is the study of the moving boundary problem modeling the growth of a spherical solid inclusion in an infinite solid matrix. The displacement in bulk is assumed infinitesimal, while the phases are modeled as isotropic elastic bodies, and the interface structure is described by a constant surface tension. Existence of solutions is proved, and their asymptotic behavior in time is studied, with particular attention to the competition between surface tension and bulk deformation.     

Thin interphase/imperfect interface in elasticity with application to coated fiber composites

The imperfect interface conditions which are equivalent to the effect of a thin elastic interphase are derived by a Taylor expansion method in terms of interface displacement and traction jumps. Plane and cylindrical interfaces are analyzed as special cases. The effective elastic moduli of a unidirectional coated fiber composite are obtained on the basis of the derived imperfect interface conditions. High accuracy of the method is demonstrated by comparison of solutions of several problems in terms of the imperfect interface conditions or explicit presence of interphase as a third phase. The problems considered are transverse shear of a coated infinite fiber in infinite matrix and effective transverse bulk and shear moduli and effective axial shear modulus of a coated fiber composite. Unlike previous elastic imperfect interface conditions in the literature, the present ones are valid for the entire range of interphase stiffness, from very small to very large.     

Equilibria, Connecting Orbits and a Priori Bounds for Semilinear Parabolic Equations with Nonlinear Boundary Conditions

We consider a semilinear parabolic equation with a nonlinear non-dissipative boundary condition. In the one-dimensional case we describe bifurcation diagrams for positive and sign-changing equilibria and connecting orbits between these equilibria. We also show that the number of sign-changing stationary solutions strongly depends on the spatial dimension. The results are based on new a priori estimates of global solutions.     

The role of constraint in the fields of a crack normal to the interface between elastically and plastically mismatched solids

Asymptotic solutions are presented for a stationary crack normal to the boundary between two elastically mismatched solids such that the crack tip is located at the interface. The second-order term in the elastic asymptotic expansion was determined as a function of elastic mismatch for a thin cracked film on a substrate and for a thin cracked lamina between two substrates. Elastic-plastic analysis was performed using both modified boundary layer formulations and full field analyses. Analytic and numerical solutions in small strain yielding identify elastic mismatch and the T-stress as the determinants of crack tip constraint. The effect of constraint on the competition between interface failure and penetration is discussed.     

Damage cascade in a softening interface

A model describing the damage at an interface which is coupled to an elastic homogeneous block is introduced. Resorting to a real-space renormalization analysis, we show that in the absence of heterogeneity localization proceeds through a cascade of bifurcations which progressively concentrates the damage from the global interface to a narrow region leading to a crack nucleation. The equivalent homogeneous interface behaviour is obtained through this entire cascade, allowing for the analysis of size effects. When random heterogeneities are introduced in the interface, prior to the onset of localization damage proceeds by a sequence of avalanches whose mean size diverges at the first bifurcation point of the homogeneous interface. The large scale features of the bifurcation cascade are preserved, while the details of the late stage are smeared out by the randomness.     

The “resultant bifurcation diagram” method and its application to bifurcation behaviors of a symmetric railway bogie system

The concept of symmetric bifurcation for a symmetric wheel-rail system is defined. After that, the time response of the system can be achieved by the numerical integration method, and an unfixed and dynamic Poincaré section and its symmetric section for the symmetric wheel-rail system are established. Then the ??resultant bifurcation diagram?? method is constructed. The method is used to study the symmetric/asymmetric bifurcation behaviors and chaotic motions of a two-axle railway bogie running on an ideal straight and perfect track, and a variety of characteristics and dynamic processes can be obtained in the results. It is indicated that, for the possible sub-critical Hopf bifurcation in the railway bogie system, the stable stationary solutions and the stable periodic solutions coexist. When the speed is in the speed range of Hopf bifurcation point and saddle-node bifurcation point, the coexistence of multiple solutions can cause the oscillating amplitude change for different kinds of disturbance. Furthermore, it is found that there are symmetric motions for lower speeds, and then the system passes to the asymmetric ones for wide ranges of the speed, and returns again to the symmetric motions with narrow speed ranges. The rule of symmetry breaking in the system is through a blue sky catastrophe in the beginning.     

Effects of material anisotropy and inhomogeneity on cavitation for composite incompressible anisotropic nonlinearly elastic spheres

The effects of material anisotropy and inhomogeneity on void nucleation and growth in incompressible anisotropic nonlinearly elastic solids are examined. A bifurcation problem is considered for a composite sphere composed of two arbitrary homogeneous incompressible nonlinearly elastic materials which are transversely isotropic about the radial direction, and perfectly bonded across a spherical interface. Under a uniform radial tensile dead-load, a branch of radially symmetric configurations involving a traction-free internal cavity bifurcates from the undeformed configuration at sufficiently large loads. Several types of bifurcation are found to occur. Explicit conditions determining the type of bifurcation are established for the general transversely isotropic composite sphere. In particular, if each phase is described by an explicit material model which may be viewed as a generalization of the classic neo-Hookean model to anisotropic materials, phenomena which were not observed for the homogeneous anisotropic sphere nor for the composite neo-Hookean sphere may occur. The stress distribution as well as the possible role of cavitation in preventing interface debonding are also examined for the general composite sphere.     

Homoclinic bifurcation at resonant eigenvalues

We consider a bifurcation of homoclinic orbits, which is an analogue of period doubling in the limit of infinite period. This bifurcation can occur in generic two parameter vector fields when a homoclinic orbit is attached to a stationary point with resonant eigenvalues. The resonance condition requires the eigenvalues with positive/negative real part closest to zero to be real, simple, and equidistant to zero. Under an additional global twist condition, an exponentially flat bifurcation of double homoclinic orbits from the primary homoclinic branch is established rigorously. Moreover, associated period doublings of periodic orbits with almost infinite period are detected. If the global twist condition is violated, a resonant side switching occurs. This corresponds to an exponentially flat bifurcation of periodic saddle-node orbits from the homoclinic branch.Partially supported by DARPA and NSF.Partially supported by the Deutsche Forschungsgemeinschaft and by Konrad-Zuse-Zentrum für Informationstechnik Berlin.     

Magma flow through elastic-walled dikes

A convection–diffusion model for the averaged flow of a viscous, incompressible magma through an elastic medium is considered. The magma flows through a dike from a magma reservoir to the Earth’s surface; only changes in dike width and velocity over large vertical length scales relative to the characteristic dike width are considered. The model emerges when nonlinear inertia terms in the momentum equation are neglected in a viscous, low-speed approximation of a magma flow model coupled to the elastic response of the rock.Stationary- and traveling-wave solutions are presented in which a Dirichlet condition is used at the magma chamber; and either a (i) free-boundary condition, (ii) Dirichlet condition, or (iii) choked-flow condition is used at the moving free or fixed-top boundary. A choked-flow boundary condition, generally used in the coupled elastic wave and magma flow model, is also used in the convection–diffusion model. The validity of this choked-flow condition is illustrated by comparing stationary flow solutions of the convection–diffusion and coupled elastic wave and magma flow model for parameter values estimated for the Tolbachik volcano region in Kamchatka, Russia. These free- and fixed-boundary solutions are subsequently explored in a conservative, local discontinuous Galerkin finite-element discretization. This method is advantageous for the accurate implementation of the choked flow and free-boundary conditions. It uses a mixed Eulerian–Lagrangian finite element with special infinite curvature basis function near the free boundary and ensures positivity of the mean aperture subject to a time-step restriction. We illustrate the model further by simulating magma flow through host rock of variable density, and magma flow that is quasi-periodic due to the growth and collapse of a lava dome.     

粘弹性界面裂纹奇异场

对于许多粘弹性问题,通常可以利用对应性原理,即由弹性问题的结果得到对应的粘弹性问题在拉普拉斯变换域内的解,再通过反演变换求得最终时域中的解.但是,由于界面裂纹场存在着振荡奇异性,弹性问题解的形式就已经非常复杂,对应的粘弹性问题要通过反演变换直接求得准确的解析解几乎是不可能的.本文在利用对应性原理时做了更简单的准静态处理,即将弹性结果中的材料参数用粘弹性材料参数做对应替代,得到了粘弹性界面裂纹场近似的经典解,并与有限元分析结果作了比较.同时,利用Comninou接触模型,对粘弹性界面裂纹在远场拉剪混合加载情况下的裂尖应力场和接触区做了考察,并与经典解作了比较.     

A Finite Strain Analysis of Cavity Formation at a Rigid Inhomogeneity

The formation of a cavity by inclusion-matrix interfacial separation is examined by analyzing the response of a plane rigid inclusion embedded in an unbounded incompressible matrix subject to remote equibiaxial dead load traction. A vanishingly thin interfacial cohesive zone, characterized by normal and tangential interface force-separation constitutive relations, is assumed to govern separation behavior. Rotationally symmetric cavity shapes (circles) are shown to be solutions of an interfacial integral equation depending on the strain energy density of the matrix, the interface force constitutive relation and the remote loading. Nonsymmetrical cavity formation, under rotationally symmetric conditions of geometry and loading, is treated within the theory of infinitesimal strain superimposed on a given finite strain state. Rotationally symmetric and nonsymmetric bifurcations are analyzed and detailed results, for the Mooney–Rivlin strain energy density and for an exponential interface force-separation law, are presented. For the nonsymmetric rigid body displacement mode, a simple formula for the critical load is presented. The effect on bifurcation behavior of interfacial shear stiffness and other interface parameters is treated as well. In particular we demonstrate that (i) for the smooth interface nonsymmetric bifurcation always precedes rotationally symmetric bifurcation, (ii) unlike rotationally symmetric bifurcation, there is no threshold value of interface parameter for which nonsymmetric bifurcation will not occur and (iii) interfacial shear may significantly delay the onset of nonsymmetric bifurcation. Also discussed is the range of validity of a nonlinear infinitesimal strain theory previously presented by the author (Levy [1]). This revised version was published online in July 2006 with corrections to the Cover Date.     

Mechanics and fracture of crack tip deformable bi-material interface

Novel interface deformable bi-layer beam theory is developed to account for local effects at crack tip of bi-material interface by modeling a bi-layer composite beam as two separate shear deformable sub-layers with consideration of crack tip deformation. Unlike the sub-layer model in the literature in which the crack tip deformations under the interface peel and shear stresses are ignored and thus a “rigid” joint is used, the present study introduces two interface compliances to account for the effect of interface stresses on the crack tip deformation which is referred to as the elastic foundation effect; thus a flexible condition along the interface is considered. Closed-form solutions of resultant forces, deformations, and interface stresses are obtained for each sub-layer in the bi-layer beam, of which the local effects at the crack tip are demonstrated. In this study, an elastic deformable crack tip model is presented for the first time which can improve the split beam solution. The present model is in excellent agreements with analytical 2-D continuum solutions and finite element analyses. The resulting crack tip rotation is then used to calculate the energy release rate (ERR) and stress intensity factor (SIF) of interface fracture in bi-layer materials. Explicit closed-form solutions for ERR and SIF are obtained for which both the transverse shear and crack tip deformation effects are accounted. Compared to the full continuum elasticity analysis, such as finite element analysis, the present solutions are much explicit, more applicable, while comparable in accuracy. Further, the concept of deformable crack tip model can be applied to other bi-layer beam analyses (e.g., delamination buckling and vibration, etc.).     

Well-Posedness, Instabilities, and Bifurcation Results for the Flow in a Rotating Hele�CShaw Cell

We study the radial movement of an incompressible fluid located in a Hele–Shaw cell rotating at a constant angular velocity in the horizontal plane. Within an analytic framework, local existence and uniqueness of solutions is proved, and it is shown that the unique rotationally invariant equilibrium of the flow is unstable. There are, however, other time-independent solutions: using surface tension as a bifurcation parameter we establish the existence of global bifurcation branches consisting of stationary fingering patterns. The same results can be obtained by fixing the surface tension while varying the angular velocity. Finally, it is shown that the equilibria on a global bifurcation branch converge to a circle as the surface tension tends to infinity, provided they stay suitably bounded.     

Modelling of progressive interface failure under combined normal compression and shear stress

The present work is concerned with an analysis of progressive interface failure under normal compressive stress and varying shear stress using the cohesive crack model. The softening model is assumed and frictional linear stress at contact is accounted for. A monotonic loading in anti-plane shear of an elastic plate bonded to a rigid substrate is considered. An analytical solution is obtained by neglecting the effect of minor shear stress component in the plate. The elastic and plate interface compliances are included into the analysis. Three types of solutions are distinguished in the progressive delamination analysis, namely short, medium and long plate solutions. The analysis of quasi-static progressive delamination process clarifies the character of critical points and post-critical response of the plate. The analytical solution provides a reference benchmark for numerical algorithms of analysis of progressive interface delamination. The case of a rigid–softening interface was treated in a companion paper, where also cyclic loading was considered.     

Pulse Dynamics in a Three-Component System: Existence Analysis

In this article, we analyze the three-component reaction-diffusion system originally developed by Schenk et al. (PRL 78:3781–3784, 1997). The system consists of bistable activator-inhibitor equations with an additional inhibitor that diffuses more rapidly than the standard inhibitor (or recovery variable). It has been used by several authors as a prototype three-component system that generates rich pulse dynamics and interactions, and this richness is the main motivation for the analysis we present. We demonstrate the existence of stationary one-pulse and two-pulse solutions, and travelling one-pulse solutions, on the real line, and we determine the parameter regimes in which they exist. Also, for one-pulse solutions, we analyze various bifurcations, including the saddle-node bifurcation in which they are created, as well as the bifurcation from a stationary to a travelling pulse, which we show can be either subcritical or supercritical. For two-pulse solutions, we show that the third component is essential, since the reduced bistable two-component system does not support them. We also analyze the saddle-node bifurcation in which two-pulse solutions are created. The analytical method used to construct all of these pulse solutions is geometric singular perturbation theory, which allows us to show that these solutions lie in the transverse intersections of invariant manifolds in the phase space of the associated six-dimensional travelling wave system. Finally, as we illustrate with numerical simulations, these solutions form the backbone of the rich pulse dynamics this system exhibits, including pulse replication, pulse annihilation, breathing pulses, and pulse scattering, among others.      

Interface damage modeled by spring boundary conditions for in-plane elastic waves

In-plane elastic wave propagation in the presence of a damaged interface is investigated. The damage is modeled as a distribution of small cracks and this is transformed into a spring boundary condition. First the scattering by a single interface crack is determined explicitly in the low frequency limit for the case of a plane wave normally incident to the interface. The transmission at an interface with a random distribution of small cracks is then determined and is compared to periodically distributed cracks. The cracked interface is then described by a distributed spring boundary condition. As an illustration the dispersion relation of the first modes in a thick plate with a damaged interface in the middle is given.     

On the bifurcations of the Lamé solutions in plane-strain elasticity

We consider the in-plane bifurcations experienced by the Lamé solutions corresponding to an elastic annulus subjected to radial tension on the curved boundaries. Numerical investigations of the relevant incremental problem reveal two main bifurcation modes: a long-wave local deformation around the central hole of the domain, or a material wrinkling-type instability along the same boundary. Strictly speaking, the latter scenario is related to the violation of the Shapiro–Lopatinskij condition in an appropriate traction boundary-value problem. It is further shown that the main features of this material instability mode can be found by using a singular-perturbation strategy.