Fundamental elastic field in an infinite plane of two-dimensional piezoelectric quasicrystal subjected to multi-physics loads

By utilizing the extended Stroh formalism, the Green's function of infinite plane is obtained for the problem of two-dimensional decagonal quasicrystals with the piezoelectric effect subjected to multi-physics loads. By numerical computations, the piezoelectric effect of the two-dimensional decagonal quasicrystals is revealed; the changes of the stress and displacement fields with multi-physics loads are discussed. The variation laws of material constants in stress and displacement fields are investigated. The results show that the effect of the phason field on the generalized displacement is larger than that on the generalized stress; and the effects of material parameters are different in diverse field.     

Application of a Viscoelastic Model to simulate the Interface of a Metal/Fibre-Reinforced Polymer Joint

A viscoelastic model with the Lemaitre-type damage is applied to simulate an interfacial adhesive zone in light weight engineering structures, like aluminum/fiber-reinforced polymer specimens. The evolution of irreversible deformation and damage progression are investigated. The joint of aluminium alloy 5754 (AA5754) and carbon fibre reinforced thermoplastic composite CF-PA66 is manufactured by means of an epoxy (1K) adhesive. The adhesive zone is considered as an interface material. The aim of the research is to study the influence of the interface geometry on the mechanical characteristics of the structure. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stress and free vibration analysis of functionally graded beams using static Green's functions

In this paper static Green's functions for functionally graded Euler-Bernoulli and Timoshenko beams are presented. All material properties are arbitrary functions along the beam thickness direction. The closed-form solutions of static Green's functions are derived from a fourth-order partial differential equation presented in [2]. In combination with Betti's reciprocal theorem the Green's functions are applied to calculate internal forces and stress analysis of functionally graded beams (FGBs) under static loadings. For symmetrical material properties along the beam thickness direction and symmetric cross-sections, the resulting stress distributions are also symmetric. For unsymmetrical material properties the neutral axis and the center of gravity axis are located at different positions. Free vibrations of functionally graded Timoshenko beams are also analyzed [3]. Analytical solutions of eigenfunctions and eigenfrequencies in closed-forms are obtained based on reference [2]. Alternatively it is also possible to use static Green's functions and Fredholm's integral equations to obtain approximate eigenfunctions and eigenfrequencies by an iterative procedure as shown in [1]. Applying the Sensitivity Analysis with Green's Functions (SAGF) [1] to derive closed-form analytical solutions of functionally graded beams, it is possible to modify the derived static Green's functions and include terms taking cracks into account, which are modeled by translational or rotational springs. Furthermore the SAGF approach in combination with the superposition principle can be used to take stiffness jumps into account and to extend static Green's functions of simple beams to that of discontinuous beams by adding new supports. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Linking macroscopic deformation processes to microstructure evolution using an FE-FFT-based micro-macro transition and non-conserved phase-fields

The purpose of this work is the multiscale FE-FFT-based prediction of macroscopic material behavior, micromechanical fields and bulk microstructure evolution in polycrystalline materials subjected to macroscopic mechanical loading. The macroscopic boundary value problem (BVP) is solved using implicit finite element (FE) methods. In each macroscopic integration point, the microscopic BVP is embedded, the solution of which is found employing fast Fourier transform (FFT), fixed-point and Green's function methods. The mean material response is determined by the stress-strain relation at the micro scale or rather the volume average of the micromechanical fields. The evolution of the microstructure is modeled by means of non-conserved phase-fields. As an example, the proposed methodology is applied to the modeling of stress-induced martensitic phase transformations in metal alloys. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Damageability of a material reinforced with unidirectional orthotropic fibers for an exponential function of long-term microstrength

The theory of long-term damageability of a homogeneous material is generalized to the case of an orthotropic fibrous composite material with a stochastic structure. Equations of mechanics of microinhomogeneous media of this structure form the base of the theory. The process of damage of components of a composite is modeled by the formation of stochastically located micropores. The criterion of fracture of a unit microvolume is characterized by its long-term strength determined by the dependence of the time of brittle fracture on the degree of closeness of the equivalent stress to its limit value, which characterizes the short-term strength on the basis of the Huber–von Mises criterion accepted as an arbitrary function of coordinates. Efficient deformation properties and the stress-strain state of an orthotropic fibrous composite with microdamages in components are determined on the base of stochastic equations of elasticity of orthotropic media. For given macrostresses and macrostrains and an arbitrary moment of time, balance equations of damage (porosity) of components are formulated. On the basis of the iteration method, we construct algorithms for calculating dependences of microdamage of components of an orthotropic fibrous material on time and dependences of macrostresses or macrostrains on time and obtain the corresponding curves for the case of a bounded function of the long-term microstrength, which is approximated by an exponential law.     

