Next Article in Journal
SERS Effect on Spin-Coated Seeding of Tilted Au-ZnO Nanorods for Low-Cost Diagnosis
Next Article in Special Issue
Experimental Study on Influence of Curing Time on Strength Behavior of SLA-Printed Samples Loaded with Different Strain Rates
Previous Article in Journal
Correction: Lara-Banda, M., et al. Alternative to Nitric Acid Passivation of 15-5 and 17-4PH Stainless Steel Using Electrochemical Techniques. Materials 2020, 13, 2836
Previous Article in Special Issue
Magnetorheological Elastomer Stress Relaxation Behaviour during Compression: Experiment and Modelling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Different Strain Hardening Models on the Behavior of Materials in the Elastic–Plastic Regime under Cyclic Loading

Department of Applied Mechanics and Mechanical Engineering, Faculty of Mechanical Engineering, Technical University of Košice, Letná 9, 042 00 Košice, Slovakia
*
Author to whom correspondence should be addressed.
Materials 2020, 13(23), 5323; https://doi.org/10.3390/ma13235323
Submission received: 5 October 2020 / Revised: 11 November 2020 / Accepted: 19 November 2020 / Published: 24 November 2020
(This article belongs to the Special Issue Manufacturing and Mechanics of Materials)

Abstract

:
In exact analyses of bodies in the elastic–plastic regime, the behavior of the material above critical stress values plays a key role. In addition, under cyclic stress, important phenomena to be taken into account are the various types of hardening and the design of the material or structure. In this process, it is important to define several groups of characteristics. These include, for instance, the initial area of plasticity or load which defines the interface between elastic and plastic deformation area. The characteristics also include the relevant law of plastic deformation which specifies the velocity direction of plastic deformation during plastic deformation. In the hardening condition, it is also important to determine the position, size and shape of the subsequent loading area. The elasto-plastic theory was used for the analysis of special compliant mechanisms that are applied for positioning of extremely precise members of the Compact Linear Collider (CLIC), e.g., cryomagnets, laser equipment, etc. Different types of deformation hardening were used to simulate the behavior of particular structural elements in the elastic–plastic regime. Obtained values of stresses and deformations may be used in further practical applications or as default values in other strain hardening model simulations.

1. Introduction

Plastic deformation of a structural element is an irreversible change in the shape and dimensions of the body due to external forces. The stress state in which the body is permanently deformed is the state of elastic–plastic stress [1].
In practical engineering, it is generally inadmissible to create plastic deformation, plastic deformation of a critical size or plastic stability failure. Under uniaxial stress, i.e., under tension, compression or bending, the plastic deformation starts significantly at the yield stress. In other stress modes and in a multiaxial state of stress, the limit state of plastic deformation is determined by a more complex procedure using different plasticity conditions [2,3].
Due to the plastic deformation of components, the geometry is permanently changed. As long as there is a functional relation between these components, their partial or complete malfunction or the overall load-bearing loss might or might not occur [4,5].
However, in certain justified or specific cases, it is possible to consider or allow plastic deformation to occur in the whole structure or in its individual elements. The possibility of allowing plastic deformation in statically certain and especially statically indeterminate tasks (structures, systems) can significantly increase the load-bearing capacity of the structure, in some cases even by several tens of percent.
The knowledge of the materials’ behavior in elastic–plastic regime plays a fundamental role in the analysis, particularly in the case of cyclic stresses where the mutual dependence of stress and deformation is formed by a hysteresis loop. Another phenomenon, which complicates the analyses, is the behavior of structural elements when, after reaching the ultimate limit of plastic load-bearing and relieving the element, repeated loading follows [4]. If the repetitive load is in its original orientation, there is an apparent increase in the yield point, the so-called Bauschinger effect or the phenomenon of deformation hardening.
Many authors are currently dealing with the issue of assessment of structural elements in the elastic–plastic regime, considering some form of deformation hardening of classical or modern materials. A few works can be mentioned, e.g., [5,6,7,8,9,10], in which the authors deal with the influence of isotropic and kinematic reinforced material in the process of sheet metal forming using simulation methods. The influence of kinematic hardening on the properties of sheets during their drawing is discussed, for instance in [6,7]. Propagation of cracks in such strengthened materials is described, for instance in [8,9,10]. Isotropic and kinematic hardening of modern, e.g., composite materials are the topic of, for instance, [11].
In engineering practice, when designing and simulating structural elements loaded under elastic–plastic deformation conditions, a designer must also rely on methodologies and algorithms in the form offered by commercial software systems for the numerical analysis using finite element methods (FEM).
One of the basic analyses is simulating the specific properties and behavior of materials in the elastic–plastic regime under cyclic loading to capture the Bauschinger effect and strengthen the materials by it. The basic mechanisms of hardening used in numerical analyses and simulations include isotropic, kinematic, combined isotropic–kinematic and others. The effect of applying different mechanisms of deformation hardening can be obtained by the comparison of the results of stress and strain analyses of mutually identical or similar tasks (practice).
The analyzed structures are used for very precise placement of cryomagnets and other devices in the Compact Linear Collider (CLIC) in CERN, Figure 1. They are part of approximately 4 m long beams carrying a mass of several hundred kilograms (up to 500 kg). Altogether, several thousand positioning systems with high positioning accuracy are mounted on a length of 2 × 21 km.
Each such positioning system consists of three actuators and other parts of the positioning system. For the active alignment itself, it is necessary to pre-align the positioning systems with an accuracy of 14 µm at a distance of 200 m, relative to a straight line. To meet these strict tolerances, active alignment must be used, where sensors from the alignment systems and control actuators adjust these components to their theoretical position. The minimum effective displacement during adjustment (resolution) is less than 0.5 µm, and the repeatability of the displacement must be less than 1 µm. Due to the above-mentioned facts, a statically indeterminate system forming a special compliance mechanism was used for positioning. The subject of the present work was to analyze the deformation and strength properties of three members of this mechanism.
The components are arranged on the beams. Their external reference surfaces are mounted on V-shaped supports. Their central axis (corresponding to the theoretical axis of movement of the particle beam) is included in a cylinder with a radius of 5 µm for each beam.
The beams are connected by a joint point geometrically defined as the intersection of the central axis of the V-shaped support of two adjacent beams in the median vertical plane. The distances between the two intersections in this plane should remain less than 10 µm, Figure 2.

