MECHANICAL CONDITIONS FOR FRACTURE OF ANISOTROPIC BODIES AND THEIR GEOMETRY IN STRESS SPACE

The conditions for fracture of anisotropic bodies and theirgeometry in stress space are proposed in this paper.The an-alytical formulae expressing the fracture conditions are es-tablished from the viewpoint of energy theory for crack pro-pagation.In stress space the limiting surface corresponding to thefracture conditions derived for anisotropic solids is quadra-tic.It is an ellipsoid in case the mean stress is greaterthan zero and it is hyperboloid in case the mean stress is smal-ler than zero.The conclusions formed by the author in the present paperhave certain generality.Some results obtained by predeces-sors appear to be special cases with respect to the presenttheory.     

Three-to-one internal resonance in multiple stepped beam systems

In this study, the vibrations of multiple stepped beams with cubic nonlinearities are considered. A three-to-one internal resonance case is investigated for the system. A general approximate solution to the problem is found using the method of multiple scales （a perturbation technique）. The modulation equations of the amplitudes and the phases are derived for two modes. These equations are utilized to determine steady state solutions and their stabilities. It is assumed that the external forcing frequency is close to the lower frequency. For the numeric part of the study, the three-to-one ratio in natural frequencies is investigated. These values are observed to be between the first and second natural frequencies in the cases of the clamped-clamped and clamped-pinned supports, and between the second and third natural frequencies in the case of the pinned-pinned support. Finally, a numeric algorithm is used to solve the three-to-one internal resonance. The first mode is externally excited for the clamped-clamped and clamped-pinned supports, and the second mode is externally excited for the pinned-pinned support. Then, the amplitudes of the first and second modes are investigated when the first mode is externally excited. The amplitudes of the second and third modes are investigated when the second mode is externally excited. The force-response, damping-response, and .frequency- response curves are plotted for the internal resonance modes of vibrations. The stability analysis is carried out for these plots.     

A continuous-discontinuous cellular automaton method for cracks growth and coalescence in brittle material简

A method of continuous-discontinuous cellular automaton for modeling the growth and coalescence of multiple cracks in brittle material is presented. The method uses the level set to track arbitrary discontinuities, and calculation grids are independent of the discontinuities and no remeshing are required with the crack growing. Based on Grif- fith fracture theory and Mohr-Coulumb criterion, a mixed fracture criterion for multiple cracks growth in brittle mate- rial is proposed. The method treats the junction and coales- cence of multiple cracks, and junction criterion and coales- cence criterion for brittle material are given, too. Besides, in order to overcome the tracking error in the level set ap- proximation for crack junction and coalescence, a dichotomy searching algorithm is proposed. Introduced the above the- ories into continuous-discontinuous cellular automaton, the present method can be applied to solving multiple crack growth in brittle material, and only cell stiffness is needed and no assembled global stiffness is needed. Some numerical examples are given to shown that the present method is efficient and accurate for crack junction, coalescence and percolation problems.     

Optimal placement of piezoelectric active bars in vibration control by topological optimization

A continuous variable optimization method and a topological optimization method are proposed for the vibration control of piezoelectric truss structures by means of the optimal placements of active bars. In this optimization model, a zero-one discrete variable is defined in order to solve the optimal placement of piezoelectric active bars. At the same time, the feedback gains are also optimized as continuous design variables. A two-phase procedure is proposed to solve the optimization problem. The sequential linear programming algorithm is used to solve optimization problem and the sensitivity analysis is carried out for objective and constraint functions to make linear approximations. On the basis of the Newmark time integration of structural transient dynamic responses, a new sensitivity analysis method is developed in this paper for the vibration control problem of piezoelectric truss structures with respect to various kinds of design variables. Numerical examples are given in the paper to demonstrate the effectiveness of the methods.     

Vibrations of FGM thin cylindrical shells with exponential volume fraction law

In this paper, the influence of an exponential volume fraction law on the vibration frequencies of thin functionally graded cylindrical shells is studied. Material properties in the shell thickness direction are graded in accordance with the exponential law. Expressions for the strain-displacement and curvature-displacement relationships are taken from Love＇s thin shell theory. The Rayleigh-Ritz approach is used to derive the shell eigenfrequency equation. Axial modal dependence is assumed in the characteristic beam functions. Natural frequencies of the shells are observed to be dependent on the constituent volume fractions. The results are compared with those available in the literature for the validity of the present methodology.     

Optimal convergence rates for three-dimensional turbulent flow equations

In this paper,the convergence rates of solutions to the three-dimensional turbulent flow equations are considered.By combining the Lp-Lq estimate for the linearized equations and an elaborate energy method,the convergence rates are obtained in various norms for the solution to the equilibrium state in the whole space when the initial perturbation of the equilibrium state is small in the H3-framework.More precisely,the optimal convergence rates of the solutions and their first-order derivatives in the L2-norm are obtained when the Lp-norm of the perturbation is bounded for some p ∈[1,6/5).     

