Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Selecting the Optimal Mine Ventilation System via a Decision Making Framework under Hesitant Linguistic Environment
Next Article in Special Issue
Investigation on the Effect of Type of Cooling on the Properties of Aluminum Alloy during Warm/Hot Hydromechanical Deep Drawing
Previous Article in Journal
Three-Way Decisions with Interval-Valued Intuitionistic Fuzzy Decision-Theoretic Rough Sets in Group Decision-Making
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Study of Dynamic Properties of Fractional Viscoplasticity Model

1
Institute of Structural Engineering, Poznań University of Technology, Piotrowo 5 street, 60-965 Poznań, Poland
2
Institute of Fundamental Technological Research, Polish Academy of Sciences, Pawińskiego 5B street, 02-106 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Symmetry 2018, 10(7), 282; https://doi.org/10.3390/sym10070282
Submission received: 17 June 2018 / Accepted: 10 July 2018 / Published: 13 July 2018

Abstract

:
The fractional viscoplasticity (FV) concept combines the Perzyna type viscoplastic model and fractional calculus. This formulation includes: (i) rate-dependence; (ii) plastic anisotropy; (iii) non-normality; (iv) directional viscosity; (v) implicit/time non-locality; and (vi) explicit/stress-fractional non-locality. This paper presents a comprehensive analysis of the above mentioned FV properties, together with a detailed discussion on a general 3D numerical implementation for the explicit time integration scheme.

1. Introduction

In search of a generalization of existing models describing experimentally observed phenomena, the concept of fractional calculus [1,2] emerged as a tool that in the recent years became widely applied. Among the areas in which this theory has found application, it is worth mentioning mechanics where one can distinguish: (i) time-fractional models; (ii) space-fractional models; and (iii) stress-fractional models. For example, in [3], the time-fractional model was used to describe the time-dependent mechanical property evolution in ductile metals. The fractional oscillators were analyzed in [4], whereas the heat and mass transfer analysis in the framework of fractional calculus was presented in [5,6]. Furthermore, the analysis and modeling of turbulent flow in a porous medium [7], fluid transport induced by the osmotic pressure of glucose and albumin [8], wave propagation in the viscoelastic material [9], non-local boundary value problem [10], and evolution for the damage variable for hyperelastic materials [11] with an application of time-fractional derivative suggests great versatility of this approach. On the other hand, the space-fractional models are successfully used in mechanics to describe the deformation of a harmonic oscillator [12], deformation of an infinite bar subjected to a self-equilibrated load distribution [13], modeling plane strain and plane stress elasticity [14], Euler–Bernoulli beam [15], Darcy’s flow in porous media [16] and fractional strain formulation [17]. Finally, the stress-fractional models [18,19] and their finite element implementations were used to study the granular soils under drained cyclic loading [20], and monotonic triaxial compression [21]. Concluding, one should emphasize that regardless of the specific formulation, the fractional operators have one common feature, namely, the ‘change’ of a selected variable is based on integration over a closed interval, thus extending the definition of integer order derivative (defined in a point) and simultaneously introducing a non-locality in a given space.
It is commonly accepted that the Theory of Thermo-Viscoplasticity (TTV), which plays a central role in the following considerations, began with the publication of Perzyna [22], which, until the present day, serves a basis for many efforts in linking experimental and numerical results for different types of materials. The main results of this theory were discussed in a great number of papers that focused on phenomena such as propagation of mechanical and thermal waves [23,24], viscosity controlled by material parameter [25,26], dispersion [26], or implicit non-locality in the time variable [27]. Nonetheless, the classical TTV formulation does not include directional viscosity, and to include the non-normality extension needs, as all classical plasticity theories, postulation of an additional potential, which is not straightforward and causes the increase of material parameters. Furthermore, the same concern is relevant to the plastic anisotropy effects in terms of the original Perzyna model; to include this effect additional variables and evolution equations for them are needed to be postulated. This limitations were resolved by the generalization of the Perzyna formulation by definition of the fractional flow rule, first proposed in [18] and later developed in [19,28,29,30].
The implementation of the fractional plastic (rate independent) rule, for the Huber-Mises-Hencky (HMH) yield criterion, in the framework of implicit and explicit procedures and with examples on material point level, was presented in [28]. This was further developed in the subsequent article [29] to any smooth and convex yield criterion but still focusing on rate independent plastic flow. Concluding, in both these articles the non-locality in the stress state was present, however the implicit time non-locality common for the viscoplastic flow was not included in them.
This paper extends the concept of FV, which was first reported in [18], for the Initial Boundary Value Problem (IBVP), and provides a detailed discussion of the model material parameters. The parametric study includes the influence of the overstress power and the relaxation time (which is understood as implicit length scale parameter, as mentioned in [27]) on the dynamic properties of the FV model. Moreover, additional fractional material parameters, which induce the directional viscosity, the non-associative, and the anisotropic plastic flow, were also discussed.

2. Fractional Viscoplasticity

2.1. Remarks on Fractional Calculus

