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]:
where:
is the stress vector with three normal and three shear stress components and
are constants sufficiently describing properties of the material.
The stress vector can be generally also determined by the three principal stresses
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
. In such a case, the number of stress vector variables can be reduced to three by using stress invariants
[
13].
The plasticity function can be then written in the form:
In the spatial Cartesian coordinate system described by the principal normal stresses
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:
. For the outer area, which is the area of plastic deformations, the following applies:
. For the area of plasticity itself, the interface between elastic and plastic regions, the following applies
.
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:
The plasticity function can be expressed as a dependent variable on the components of deviatoric stress in the following form:
Assuming isotropy of the material, this function is then expressed using the principal values of deviatoric stress:
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).
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
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
and
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]:
The parameter is a proportionality factor that is always non-negative and of indeterminate size. Therefore, vector is also of indeterminate size.
The last relation shows that in the incremental theory of plasticity, the rate (increment) of plastic deformation
is expressed as a function of instantaneous stress, strain rate and stress rate [
40,
41,
42,
43,
44,
45]:
The plastic potential theory is equivalent to the principle of maximum dissipation rate. According to this principle, the actual dissipation rate
corresponding to a given plastic deformation rate
is greater than the fictive dissipation rate
, which is expressed from the actual plastic deformation rate
and a certain stress
on or within the surface (
Figure 4). Thus, the principle of maximum dissipation rate implies [
46,
47,
48,
49]:
The sign of equality applies only in the special case of neutral loading.
The scalar product of vectors
and
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
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
and
. History parameters can be expressed in simpler cases either by dissipation of plastic energy
according to [
13]:
or by an equivalent plastic deformation that is proportional to the second invariant of the plastic deformation velocity deviator
according to [
13]:
The symbol represents deviator.
The hardening condition then defines the hardening function
which determines the subsequent loading area during plastic deformation. The hardening function is expressed as [
13,
55]:
The full time differential of the loading function is expressed by the relation [
13,
14]:
or using Equation (12) according to [
13]:
Components of the deformation vector
or the deformation rate occurring in the elastic–plastic body can be additively decomposed into the elastic portion
and the plastic portion
. Then, the following applies [
13]:
The elastic components of the deformation vector
are coupled to stress vector components
and the flexibility matrix
by the linear Hooke’s law according to [
13]:
At the same time, when using Equation (7) [
13]:
Using Equation (14), the following applies [
13]:
From the relation Equation (17) dependencies (relationships) can be expressed [
13,
14]:
and
From the last relation, it is possible to express the proportionality factor [
13]:
Using Equations (19) and (21) and the elastic–plastic symmetrical matrix
, the stress velocity vector is expressed [
13]:
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]:
where
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]:
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
, the Poisson’s ratio 0.3, the yield stress
, the tensile strength
. 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
in the longitudinal
y-axis direction for isotropic hardening is shown in
Figure 12a, and the stress field
for kinematic hardening is shown in
Figure 12b. At the same time, the locations and magnitudes of extreme values
,
,
and
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
in the longitudinal
-axis direction for isotropic hardening is shown in
Figure 14a, and the field of strains
for kinematic hardening in
Figure 14b. At the same time, the locations and magnitudes of extreme values
,
,
and
of these deformations are identified.
The dependence of the instantaneous maximum normal stress
on the corresponding instantaneous strain
at the respective isotropic hardening point is shown in
Figure 15a, and the dependence
on
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
,
and residual strains
,
due to plastic deformation after the end of loading.
When comparing all corresponding variables in the tensile simulation for both hardening models, their relational dependence was found, namely: , , , , , .
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
according to von Mises theory for isotropic hardening is shown in
Figure 17a, and the field of equivalent stresses
for kinematic hardening in
Figure 17b. At the same time, locations and magnitudes of extreme values
and
of these stresses are identified. The values
and
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
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
for isotropic hardening is shown in
Figure 18a, and the field of equivalent von Mises strains
for kinematic hardening in
Figure 18b. At the same time, the locations and magnitudes of extreme values
and
of these deformations are identified.
As in the case of stress analysis, the values and 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
on the corresponding instantaneous strain
at the respective isotropic hardening point is shown for
in
Figure 19a,c, and the dependence
on
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
,
and residual strains
,
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 in the torsional simulation for both hardening models, their relational dependence was found, which is opposite as in the tensile simulation, namely: , , , .
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
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- Desmorat, R. Non-saturating nonlinear kinematic hardening laws. Comptes Rendus Mécanique 2010, 338, 146–151. [Google Scholar] [CrossRef]
- 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]
- 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]
- Servít, R.; Crha, M.; Doležalová, E. Teorie Pružnosti a Plasticity; SNTL: Praha, Czech Republic, 1981. [Google Scholar]
- Trebuňa, F.; Šimčák, F. Odolnosť Prvkov Mechanických Sústav; TUKE: Košice, Slovakia, 2004. [Google Scholar]
- 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]
- 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]
- 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]
- 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]
- 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]
- Bocko, J.; Segľa, Š. Numerické Metódy Mechaniky Tuhých a Poddajných Telies; TUKE: Košice, Slovakia, 2016. [Google Scholar]
- 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]
- Wang, H.; Xu, T.; Shou, B. Determination of material strengths by hydraulic bulge test. Materials 2017, 10, 23. [Google Scholar] [CrossRef] [Green Version]
- 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]
- 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]
- Božek, P.; Pivarčiová, E. Registration of holographic images based on integral transformation. Comput. Inform. 2012, 31, 1369–1383. [Google Scholar]
- 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]
- 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]
- Kurylo, P. Modernization and optimization of phosphoric cast iron casting. Metals 2019, 9, 1060. [Google Scholar] [CrossRef] [Green Version]
- Tran, T.M.; Le, K.C. Kinematic hardening and Bauschinger effect. PAMM 2019, 19, e201900011. [Google Scholar] [CrossRef]
- Puzrin, A.M.; Houlsby, G.T. Fundamentals of kinematic hardening hyperplasticity. Int. J. Solids Struct. 2001, 38, 3771–3794. [Google Scholar] [CrossRef]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- Lewandowski, D. Modeling of magnetorheological elastomers using the elastic–plastic model with kinematic hardening. Materials 2019, 12, 892. [Google Scholar] [CrossRef] [Green Version]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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].
Figure 2.
Scheme of position adjusting.
Figure 2.
Scheme of position adjusting.
Figure 3.
Limit prism of plasticity.
Figure 3.
Limit prism of plasticity.
Figure 4.
Convexity of the loading area shown on its planar section.
Figure 4.
Convexity of the loading area shown on its planar section.
Figure 5.
Loading areas for an isotropically hardened material.
Figure 5.
Loading areas for an isotropically hardened material.
Figure 6.
Loading areas for a kinematically hardened material.
Figure 6.
Loading areas for a kinematically hardened material.
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.
Figure 8.
Shape and dimensions of the element α.
Figure 8.
Shape and dimensions of the element α.
Figure 9.
Shape and dimensions of the element β.
Figure 9.
Shape and dimensions of the element β.
Figure 10.
Shape and dimensions of the element γ.
Figure 10.
Shape and dimensions of the element γ.
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.
Figure 12.
Field and extreme values of normal stresses. (a) for isotropic hardening. (b) for kinematic hardening.
Figure 12.
Field and extreme values of normal stresses. (a) for isotropic hardening. (b) for kinematic hardening.
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.
Figure 14.
Field and extreme values of strains. (a) for isotropic hardening. (b) for kinematic hardening
Figure 14.
Field and extreme values of strains. (a) for isotropic hardening. (b) for kinematic hardening
Figure 15.
Graphical dependencies. (a) of stress on deformation for isotropic hardening; (b) of stress on deformation for kinematic hardening.
Figure 15.
Graphical dependencies. (a) of stress on deformation for isotropic hardening; (b) of stress on deformation for kinematic hardening.
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.
Figure 17.
Field and extreme values of equivalent stresses. (a) for isotropic hardening. (b) for kinematic hardening.
Figure 17.
Field and extreme values of equivalent stresses. (a) for isotropic hardening. (b) for kinematic hardening.
Figure 18.
Field and extreme values of strains. (a) for isotropic hardening. (b) for kinematic hardening.
Figure 18.
Field and extreme values of strains. (a) for isotropic hardening. (b) for kinematic hardening.
Figure 19.
Graphical dependencies. (a) Of stress on deformation for isotropic hardening of the element α. (b) Of stress on the deformation for kinematic hardening of the element α. (c) Of stress on the deformation for isotropic hardening of the element β. (d) Of stress on the deformation for kinematic hardening of the element β.
Figure 19.
Graphical dependencies. (a) Of stress on deformation for isotropic hardening of the element α. (b) Of stress on the deformation for kinematic hardening of the element α. (c) Of stress on the deformation for isotropic hardening of the element β. (d) Of stress on the deformation for kinematic hardening of the element β.
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).