2. An Analytical Description of Elastic–Plastic Material Behavior

An analytical description of the mechanical behavior is necessary in the analysis of structural elements made of elastic–plastic material with reinforcement [13,14,15,16,17,18]. In particular, it is necessary to define:
  • Initial area of plasticity, which specifies the interface between the elastic and plastic deformation area.
  • Law of plastic transformation, which specifies the direction of plastic deformation velocity during plastic deformation (transformation).
  • Hardening condition, which specifies the position, size and shape of the subsequent load area.

2.1. Initial Area of Plasticity

For the three-dimensional stress state, the criterion characterizing the transition of a material from the elastic to the plastic state can be described by the condition of plasticity, expressed by the function [15,19,20]:
f σ , k i = 0
where:
σ = σ x , σ y , σ z , τ x y , τ x z , τ y z T
is the stress vector with three normal and three shear stress components and k i   are constants sufficiently describing properties of the material.
The stress vector can be generally also determined by the three principal stresses σ 1 , σ 2 , σ 3 and the corresponding three direction angles. When considering the isotropic material in both the elastic and the plastic region, these stresses do not depend on the choice of coordinate axes x , y , z . In such a case, the number of stress vector variables can be reduced to three by using stress invariants I 1 , I 2 , I 3 [13].
The plasticity function can be then written in the form:
f I 1 , I 2 , I 3 , k i   =   0
In the spatial Cartesian coordinate system described by the principal normal stresses σ 1 , σ 2 , σ 3 ,   the plasticity condition can be represented graphically by the so-called flat plasticity [16]. The plasticity area divides the stress space into two areas. For the inner area, which is the area of elastic deformations, the following applies: f < 0 . For the outer area, which is the area of plastic deformations, the following applies: f > 0 . For the area of plasticity itself, the interface between elastic and plastic regions, the following applies f   =   0 .
If this area is fixed, i.e., unchanged relative to the coordinate system, and there is no change in the stresses in space when loading, unloading and reloading, this area is referred to as the initial surface of plasticity [21,22,23].
When considering a material whose plastic state is not influenced by the mean normal stress:
σ m =   σ 1 + σ 2 + σ 3 / 3
The plasticity function can be expressed as a dependent variable on the components of deviatoric stress in the following form:
f s x , s y , s z , τ x y , τ x z , τ y z , k i   =   0
Assuming isotropy of the material, this function is then expressed using the principal values of deviatoric stress:
s 1 , s 2 , s 3 , k i f σ 1 σ m , σ 2 σ m , σ 3 σ m , k i = 0
The particular geometric shape of the plasticity area in the three-dimensional space of principal stresses depends on the theory applied to the plasticity condition. In the case of shear stress stability (von Mises), the plasticity condition is expressed by the limit cylindrical area of plasticity. In the case of maximum shear stress stability (Tresca), the plasticity condition is expressed by the boundary area of the regular hexagonal prism of plasticity (see Figure 3). σ 0 is the relevant critical stress (yield strength of the material). Three spatial limit hexagons of plasticity 1-2, 2-3 and 1-3 as the intersections of the limit prism of plasticity and three coordinate planes are also marked in Figure 3. The axis of the respective area is identical to the so-called hydrostatic axis which passes through the origin of the coordinate system and is equidistant from all three major axes [24,25,26,27,28,29].

2.2. Law of Plastic Transformation

The law of plastic deformation is closely related to the condition of plasticity on the basis of the so-called plastic potential theory. The relevant theory which applies to strengthened materials can be formulated in terms of convexity and normality [30,31,32,33,34].
Depending on the convexity condition, the loading or plasticity area decomposes on one side of its tangent plane and is therefore convex. Depending on the normality condition at the regular points of the loading area, the plastic deformation velocity vector ε ˙ p is normal to the loading area [13,35,36,37,38,39].
The law of plastic transformation is an associative law directly formulated from the condition of normality. It assumes functional identity of loading areas with the function of plastic potential. If vectors ε ˙ p and F f / σ act in the direction of the outer normal to the loading area and at the end point of the vector σ , i.e., at the loading point, then both vectors act on the same beam and the rate of plastic deformation is expressed in the form [13]:
ε ˙ p = F λ ˙ f σ λ ˙
The parameter λ ˙ is a proportionality factor that is always non-negative and of indeterminate size. Therefore, vector ε ˙ p is also of indeterminate size.
The last relation shows that in the incremental theory of plasticity, the rate (increment) of plastic deformation ε ˙ p is expressed as a function of instantaneous stress, strain rate and stress rate [40,41,42,43,44,45]:
ε ˙ p = ε ˙ p σ , ε ˙ , σ ˙
The plastic potential theory is equivalent to the principle of maximum dissipation rate. According to this principle, the actual dissipation rate σ T ε ˙ p corresponding to a given plastic deformation rate ε ˙ p   is greater than the fictive dissipation rate σ + T ε ˙ p , which is expressed from the actual plastic deformation rate ε ˙ p and a certain stress σ + on or within the surface (Figure 4). Thus, the principle of maximum dissipation rate implies [46,47,48,49]:
σ σ + T ε ˙ p 0
The sign of equality applies only in the special case of neutral loading.
The scalar product of vectors σ σ + T and ε ˙ p defines a blunt, i.e., a right or an acute angle between these vectors for any σ + . At the same time, the plastic deformation velocity vector ε ˙ p must be normal to the loading area and the loading area itself must be convex [50,51,52,53,54].

2.3. Condition and Function of Hardening