Fractional calculus (FC) introduces a new, universal method for calculating the intensity of changes of various quantities in mathematical models describing experimentally observed phenomena. FC implies a generalization of integer order derivatives, by fractional derivatives (FD). The selection of the FD definition (from an infinite number) can use a type of material as a criterion to obtain the best fitting of the constitutive model to a given experimental evidence. All definitions of the FD have a common property, namely they include summation over an interval abandoning the integer order derivative definition given at a single point; therefore they are called non-local. The classical derivative can be regarded a special case of the FD when its order becomes integer.
In order to explain the FD concept, let us consider a generalized fractional differential operator B P α as a composition of fractional integral K P α with classical integer (n-th) differential operator [31]
B P α = K P n - α d n d t n ,
where α is the order of the derivative, n = α + 1 , · denotes the floor function, P is a parameter set (described below) and ∘ denotes the composition operator. B P α is referred to as the fractional differential operator B (B-op) of order α and p-set P, and analogously K P α identifies the K (K-op) fractional integer operator of order α and p-set P.
The definition of K for the parameter set P = a , t , b , p , q can be given as
K P α f t = p a t k α ( t , τ ) f ( τ ) d τ + q t b k α ( τ , t ) f ( τ ) d τ ,
where t a , b and a < t < b , p , q are real numbers, and k α ( t , τ ) is a kernel that depends on the order of the derivative α . It can be shown that if k α is a difference kernel, i.e., k α ( t , τ ) = k α ( t - τ ) and k α L 1 0 , b - a then L 1 b , a L 1 b , a is well defined, bounded and linear. For explicit definition, the special form of the kernel function can be assumed
k α ( t - τ ) = 1 Γ ( α ) t - τ α - 1 ,
then for P = a , t , b , 1 , 0
K P α f t = 1 Γ ( α ) a t ( t - τ ) α - 1 f ( τ ) d τ = a I t α f ( t ) ,
is obtained or, if P = a , t , b , 0 , 1 then
K P α f t = 1 Γ ( α ) t b ( τ - t ) α - 1 f ( τ ) d τ = t I b α f ( t ) ,
where Γ is the Euler gamma function. Equations (4) and (5) describe the left and right Riemann-Liouville fractional integrals of the order α , respectively. The application of these operators in Equation (1) leads to the following fractional derivative definitions:
B P α f t = a C D t α f t = 1 Γ ( n - α ) a t f ( n ) ( τ ) ( t - τ ) α - n + 1 d τ ,
for t > a , and
- B P α f t = t C D b α f t = ( - 1 ) n Γ ( n - α ) t b f ( n ) ( τ ) ( τ - t ) α - n + 1 d τ ,
for t < b . The FD operators a C D t α f t and t C D b α f t are known as the left- and right-sided Caupto fractional integrals.
Finally, for the purpose of further definition of the FV, the Riesz-Caputo (RC) derivative can be expressed as a linear combination of previously given left and right Caputo derivatives
a R C D b α f t = 1 2 a C D t α f t + ( - 1 ) n t C D b α f t .
It can be shown that for the RC derivative the fundamental property of integer order derivatives is preserved, that is, the derivative of a constant is zero.

2.2. Basic Concepts

In the following section Voigt notation is applied, thus the second rank tensors are ordered as ( 6 × 1 column matrix)
t = t 11 t 22 t 33 t 23 t 13 t 12 T = t 1 t 2 t 3 t 4 t 5 t 6 T ,
whereas the fourth order tensors are represented by 6 × 6 matrices ordered in accordance with the rule used in Equation (9).
Deformation assumes the additive decomposition of total strain, therefore
ε = ε e + ε v p ,
or in a rate form
ε ˙ = ε ˙ e + ε ˙ v p ,
where is the total strain, ε e is the elastic strain and ε v p is the viscoplastic strain. Next, due to thermodynamic restrictions, the elastic strain is related to elastic stress through Hooke’s law
σ e = L e ε e ,
where σe denotes the Cauchy stress tensor and L e denotes the elastic constitutive tensor. The rate of viscoplastic strain is analogous to the classical viscoplastic definition, namely
ε ˙ v p = Λ p ,
where Λ is a scalar multiplier and p is the second order unit tensor which governs the direction of viscoplastic flow. As the p tensor is normalized, the magnitude of ε ˙ v p depends solely on the Λ parameter.
Following the concept introduced by Perzyna [22], this parameter is expressed as
Λ = γ Φ ( F ) ,
where γ = 1 T m is the viscosity parameter, Φ is the overstress function that depends on the rate-independent yield surface F, and · is Macaulay brackets. It is well-known that γ introduces implicit time non-locality in the viscoplastic model [27]. Furthermore, the function Φ has the following form
Φ ( F ) = F m v p = J 2 κ - 1 m v p ,
where J 2 denotes the second invariant of stress deviator and κ is the static yield stress in simple shear.
Finally, the remaining object needed to be defined is the tensor p . In this place, the difference between the classical theory of viscoplasticity and the new approach is most evident. Let us recall, that in the classical formulation the direction of yield is normal to yield surface and p can be written as
p = F F - 1 .
It is also well known, that Equations (16) and (15) indicate that the viscoplastic strain is coaxial with the deviatoric stress tensor (associated flow). As a result, the volume change can occur in the range of elastic deformations only. For modern materials such as metal-matrix composites, this assumption is no longer valid. Therefore, the constitutive model should be modified to capture this phenomenon.
The fractional approach assumes the application of the RC operator to p definition [18]. In such a case, Equation (16) is generalized to the form
p = D α F D α F - 1 ,
where D α stands for the RC operator (see Equation (8)). It is worth noting that the proposed formulation of p introduces the anisotropy of viscoplastic flow and furthermore (due to non-associativity) develops a tool to control the volume change in the plastic range of material behaviour [18]. Another essential remark is that Equation (17) introduces explicit stress-fractional non-locality in the overall model. It is important that the thermodynamic restrictions are formulated in a standard manner, and because of complicated structure of Equation (17) they are checked incrementally in the numerical procedure (see [29] for a detailed discussion).