Displacements and stresses in composite multi-layered media due to varying temperature and concentrated load

This paper deals with the determination of the thermo-elastic displacements and stresses in a multi-layered body set up in different layers of different thickness having different elastic properties due to the application of heat and a concentrated load in the uppermost surface of the medium. Each layer is assumed to be made of homogeneous and isotropic elastic material. The relevant displacement components for each layer are taken to be axisymmetric about a line, which is perpendicular to the plane surfaces of all layers. The stress function for each layer, therefore, satisfies a single equation in absence of any body forces. The equation is then solved by integral transform technique. Analytical expressions for thermo-elastic displacements and stresses in the underlying mass and the corresponding numerical codes are constructed for any number of layers. However, the numerical comparison is made for three and four layers.     

Stochastic Hopf bifurcation in quasi-integrable-Hamiltonian systems

A new procedure is developed to study the stochastic Hopf bifurcation in quasiintegrable-Hamiltonian systems under the Gaussian white noise excitation. Firstly, the singular boundaries of the first-class and their asymptotic stable conditions in probability are given for the averaged Ito differential equations about all the sub-system‘s energy levels with respect to the stochastic averaging method. Secondly, the stochastic Hopf bifurcation for the coupled sub-systems are discussed by defining a suitable bounded torus region in the space of the energy levels and employing the theory of the torus region when the singular boundaries turn into the unstable ones. Lastly, a quasi-integrable-Hamiltonian system with two degrees of freedom is studied in detail to illustrate the above procedure.Moreover, simulations by the Monte-Carlo method are performed for the illustrative example to verify the proposed procedure. It is shown that the attenuation motions and the stochastic Hopf bifurcation of two oscillators and the stochastic Hopf bifurcation of a single oscillator may occur in the system for some system‘s parameters. Therefore, one can see that the numerical results are consistent with the theoretical predictions.     

Frequency domain fundamental solutions for a poroelastic half-space

In frequency domain, the fundamental solutions for a poroelastic half-space are re-derived in the context of Biot＇s theory. Based on Biot＇s theory, the governing field equations for the dynamic poroelasicity are established in terms of solid displacement and pore pressure. A method of potentials in cylindrical coordinate system is proposed to decouple the homogeneous Biot＇s wave equations into four scalar Helmholtz equations, and the general solutions to these scalar wave equations are obtained. After that, spectral Green＇s functions for a poroelastic full-space are found through a decomposition of solid displacement, pore pressure, and body force fields. Mirror-image technique is then applied to construct the half-space fundamental solutions.Finally, transient responses of the half-space to buried point forces are examined.     

REFLECTION CHARACTERISTICS OF LONGITUDINAL WAVES IN A SEMI-INFINITE CYLINDRICAL ROD CONNECTED TO AN INFINITE VISCOELASTIC STRATUM

Reflection characteristics of longitudinal strain waves in a semi-infinite elastic rod con-nected to a viscoelastic stratum are investigated analytically.The three-dimensional viscoelasticity the-ory is applied to the stratum,and the Laplace transform with respect to time and the numerical inverseLaplace transform by means of Laguerre function are used.The time histories for the longitudinalstrain of an arbitrary point of the rod are presented.Two typical viscoelastic models are considered,one is the usual Maxwell-Voigt model,the other is whose relaxation function is given by a power law.The numerical results for the two models are presented and compared each other and also with previ-ously published results for the elastic stratum.     

Hybrid algorithm for modeling of fluid-structure interaction in incompressible, viscous flows

The objective of this paper is to present and to validate a new hybrid coupling (HC) algorithm for modeling of fluid-structure interaction (FSI) in incompressible, viscous flows. The HC algorithm is able to avoid numerical instability issues associated with artificial added mass effects, which are often encountered by standard loosely coupled (LC) and tightly coupled (TC) algorithms, when modeling the FSI response of flexible structures in incompressible flow. The artificial added mass effect is caused by the lag in exchange of interfacial displacements and forces between the fluid and solid solvers in partitioned algorithms. The artificial added mass effect is much more prominent for light/flexible structures moving in water, because the fluid forces are in the same order of magnitude as the solid forces, and because the speed at which numerical errors propagate in an incompressible fluid. The new HC algorithm avoids numerical instability issues associated with artificial added mass effects by embedding Theodorsen’s analytical approximation of the hydroelastic forces in the solution process to obtain better initial estimates of the displacements. Details of the new HC algorithm are presented. Numerical validation studies are shown for the forced pitching response of a steel and a plastic hydrofoil. The results show that the HC algorithm is able to converge faster, and is able to avoid numerical instability issues, compared to standard LC and TC algorithms, when modeling the transient FSI response of a plastic hydrofoil. Although the HC algorithm is only demonstrated for a NACA0009 hydrofoil subject to pure pitching motion, the method can be easily extended to model general 3-D FSI response and stability of complex, flexible structures in turbulent, incompressible, multiphase flows.     

Micromechanical analysis of the behavior of stiff clay

Cementations formed in geological timescale are observed in various stiff clays.A micromechanical stress strain model is developed for modeling the effect of cementation on the deformation behavior of stiff clay.The proposed approach considers explicitly cementations at intercluster contacts,which is different from conventional model.The concept of inter-cluster bonding is introduced to account for an additional cohesion in shear sliding and a higher yield stress in normal compression.A damage law for inter-cluster bonding is proposed at cluster contacts for the debonding process during mechanical loading.The model is used to simulate numerous stress-path tests on Vallericca stiff clay.The applicability of the present model is evaluated through comparisons between the predicted and the measured results.In order to explain the stress-induced anisotropy arising from externally applied load,the evolution of local stresses and local strains at inter-cluster planes are discussed.     

Study of modeling unsteady blade row interaction in a transonic compressor stage part 1: code development and deterministic correlation analysis

The average-passage equation system (APES) provides a rigorous mathematical framework for accounting for the unsteady blade row interaction through multistage compressors in steady state environment by introducing deterministic correlations (DC) that need to be modeled to close the equation system.The primary purpose of this study is to provide insight into the DC characteristics and the influence of DC on the time-averaged flow field of the APES.In Part 1 of this two-part paper,firstly a 3D viscous unsteady and time-averaging flow CFD solver is developed to investigate the APES technique.Then steady and unsteady simulations are conducted in a transonic compressor stage.The results from both simulations are compared to highlight the significance of the unsteady interactions.Furthermore,the distribution characteristics of DC are studied and the DC at the rotor/stator interface are compared with their spatial correlations (SC).Lastly,steady and time-averaging (employing APES with DC) simulations for the downstream stator alone are conducted employing DC derived from the unsteady results.The results from steady and time-averaging simulations are compared with the time-averaged unsteady results.The comparisons demonstrate that the simulation employing APES with DC can reproduce the time-averaged field and the 3D viscous time-averaging flow solver is validated.     

Effect of induced magnetic field on natural convection in vertical concentric annuli

In the present paper,we have considered the steady fully developed laminar natural convective flow in open ended vertical concentric annuli in the presence of a radial magnetic field.The induced magnetic field produced by the motion of an electrically conducting fluid is taken into account.The transport equations concerned with the considered model are first recast in the non-dimensional form and then unified analytical solutions for the velocity,induced magnetic field and temperature field are obtained for the cases of isothermal and constant heat flux on the inner cylinder of concentric annuli.The effects of the various physical parameters appearing into the model are demonstrated through graphs and tables.It is found that the magnitude of maximum value of the fluid velocity as well as induced magnetic field is greater in the case of isothermal condition compared with the constant heat flux case when the gap between the cylinders is less or equal to 1.70 times the radius of inner cylinder,while reverse trend occurs when the gap between the cylinders is greater than 1.71 times the radius of inner cylinder.These fields are almost the same when the gap between the cylinders is equal to 1.71 times the radius of inner cylinder for both the cases.It is also found that as the Hartmann number increases,there is a flattening tendency for both the velocity and the induced magnetic field.The influence of the induced magnetic field is to increase the velocity profiles.     

Indentation of a compressible soft electroactive half-space: Some theoretical aspects

We theoretically study the indentation response of a compressible soft electroactive material by a rigid punch. The half-space material is assumed to be initially subjected to a finite deformation and an electric biasing field. By adopting the linearized theory for incremental fields, which is established on the basis of a general nonlinear theory for electroelasticity, the appropriate equations governing the perturbed infinitesimal elastic and electric fields are derived particularly when the material is subjected to a uniform equibiaxial stretch and a uniform electric displacement. A general solution to the governing equations is presented, which is concisely expressed in terms of four quasi-harmonic functions. By adopting the potential theory method, exact contact solutions for three common perfectly conducting rigid indenters of flat-ended circular, conical and spherical geometries can be derived, and some explicit relations that are of practical importance are outlined.     

On the perturbation solution of interface-reaction controlled diffusion in solids