Inplane Stress Singularities at the Corner of an Anisotropic Material Discontinuity

In the mechanics of composite laminates the local mechanical inplane fields at corners of anisotropic material discontinuities are of particular interest since they can have singular behavior. In the present study, the stress and strain fields in the local near field of such corners are investigated by an asymptotic analysis. The order of the singularity of these mechanical inplane fields are determined in closed‐form manner by use of the complex potential method based on Lekhnitskii's approach. Various different geometrical setups and material combinations of corners with material discontinuities are investigated with regard to their effect on the singular behavior of the mechanical fields present. These examples show that the order of singularity considered is clearly weaker than the typical crack tip singularity in fracture mechanics. Nevertheless, it may render the corner a critical location for the onset of failure. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Shear of orthogonally reinforced spherofibrous composites

The problem of determining the shear characteristics and interphase stress concentration of fibrous composites with spherical inclusions is examined on the basis of a three-phase model. Stress fields caused by diffusion interaction of phases are neglected. The elastic moduli of the composite are investigated and compared with those obtained from a two-phase model. The general formula for determination of the shear modulus of triorthogonally reinforced compsites is derived using previously investigated relationships for averaged stress fields. The matrix of these compsites contained spherical cavities. The dependence of integral characteristics of three-phase composites on their bulk phase concentration was investigated. The stresses between phases were studied as a function of composite structure.A. A. Blagonravov Machine-Science Institute, Russian Academy of Sciences, Moscow, Russia. Translated from Mekhanika Kompozitnykh Materialov, No. 1, 104–111, January–February, 1997.     

Dynamic modeling for rotating composite Timoshenko beam and analysis on its bending-torsion coupled vibration

A formulation is presented for steady-state dynamic responses of rotating bending-torsion coupled composite Timoshenko beams (CTBs) subjected to distributed and/or concentrated harmonic loadings. The separation of cross section's mass center from its shear center and the introduced coupled rigidity of composite material lead to the bending-torsion coupled vibration of the beams. Considering those two coupling factors and based on Hamilton's principle, three partial differential non-homogeneous governing equations of vibration with arbitrary boundary conditions are formulated in terms of the flexural translation, torsional rotation and angle rotation of cross section of the beams. The parameters for the damping, axial load, shear deformation, rotation speed, hub radius and so forth are incorporated into those equations of motion. Subsequently, the Green's function element method (GFEM) is developed to solve these equations in matrix form, and the analytical Green's functions of the beams are given in terms of piecewise functions. Using the superposition principle, the explicit expressions of dynamic responses of the beams under various harmonic loadings are obtained. The present solving procedure for Timoshenko beams can be degenerated to deal with for Rayleigh and Euler beams by specifying the values of shear rigidity and rotational inertia. Cantilevers with bending-torsion coupled vibration are given as examples to verify the present theory and to illustrate the use of the present formulation. The influences of rotation speed, bending-torsion couplings and damping on the natural frequencies and/or shape functions of the beams are performed. The steady-state responses of the beam subjected to external harmonic excitation are given through numerical simulations. Remarkably, the symmetric property of the Green's functions is maintained for rotating bending-torsion coupled CTBs, but there will be a slight deviation in the numerical calculations.     

Static and dynamic analysis of cracked or weakened structures