3. Implementation

Introduction of the three-dimensional fractional viscoplastic model requires a numerical procedure that governs the solution. Since our considerations are focused on extreme dynamic processes, the explicit time integration was chosen for finite element method. Therefore, the ABAQUS/Explicit code was utilized together with the user subroutine VUMAT. The critical steps of the implementation are presented below.
In the first step, Hooke’s law is written as
σ 11 σ 22 σ 33 σ 23 σ 13 σ 12 = 2 G + λ λ λ 0 0 0 λ 2 G + λ λ 0 0 0 λ λ 2 G + λ 0 0 0 0 0 0 G 0 0 0 0 0 0 G 0 0 0 0 0 0 G ε 11 ε 22 ε 33 2 ε 23 2 ε 13 2 ε 12 ,
where G = E / 2 ( 1 + ν ) and λ = E ν / ( 1 + ν ) ( 1 - 2 ν ) are elastic constants and E and ν denote Young’s modulus and Poisson’s ratio, respectively. Next, because of application of the HMH yield criterion, the yield function in a matrix form is as follows
F ( ) = σ 11 σ 22 σ 33 σ 23 σ 13 σ 12 T 1 - 1 2 - 1 2 0 0 0 - 1 2 1 - 1 2 0 0 0 - 1 2 - 1 2 1 0 0 0 0 0 0 3 0 0 0 0 0 0 3 0 0 0 0 0 0 3 σ 11 σ 22 σ 33 σ 23 σ 13 σ 12 - 3 κ 2 = 0 .
Afterwards, the increment of viscoplastic strain (see Equations (13) and (17)) can be written in the form
Δ ε v p = Δ t Λ p = Δ t Λ D σ 11 α F D σ 22 α F D σ 33 α F D σ 23 α F D σ 13 α F D σ 12 α F T D α F ,
where (cf. [28])
p 11 p 22 p 33 p 23 p 13 p 12 = k ( 11 ) M - 1 2 k ( 11 ) M - 1 2 k ( 11 ) M 0 0 0 - 1 2 k ( 22 ) M k ( 22 ) M - 1 2 k ( 22 ) M 0 0 0 - 1 2 k ( 33 ) M - 1 2 k ( 33 ) M k ( 33 ) M 0 0 0 0 0 0 3 k ( 23 ) M 0 0 0 0 0 0 3 k ( 13 ) M 0 0 0 0 0 0 3 k ( 12 ) M σ 11 σ 22 σ 33 σ 23 σ 13 σ 12 + k ( 11 ) Q k ( 22 ) Q k ( 33 ) Q 3 k ( 23 ) Q 3 k ( 13 ) Q 3 k ( 12 ) Q .
Symbols in Equation (21) denotes
k ( i j ) M = Γ ( 2 ) Γ ( 2 - α ) Δ ( i j ) L 1 - α + Δ ( i j ) R 1 - α ,
k ( i j ) Q = Γ ( 2 ) Γ ( 2 - α ) - 1 2 Γ ( 3 ) Γ ( 3 - α ) Δ ( i j ) R 2 - α - Δ ( i j ) L 2 - α ,
Δ L = Δ ( 11 ) L Δ ( 22 ) L Δ ( 33 ) L Δ ( 23 ) L Δ ( 13 ) L Δ ( 12 ) L T ,
and
Δ R = Δ ( 11 ) R Δ ( 22 ) R Δ ( 33 ) R Δ ( 23 ) R Δ ( 13 ) R Δ ( 12 ) R T ,
where
a ( i j ) = σ i j - Δ ( i j ) L , b ( i j ) = σ i j + Δ ( i j ) R .
Terminals a ( i j ) , b ( i j ) are needed to define the partial fractional derivatives in Equation (17) that enforce the directional nature of the fractional viscoplastic flow–the subscripts L and R corresponds to the left and the right Caputo derivatives, respectively. In addition, by introducing sections that extend the calculation beyond the material point, a virtual neighbourhood is obtained that results in a non-locality in a stress state.
The interpretation of the virtual surrounding in a stress state depends on the specific material (see [28]), but in general it could be understood as a (homogenized) phenomenological measure of some instability, e.g., for metallic materials it is connected with dislocation nucleation [32,33,34], nucleation of voids [35] or breakup of grains [36,37,38] (see review paper [39]). By way of illustration, Figure 1 shows the cross-section of this virtual neighbourhood in the σ 2 - σ 3 plane.
The analysis of Equation (21) shows differences in relation to the classical viscoplasticity where the change in volume may only occur in the elastic range—in the classical case, the trace of the p tensor equals 0. This condition is abandoned in the fractional formulation (when α ( 0 , 1 ) ), thus explicitly providing a tool to control the evolution of volume in the plastic range through α and parametric vectors Δ L and Δ R . Moreover, the versatility of the fractional approach is proven for α = 1 , for which the associated plastic flow as a special case is obtained.
Finally, the flowchart was formulated that presents the general calculation scheme for the elasto-viscoplastic material in the framework of the fractional viscoplastic flow rule for explicit time integration in VUAMT subroutine (see Figure 2). The VUAMT subroutine aims at determination of the values of Cauchy stresses and updating strains and internal variables at time t n + 1 based on the knowledge of these parameters at the previous moment t n . The procedure starts with the calculation of the elastic trial stress, which is later used to establish the value of the yield criterion. If this criterion is fulfilled, the plastic multiplier Λ and the direction of plastic flow p are calculated according to the flow rule; otherwise the elastic step is conducted.