The reinforcement condition, specifying the position, size and shape of the subsequent loading area, is determined by the dependence on the plastic history π which corresponds to ε p   and k i . History parameters can be expressed in simpler cases either by dissipation of plastic energy W ˙ p l = σ T   ε ˙ p according to [13]:
ϰ W ˙ p l d t = σ T   ε ˙ p d t
or by an equivalent plastic deformation that is proportional to the second invariant of the plastic deformation velocity deviator J 2   d e v s ˙ p   according to [13]:
ϰ 2 3 J 2   d e v s ˙ p 1 / 2 d t
The symbol d e v · represents deviator.
The hardening condition then defines the hardening function H σ , π which determines the subsequent loading area during plastic deformation. The hardening function is expressed as [13,55]:
H σ , π = 1 λ ˙ f ε p T ε ˙ p + f κ i κ ˙ i
The full time differential of the loading function is expressed by the relation [13,14]:
f ˙ = F T σ ˙ + f ε p T   ε ˙ p + f κ i κ ˙ i
or using Equation (12) according to [13]:
f ˙ = F T σ ˙ H λ ˙
Components of the deformation vector ε or the deformation rate occurring in the elastic–plastic body can be additively decomposed into the elastic portion ε e and the plastic portion ε p . Then, the following applies [13]:
ε = ε e + ε p ,   or   ε ˙ = ε ˙ e + ε ˙ p
The elastic components of the deformation vector ε e are coupled to stress vector components σ and the flexibility matrix D by the linear Hooke’s law according to [13]:
ε e = D 1 σ   or     ε ˙ e = D 1 σ ˙
At the same time, when using Equation (7) [13]:
ε ˙ = ε ˙ e + ε ˙ p = D 1 σ ˙ + F λ ˙
Using Equation (14), the following applies [13]:
ε ˙ = D 1 + H 1 F F T σ ˙
From the relation Equation (17) dependencies (relationships) can be expressed [13,14]:
D ε ˙ = σ ˙ + D F λ ˙
and
F T D ε ˙ = H λ ˙ + F T D F λ ˙
From the last relation, it is possible to express the proportionality factor [13]:
λ ˙ = F T D ε ˙ H + F T D F
Using Equations (19) and (21) and the elastic–plastic symmetrical matrix D e p , the stress velocity vector is expressed [13]:
σ ˙ = D D F F T D H + F T D F ε ˙ = D e p ε ˙
The previous relationships determine the general reinforcement condition and the function of hardening which determines the subsequent loading area during plastic deformation. They can be further modified for specific cases to accurately describe individual modifications of material models, e.g., a thermoplastic or elastic–plastic material, either with or without deformation hardening, etc. [56,57,58,59,60]

2.4. Isotropic Hardening

In the case of isotropic hardening, the loading or plasticity area increases in all directions equally and it is geometrically similar to the initial area. For isotropic hardening, e.g., in the case of the von Mises plasticity condition, only the radius of the plasticity area changes (see Figure 5). Generally, for the isotropic state of the loading area, Equation (23) applies [13]:
f σ , Y = f ¯ σ Y = 0
where Y is a scalar variable related to the development of dislocation structures during loading.
The material model with pure isotropic hardening can be used with sufficient accuracy only in the case of monotonic loading in technological operations such as forming. In the case of cyclic loading, this principle of hardening is not sufficiently precise, as it fails to describe the Bauschinger effect (see Figure 5). It also does not correspond to the typical behavior of ductile materials [13]:

2.5. Kinematic Hardening

Kinematic hardening is applied in the cases where the center of the loading or plasticity area is displaced non-symmetrically with respect to the origin of the coordinate system. The size and the shape of the area do not change at all in this case but move as a rigid unit beyond the loading point (see Figure 6) [13].
The internal control variable is the kinematic stress tensor α which defines the actual position of the center of the plasticity area and it is related to the internal stress. Then, in purely kinematic hardening, the plasticity condition is expressed in the form of [36,37,38,39,40,41,42,43,44,45]:
f σ , α = f ¯ σ α σ Y = 0
Purely kinematic hardening can be used for steel whose yield strength is similar to both the static and cyclic deformation curve [35,36,37,38,39].

3. Numerical Simulations of the Effects of Isotropic and Kinematic Hardening

A series of analyses were performed in the SolidWorks programming environment (2012; Dassault Systèmes SolidWorks Corporation, Waltham, MA, USA) to verify the degree of consistency of the simulation results of different methods of material hardening. These analyses were based on the methodology of applying different material models in the conditions of uniaxial or multiaxial stress for basic stress modes [61,62,63,64].
Virtual objects of these analyses were a group of structural elements—universal joints in the form of elastic elements, internally designated as α, β, γ, Figure 7, Figure 8, Figure 9 and Figure 10. Real models of these elements were developed by ZTS VVÚ Košice (Košice, Slovakia) and they are part of a positioning device on the Compact Linear Collider of the European Organization for Nuclear Research (CERN).
Two methods of virtual loading were applied, identical to the assumed stress of real elements using the isotropic and kinematic model of material hardening. The first one included tension above the yield strength followed by relief and the second one included torsion above the yield strength followed by relief. The torsion was done in both directions.
The following study parameters were used:
  • control method: controlled by force increments
  • iterative method: Newton-Raphson method
  • integration method: Newmark method.
Individual simulations were done using isotropic and kinematic models of material strengthening.
In this paper, the elements α and β were chosen for a non-linear static study in the FEM programming environment. The results of some computations were compared.
The elastic members of the positioning system with extremely high positioning accuracy were subjected to numerical simulation. Such positioning systems are intended for positioning cryomagnets, laser devices and other special technologies and devices used in particle accelerators used in particle physics experiments. The solution of the positioning system for positioning compact linear cryomagnets is adapted to the new original concept of two-beam acceleration of the accelerator itself, in which high-frequency power is obtained from a low-energy particle beam, which has a high intensity and is transmitted to a parallel high-energy accelerating beam.
For the individual basic parts of the positioning system, such as the supporting members, the actuating members of the joint, their conceptual structural solutions were designed. The positioning frame includes elastic members, Figure 7, Figure 8, Figure 9 and Figure 10. To verify the required strength and stiffness of the proposed parts, as well as of the whole support system for positioning the cryomagnets, control computations were realized by the finite element method.
Numerical analyses included monitoring of differences or, on the contrary, the degree of conformity of relevant simulation results, which included information on the distribution and magnitude of the maximum equivalent stresses and deformations, the course of mutual dependence of stresses and deformations for the individual hardening models used. The plasticity condition—shear stress intensity constancy (von Mises)—was chosen according to the conclusions of the analyses described in articles [45,46,47,48,49,50,51,52,53].

3.1. Numerical Simulation with Tensile Stress Followed by Unloading

The elements are made of material 34CrNiMo6 with basic material properties: the Young’s modulus 2.1 · 10 6   MPa , the Poisson’s ratio 0.3, the yield stress 900   MPa , the tensile strength 1100   MPa . For this group of simulations, meshing, boundary conditions and load of the virtual object were chosen (Figure 11a,b) with the time-dependent chart of loading according to Figure 11c. In all computations, the size of elements in critical points (e.g., notches) was approximately 0.1 mm. For meshing, parabolic tetrahedral solid elements were used.
This virtual object was then subjected to a numerical simulation using an isotropic and kinematic hardening model. In this way, fields and extremes of stresses and deformations were obtained.
The field of normal stresses σ yI in the longitudinal y-axis direction for isotropic hardening is shown in Figure 12a, and the stress field σ yK for kinematic hardening is shown in Figure 12b. At the same time, the locations and magnitudes of extreme values σ yI max , σ yI min , σ yK max and σ yK min of these stresses are identified.
Another interesting feature is the occurrence of compressive normal stresses. These probably occurred due to complex shape configurations and stress concentrators. Figure 13 identifies stress concentration points of these compressive stresses using the isotropic material hardening model.
The field of strains ε yI in the longitudinal y -axis direction for isotropic hardening is shown in Figure 14a, and the field of strains ε yK for kinematic hardening in Figure 14b. At the same time, the locations and magnitudes of extreme values ε yI max , ε yI min , ε yK max and ε yK min of these deformations are identified.
The dependence of the instantaneous maximum normal stress σ yI max on the corresponding instantaneous strain ε yI   at the respective isotropic hardening point is shown in Figure 15a, and the dependence σ yK max on ε yK   for kinematic hardening in Figure 15b, in both cases for the entire loading cycle. In both dependencies, it is also possible to clearly identify residual normal stresses σ yI res , σ yK res and residual strains ε yI res , ε yK res due to plastic deformation after the end of loading.
When comparing all corresponding variables Q in the tensile simulation for both hardening models, their relational dependence Q I > Q K was found, namely: σ yI max > σ yK max , σ yI min > σ yK min , ε yI max > ε yK max , ε yI min > ε yK min , σ yI res > σ yK res , ε yI res > ε yK res .

3.2. Numerical Simulation with Torsional Stress Followed by Unloading

For this group of simulations, meshing and boundary conditions of the virtual object were chosen according to Figure 16a. A controlled loading process with a loading torque acting around the y-axis (Figure 16a) with both orientations and time course according to Figure 16b was applied.
This virtual object was then subjected to the numerical simulation using an isotropic and kinematic hardening model. In this way, fields and extremes of stresses and deformations were obtained.
The field of equivalent stresses σ VM I according to von Mises theory for isotropic hardening is shown in Figure 17a, and the field of equivalent stresses σ VM K for kinematic hardening in Figure 17b. At the same time, locations and magnitudes of extreme values σ VM I max and σ VM K   max of these stresses are identified. The values σ VM I max and σ VM K max for both isotropic and kinematic hardening were achieved in the second loading cycle, i.e., in the torsion in the opposite direction. However, unlike the isotropic hardening simulation, in the latter case, the maximum values of the equivalent stresses σ VM K max exceeded the (strength) limit of the material. Simulation conditions were the same in both cases. Thus, it can be assumed that if a material model with kinematic hardening is applied, the structural element would be damaged.
The field of equivalent von Mises strains ε eq I for isotropic hardening is shown in Figure 18a, and the field of equivalent von Mises strains ε eq K for kinematic hardening in Figure 18b. At the same time, the locations and magnitudes of extreme values ε eq I max and ε eq K max of these deformations are identified.
As in the case of stress analysis, the values ε eq I max and ε eq K max for both isotropic and kinematic hardening were achieved in the second loading cycle, i.e., in the torsion in the opposite direction.
The dependence of the instantaneous maximum equivalent stress σ VM I max on the corresponding instantaneous strain ε eq I at the respective isotropic hardening point is shown for α , β in Figure 19a,c, and the dependence σ VM K max on ε eq K for kinematic hardening for α, β in Figure 19b,d, in both cases for two complete loading cycles. In both dependencies, it is also possible to clearly identify residual normal stresses σ VM I res , σ VM K res and residual strains ε eq I res , ε eq K res due to the plastic deformation after the end of loading. In addition, at the end of each of the two load cycles, the stresses began to increase slightly again.
When comparing all corresponding variables Q in the torsional simulation for both hardening models, their relational dependence Q I < Q K was found, which is opposite as in the tensile simulation, namely: σ VM I max < σ VM K max , ε eq I max < ε eq K max , σ VM I res < σ VM K res , ε eq I res < ε eq K res .

4. Conclusions

When designing a positioning system for cryomagnet positioning, all proposed parts must be fully functional. To verify the required strength and stiffness of the proposed parts of the supporting system of cryomagnet positioning, it was necessary to perform control computations and analyses using FEM. In the present article, attention was paid to the flexible members (special compliance mechanisms), which are part of the positioning device of the collider magnets at CERN. Simulations were performed for isotropic and kinematic hardening of individual elastic members.
Based on the thorough and detailed examination and comparison of the obtained data, several preliminary partial conclusions can be stated.
In the first case, each member was stressed to pull above the yield point with consequent relief. In the second case, the members were subjected to torsional stresses equally above the yield point with consequent relief, the torsion being applied first to one side and then to the other. Differences between the maximum stresses, displacements and stress-strain curves were observed using individual hardening models. Their difference was on the order of a few percent, in extreme cases even several tens of percent, depending on the quantity considered.
For the individual elastic members, the resulting values of stresses, deformations and their dependence curves for individual types of hardenings were the same, possibly comparable. From the curves of the dependence of the stress on the deformation under torsional load for the elastic element α, a slight difference in the second cycle during unloading is evident, Figure 19a,b. For the elastic element β, it is clear from the dependence that the stresses at the end of both cycles started to rise again, Figure 19c,d. Under torsional stress, the results were specific to the individual elastic members, so it was not possible to draw general conclusions.
The results of the calculations show that the elastic members are suitable for the given load. Moreover, the element stiffnesses were used for further analysis of structural deformations.
However, the above-stated conclusions can be considered preliminary due to the limited number of loading cases. At the same time, they were influenced by specific configurations of design and stress conditions of structural elements without any experimental verification.