Stiffness modifications in engineering structures, for example due to damage and cracking, will inevitably also lead to changes in deformations, internal forces, natural frequencies and mode shapes of the structures. In this paper, an efficient and simple method for sensitivity analysis of cracked or weakened structures under time-harmonic loading is presented. The method is based on a comparison between the strain energy and the kinetic energy of an uncracked structure and that of a cracked structure in conjunction with the application of exact or approximate Green's functions as described in [3] for the static case. The present analysis enables the prediction of any changes in the displacements and stresses and has a lower computational effort as compared to available classical methods, because only the damaged region has to be re-considered in the method. Green's functions are taken as a basis of the approach, which have the ability to weight the influence of the stiffness modifications in a region of a structure and show how sensitive other regions respond to the stiffness modifications. Based on linear elastic fracture mechanics, cracked or damaged regions are approximated by spring models in the analytical solution of some simple beam problems, while cracked finite elements are used for complicated cases where analytical solutions cannot be obtained. Sensitivity analysis with Green's functions (SAGF) approach is applied to static and dynamic analysis of cracked and weakened structures, which consist of homogeneous materials or fiber reinforced composites like reinforced concretes. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Approximate Separability of the Green's Function of the Helmholtz Equation in the High Frequency Limit

The minimum number of terms that are needed in a separable approximation for a Green's function reveals the intrinsic complexity of the solution space of the underlying differential equation. It also has implications for whether low‐rank structures exist in the linear system after numerical discretization. The Green's function for a coercive elliptic differential operator in divergence form was shown to be highly separable [2], and efficient numerical algorithms exploiting low‐rank structures of the discretized systems were developed. In this work, a new approach to study the approximate separability of the Green's function of the Helmholtz equation in the high‐frequency limit is developed. We show (1) lower bounds based on an explicit characterization of the correlation between two Green's functions and a tight dimension estimate for the best linear subspace to approximate a set of decorrelated Green's functions, (2) upper bounds based on constructing specific separable approximations, and (3) sharpness of these bounds for a few case studies of practical interest. © 2018 Wiley Periodicals, Inc.     

Wavelet compression of integral operators arising from boudary-domain integral method

The boundary-domain integral method uses Green's functions to write integral representations of partial differential equations. Since Green's functions are non-local, the systems of linear equations arising from the discretization of integral representations are fully populated. Several algorithms have been proposed, which yield a data-sparse approximation of these systems, such as the fast multipole method, adaptive cross approximation and among others, wavelet compression. In the framework of solving the Navier-Stokes equations in velocity-vorticity form one may utilize the boundary-domain integral method for the solution of the kinematics equation to calculate the boundary vorticity values. Since the kinematics equation is a Poisson type equation, usually its integral representation is written with the Green's function for the Laplace operator. In this work, we introduce a false time into the equation and parabolize its nature. Thus, a time-dependent Green's function may be used. This provides a new parameter, the time step, which can be set to control the Green's function. The time-dependent Green's function is a global function, but by carefully choosing the time step, its behaviour is almost local. This makes it a good candidate for wavelet compression, yielding much better compression ratios at a given accuracy than when using the Green's function for the Laplace operator. However, as false time is introduced, several time steps must be solved in order to reach a converged solution. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Simulation of Uncertain Inelastic Material Behaviour using the Stochastic Finite Element Method

In this work, an approach for computing three-dimensional structures with random material properties, such as the yield stress, Young's modulus and hardening parameters is proposed. The random material properties are represented as random fields which are realized with the Spectral Representation Method (SPRM). The proposed approach is coupled with Monte Carlo Simulation (MCS) to determine the response statistics of a simple mechanical structure. The numerical results are compared with those obtained from classical Latin Hypercube Sampling (LHS). (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Modelling of Lamb wave propagation in composite plate excited by surface bonded piezoelectrical actuators

Transient propagation of Lamb waves in a multilayered infinite composite plate with the arbitrary elastic anisotropy of each layer is considered in this paper. Integral representations of wave fields, an algorithm of computing the Green's matrix of the multilayered medium in Fourier domain and some methods of the Fourier transform inversion are used as basic analytical research instruments. The piezoelectrical transducers frequently used for Lamb wave excitation are modelled in this paper as surface sources with stresses concentrated on the boundaries of actuators. The contribution of each wave mode to a whole amplitude, focusing of Lamb waves and influence of the geometry of the surface source are discussed in this paper by numerical examples. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A closed‐form analysis of material and geometry effects on stress singularities at unsymmetric bimaterial notches