4. Parametric Study: Uniaxial Tension

4.1. Description of the Numerical Experiment

The conducted parametric study is focused on the material point level represented by a unit cube with dimensions of 1 × 1 × 1 mm discretized by a single finite element C38DR (linear, eight-node brick with reduce integration). The boundary conditions required to achieve uniaxial constraints are shown in Figure 3. Basic mechanical properties were assumed as for the carbon steel, therefore the elastic range was characterized by Young’s modulus E = 205 GPa and Poisson’s ratio ν = 0.27. The fractional flow rule presented in Section 2 was applied in the plastic range. The static yield stress in simple shear for the selected material was κ = 605 MPa. Other material parameters, depending on the studied case, were chosen as described below.
The analyzed cases of fractional flow were divided into two groups to show how various combinations of model parameters influence plastic deformation. The first group is focused on the value of the stress-fractional non-locality spread ( Δ L , R ) and the order ( α ) of the fractional flow. In the second group, the influence of the material parameters T m and m under various speeds of the imposed displacement is closely studied. Anticipating the anisotropic behaviour of the fractional material model, these two groups were further subdivided according to the direction where the dominant viscoplastic flow was expected. Hence, two cases were formed for tension direction ( Δ 22 L , R = 0.005 κ 3.0 MPa ) and direction perpendicular to tension ( Δ 11 L , R = 0.005 κ 3.0 MPa ). In each of those cases, other values of the Δ were set to 0.0017 κ 1.0 MPa .
Two kinds of plots were used to present the results of the parametric study. The first kind exemplify the relation between three normal strains ε 11 , ε 22 , ε 33 . The second type shows the stress–strain relation in the tension (2) direction (it should be pointed out that for this kind of plots, a ‘softening’ is observed, especially for highest tension velocities, however, this effect is due to the lateral stresses induced by the inertia effects and is not due to constitutive model κ = c o n s t . See Figure 4, where this effect is negligible due to relatively small tension velocity v = 1 m s ). The research on influence of the fractional derivative order was performed for a set α { 0 . 1 , 0 . 25 , 0 . 5 , 0 . 75 , 0 . 99 , 1 . 0 } —as mentioned earlier, fractional generalization of the viscoplasticity reduces to the classical solution for α = 1 . The study of T m and m was conducted for three different velocities of tension, i.e., v = 1 , 25 and 50 m s .

4.2. Influence of the Order of FV and Non-Locality in a Stress State on Plastic Flow

4.2.1. Study of Intensified Plastic Flow in Tension Direction for Different Orders of Flow

Figure 4 presents the material response to the applied tension velocity of v = 1 m s for different flow intensities in tension direction and flow orders. Increasing the flow intensity parameter ( Δ L , R ) causes higher evolution of the plastic flow in the chosen direction but in this case the velocity of the load is not sufficient to reveal different behaviour in the stress–strain relation for different values of α (see Figure 4). Next, for the same configuration of material parameters higher velocity is applied, namely v = 25 m s . At this speed (Figure 5), a slight waveform begins to be visible for Δ 22 L , R = 1 . 0 MPa . The amplitude of the stress signal increases with the increases Δ 22 L , R . Additionally, the influence of α is shown because the decrease in its value translates into greater amplitude of oscillations. So, we conclude that both fractional parameters, control the dynamic properties of a fractional model.

4.2.2. Study of Intensified Plastic Flow Perpendicular to the Tension Direction for Different Orders of Flow

Results presented in this section were obtained for a parameter set similar to this in Section 4.2.1 with the difference that flow intensity is increased in the direction perpendicular to tension, namely Δ 11 L , R = 3 . 0 MPa . Others components of vectors in Equations (24) and (25) equal 1. As in the discussion in the previous section, the tension velocity of v = 1 m s is insufficient to reveal the influence of α on the stress–strain relation (see Figure 6). However, Figure 7 shows that the material prefers to deform in (1) direction when the magnitude of Δ 11 L , R growths, hence the ε 11 / ε 33 ratio is greater then 1. Moreover, the intensity of the flow in the preferred direction increases as the value of α diminishes to 0. Next, as before, the velocity is increased to v = 25 m s . Figure 8 shows that when the value of Δ L , R rises, greater amplitude of the oscillation and hardening of the material can be observed. This last effect is inversely proportional to the order of the fractional flow. The relation between ε 11 , ε 22 and ε 33 is presented in Figure 9 and is very similar to Figure 7 with the only distinction that slight oscillation occurs as a result of higher velocity. So, we conclude that both fractional parameters control the anisotropic properties of a fractional model in the plastic range.

4.3. Influence of the Relaxation Time and the Overstress Power

4.3.1. Study of the Fractional Flow Under Different Dynamic Loading Rates for Intensified Plastic Flow in Tension Direction