Insertion of species A into species B forms a product P through two kinetic processes, namely, (1) the chemical reaction between A and B that occurs at the B-P interface, and (2) the diffusion of species A in product P. These two processes are symbiotic in that the chemical reaction provides the driving force for the diffusion, while the diffusion sustains the chemical reaction by providing sufficient reactant to the reactive interface. In this paper, a mathematical framework is developed for the coupled reaction-diffusion processes. The resulting system of boundary and initial value problem is solved analytically for the case of interface-reaction controlled diffusion, i.e., the rate of diffusion is much faster than the rate of chemical reaction at the interface so that the final kinetics are limited by the interface chemical reaction. Asymptotic expressions are given for the velocity of the reactive interface and the concentration of diffusing species under two different boundary conditions.     

Assessment of composite failure and ultimate strength without experiment on composite

This paper attempts to estimate the ultimate strength of a laminated composite only based on its con- stituent properties measured independently. Three important issues involved have been systematically addressed, i.e., stress calculation for the constituent fiber and matrix materials, failure detection for the lamina and laminate upon the internal stresses in their constituents, and input data determination of the constituents from monolithic measurements. There are three important factors to influence the accuracy of the strength prediction. One is the stress concentration factor （SCF） in the matrix. Another is matrix plasticity. The third is thermal residual stresses in the constituents. It is these three factors, however, that have not been sufficiently well realized in the composite community. One can easily find out the elastic and strength parameters of a great many laminae and laminates in the current literature. Unfortunately, necessary information to determine the SCF, the matrix plasticity, and the thermal residual stresses of the composites is rare or incomplete. A useful design methodology is demonstrated in the paper.     

Analysis of axial strain in one-dimensional loading by different models

Different phenomenological equations based on plasticity, primary creep （as a viscoplastic mechanism）, secondary creep （as another viscoplastic mechanism） and different combinations of these equations are presented and used to describe the material inelastic deformation in uniaxial test. Agreement of the models with experimental results and with the theoretical concepts and physical realities is the criterion of choosing the most appropriate formulation for uniaxial test. A model is thus proposed in which plastic deformation, primary creep and secondary creep contribute to the inelastic deformation. However, it is believed that the hardening parameter is composed of plastic and primary creep parts. Accordingly, the axial plastic strain in a uniaxial test may no longer be considered as the hardening parameter. Therefore, a proportionality concept is proposed to calculate the plastic contribution of deformation.     

Stability and Hopf bifurcation of a delayed ratio-dependent predator-prey system

Since the ratio-dependent theory reflects the fact that predators must share and compete for food, it is suitable for describing the relationship between predators and their preys and has recently become a very important theory put forward by biologists. In order to investigate the dynamical relationship between predators and their preys, a so-called Michaelis-Menten ratio-dependent predator-prey model is studied in this paper with gestation time delays of predators and preys taken into consideration. The stability of the positive equilibrium is investigated by the Nyquist criteria, and the existence of the local Hopf bifurcation is analyzed by employing the theory of Hopf bifurcation. By means of the center manifold and the normal form theories, explicit formulae are derived to determine the stability, direction and other properties of bifurcating periodic solutions. The above theoretical results are validated by numerical simulations with the help of dynamical software WinPP. The results show that if both the gestation delays are small enough, their sizes will keep stable in the long run, but if the gestation delays of predators are big enough, their sizes will periodically fluc-tuate in the long term. In order to reveal the effects of time delays on the ratio-dependent predator-prey model, a ratiodependent predator-prey model without time delays is considered. By Hurwitz criteria, the local stability of positive equilibrium of this model is investigated. The conditions under which the positive equilibrium is locally asymptotically stable are obtained. By comparing the results with those of the model with time delays, it shows that the dynamical behaviors of ratio-dependent predator-prey model with time delays are more complicated. Under the same conditions, namely, with the same parameters, the stability of positive equilibrium of ratio-dependent predator-prey model would change due to the introduction of gestation time delays for predators and preys. Moreover, with the variation of time delays, the positive equilibrium of the ratio-dependent predator-prey model subjects to Hopf bifurcation.     

A class of exact solutions for the incompressible viscous magnetohydrodynamic flow over a porous rotating disko

The present paper is concerned with a class of exact solutions to the steady Navier-Stokes equations for the incompressible Newtonian viscous electrically conducting fluid flow due to a porous disk rotating with a constant angular speed.The three-dimensional hydromagnetic equations of motion are treated analytically to obtained exact solutions with the inclusion of suction and injection.The well-known thinning/thickening flow field effect of the suction/injection is better understood from the constructed closed form velocity equations.Making use of this solution,analytical formulas for the angular velocity components as well as for the permeable wall shear stresses are derived.Interaction of the resolved flow field with the surrounding temperature is further analyzed via the energy equation.The temperature field is shown to accord with the dissipation and the Joule heating.As a result,exact formulas are obtained for the temperature field which take different forms corresponding to the condition of suction or injection imposed on the wall.