An important issue in the mechanics of adhesive bonds is the knowledge of local mechanical fields. In the present study, an asymptotic analysis of the stress fields near an unsymmetric bimaterial notch with arbitrary opening angle is performed. Using the complex potential method, the order of the singularity of the stress fields at a notch tip can be determined in closed‐form analytical manner, so that the dependency of the occurring singularity exponents on geometry and material properties can be studied systematically.     

Green's functions for elliptic and parabolic systems with Robin-type boundary conditions

The aim of this paper is to investigate Green's function for parabolic and elliptic systems satisfying a possibly nonlocal Robin-type boundary condition. We construct Green's function for parabolic systems with time-dependent coefficients satisfying a possibly nonlocal Robin-type boundary condition assuming that weak solutions of the system are locally Hölder continuous in the interior of the domain, and as a corollary we construct Green's function for elliptic system with a Robin-type condition. Also, we obtain Gaussian bound for Robin Green's function under an additional assumption that weak solutions of Robin problem are locally bounded up to the boundary. We provide some examples satisfying such a local boundedness property, and thus have Gaussian bounds for their Green's functions.     

Multiscale Modeling of the Interphase Interactions in the Mechanics of Materials

The consistent and correct model of media taking into account scale effects (cohesion and adhesive interactions) is constructed as a special case of the Cosserat's pseudocontinuum model. The variant of the interphase layer theory is elaborated, which includes the following moments: formal mathematical statement, physical constitutive equations, numerical estimations of an interphase layer influence on the stress state and energy density distribution in a composite. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Minimization principle and generalized Fourier series for discontinuous Sturm-Liouville systems in direct sum spaces

By modifing the Green''s function method we study certain spectral aspects of discontinuous Sturm-Liouville problems with interior singularities. Firstly, we define four eigen-solutions and construct the Green''s function in terms of them. Based on the Green''s function we establish the uniform convergeness of generalized Fourier series as eigenfunction expansion in the direct sum of Lebesgue spaces $L_2$ where the usual inner product replaced by new inner product. Finally, we extend and generalize such important spectral properties as Parseval equation, Rayleigh quotient and Rayleigh-Ritz formula (minimization principle) for the considered problem.     

A modified Green's function for the internal gravitational wave equation in a layer of a stratified medium with a constant shear flow

The construction of a modified Green's function for the internal gravitational wave (IGW) equation in a layer of a stratified medium when there are constant mean shear flows is considered and the basic properties of the corresponding eigenvalue problems and the modified eigenfunctions and eigenvalues are investigated. It is shown that each mode of the modified Green's function consists of a sum of three terms describing (1) the IGWs that propagate from the source, (2) the effects of a time varying source, localized in a certain neighbourhood of it, and (3) the effects of the displacement of the fluid (an internal discontinuity) caused by the source. The resulting expressions are analysed out for a constant and oscillating source of the generation of IGWs in which each of the terms of Green's function is represented in the form of simple quadratures.     

Effective bulk moduli of materials containing stochastically distributed nano-inhomogeneities with surface stresses

In the present contribution, a mathematical model for the investigation of the effective properties of a material with randomly distributed nano-particles is proposed. The surface effect is introduced via Gurtin-Murdoch equations describing properties of the matrix/nano-particle interface. They are added to the system of stochastic differential equations formulated within the framework of linear elasticity. The homogenization problem is reduced to finding a statistically averaged solution of the system of stochastic differential equations. These equations are based on the fundamental equations of linear elasticity, which are coupled with surface/interface elasticity accounting for the presence of surface tension. Using Green's function this system is transformed to a system of statistically non-linear integral equations. It is solved by the method of conditional moments. Closed-form expressions are derived for the effective moduli of a composite consisting of a matrix with randomly distributed spherical inhomogeneities. The radius of the nano-particles is included in the expression for the bulk moduli. As numerical examples, nano-porous aluminum and nano-porous gold are investigated assuming that only the influence of the interface effects on the effective bulk modulus is of interest. The dependence of the normalized bulk moduli of nano-porous aluminum on the pore volume fraction (for certain radii of nano-pores) are compared to and discussed in the context of other theoretical predictions. The effective Young's modulus of nano-porous gold as a function of pore radius (for fixed void volume fraction) and the normalized Young's modulus vs. the pore volume fraction for different pore radii are analyzed. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)