Here we assume that the intensified plastic flow, determined by Δ 22 L , R = 3.0 MPa , is in tension direction. Figure 10 presents the effect of different relaxation times for various velocities of tension. It should be pointed out that in order to increase clarity of interpretation both the classical ( α = 1 ) and the fractional ( α = 0 . 75 ) solutions are compared on each graph. As can be seen, when the relaxation time grows, the hardening of the material as well as the stress waves oscillations increase. The latter is especially pronounced for the relaxation time T m = 2.5 e - 5 s .
Figure 11 presents the result of increasing the value of the overstress parameter m. For the fractional ( α = 0.75 ) viscoplastic material the stress level is smaller in relation to the classical ( α = 1.0 ) viscoplastic solution (Figure 10 and Figure 11).
Results discussed above indicate that the relaxation time and overstress power, together with fractional parameters, control the level of strain rate hardening and stress waves oscillation amplitude.

4.3.2. Study of the Fractional Flow Under Different Dynamic Loading for the Intensified Plastic Flow Perpendicular to the Tension Direction

In this section, it is assumed that fractional flow is intensified in the direction perpendicular to tension load ( Δ 11 L , R = 3 . 0 M P a ). Figure 12 shows that raising relaxation time causes similar effects to those discussed in the previous section. The biggest change occurs for the relaxation time T m = 2 . 5 e - 5 s both in the stress level and the stress wave oscillations frequency. In Figure 13 a clear anisotropy of the plastic deformation for α = 0 . 75 can be noticed. As before, it is observed that increasing the value of m causes a growth in the material strain hardening without any apparent influence on the frequency of the stress wave (see Figure 14). By analogy to what is presented for T m , in Figure 15 the dominant nature of ε 11 can be observed, in addition to more pronounced oscillations for greater values of m. The analysis of the σ - ε relation for α = 1 and α = 0 . 75 also revealed that the stress levels of the fractional model are generally greater than for the classical solution (see Figure 12 and Figure 14). This last observation is different from what was discovered in the previous section.
In conclusion, as before, a prominent influence of the relaxation time and overstress power, together with fractional parameters, on the level of strain rate hardening and stress waves oscillation amplitude is observed. Additionally, the impact on dynamic properties of deformation anisotropy for fractional material was confirmed.

4.4. Study of the Disperse Character of the Fractional Viscoplastic Stress Waves

The last study performed was the analysis of the disperse character of the fractional viscoplastic stress waves. As shown in the previous sections, the uniaxial dynamic deformation induces a stress waves. The regularity of the oscillations was determined by averaging intervals between peaks and then calculating the frequency. Table 1 lists the stress wave frequencies for Δ 22 L , R = 3 . 0 MPa and Δ 11 L , R = 3 . 0 MPa , which were depicted in the middle graphs of Figure 10 and Figure 12. The material parameter Δ L , R does not appear in the classical viscoplasticity, so there is no change in the stress wave frequency for various values of this parameter when α = 1 , thus tension velocity is only important. However, for fractional material, i.e., when α = 0 . 75 , both flow intensity parameters and tension velocities modulate the frequency of stress waves. For α = 0 . 75 the frequencies are higher when the distinguished direction is co-linear with tension ( Δ 22 = 3 . 0 MPa ).
The results presented in Table 2 correspond to the investigation of the role of the relaxation time and overstress power discussed in Section 4.3. Regardless of the value of α , the largest change in frequency can be observed between columns 2 and 3, that is for T m = 2 . 5 e - 6 s and T m = 2 . 5 e - 5 s . The comparison of the values in columns 1 ( T m = 2 . 5 e - 7 s ) and 2 ( T m = 2 . 5 e - 6 s ) shows that the stress wave frequencies are the same or very similar. The impact of m shows that the significant influence of this parameter was only recorded for T m = 2 . 5 e - 5 s . For both α = 0 . 75 and α = 1 the increase of the overstress power results in the increase of the frequency of the stress wave.
One concludes that the relaxation time and overstress power, together with fractional parameters, control the dispersive character of stress waves, and even more, makes this attribute directional.

5. Conclusions

The analysis of the dynamic properties of the Perzyna model of viscoplasticity (implicit time, non-local) generalized using fractional calculus (explicit, stress-fractional, non-local) leads to the following conclusions:
  • Fractional viscoplasticity introduces an additional set of material parameters, namely flow order α and virtual stress state surrounding Δ .
  • Fractional parameters α and Δ control the dynamic properties of the fractional model, especially hardening, the character of the stress waves, and plastic anisotropy.
  • The direction of the flow vector is controlled by Δ , which in general leads to non-normality of plastic flow.
  • As in the classical Perzyna model, the relaxation time T m and the overstress power m affect the strain rate hardening and the character of the stress waves.
  • Induced plastic anisotropy of the fractional model should be regarded not only in the classical sense as directional deformation but also as directional viscosity, which results in directional dispersive character.
The above results are fundamental from the point of view of modeling strain localization and damage phenomena. Both these aspects will serve as a base for future studies.

Author Contributions

All of the authors have contributed to the writing of this paper. They read and approved the manuscript.

Funding

The National Science Centre, Poland under Grant No. 2017/27/B/ST8/00351.

Acknowledgments