Author Contributions

Conceptualization, data curation, methodology, P.S.; conceptualization, formal analysis, project administration, P.F.; formal analysis, writing—original draft, writing—review and editing, I.D.; formal analysis, project administration, validation, J.B.; formal analysis, software, visualization, J.K.; formal analysis, writing—original draft, B.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Ministry of Education of Slovakia Foundation, Grant project VEGA No. 1/0290/18, VEGA No. 1/0500/20 and project KEGA No. 027TUKE-4/2020.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bouhamed, A.; Jrad, H.; Said, L.B.; Wali, M.; Dammak, F. A non-associated anisotropic plasticity model with mixed isotropic–kinematic hardening for finite element simulation of incremental sheet metal forming process. Int. J. Adv. Manuf. Technol. 2019, 100, 929–940. [Google Scholar] [CrossRef]
  2. Vladimirov, I.N.; Pietryga, M.P.; Reese, S. Anisotropic finite elastoplasticity with nonlinear kinematic and isotropic hardening and application to sheet metal forming. Int. J. Plast. 2010, 26, 659–687. [Google Scholar] [CrossRef]
  3. Sága, M.; Kopas, P.; Uhríčík, M. Modeling and experimental analysis of the aluminium alloy fatigue damage in the case of bending–torsion loading. Procedia Eng. 2012, 48, 599–606. [Google Scholar] [CrossRef] [Green Version]
  4. Karvan, P.; Varvani-Farahani, A. Isotropic-kinematic hardening framework to assess ratcheting response of steel samples undergoing asymmetric loading cycles. Fatigue Fract. Eng. Mater. Struct. 2018, 42, 295–306. [Google Scholar] [CrossRef] [Green Version]
  5. Stache, J.M.; Peters, J.F.; Hammi, Y.; Vahedifard, F. A kinematic hardening model based on endochronic theory for complex stress histories. Comput. Geotech. 2019, 114, 103–117. [Google Scholar] [CrossRef]
  6. Varvani-Farahani, A. A comparative study in descriptions of coupled kinematic hardening rules and ratcheting assessment over asymmetric stress cycles. Fatigue Fract. Eng. Mater. Struct. 2017, 40, 882–893. [Google Scholar] [CrossRef]
  7. Juul, K.J.; Martínez-Pañeda, E.; Nielsen, K.L.; Niordson, C.F. Steady-state fracture toughness of elastic-plastic solids: Isotropic versus kinematic hardening. Eng. Fract. Mech. 2019, 207, 254–268. [Google Scholar] [CrossRef] [Green Version]
  8. Martínez-Pañeda, E.; Fleck, N.A. Crack growth resistance in metallic alloys: The role of isotropic versus kinematic hardening. J. Appl. Mech. 2018, 85, 111002. [Google Scholar] [CrossRef] [Green Version]
  9. Badreddine, H.; Saanouni, K.; Labergère, C.; Duval, J.L. Effect of the kinematic hardening on the plastic anisotropy parameters for metallic sheets. Comptes Rendus Mécanique 2018, 346, 678–700. [Google Scholar] [CrossRef]
  10. Desmorat, R. Non-saturating nonlinear kinematic hardening laws. Comptes Rendus Mécanique 2010, 338, 146–151. [Google Scholar] [CrossRef]
  11. Lucchetta, A.; Auslender, F.; Bornert, M.; Kondo, D. A double incremental variational procedure for elastoplastic composites with combined isotropic and linear kinematic hardening. Int. J. Solids Struct. 2019, 158, 243–267. [Google Scholar] [CrossRef]
  12. Kmec, J.; Čáni, D.; Kubák, V. Positioning System for the Positioning of Cryomagnets; Report Req-00169-0003/4/2016; ZTS VVU Košice: Košice, Slovakia, 2016. (In Slovak) [Google Scholar]
  13. Servít, R.; Crha, M.; Doležalová, E. Teorie Pružnosti a Plasticity; SNTL: Praha, Czech Republic, 1981. [Google Scholar]
  14. Trebuňa, F.; Šimčák, F. Odolnosť Prvkov Mechanických Sústav; TUKE: Košice, Slovakia, 2004. [Google Scholar]
  15. Bocko, J.; Delyová, I.; Sivák, P.; Tomko, M. Selection of a significant numerical model of plasticity for the purpose of numerical analysis of plastic reinforcement. Am. J. Mech. Eng. 2017, 5, 334–340. [Google Scholar]
  16. Halama, R.; Markopoulos, A.; Fojtík, F.; Fusek, M.; Poruba, Z.; Famfulík, J. Effect of stress amplitude on uniaxial ratcheting of aluminum alloy 2124-T851: Einfluss der Mittelspannung auf das Ratcheting-Verhalten der Aluminiumlegierung 2124-T851. Mater. Und Werkst. 2017, 48, 814–819. [Google Scholar] [CrossRef]
  17. Halama, R.; Fojtik, F.; Markopoulos, A. Memorization and other transient effects of ST52 steel and its FE description. In Applied Mechanics and Materials; Trans Tech Publications Ltd.: Bäch SZ, Switzerland, 2014; pp. 48–53. [Google Scholar]
  18. Paška, Z.; Halama, R.; Fojtík, F.; Gajdoš, P. Strain response determination in notched specimens under multiaxial cyclic loading by DICM, Experimental Stress Analysis 2017. In Proceedings of the 55th International Scientific Conference on Experimental Stress Analysis 2017, Novy Smokovec, Slovakia, 30 May–1 June 2017; pp. 143–149. [Google Scholar]
  19. Lozhkin, A.; Božek, P.; Lyalin, V.; Tarasov, V.; Tothova, M.; Sultanov, R. Reverse and direct methods for solving the characteristic equation. AIP Conf. Proc. 2016, 1738, 480140. [Google Scholar]
  20. Bocko, J.; Segľa, Š. Numerické Metódy Mechaniky Tuhých a Poddajných Telies; TUKE: Košice, Slovakia, 2016. [Google Scholar]
  21. Borges, M.F.; Antunes, F.V.; Prates, P.A.; Branco, R.; Cruces, A.S.; Lopez-Crespo, P. Effect of kinematic hardening parameters on fatigue crack growth. Theor. Appl. Fract. Mech. 2020, 106, 102501. [Google Scholar] [CrossRef]
  22. Wang, H.; Xu, T.; Shou, B. Determination of material strengths by hydraulic bulge test. Materials 2017, 10, 23. [Google Scholar] [CrossRef] [Green Version]
  23. Reis, L.C.; Rodrigues, C.A.; Oliveira, M.C.; Sakharova, N.A.; Fernandes, J.V. Characterization of the plastic behaviour of sheet metal using the hydraulic bulge test. In Proceedings of the First European Community on Computational Methods in Applied Sciences (ECCOMAS) Young Investigators Conference, Aveiro, Portugal, 24–27 April 2012; Volume 67. [Google Scholar]
  24. Li, Z.; Liu, H. An isotropic-kinematic hardening model for cyclic shakedown and ratcheting of sand. Soil Dyn. Earthq. Eng. 2020, 138, 106329. [Google Scholar] [CrossRef]
  25. Božek, P.; Pivarčiová, E. Registration of holographic images based on integral transformation. Comput. Inform. 2012, 31, 1369–1383. [Google Scholar]
  26. Mercier, S.; Kowalczyk-Gajewska, K.; Czarnota, C. Effective behavior of composites with combined kinematic and isotropic hardening based on additive tangent Mori–Tanaka scheme. Compos. Part B Eng. 2019, 174, 107052. [Google Scholar] [CrossRef] [Green Version]
  27. Garcia, A.; Trinidad, J.; Otegi, N.; Mendiguren, J.; de Argandoña, E.S.; Silvestre, E.; Galdos, L. The influence of the kinematic hardening on the FEM simulation of tension levelling process. Procedia Manuf. 2020, 47, 1381–1386. [Google Scholar] [CrossRef]
  28. Kurylo, P. Modernization and optimization of phosphoric cast iron casting. Metals 2019, 9, 1060. [Google Scholar] [CrossRef] [Green Version]
  29. Tran, T.M.; Le, K.C. Kinematic hardening and Bauschinger effect. PAMM 2019, 19, e201900011. [Google Scholar] [CrossRef]
  30. Puzrin, A.M.; Houlsby, G.T. Fundamentals of kinematic hardening hyperplasticity. Int. J. Solids Struct. 2001, 38, 3771–3794. [Google Scholar] [CrossRef]
  31. Božek, P.; Chmelikova, G. Virtual Technology Utilization in Teaching. In Proceedings of the 2011 14th International Conference on Interactive Collaborative Learning, Piestany, Slovakia, 21–23 September 2011; IEEE: Piscataway, NJ, USA, 2011; pp. 409–413. [Google Scholar]
  32. Younas, N.; Chalal, H.; Abed-Meraim, F. Finite Element Simulation of Sheet Metal Forming Processes using Non-Quadratic Anisotropic Plasticity Models and Solid—Shell Finite Elements. Procedia Manuf. 2020, 47, 1416–1423. [Google Scholar] [CrossRef]
  33. Shchenyatsky, A.; Kotelnikov, M.; Božek, P.; Romanov, A. Analysis of the effect, the conditions of sensing element mounting in the housing have on the stress-strain state and frequency-domain characteristics of products. Mater. Today Proc. 2019, 19, 2568–2570. [Google Scholar] [CrossRef]
  34. Lee, E.H.; Stoughton, T.B.; Yoon, J.W. Kinematic hardening model considering directional hardening response. Int. J. Plast. 2018, 110, 145–165. [Google Scholar] [CrossRef]
  35. Kilin, A.; Božek, P.; Karavaev, Y.; Klekovkin, A.; Shestakov, M. Experimental investigations of a highly maneuverable mobile omniwheel robot. Int. J. Adv. Robot. Syst. 2017, 14, 1–9. [Google Scholar] [CrossRef]
  36. Zhang, Y.; Chen, H.; Jia, Y.F.; Li, D.F.; Yuan, G.J.; Zhang, X.C.; Tu, S.T. A modified kinematic hardening model considering hetero-deformation induced hardening for bimodal structure based on crystal plasticity. Int. J. Mech. Sci. 2020, 191, 106068. [Google Scholar] [CrossRef]
  37. Božek, P.; Lozkin, A.; Gorbushin, A. Geometrical method for increasing precision of machine building parts. Procedia Eng. 2016, 149, 576–580. [Google Scholar] [CrossRef] [Green Version]
  38. Maleki, M.; Pouyan, H. A kinematic hardening based model for unsaturated soils considering different hydraulic conditions. Int. J. Numer. Anal. Methods Geomech. 2016, 40, 2271–2290. [Google Scholar] [CrossRef]
  39. Fan, Q.L.; Xiao, G.F.; Chen, X.D. Analysis of Cyclic Response of Bucket Foundations Based on Simplified Kinematic Hardening Model. Adv. Civ. Eng. 2020, 2020, 9324935. [Google Scholar]
  40. Pivarčiová, E.; Božek, P.; Domnina, K.; Škultéty, E.; Fedosov, S. Interferometric Measurement of Heat Transfer above New Generation Foam Concrete. Meas. Sci. Rev. 2019, 19, 153–160. [Google Scholar] [CrossRef] [Green Version]
  41. Kalentev, E.; Vaclav, S.; Božek, P.; Tarasov, V.; Korshunov, A. Numerical analysis of the stress-strain state of a rope strand with linear contact under tension and torsion loading conditions. Adv. Sci. Technol. Res. J. 2017, 11, 231–239. [Google Scholar] [CrossRef] [Green Version]
  42. Fedosov, S.; Domnina, K.; Pivarčiová, E. The construction of the function of the ultimate goal of the technological process of non-autoclaved foam concrete obtaining. Acta Tecnol. 2019, 5, 59–62. [Google Scholar] [CrossRef]
  43. Frankovský, P.; Delyová, I.; Sivák, P.; Kurylo, P.; Pivarčiová, E.; Neumann, V. Experimental assessment of time-limited operation and rectification of a bridge crane. Materials 2020, 13, 2708. [Google Scholar] [CrossRef]
  44. Meng, L.; Chen, W.; Yan, Y.; Kitamura, T.; Feng, M. Modelling of creep and plasticity deformation considering creep damage and kinematic hardening. Eng. Fract. Mech. 2019, 218, 106582. [Google Scholar] [CrossRef]
  45. Han, S.; Shi, D.; Yang, X.; Huang, J.; Sun, Y. A hypothetical dislocation well model for kinematic hardening in cyclic plasticity. Int. J. Plast. 2018, 110, 220–247. [Google Scholar] [CrossRef]
  46. Lewandowski, D. Modeling of magnetorheological elastomers using the elastic–plastic model with kinematic hardening. Materials 2019, 12, 892. [Google Scholar] [CrossRef] [Green Version]
  47. Kopas, P.; Blatnický, M.; Sága, M.; Vaško, M. Identification of mechanical properties of weld joints of AlMgSi07.F25 aluminium alloy. Metalurgija 2017, 56, 99–102. [Google Scholar]
  48. Bocko, J.; Delyová, I.; Frankovský, P.; Neumann, V. Lifetime assesment of the technological equipment for robotic workplace. Int. J. Appl. Mech. 2020, 12. [Google Scholar] [CrossRef]
  49. Valiollahi, A.; Shojaeifard, M.; Baghani, M. Implementing stretch-based strain energy functions in large coupled axial and torsional deformations of functionally graded cylinder. Int. J. Appl. Mech. 2019, 11, 1950039. [Google Scholar] [CrossRef]
  50. Kráčalík, M. Selection of the crack driving force concept in context of linear-elastic and elastic-plastic fracture mechanics. Acta Mechatron. 2018, 4, 1–6. [Google Scholar] [CrossRef]
  51. Cicero, S.; Fuentes, J.D.; Torabi, A.R. Using the equivalent material concept and the average strain energy density to analyse the fracture behaviour of structural materials. Appl. Sci. 2020, 10, 1601. [Google Scholar] [CrossRef] [Green Version]
  52. Domanski, T.; Sapietova, A.; Saga, M. Application of Abaqus software for the modeling of surface progressive hardening. Procedia Eng. 2017, 177, 64–69. [Google Scholar] [CrossRef]
  53. Wu, S.J.; Chin, P.C.; Liu, H. Measurement of elastic properties of brittle materials by ultrasonic and indentation methods. Appl. Sci. 2019, 9, 2067. [Google Scholar] [CrossRef] [Green Version]
  54. Klačková, I.; Zajačko, I.; Lenhard, R.; Gritsuk, I.; Wiecek, D. Simulation of wood biomass combustion in hot water boiler. Mater. Sci. Eng. 2019, 776, 5. [Google Scholar] [CrossRef]
  55. Kopas, P.; Saga, M.; Baniari, V.; Vasko, M.; Handrik, M. A plastic strain and stress analysis of bending and torsion fatigue specimens in the low-cycle fatigue region using the finite element methods. Procedia Eng. 2017, 177, 526–531. [Google Scholar] [CrossRef]
  56. Akatov, N.; Klačková, I.; Mingaleva, Z.; Galieva, G.; Shaidurova, N. Expert technology for risk management in the implementation of QRM in a high-tech industrial enterprise. Manag. Syst. Prod. Eng. 2019, 27, 250–254. [Google Scholar] [CrossRef] [Green Version]
  57. Zapoměl, J.; Dekýš, V.; Ferfecki, P.; Sapietová, A.; Sága, M.; Žmindák, M. Identification of material damping of a carbon composite bar and study of its effect on attenuation of its transient lateral vibrations. Int. J. Appl. Mech. 2015, 7, 1550081. [Google Scholar] [CrossRef]
  58. Kuryło, P.; Klekiel, T.; Pruszyński, P. A mechanical study of novel additive manufactured modular mandible fracture fixation plates-preliminary study with finite element analysis. Injury 2020, 51, 1527–1535. [Google Scholar]
  59. Beno, M.; Zvoncan, M.; Kovac, M.; Peterka, J. Circular interpolation and positioning accuracy deviation measurement on five axis machine tools with different structures. Tehnicki Vjestnik Technical Gazette 2013, 20, 479–484. [Google Scholar]
  60. Peterka, J.; Pokorny, P.; Vaclav, S. CAM strategies and surface accuracy. In Annals of DAAAM & Proceedings; DAAAM International Vienna: Wolkersdorf, Austria, 2008; pp. 1061–1062. [Google Scholar]
  61. Zajačko, I.; Gál, T.; Ságová, Z.; Mateichyk, V.; Wiecek, D. Application of artificial intelligence principles in mechanical engineering. MATEC Web Conf. 2018, 244, 01027. [Google Scholar] [CrossRef]
  62. Kuric, I. New methods and trends in product development and planning. In Proceedings of the 1st International Conference on Quality and Innovation in Engineering and Management (QIEM), Cluj Napoca, Romania, 17–19 March 2011; pp. 453–456. [Google Scholar]
  63. Kuric, I.; Cisar, M.; Tlach, V. Technical diagnostics at the department of automation and production systems. In Advances in Intelligent Systems and Computing; Springer: Cham, Switzerland, 2019; Volume 835, pp. 474–484. [Google Scholar]
  64. Roszkos, C.S.; Bocko, J.; Kula, T.; Šarloši, J. Static and dynamic analyses of aluminum foam geometric models using the homogenization procedure and the FEA. Compos. Part B Eng. 2019, 171, 361–374. [Google Scholar] [CrossRef]
