Fractional WKB approximation

Wentzel–Kramer–Brillouin (WKB) approximation for fractional systems is investigated in this paper using the fractional calculus. In the fractional case, the wave function is constructed such that the phase factor is the same as the Hamilton’s principle function S. To demonstrate our proposed approach, two examples are investigated in detail.     

A hybrid technique for the solution of unsteady Maxwell fluid with fractional derivatives due to tangential shear stress

The article describes the unsteady motion of viscoelastic fluid for a Maxwell model with fractional derivatives. The flow is produced by cylinder, considering time dependent quadratic shear stress ft² on Maxwell fluid with fractional derivatives. The fractional calculus approach is used in the constitutive relationship of Maxwell model. By applying Laplace transform with respect to time t and modified Bessel functions, semianalytical solutions for velocity function and tangential shear stress are obtained. The obtained semianalytical results are presented in transform domain, satisfy both initial and boundary conditions. Our solutions particularized to Newtonian and Maxwell fluids having typical derivatives. The inverse Laplace transform has been calculated numerically. The numerical results for velocity function are shown in Table by using MATLAB program and compared them with two other algorithms in order to provide validation of obtained results. The influence of fractional parameters and material constants on the velocity field and tangential stress is analyzed by graphs.     

Strain gradient plasticity,strengthening effects and plastic limit analysis

Within the framework of isotropic strain gradient plasticity, a rate-independent constitutive model exhibiting size dependent hardening is formulated and discussed with particular concern to its strengthening behavior. The latter is modelled as a (fictitious) isotropic hardening featured by a potential which is a positively degree-one homogeneous function of the effective plastic strain and its gradient. This potential leads to a strengthening law in which the strengthening stress, i.e. the increase of the plastically undeformed material initial yield stress, is related to the effective plastic strain through a second order PDE and related higher order boundary conditions. The plasticity flow laws, with the role there played by the strengthening stress, are addressed and shown to admit a maximum dissipation principle. For an idealized elastic perfectly plastic material with strengthening effects, the plastic collapse load problem of a micro/nano scale structure is addressed and its basic features under the light of classical plastic limit analysis are pointed out. It is found that the conceptual framework of classical limit analysis, including the notion of rigid-plastic behavior, remains valid. The lower bound and upper bound theorems of classical limit analysis are extended to strengthening materials. A static-type maximum principle and a kinematic-type minimum principle, consequences of the lower and upper bound theorems, respectively, are each independently shown to solve the collapse load problem. These principles coincide with their respective classical counterparts in the case of simple material. Comparisons with existing theories are provided. An application of this nonclassical plastic limit analysis to a simple shear model is also presented, in which the plastic collapse load is shown to increase with the decreasing sample size (Hall–Petch size effects).     

A Fractional Model of Continuum Mechanics

Although there has been renewed interest in the use of fractional models in many application areas, in reality fractional analysis has a long and distinguished history and can be traced back to the likes of Leibniz (Letter to L’Hospital, 1695), Liouville (J. éc. Polytech. 13:71, 1832), and Riemann (Gesammelte Werke, p. 62, 1876). Recent publications (Podlubny in Math. Sci. Eng. 198, 1999; Sabatier et al. in Advances in fractional calculus: theoretical developments and applications in physics and engineering, Springer, Berlin, 2007; Das in Functional fractional calculus for system identification and controls, Springer, Berlin, 2007) demonstrate that fractional derivative models have found widespread applications in science and engineering. Late fundamental considerations have led to the introduction of fractional calculus in continuum mechanics in an attempt to develop non-local constitutive relations (Lazopoulos in Mech. Res. Commun. 33:753–757, 2006). Attempts have also been made to model microscopic forces using fractional derivatives (Vazquez in Nonlinear waves: classical and quantum aspects, pp. 129–133, 2004). Our approach in this paper differs from previous theoretical work, in that we develop a general framework directly from the classical continuum mechanics, by defining the laws of motion and the stresses using fractional derivatives. The timeliness and relevance of this work is justified by the surge in interest in applications of fractional order models to biological, physical and economic systems. The aim of the present paper is to lay the foundations for a new non-local model of continuum mechanics based on fractional order derivatives which we will refer to as the fractional model of continuum mechanics. Following the theoretical development, we apply this framework to two one-dimensional model problems: the deformation of an infinite bar subjected to a self-equilibrated load distribution, and the propagation of longitudinal waves in a thin finite bar.     

Nonlocal elasticity: an approach based on fractional calculus

Fractional calculus is the mathematical subject dealing with integrals and derivatives of non-integer order. Although its age approaches that of classical calculus, its applications in mechanics are relatively recent and mainly related to fractional damping. Investigations using fractional spatial derivatives are even newer. In the present paper spatial fractional calculus is exploited to investigate a material whose nonlocal stress is defined as the fractional integral of the strain field. The developed fractional nonlocal elastic model is compared with standard integral nonlocal elasticity, which dates back to Eringen’s works. Analogies and differences are highlighted. The long tails of the power law kernel of fractional integrals make the mechanical behaviour of fractional nonlocal elastic materials peculiar. Peculiar are also the power law size effects yielded by the anomalous physical dimension of fractional operators. Furthermore we prove that the fractional nonlocal elastic medium can be seen as the continuum limit of a lattice model whose points are connected by three levels of springs with stiffness decaying with the power law of the distance between the connected points. Interestingly, interactions between bulk and surface material points are taken distinctly into account by the fractional model. Finally, the fractional differential equation in terms of the displacement function along with the proper static and kinematic boundary conditions are derived and solved implementing a suitable numerical algorithm. Applications to some example problems conclude the paper.     

A non-associated flow rule for sheet metal forming

An improved model of material behavior is proposed that shows good agreement with experimental data for both yield and plastic strain ratios in uniaxial, equi-biaxial, and plane-strain tension under proportional loading for steel, aluminum and possibly other alloys. This model is based on a non-associated flow rule in which the plastic potential and yield surface functions are defined by quadratic functions of the stress tensor. The plastic potential aspect of the model is identical to that proposed by Hill for a quadratic anisotropic plastic potential defined in terms of measured r values. The new model differs in that the yield surface, although also defined by a quadratic function of the stress tensor, is defined independently of the plastic potential in terms of measured yield stresses. The model is developed and implemented in an FEM code that is based on a convected coordinate system. Since the associated flow rule, which assumes equivalency between the plastic potential and yield functions, is commonly accepted as a valid law in the theory of plastic deformation of most metals, the arguments for the associated flow rule are also discussed.     

A note on the definition of fractional derivatives applied in rheology

It is known that there exist obivious differences between the two most commonly used definitions of fractional derivatives—Riemann–Liouville (R–L) definition and Caputo definition. The multiple definitions of fractional derivatives in fractional calculus have hindered the application of fractional calculus in rheology. In this paper, we clarify that the R–L definition and Caputo definition are both rheologically imperfect with the help of mechanical analogues of the fractional element model (Scott–Blair model). We also clarify that to make them perfect rheologically, the lower terminals of both definitions should be put to ∞. We further prove that the R–L definition with lower terminal a →∞ and the Caputo definition with lower terminal a →∞ are equivalent in the differentiation of functions that are smooth enough and functions that have finite number of singular points. Thus we can define the fractional derivatives in rheology as the R–L derivatives with lower terminal a →∞ (or, equivalently, the Caputo derivatives with lower terminal a →∞) not only for steady-state processes, but also for transient processes. Based on the above definition, the problems of composition rules of fractional operators and the initial conditions for fractional differential equations are discussed, respectively. As an example we study a fractional oscillator with Scott–Blair model and give an exact solution of this equation under given initial conditions.     

Fractional plasticity for over-consolidated soft soil

The stress–strain response of over-consolidated soft soil, e.g., clay, is dependent on its pre-consolidation history and material state. In this study, a state-dependent constitutive model for over-consolidated soft soils is proposed by extending the fractional plasticity originally developed for granular soil. The state-dependent stress-dilatancy and peak failure behaviour of over-consolidated soft soil are analytically captured through stress-fractional gradient of the current yielding surface. In addition, a reference yielding surface describing the pre-consolidation history of soft soil is proposed. A combined hardening rule expressed as a function of both the incremental plastic volumetric and shear strains is suggested. To validate the proposed model, a series of drained and undrained test results of different soft soils with a wide range of over-consolidation ratios are simulated and compared. It is found that without using additional plastic potentials and state parameters, the developed fractional model can capture the state-dependent hardening and softening responses of soft soils subjected to different over-consolidation states.

     

An exact solution for a model of pressure-dependent plasticity in an un-steady plane strain process

An incompressible material obeying a pressure-dependent yield condition is confined between two planar plates which are inclined at an angle 2α. The plates intersect in a hinged line and the angle α slowly decreases from an initial value. An initial/boundary value problem for the flow of the material is formulated and solved for the stress and the velocity fields, the solution being in closed form. The material is assumed to obey a special case of the double slip and rotation model, which generalises the classical plastic potential model and is also a variant of the double shearing model. The solution for the velocity field may exhibit sliding or sticking at the plates. Solutions which exhibit sticking may have a rigidly rotating zone in the region adjacent to the plates. It is shown that sliding occurs when the value of α is less than a certain critical value α_c ; that sticking occurs without a rigid zone if α exceeds or equals α_c but is less than a second critical value α₀; and that sticking with a rigid zone adjacent to the plates occurs if α exceeds α₀. The values of α_c and α₀ coincide for a certain range of model parameters. Solutions which exhibit sliding are singular. Qualitative features of the solution found are compared with those of the solution for the classical plastic potential model.     

A non-associated constitutive model with mixed iso-kinematic hardening for finite element simulation of sheet metal forming

In this paper an anisotropic material model based on non-associated flow rule and mixed isotropic–kinematic hardening was developed and implemented into a user-defined material (UMAT) subroutine for the commercial finite element code ABAQUS. Both yield function and plastic potential were defined in the form of Hill’s [Hill, R., 1948. A theory of the yielding and plastic flow of anisotropic metals. Proc. R. Soc. Lond. A 193, 281–297] quadratic anisotropic function, where the coefficients for the yield function were determined from the yield stresses in different material orientations, and those of the plastic potential were determined from the r-values in different directions. Isotropic hardening follows a nonlinear behavior, generally in the power law form for most grades of steel and the exponential law form for aluminum alloys. Also, a kinematic hardening law was implemented to account for cyclic loading effects. The evolution of the backstress tensor was modeled based on the nonlinear kinematic hardening theory (Armstrong–Frederick formulation). Computational plasticity equations were then formulated by using a return-mapping algorithm to integrate the stress over each time increment. Either explicit or implicit time integration schemes can be used for this model. Finally, the implemented material model was utilized to simulate two sheet metal forming processes: the cup drawing of AA2090-T3, and the springback of the channel drawing of two sheet materials (DP600 and AA6022-T43). Experimental cyclic shear tests were carried out in order to determine the cyclic stress–strain behavior and the Bauschinger ratio. The in-plane anisotropy (r-value and yield stress directionalities) of these sheet materials was also compared with the results of numerical simulations using the non-associated model. These results showed that this non-associated, mixed hardening model significantly improves the prediction of earing in the cup drawing process and the prediction of springback in the sidewall of drawn channel sections, even when a simple quadratic constitutive model is used.     

Phenomenological study of parabolic and spherical indentation of elastic-ideally plastic material

A phenomenological study of parabolic and spherical indentation of elastic ideally plastic materials was carried out by using precise results of finite elements calculations. The study shows that no “pseudo-Hertzian” regime occurs during spherical indentation. As soon as the yield stress of the indented material is exceeded, a deviation from the, purely elastic Hertzian contact behaviour is found. Two elastic–plastic regimes and two plastic regimes are observed for materials of very large Young modulus to Yield stress ratio, E/σ_y. The first elastic–plastic regime corresponds to a strong evolution of the indented plastic zone. The first plastic regime corresponds to the commonly called “fully plastic regime”, in which the average indentation pressure is constant and equal to about three times the yield stress of the indented material. In this regime, the contact depth to penetration depth ratio tends toward a constant value, i.e. h_c/h = 1.47. h_c/h is only constant for very low values of yield strain (σ_y/E lower than 5 × 10^?6) when aE¹/Rσ_y is higher than 10,000. The second plastic regime corresponds to a decrease in the average indentation pressure and to a steeper increase in the pile-up. For materials with very large E/σ_y ratio, the second plastic regime appears when the value of the non-dimensional contact radius a/R is lower than 0.01. In the case of spherical and parabolic indentation, results show that the first plastic regime exists only for elastic-ideally plastic materials having an E/σ_y ratio higher than approximately 2.000.     

Model for charge/discharge-rate-dependent plastic flow in amorphous battery materials

Plastic flow is an important mechanism for relaxing stresses that develop due to swelling/shrinkage during charging/discharging of battery materials. Amorphous high-storage-capacity Li–Si has lower flow stresses than crystalline materials but there is evidence that the plastic flow stress depends on the conditions of charging and discharging, indicating important non-equilibrium aspects to the flow behavior. Here, a mechanistically-based constitutive model for rate-dependent plastic flow in amorphous materials, such as Li_xSi alloys, during charging and discharging is developed based on two physical concepts: (i) excess energy is stored in the material during electrochemical charging and discharging due to the inability of the amorphous material to fully relax during the charging/discharging process and (ii) this excess energy reduces the barriers for plastic flow processes and thus reduces the applied stresses necessary to cause plastic flow. The plastic flow stress is thus a competition between the time scales of charging/discharging and the time scales of glassy relaxation. The two concepts, as well as other aspects of the model, are validated using molecular simulations on a model Li–Si system. The model is applied to examine the plastic flow behavior of typical specimen geometries due to combined charging/discharging and stress history, and the results generally rationalize experimental observations.     

Non-linear viscoelastic behavior of polymer melts interpreted by fractional viscoelastic model

Very recently, researchers dealing with constitutive law pertinent viscoelastic materials put forward the successful idea to introduce viscoelastic laws embedded with fractional calculus, relating the stress function to a real order derivative of the strain function. The latter consideration leads to represent both, relaxation and creep functions, through a power law function. In literature there are many papers in which the best fitting of the peculiar viscoelastic functions using a fractional model is performed. However there are not present studies about best fitting of relaxation function and/or creep function of materials that exhibit a non-linear viscoelastic behavior, as polymer melts, using a fractional model. In this paper the authors propose an advanced model for capturing the non-linear trend of the shear viscosity of polymer melts as function of the shear rate. Results obtained with the fractional model are compared with those obtained using a classical model which involves classical Maxwell elements. The comparison between experimental data and the theoretical model shows a good agreement, emphasizing that fractional model is proper for studying viscoelasticity, even if the material exhibits a non-linear behavior.     

A new interpretation for the dynamic behaviour of complex fluids at the sol–gel transition using the fractional calculus

We propose to analyse power law shear stress relaxation modulus observed at the sol–gel transition (SGT) in many gelling systems in terms of fractional calculus. We show that the critical gel (gel at SGT) can be associated to a single fractional element and the gel in the post-SGT state to a fractional Kelvin–Voigt model. In this case, it is possible to give a physical interpretation to the fractional derivative order. It is associated to the power law exponent of the shear modulus related to the fractal dimension of the critical gel. A preliminary experimental application to silica alkoxide-based systems is given.

A fractional differential constitutive model for dynamic stress intensity factors of an anti-plane crack in viscoelastic materials

Fractional differential constitutive relationships are introduced to depict the history of dynamic stress inten- sity factors （DSIFs） for a semi-infinite crack in infinite viscoelastic material subjected to anti-plane shear impact load. The basic equations which govern the anti-plane deformation behavior are converted to a fractional wave-like equation. By utilizing Laplace and Fourier integral transforms, the fractional wave-like equation is cast into an ordinary differential equation （ODE）. The unknown function in the solution of ODE is obtained by applying Fourier transform directly to the boundary conditions of fractional wave-like equation in Laplace domain instead of solving dual integral equations. Analytical solutions of DSIFs in Laplace domain are derived by Wiener-Hopf technique and the numerical solutions of DSIFs in time domain are obtained by Talbot algorithm. The effects of four parameters α, β, b1, b2 of the fractional dif- ferential constitutive model on DSIFs are discussed. The numerical results show that the present fractional differential constitutive model can well describe the behavior of DSIFs of anti-plane fracture in viscoelastic materials, and the model is also compatible with solutions of DSIFs of anti-plane fracture in elastic materials.     

How to compute plastic zones of heterogeneous materials: A simple approach using classical continuum and fracture mechanics

This is a short technical paper on how to use classical continuum and fracture mechanics to calculate the plastic zones caused by cracks on heterogeneous or composite materials. As an example, a sample consisting of an α-phase and β-phase is used. A crack is introduced to the sample, and stress is then applied. The plastic zone in front of the crack resulting from the applied stress is then calculated using commercial software. The concept uses two-level modeling: a global model using homogenized stiffness from a unit cell of heterogeneous material and a local model for the α-phase and β-phase. While this paper is written for general purposes, a concrete example using ferrite and martensite is also presented along with the experimental data. General agreement between the model and the experiment is observed. This method eliminates the need for a cumbersome analytical approach.     

A comparison between crystal and isotropic strain gradient plasticity theories with accent on the role of the plastic spin

In the small deformation range, we consider crystal and isotropic “higher-order” theories of strain gradient plasticity, in which two different types of size effects are accounted for: (i) that dissipative, entering the model through the definition of an effective measure of plastic deformation peculiar of the isotropic hardening function and (ii) that energetic, included by defining the defect energy (i.e., a function of Nye's dislocation density tensor added to the free energy; see, e.g., [Gurtin, M.E., 2002. A gradient theory of single-crystal viscoplasticity that accounts for geometrically necessary dislocations. J. Mech. Phys. Solids 50, 5–32]). In order to compare the two modellings, we recast both of them into a unified deformation theory framework and apply them to a simple boundary value problem for which we can exploit the Γ-convergence results of [Bardella, L., Giacomini, A., 2008. Influence of material parameters and crystallography on the size effects describable by means of strain gradient plasticity. J. Mech. Phys. Solids 56 (9), 2906–2934], in which the crystal model is made isotropic by imposing that any direction be a possible slip system. We show that the isotropic modelling can satisfactorily approximate the behaviour described by the isotropic limit obtained from the crystal modelling if the former constitutively involves the plastic spin, as in the theory put forward in Section 12 of [Gurtin, M.E., 2004. A gradient theory of small-deformation isotropic plasticity that accounts for the Burgers vector and for dissipation due to plastic spin. J. Mech. Phys. Solids 52, 2545–2568]. The analysis suggests a criterium for choosing the material parameter governing the plastic spin dependence into the relevant Gurtin model.     

Basic equations for the microscopic theory of plasticity of polycrystalline materials with given texture

An expression for the yield stress of anisotropic materials is applied to the anisotropic strength of hard rolled copper foils whose crystallographic texture is known. We assume that this crystallographic texture is the only cause of the anisotropic plastic behaviour of the material. The model used for the yield stress is also used to deduce:

1. Stress-strain relations for isotropic polycrystalline materials;
2. A formula for the fully plastic strain tensor, applied to anisotropic hard rolled copper foils.

For the anisotropic copper foils considered the calculated curves of the yield stress and of the strain tensor as a function of the angle x between rolling and tensile direction agree qualitatively with the measured values. However, the theory is not complete, since the yield stress and the plastic strain tensor are both a function of a parameter Q, the fraction of the number of available crystallographic slip planes on which the maximum shear stress has reached the critical value τ_a. We assume that for “fully” plastic deformation a certain critical fraction Q _e of the total number of slip planes has to be active. The fraction Q _e is called the critical active quantity. With the parameter Q _e we adjust the calculated curves to the measured ones. The dependence of Q _e on the properties of the material (e.g. the crystallographic texture) is discussed in Appendix I.