This work is supported by the National Science Centre, Poland under Grant No. 2017/27/B/ST8/00351.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Podlubny, I. Fractional Differential Equations. In Mathematics in Science and Engineering; Academin Press: Cambridge, MA, USA, 1999; Volume 198. [Google Scholar]
  2. Samko, S.; Kilbas, A.; Marichev, O. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach: Amsterdam, The Netherlands, 1993. [Google Scholar]
  3. Meng, R.; Yin, D.; Zhou, C.; Wu, H. Fractional description of time-dependent mechanical property evolution in materials with strain softening behavior. Appl. Math. Modell. 2016, 40, 398–406. [Google Scholar] [CrossRef]
  4. Li, M. Three Classes of Fractional Oscillators. Symmetry 2018, 10, 40. [Google Scholar] [CrossRef]
  5. Zhao, J.; Zheng, L.; Zhang, X.; Liu, F. Convection heat and mass transfer of fractional MHD Maxwell fluid in a porous medium with Soret and Dufour effects. Int. J. Heat Mass Trans. 2016, 103, 203–210. [Google Scholar] [CrossRef]
  6. Jinhu, Z.; Liancun, Z.; Xinxin, Z.; Fawang, L. Mixed convection heat transfer of viscoelastic fluid along an inclined plate obeying the fractional constitutive laws. Heat Trans. Res. 2017, 48, 1165–1178. [Google Scholar]
  7. Zhang, X.; Liu, L.; Wu, Y. The uniqueness of positive solution for a fractional order model of turbulent flow in a porous medium. Appl. Math. Lett. 2014, 37, 26–33. [Google Scholar] [CrossRef]
  8. Cherniha, R.; Gozak, K.; Waniewski, J. Exact and Numerical Solutions of a Spatially-Distributed Mathematical Model for Fluid and Solute Transport in Peritoneal Dialysis. Symmetry 2016, 8, 50. [Google Scholar] [CrossRef]
  9. Atanackovic, T.; Janev, M.; Oparnica, L.; Pilipovic, S.; Zorica, D. Space–time fractional Zener wave equation. Proc. Math. Phys. Eng. Sci. 2015, 471, 20140614. [Google Scholar]
  10. Ren, T.; Li, S.; Zhang, X.; Liu, L. Maximum and minimum solutions for a nonlocal p-Laplacian fractional differential system from eco-economical processes. Bound. Value Probl. 2017, 2017, 118. [Google Scholar] [CrossRef]
  11. Sumelka, W.; Voyiadjis, G. A hyperelastic fractional damage material model with memory. Int. J. Solids Struct. 2017, 124, 151–160. [Google Scholar] [CrossRef]
  12. Klimek, M. Fractional sequential mechanics—Models with symmetric fractional derivative. Czechoslov. J. Phys. 2001, 51, 1348–1354. [Google Scholar] [CrossRef]
  13. Drapaca, C.; Sivaloganathan, S. A Fractional Model of Continuum Mechanics. J. Elast. 2012, 107, 107–123. [Google Scholar] [CrossRef]
  14. Sumelka, W.; Szajek, K.; Łodygowski, T. Plane strain and plane stress elasticity under fractional continuum mechanics. Arch. Appl. Mech. 2015, 89, 1527–1544. [Google Scholar] [CrossRef]
  15. Tomasz, B. Analytical and numerical solution of the fractional Euler–Bernoulli beam equation. J. Mech. Mater. Struct. 2017, 12, 23–34. [Google Scholar]
  16. Lazopoulos, K.; Lazopoulos, A. Fractional vector calculus and fluid mechanics. J. Mech. Behav. Mater. 2017, 26, 43–54. [Google Scholar] [CrossRef]
  17. Peter, B. Dynamical Systems Approach of Internal Length in Fractional Calculus. Eng. Trans. 2017, 65, 209–215. [Google Scholar]
  18. Sumelka, W. Fractional viscoplasticity. Mech. Res. Commun. 2014, 56, 31–36. [Google Scholar] [CrossRef]
  19. Sun, Y.; Shen, Y. Constitutive model of granular soils using fractional-order plastic-flow rule. Int. J. Geomech. 2017, 17, 04017025. [Google Scholar] [CrossRef]
  20. Sun, Y.; Xiao, Y. Fractional order model for granular soils under drained cyclic loading. Int. J. Numer. Anal. Meth. Geomech. 2017, 41, 555–577. [Google Scholar] [CrossRef]
  21. Sun, Y.; Xiao, Y. Fractional order plasticity model for granular soils subjected to monotonic triaxial compression. Int. J. Solids Struct. 2017, 118-119, 224–234. [Google Scholar] [CrossRef]
  22. Perzyna, P. The constitutive equations for rate sensitive plastic materials. Q. Appl. Math. 1963, 20, 321–332. [Google Scholar] [CrossRef] [Green Version]
  23. Glema, A.; Łodygowski, T. On importance of imperfections in plastic strain localization problems in materials under impact loading. Arch. Mech. 2002, 54, 411–423. [Google Scholar]
  24. Glema, A. Analysis of wave nature in plastic strain localization in solids. In Rozprawy; Publishing House of Poznan University of Technology: Poznan, Poland, 2004; Volume 379. (In Polish) [Google Scholar]
  25. Glema, A.; Łodygowski, T.; Perzyna, P. Interaction of deformation waves and localization phenomena in inelastic solids. Comput. Meth. Appl. Mech. Eng. 2000, 183, 123–140. [Google Scholar] [CrossRef] [Green Version]
  26. Glema, A.; Łodygowski, T.; Perzyna, P. Localization of plastic deformations as a result of wave interaction. Comput. Assist. Mech. Eng. Sci. 2003, 10, 81–91. [Google Scholar]
  27. Perzyna, P. The Thermodynamical Theory of Elasto-Viscoplasticity. Eng. Trans. 2005, 53, 235–316. [Google Scholar]
  28. Sumelka, W.; Nowak, M. Non-normality and induced plastic anisotropy under fractional plastic flow rule: A numerical study. Int. J. Numer. Anal. Meth. Geomech. 2016, 40, 651–675. [Google Scholar] [CrossRef]
  29. Sumelka, W.; Nowak, M. On a general numerical scheme for the fractional plastic flow rule. Mech. Mater. 2017, 116, 120–129. [Google Scholar] [CrossRef]
  30. Xiao, R.; Sun, H.; Chen, W. A finite deformation fractional viscoplastic model for the glass transition behavior of amorphous polymers. Int. J. Non-Linear Mech. 2017, 93, 7–14. [Google Scholar] [CrossRef]
  31. Odzijewicz, T.; Malinowska, A.; Torres, D. Green’s theorem for generalized fractional derivatives. Fract. Calc. Appl. Anal. 2013, 16, 64–75. [Google Scholar] [CrossRef] [Green Version]
  32. Ziegler, H. An Introduction to Thermomechanics; North Holland Series in Applied Mathematics and Mechanics; North-Holland Publishing Company: North Holland, The Netherlands, 1983; Volume 21. [Google Scholar]
  33. Dao, M.; Asaro, R. Non-Schmid effects and localized plastic flow in intermetallic alloys. Mater. Sci. Eng. A 1993, 170, 143–160. [Google Scholar] [CrossRef]
  34. Racherla, V.; Bassani, J. Strain burst phenomena in the necking of a sheet that deforms by non-associated plastic flow. Modell. Simul. Mater. Sci. Eng. 2007, 15, S297–S311. [Google Scholar] [CrossRef]
  35. Marin, E.; McDowell, D. Models for Compressible Elasto-Plasticity Based on Internal State Variables. Int. J. Damage Mech. 1998, 7, 47–83. [Google Scholar] [CrossRef] [Green Version]
  36. Leffers, T. Lattice rotations during plastic deformation with grain subdivision. Mater. Sci. Forum 1994, 157–162, 1815–1820. [Google Scholar] [CrossRef]
  37. Hughes, D.; Liu, Q.; Chrzan, D.; Hansen, N. Scaling of microstructural parameters: Misorientations of deformation induced boundaries. Acta Mater. 1997, 45, 105–112. [Google Scholar] [CrossRef]
  38. Steinmann, P.; Kuhl, E.; Stein, E. Aspects of non-associated single crystal plasticity: Influence of non-schmid effects and localization analysis. Int. J. Solids Struct. 1998, 35, 4437–4456. [Google Scholar] [CrossRef]
  39. McDowell, D. Viscoplasticity of heterogeneous metallic materials. Mater. Sci. Eng. R 2008, 62, 67–123. [Google Scholar] [CrossRef]