Figure 1. Detail of the collider assembly [12].
Figure 1. Detail of the collider assembly [12].
Materials 13 05323 g001
Figure 2. Scheme of position adjusting.
Figure 2. Scheme of position adjusting.
Materials 13 05323 g002
Figure 3. Limit prism of plasticity.
Figure 3. Limit prism of plasticity.
Materials 13 05323 g003
Figure 4. Convexity of the loading area shown on its planar section.
Figure 4. Convexity of the loading area shown on its planar section.
Materials 13 05323 g004
Figure 5. Loading areas for an isotropically hardened material.
Figure 5. Loading areas for an isotropically hardened material.
Materials 13 05323 g005
Figure 6. Loading areas for a kinematically hardened material.
Figure 6. Loading areas for a kinematically hardened material.
Materials 13 05323 g006
Figure 7. Detail of the location of the elastic members in the positioning frame.
Figure 7. Detail of the location of the elastic members in the positioning frame.
Materials 13 05323 g007
Figure 8. Shape and dimensions of the element α.
Figure 8. Shape and dimensions of the element α.
Materials 13 05323 g008
Figure 9. Shape and dimensions of the element β.
Figure 9. Shape and dimensions of the element β.
Materials 13 05323 g009
Figure 10. Shape and dimensions of the element γ.
Figure 10. Shape and dimensions of the element γ.
Materials 13 05323 g010
Figure 11. Tensile simulations. (a) Meshing, boundary conditions and load for the element α. (b) Meshing, boundary conditions and load for the element β. (c) Time dependence of the loading force.
Figure 11. Tensile simulations. (a) Meshing, boundary conditions and load for the element α. (b) Meshing, boundary conditions and load for the element β. (c) Time dependence of the loading force.
Materials 13 05323 g011
Figure 12. Field and extreme values of normal stresses. (a) σ yI   for isotropic hardening. (b) σ yK for kinematic hardening.
Figure 12. Field and extreme values of normal stresses. (a) σ yI   for isotropic hardening. (b) σ yK for kinematic hardening.
Materials 13 05323 g012
Figure 13. Identification of concentration points of compressive stresses using the isotropic material hardening model.
Figure 13. Identification of concentration points of compressive stresses using the isotropic material hardening model.
Materials 13 05323 g013
Figure 14. Field and extreme values of strains. (a) ε yI   for isotropic hardening. (b) ε yK for kinematic hardening
Figure 14. Field and extreme values of strains. (a) ε yI   for isotropic hardening. (b) ε yK for kinematic hardening
Materials 13 05323 g014
Figure 15. Graphical dependencies. (a) of stress σ yI max on deformation ε yI for isotropic hardening; (b) of stress σ yK max on deformation ε yK for kinematic hardening.
Figure 15. Graphical dependencies. (a) of stress σ yI max on deformation ε yI for isotropic hardening; (b) of stress σ yK max on deformation ε yK for kinematic hardening.
Materials 13 05323 g015
Figure 16. Torsional simulations. (a) Meshing, boundary conditions and load. (b) Time dependence of the loading torque.
Figure 16. Torsional simulations. (a) Meshing, boundary conditions and load. (b) Time dependence of the loading torque.
Materials 13 05323 g016
Figure 17. Field and extreme values of equivalent stresses. (a) σ VM I for isotropic hardening. (b) σ VM K for kinematic hardening.
Figure 17. Field and extreme values of equivalent stresses. (a) σ VM I for isotropic hardening. (b) σ VM K for kinematic hardening.
Materials 13 05323 g017
Figure 18. Field and extreme values of strains. (a) ε eq I for isotropic hardening. (b) ε eq K for kinematic hardening.
Figure 18. Field and extreme values of strains. (a) ε eq I for isotropic hardening. (b) ε eq K for kinematic hardening.
Materials 13 05323 g018
Figure 19. Graphical dependencies. (a) Of stress σ VM I max on deformation ε eq I for isotropic hardening of the element α. (b) Of stress σ VM K max on the deformation ε eq K for kinematic hardening of the element α. (c) Of stress σ VM I max on the deformation ε eq I for isotropic hardening of the element β. (d) Of stress σ VM K max on the deformation ε eq K for kinematic hardening of the element β.
Figure 19. Graphical dependencies. (a) Of stress σ VM I max on deformation ε eq I for isotropic hardening of the element α. (b) Of stress σ VM K max on the deformation ε eq K for kinematic hardening of the element α. (c) Of stress σ VM I max on the deformation ε eq I for isotropic hardening of the element β. (d) Of stress σ VM K max on the deformation ε eq K for kinematic hardening of the element β.
Materials 13 05323 g019aMaterials 13 05323 g019b
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sivák, P.; Frankovský, P.; Delyová, I.; Bocko, J.; Kostka, J.; Schürger, B. Influence of Different Strain Hardening Models on the Behavior of Materials in the Elastic–Plastic Regime under Cyclic Loading. Materials 2020, 13, 5323. https://doi.org/10.3390/ma13235323

AMA Style

Sivák P, Frankovský P, Delyová I, Bocko J, Kostka J, Schürger B. Influence of Different Strain Hardening Models on the Behavior of Materials in the Elastic–Plastic Regime under Cyclic Loading. Materials. 2020; 13(23):5323. https://doi.org/10.3390/ma13235323

Chicago/Turabian Style

Sivák, Peter, Peter Frankovský, Ingrid Delyová, Jozef Bocko, Ján Kostka, and Barbara Schürger. 2020. "Influence of Different Strain Hardening Models on the Behavior of Materials in the Elastic–Plastic Regime under Cyclic Loading" Materials 13, no. 23: 5323. https://doi.org/10.3390/ma13235323

APA Style

Sivák, P., Frankovský, P., Delyová, I., Bocko, J., Kostka, J., & Schürger, B. (2020). Influence of Different Strain Hardening Models on the Behavior of Materials in the Elastic–Plastic Regime under Cyclic Loading. Materials, 13(23), 5323. https://doi.org/10.3390/ma13235323

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