Figure 1. Virtual surrounding of a material point.
Figure 1. Virtual surrounding of a material point.
Symmetry 10 00282 g001
Figure 2. VUMAT subroutine flowchart for the fractional viscoplastic rule.
Figure 2. VUMAT subroutine flowchart for the fractional viscoplastic rule.
Symmetry 10 00282 g002
Figure 3. Unit cubic model restricted to uniaxial tension.
Figure 3. Unit cubic model restricted to uniaxial tension.
Symmetry 10 00282 g003
Figure 4. Influence of the order α and the value of material parameter Δ 22 on the stress–strain relation, for: v = 1 m s , T m = 2 . 5 e - 6 s , m = 1 .
Figure 4. Influence of the order α and the value of material parameter Δ 22 on the stress–strain relation, for: v = 1 m s , T m = 2 . 5 e - 6 s , m = 1 .
Symmetry 10 00282 g004
Figure 5. Influence of the order α and the value of material parameter Δ 22 on the stress–strain relation, for: v = 25 m s , T m = 2 . 5 e - 6 s , m = 1 .
Figure 5. Influence of the order α and the value of material parameter Δ 22 on the stress–strain relation, for: v = 25 m s , T m = 2 . 5 e - 6 s , m = 1 .
Symmetry 10 00282 g005
Figure 6. Influence of the order α and the value of material parameter Δ 11 on the stress–strain relation, for: v = 1 m s , T m = 2 . 5 e - 6 s , m = 1 .
Figure 6. Influence of the order α and the value of material parameter Δ 11 on the stress–strain relation, for: v = 1 m s , T m = 2 . 5 e - 6 s , m = 1 .
Symmetry 10 00282 g006
Figure 7. Influence of the order α and the value of material parameter Δ 11 on the relation between three normal stresses, for: v = 1 m s , T m = 2 . 5 e - 6 s , m = 1 .
Figure 7. Influence of the order α and the value of material parameter Δ 11 on the relation between three normal stresses, for: v = 1 m s , T m = 2 . 5 e - 6 s , m = 1 .
Symmetry 10 00282 g007
Figure 8. Influence of the order α and the value of material parameter Δ 11 on the stress–strain relation, for: v = 25 m s , T m = 2 . 5 e - 6 s , m = 1 .
Figure 8. Influence of the order α and the value of material parameter Δ 11 on the stress–strain relation, for: v = 25 m s , T m = 2 . 5 e - 6 s , m = 1 .
Symmetry 10 00282 g008
Figure 9. Influence of the order α and the value of material parameter Δ 11 on the relation between three normal stresses, for: v = 25 m s , T m = 2 . 5 e - 6 s , m = 1 .
Figure 9. Influence of the order α and the value of material parameter Δ 11 on the relation between three normal stresses, for: v = 25 m s , T m = 2 . 5 e - 6 s , m = 1 .
Symmetry 10 00282 g009
Figure 10. Influence of the relaxation parameter T m and the value of applied velocity field v on the stress–strain relation, for: α = 0 . 75 , m = 1 , Δ 22 = 3 . 0 .
Figure 10. Influence of the relaxation parameter T m and the value of applied velocity field v on the stress–strain relation, for: α = 0 . 75 , m = 1 , Δ 22 = 3 . 0 .
Symmetry 10 00282 g010
Figure 11. Influence of the material parameter m and the value of applied velocity field v on the stress–strain relation, for: α = 0 . 75 , T m = 2 . 5 e - 6 s , Δ 22 = 3 . 0 .
Figure 11. Influence of the material parameter m and the value of applied velocity field v on the stress–strain relation, for: α = 0 . 75 , T m = 2 . 5 e - 6 s , Δ 22 = 3 . 0 .
Symmetry 10 00282 g011
Figure 12. Influence of the relaxation parameter T m and the value of applied velocity field v on the stress–strain relation, for: α = 0 . 75 , m = 1 , Δ 11 = 3 . 0 .
Figure 12. Influence of the relaxation parameter T m and the value of applied velocity field v on the stress–strain relation, for: α = 0 . 75 , m = 1 , Δ 11 = 3 . 0 .
Symmetry 10 00282 g012
Figure 13. Influence of the relaxation parameter T m and the value of applied velocity field v on the relation between three normal stresses, for: α = 0 . 75 , m = 1 , Δ 11 = 3 . 0 .
Figure 13. Influence of the relaxation parameter T m and the value of applied velocity field v on the relation between three normal stresses, for: α = 0 . 75 , m = 1 , Δ 11 = 3 . 0 .
Symmetry 10 00282 g013
Figure 14. Influence of the material parameter m and the value of applied velocity field v on the stress–strain relation, for: α = 0 . 75 , T m = 2 . 5 e - 6 s , Δ 11 = 3 . 0 .
Figure 14. Influence of the material parameter m and the value of applied velocity field v on the stress–strain relation, for: α = 0 . 75 , T m = 2 . 5 e - 6 s , Δ 11 = 3 . 0 .
Symmetry 10 00282 g014
Figure 15. Influence of the material parameter m and the value of applied velocity field v on the relation between three normal stresses, for: α = 0 . 75 , T m = 2 . 5 e - 6 s , Δ 11 = 3 . 0 .
Figure 15. Influence of the material parameter m and the value of applied velocity field v on the relation between three normal stresses, for: α = 0 . 75 , T m = 2 . 5 e - 6 s , Δ 11 = 3 . 0 .
Symmetry 10 00282 g015
Table 1. The stress wave frequencies for T m = 2 . 5 e - 6 s , m = 1
Table 1. The stress wave frequencies for T m = 2 . 5 e - 6 s , m = 1
Δ 22 = 3.0 Δ 11 = 3.0
v = 25 m s α = 1 2.085 MHz2.085 MHz
α = 0 . 75 2.108 MHz1.996 MHz
v = 50 m s α = 1 2.073 MHz2.073 MHz
α = 0 . 75 2.157 MHz2.028 MHz
Table 2. The stress wave frequencies for Δ 22 = 3 . 0 M P a , v = 50 m s
Table 2. The stress wave frequencies for Δ 22 = 3 . 0 M P a , v = 50 m s
T m
2.5e-72.5e-62.5e-5
m = 1 α = 1 2.073 MHz2.073MHz2.274 MHz
α = 0 . 75 2.157 MHz2.157 MHz2.288 MHz
m = 2 α = 1 2.073 MHz2.085 MHz2.182 MHz
α = 0 . 75 2.157 MHz2.157 MHz2.207 MHz
m = 3 α = 1 2.073 MHz2.085 MHz2.169 MHz
α = 0 . 75 2.157 MHz2.157 MHz2.182 MHz

Share and Cite

MDPI and ACS Style

Szymczyk, M.; Nowak, M.; Sumelka, W. Numerical Study of Dynamic Properties of Fractional Viscoplasticity Model. Symmetry 2018, 10, 282. https://doi.org/10.3390/sym10070282

AMA Style

Szymczyk M, Nowak M, Sumelka W. Numerical Study of Dynamic Properties of Fractional Viscoplasticity Model. Symmetry. 2018; 10(7):282. https://doi.org/10.3390/sym10070282

Chicago/Turabian Style

Szymczyk, Michał, Marcin Nowak, and Wojciech Sumelka. 2018. "Numerical Study of Dynamic Properties of Fractional Viscoplasticity Model" Symmetry 10, no. 7: 282. https://doi.org/10.3390/sym10070282

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop