Next Article in Journal
Rain Area Detection in South-Western Kenya by Using Multispectral Satellite Data from Meteosat Second Generation
Previous Article in Journal
A Crosstalk- and Interferent-Free Dual Electrode Amperometric Biosensor for the Simultaneous Determination of Choline and Phosphocholine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Preisach Elasto-Plastic Model for Mild Steel Hysteretic Behavior-Experimental and Theoretical Considerations

1
Department of Technical Sciences, Civil Engineering, State University of Novi Pazar, 36300 Novi Pazar, Serbia
2
College of Civil and Architecture Engineering, Chuzhou University, Chuzhou 239000, China
3
Jiangxi Province Key Laboratory of Environmental Geotechnical Engineering and Hazards Control, Jiangxi University of Science and Technology, Ganzhou 341000, China
4
Nantong Ocean and Coastal Engineering Research Institute, Hohai University, Nantong 226000, China
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(10), 3546; https://doi.org/10.3390/s21103546
Submission received: 2 April 2021 / Revised: 25 April 2021 / Accepted: 14 May 2021 / Published: 19 May 2021

Abstract

:
The Preisach model already successfully implemented for axial and bending cyclic loading is applied for modeling of the plateau problem for mild steel. It is shown that after the first cycle plateau disappears an extension of the existing Preisach model is needed. Heat dissipation and locked-in energy is calculated due to plastic deformation using the Preisach model. Theoretical results are verified by experiments performed on mild steel S275. The comparison of theoretical and experimental results is evident, showing the capability of the Presicah model in predicting behavior of structures under cyclic loading in the elastoplastic region. The purpose of this paper is to establish a theoretical background for embedded sensors like regenerated fiber Bragg gratings (RFBG) for measurement of strains and temperature in real structures. In addition, the present paper brings a theoretical base for application of nested split-ring resonator (NSRR) probes in measurements of plastic strain in real structures.

1. Introduction

The phenomenon of hysteresis has been a long-standing research interest. It occurs in various branches of physics: magnetism, adsorption, viscoelasticity, plasticity, etc. Its origin starts in 1935 after the paper published by Preisach [1] dealing with the problem of magnetism. After that the Presiach model extended to the problem of adsorption [2]. The first application was in mechanics, for solving problems of cyclic plasticity of axially loaded bars, which was elaborated in [3,4]. After that, the method was extended to elasto-plastic cyclic bending, as described in [5,6]. Low cycle fatigue due to earthquake loading was elaborated in [6]. The problem of the plateau associated with the mild steel is considered in [7].
In this paper theoretically and experimentally all phenomenon of elasto-plastic cyclic axial loading of the mild steel S275 are examined: Lüder dilatation, locked-in energy, dissipated heat, and temperature increase of the surface of specimens due to plastic deformation. For other types of mild, so called construction steels, such are S235 and S355, results can be found in [7].
The motivation of the authors for this paper is to use advantages of the Preisach model to make as simple as possible the calculation of stresses and strains in structures under extremely severe loadings. In addition, this paper has the intention to make connection with the scientists who are investigating reliable embedded sensors, which are beneficial for structural health monitoring (SHM).
To measure strains and temperature in real structures in civil and mechanical engineering and aeronautics, embedded sensors can be used, as explained in [8]. The purpose of this paper is to establish stress–strain relations and construction of cyclic stress–strain loops for arbitrary elements in structures. This is very important for structures resistant to earthquakes or bridges exposed to heavy moving loads. In addition, aircraft structures are usually in the situation that are exposed to cyclic loads during service life. Also, temperature rises due to plastic deformation can be calculated using the algorithm given in this paper. This temperature can be measured with embedded elements as shown in [8]. As stated in [8]: “For the tensile tests, a cyclic procedure was chosen, which allowed us to distinguish between the elastic and plastic deformation of the specimen. An analytical model, which described the elastic part of the tensile test, was introduced and showed good agreement with the measurements.” Readers of the paper are invited to see the explanation how to measure strains (elastic and plastic) given in reference [8] (Section 2.2 and Section 3). It is in detail explained theoretically in Section 2.2 as an application of FBG embedded in cast aluminum during a tensile test. The test was modeled as cyclic uniaxial tension with elastic and plastic strain. The expression connecting the change of the Bragg wavelength due to external strain is derived. Than this result is applied in Paragraph 3 where experimentally obtained results are discussed. During experiments elastic (reversible) and plastic (permanent) strain are measured using FBG with satisfactorily agreement with classical extensometer measurements.
In addition, the present paper establishes a theoretical background, confirmed by experimental investigation, for application of a nested split-ring resonator (NSRR) probe in measurements of plastic strain in real structures [9]. Explanation of the application of this sensor, which was used for the first time to measure elastic and plastic strains in the reinforcing bar (rebar) of a simple supported reinforced concrete beam is explained in [9]: “The system utilizes electromagnetic coupling between the transceiver antenna located outside the beam, and the sensing probes placed on the reinforcing bar (rebar) surface inside the beam. The probes were designed in the form of a nested split-ring resonator, a metamaterial-based structure chosen for its compact size and high sensitivity/resolution, which is at μm/microstrains level. Experiments were performed in both the elastic and plastic deformation cases of steel rebars, and the sensing system was demonstrated to acquire telemetric data in both cases.” In Figure 9c of the paper [9] it can be seen that a plateau in the stress–strain diagram appears. This phenomenon is thoroughly examined in our paper.

2. Theoretical Model

The Preisach model, originally developed in 1935 to describe hysteresis in ferromagnetism, was named after its author [1]. This model represents the mapping of the input data functions into the output data function. It is based on the elementary nonlinear hysteresis operator Gα,β, which is a discontinuous operator with local memory. However, by superposing these operators within the domain Γ, the Preisach hysteresis operator Γ ^ u t is formed as a continuous system of infinitely many elementary operators connected in parallel (or in series):
f t = Γ ^ u t = α β G α , β u t μ α , β d α d β
where μ α , β represents the Preisach weight function, according to which the elementary operators are arranged, and Gα,β is the elementary hysteresis operator shown in Figure 1.
Structural mild steel behavior under monotonic loading is characterized by a phenomenon called the Lüders band, whose main feature is the development of a horizontal yield plateau. Those phenomena occur during the transition from the elastic region to the region of nonlinear plastic hardening.
The values of the upper and lower yield limits, and the Lüder dilatation, depend on a large number of parameters, such as the grain size of the steel [10], the rate of deformation, and the carbon content [11].
This local instability is characteristic only of monotonic loading, while the horizontal yield plateau disappears due to reversible loading.

2.1. Single Crystal Preisach Model

In this paper, a new type of Preisach model for response of structural mild steel under constant cyclic loading will be developed. Its basis is a model that describes the behavior of this type of steel under monotonic loading [7].
The behavior of the steel S275, under monotonic loading, is characterized by an unstable elastoplastic transition which reduces their workability and ductility. This phenomenon is a result of the separation of the free atoms (usually carbon or nickel) and their incorporation at the sites of existing and newly formed dislocations within the lattice formed by atoms of iron [12].
The initial part of the stress–strain diagram of these steel is linear and proportional to the modulus of elasticity E until the stress σuy is reached This limit is called the upper yield stress. Reaching a given value is succeeded by a sudden drop in stress to the value of the lower yield point σly. This is followed by the formation of a yield plateau, with an approximate magnitude of 1–3% of the total dilatation, according to [13]. The cause of this behavior of the mentioned steel lies in the microstructure and atomic lattice. Free carbon and nitrogen atoms surround the dislocation, causing a high level of stress required to move it within the crystal lattice σuy. After initiating the dislocation, its motion continues relatively easily, forming a yield plateau, until a new regrouping of atoms within the crystal lattice occurs.
This phenomenon is called the Lüders band, and the length of the yield plateau is the Lüder dilatation εL. The Lüders band represents a local, inhomogeneous deformation, which causes the development of a yield plateau until the Lüders band spreads over the entire sample, followed by material hardening and homogeneous deformation (Figure 2).
The upper yield stress value σuy, of the mild construction steel, is very sensitive to small stress concentrations, alignment of the sample inside the machine jaws, as well as to other parameters. For this reason, the upper yield strength is neglected, and the value of the stress at which the transition from the elastic to the plastic region occurs is taken to be the value of the lower yield strength σly.
The basic principle in modeling elastoplastic material behavior is based on defining an analog mechanical model, determined by an appropriate set of algebraic and/or differential equations. A material model, describing the response of the subject steel under monotonic axial stress, is constructed by modifying the existing three-element model [3,4,5,6]. By introducing a delay element, the delay of material hardening is achieved after reaching the yield strength Y = σly.
The concavity of the σ-ε curve in the hardening zone is achieved by its approximation with three lines, which makes the working diagram quintuple linear. This is achieved by introducing additional elements into the three linear mechanical model, as shown in the Figure 3.
The material properties of the mechanical model, shown in Figure 3b, are defined by Equation (2).
E = E 0 L 0 + L 1 + L 2 / L 0 E a = E 0 × E 1 / E 0 + E 1 E 1 = h 1 L 0 + L 1 + L 2 / L 1 E h = E 0 × E 1 × E 2 / E 0 + E 1 + E 2 E 2 = h 2 L 0 + L 1 + L 2 / L 2
It is possible to define a new hysteresis mechanical model based on the working diagram shown in Figure 3, which describes the behavior of a structural mild steel single crystal under monotonic loading:
σ ( t ) = E 2 α β δ α β G α , β ε ( t ) d α d β E 2 α β δ α β 2 ε 1 G α , β ε ( t ) d α d β + E a 2 α β δ α β 2 ε 2 G α , β ε ( t ) d α d β + E h E a 2 α β δ α β 2 ε 3 G α , β ε ( t ) d α d β E h 2 α β δ α β 2 ε 4 G α , β ε ( t ) d α d β .
And its appropriate Preisach function:
μ ( α , β ) = 1 2 2 f α , β α β = E 2 × δ ( α β ) δ ( α β 2 ε 1 ) + E a 2 × δ ( α β 2 ε + E h E a 2 δ α β 2 ε 3 E h 2 × δ α β 2 ε 4 .
The formation of a horizontal yield plateau and a phenomenon called the Lüders band represent the characteristic behavior of mild steels under monotonic loading [14]. Under the first load direction switch, a vanishing of the yield plateau occurs, with a complete transformation of the material behavior into an ideal elastoplastic behavior with hardening.
The cyclic behavior of the examined steel types is characterized by the formation of regular hysteresis loops, with no flattening. A new hysteresis model is established by excluding the yield plateau and modifying the material parameters.
A satisfactory approximation of hysteresis loops is possible using a three-element model [3]. In order to better match the real response of the material, a five-element model with a trilinear working diagram, shown in Figure 4, was used in this paper.
The material properties of the mechanical model, shown in Figure 4b, are defined by Equation (2), where substitution should be made:
E a = E a c   a n d   E h = E h c
It is possible to define a new hysteresis mechanical model based on the working diagram shown in Figure 4, which describes the behavior of a structural mild steel single crystal under cyclic loading:
σ c ( t ) = E 2 α β δ α β G α , β ε ( t ) d α d β + E a c E 2 α β δ α β 2 ε 1 c G α , β ε ( t ) d α d β + E h c E a c 2 α β δ α β 2 ε 2 c G α , β ε ( t ) d α d β .
By switching to the single integral, where the integration of the first part on the right hand side of Equation (6) takes part along line β = α , second along line β = α 2 ε 1 c , and third along line β = α 2 ε 2 c Equation (6) for stress becomes:
σ c ( t ) = E 2 ε s ε s G α , α ε ( t ) d α + E a c E 2 2 ε 1 c ε s ε s G α , α 2 ε 1 c ε ( t ) d α + E h c E a c 2 2 ε 2 c ε s ε s G α , α 2 ε 2 c ε ( t ) d α
The first integral of Equation (7) represents the stress due to elastic deformation. The second and third addends describe the behavior of the material after reaching the yield strengths Y 1 c and Y 2 c , respectively.

2.2. Polycrystal Preisach Model

To model the real polycrystalline material, such as mild steel S275, parallel or series connections of infinitely many single crystal elements (Figure 5) are exploited as in [15]. In this paper we are using parallel connections of monocrystalline elements with the range of the yield stress Yimin ≤ Yi ≤ Yimax. Then the stress in the polycrystalline element reads:
σ ( t ) = Y i c , min Y i c , max p ( Y i c ) σ ( Y i c , t ) d Y i c
In Equation (8) p(Yi) is the yield strength probability density function. Assuming that the yield limits Y1 are the same in all parallel-connected individual units and defining that the distribution functions of other Yi values are uniform, as in papers [3,4,5,6,7].
p ( Y i c ) = 1 Y i c , max Y i c , min = c o n s t
The total stress, due to strain as an input, becomes:
σ c ( t ) = E 2 ε s ε s G α , α ε ( t ) d α + E a c E E p Y 1 c Y 1 c , min Y 1 c , max 2 ε 1 c ε s ε s G α , α 2 ε 1 c ε ( t ) d α d Y 1 c + E h c E a c E p Y 2 c Y 2 c , min Y 2 c , max 2 ε 2 c ε s ε s G α , α 2 ε 2 c ε ( t ) d α d Y 2 c
Since the first addend of Equation (10) does not depend on Y i c , on the basis of the second it is α β = 2 ε 1 c , and on the basis of the third α β = 2 ε 2 c . The term β can be reintroduced into the equation, with shifts d β × E / 2 = d Y 1 c , and d β × E a c / 2 = d Y 2 c , where the negative sign of the shift is lost due to the change of integration boundaries within triangles:
σ c ( t ) = E 2 ε s ε s G α , α ε ( t ) d α + E E a c E 4 p Y 1 c A G α , β ε ( t ) d α d β + E a c E h c E a c 4 p Y 2 c B G α , β ε ( t ) d α d β .
The integration domains in Equation (11) represent the areas of the bands between the corresponding lines in a bounded triangle (Figure 6), because the Preisach function exists only in these domains and otherwise is zero. Domain A′ represents the area between the lines α β = 2 ε 1 c , f u l l , and α β = 2 ε 1 c , i n i t , while domain B′ represents the area between the lines α β = 2 ε 2 c , f u l l and α β = 2 ε 2 c , i n i t .
The Preisach function for polycrystalline internal hysteresis loops is defined as:
μ c ( α , β ) = E 2 δ α β + E E a c E 4 H α β 2 ε 1 c , i n i t H α β 2 ε 1 c , f u l l + E a c E h c E a c 4 H α β 2 ε 2 c , i n i t H α β 2 ε 2 c , f u l l .
The geometric interpretation of Equation (12) represents the Preisach triangle shown in Figure 6.

2.3. Wiping Out and Congruency Properties of Preisach Model

Using the geometric interpretation of the Preisach model, the calculation of the double integral in Equation (11) can be avoided. Determination of the integral value comes down to the estimation of the area inside the Preisach triangle. Observing the elementary operators values, it can be stated that at any time, the Preisach triangle consists of points at which elementary operators are at a “switch on” position and points at which they are at a “switch off” position. It can be noticed that the value of the output, at some arbitrary moment, depends on the division of the boundary triangle into positive and negative sets A+ and A−. Increasing and decreasing of the input leads to a change in this redistribution and the formation of a staircase line L(t), which represents the boundary between these sets.
The vertices of the staircase line are the extreme values of the input data so the L(t) line represents the memory of this operator. The memory formation in Preisach model is achieved by changing the shape of the staircase line L(t).
The wiping-out property of Preisach model can be described through its geometric interpretation. This phenomenon characteristic is that each local maximum erases the vertices of the line L(t) whose α coordinates are below that maximum, and each local minimum erases vertices whose β coordinates are above that minimum (Figure 7). Only the alternative values of the dominant extremes of the input data are stored in the Preisach model, while all others are deleted.
The Preisach model, therefore, possesses selective memory. The elementary hysteresis operator (Figure 1) has local memory, but by combining a large number of these operators, non-local model memory is achieved. This memory does not depend on the entire load history, but only on the dominant values of the input extremes, due to the wiping-out property.
Congruence is another fundamental property of the Preisach model, which can be easily demonstrated through the geometric interpretation of the model (Figure 8). If we observe two states, defined by Preisach triangles T1 and T2, in which the input varies in the same range, it can be noticed that the changes in the areas of the triangles are equal in both cases (ΔT1 = ΔT2), regardless of previous load histories. It follows that the output increments for two states with different previous load histories are mutually equal Δf1 = Δf2. The consequence of this property is that all inner loops, created by varying the input between the same dominant extremes, have the same shape and size, but a different position within the main loop, which depends on the previous load history. The complete matching of these loops could be achieved by their translation in the f-axis direction.
The wiping-out and congruence properties represent necessary and sufficient conditions for hysteresis nonlinearity to be defined by the Preisach model [16]. These phenomena are essential for further consideration of energy losses achieved during random load histories.

2.4. Heat Loses

The term hysteresis is mainly related to the appearance of the hysteresis loop. However, hysteresis is also associated with loss of energy manifested through energy losses. Energy losses accomplished during the formation of hysteresis loops were first considered in electromagnetism. The hysteresis energy losses in magnetism are defined by Charles Proteus Steinmetz as the surface of the hysteresis loop [17].
Thanks to this observation, the determination of hysteresis losses at cyclic loads is based on the principle of energy conservation. The analysis of the energy consumed in the formation of hysteresis is important for the field of continuum mechanics.
Determination of energy losses during the formation of hysteresis loops, for an arbitrary history of excitation, can be performed on the basis of Preisach hysteresis operator and its weight function μ(α, β). The general solution was defined by Mayergoyz [16], while the analysis of energy dissipation in the plastic domain for a three-element model of materials is presented in the paper [3].
The overall plastic work Wp, under cyclic loading, is spent on thermal changes S, and the energy remains trapped inside the material [4], called locked-in energy WL.
W p = S + W L s
If the heat losses S are defined using the Preisach hysteresis model, the total loss can be obtained as the sum of the losses in elementary hysteresis operators. The realized loss in one elementary hysteresis operator Gα,β, at full cycle, as shown in Figure 1, is:
s c y c l e = 2 α β
while the energy loss during one input change (switch up or switch down) is:
s = α β
Since the product μ(α, β) × , β can be considered as a rectangular loop with an output value ± μ (α, β), the energy loss in such a loop, during one input change, is defined as the product μ(α, β) (αβ). Whereas the Preisach hysteresis model is a set of a large (infinite) number of elementary hysteresis operators, the total energy losses are obtained by summing the energy losses of all elementary hysteresis operators. This is achieved by integration over the domain Ω where the value of the input within the Preisach triangle has changed. Energy losses at the cross-sectional level are then defined as:
S = Ω μ α , β α β d α d β
The double integral in Equation (16) is determined as the volume with basis area Ω in the αβ plane and height in the z-direction with the value of (αβ). From Equation (16), to find the fraction of the plastic work dissipated into heat, one should integrate the Preisach function for a chosen model within the area of the limiting triangle. Following the procedure explained in [3] for a parallel connection of infinitely many slip and sliding elements connected in series one obtains:
S = 1 3 2 Y min + σ + E × ε p × ε p
For mild steel material constants that are used in the above equation are: E = 210 GPa, Ymin = 16 KN/cm2, σ = Ymax = 2 Ymin = 32 KN/cm2 and εp = 1.25%.
Then from Equation (14) it is obtained for one monotonic cycle:
S = 13,600   kNm / m 3 = 1.36   kNcm / cm 3
One can easily calculate, following the procedure given in [3], the amount of plastic work Wp as:
W p = 0 ε p σ d ε p = 1 3 Y m i n + 2 σ + 1 2 E ε p ε p
Finally, the difference is:
W L = W p S = 1 3 σ Y m i n 1 2 E ε p ε p
This represents the locked-in energy stored during the primary loading from initial, zero state stress to the state of stress Ymin ≤ σ ≤ Ymax.

3. Temperature Field in Cyclically Loaded Cylindrical Specimen in Plastic Region

The Preisach model is capable to determine heat generation in the material during cyclic loading in the plastic region. Based on Equation (16), heat generated in specimen S is given by Equation (17). S is the volumetric source of heat and it is measured in Wh/m3, with the conversation factor: 1 Wh = 3.6 kNm.
From the description of the experiment given in Section 5, it can be seen that the loading, even cyclic, is slow and the strain rate is not considered. The temperature field is in the steady state condition. Let us consider the temperature field within the cylinder. For this purpose we begin with a sketch of a cylindrical section with a radius r and of width Δr, and extended the entire length L of the cylinder as shown in Figure 9 and Figure 10. The energy balance for the cylindrical shell under steady state condition is [18,19]:
Q r + 2 π × r × Δ r × L × S = Q r + Δ r
Rearranging the above equation, dividing by Δr and taking the limit as Δr tends to zero, leads to:
l i m Δ r 0 Q r + Δ r Q r Δ r = 2 π r L S
Or:
d Q d r = 2 π r L S
The heat flow rate can be expressed as:
Q = q r A = q r 2 π r L
Taking Fourier’s law, with the k as conductivity constant,
q r = k d T d r
leads to:
Q = 2 π r L k d T d r
Substituting (23) into (20) one obtains:
1 r d d r r d T d r = S k
Multiplying by r, integrating once with respect to r and finally dividing by r leads to:
d T d r = S 2 k r + C 1 r + C 2
After one more integration the temperature field is:
T r = S 4 k r 2 + c 1 l n r + C 2
In order that the temperature should be of finite value at r = 0 leads to C1 = 0. Then with the integration constant C2 we can find from the boundary condition that for r = R, T(R) = Ts:
C 2 = T s + S 4 k R 2
Substituting (30) into (29) we obtain the final temperature distribution:
T r = T s + S 4 k R 2 r 2
The maximum temperature is at the center of the cylinder (specimen), as is expected for r = 0:
T r = T s + S 4 k R 2
The most important for us is to theoretically determine the surface temperature Ts. Total heat generated within the cylinder, which is the product of heat generated rate per unit volume S and the volume of the cylinder, will be transferred to the surrounding area. On the other hand, the total heat transferred to the outside of the cylinder is the product of the heat flux as convective heat transfer and the surface area of the cylinder:
Q = π R 2 L S = h 0 T s T 0 2 π R L
From the above equation, finally we obtain:
T s = T 0 + R 2 h 0 S
Usual values for the conduction coefficient k for steel are k = 16–24 W/mK, and the convection coefficient h0 = 10–100 W/m2K.
Substituting the value of heat generation given by (18) into (34), taking for h0 = 25 W/m2K, one obtains.
T = 22 + 3.777 = 25.777 °C

4. Preliminary Analysis of Load History

The Preisach model presents a very powerful tool in the analysis of the cyclic behavior of ductile materials and accompanying energy phenomena. However, loading in real structures in civil, mechanical and aerospace engineering is mainly characterized by a random load history. Solving complicated problems is simplified by prior analysis of the load history using some of the methods for counting cycles within a random load history.
Several methods for determining the number of cycles have been developed, where each has found its application in solving a certain type of problem. The purpose of each method is to display random load histories through a set of load histories with constant amplitudes. It should be noted that before applying any cycle counting method, load histories are filtered so that only values of local extremes are retained from real load histories. In general, all existing methods can be classified into three groups [20,21,22]:
Peak counting methods
Range counting methods
Level crossing methods.
The cyclic load’s amplitude is usually a random variable. The rainflow method is the most common method of determining the number of cycles of such a load. Its prevalence and acceptance are reflected in the physical nature of the full cycle, which is represented by a closed hysteresis loop on the stress–strain diagram.
The concept of this method was first given by Matsuishi and Endo [23], where an analogy was made with raindrops, which flow through a row of canopy roofs. This method is developed for counting half-cycles and then pairing them into full cycles (range counting method).
This algorithm requires knowledge of the entire load history, which does not allow immediate data processing and requires large memory for storing the same. To eliminate these shortcomings, several algorithms have been developed, where the most well-known are:
-
three-points algorithm (Figure 11a)
Y = P i + 1 P i X = P i + 2 P i + 1
wherein i = 1, 2, 3, … M − 2, and M is the length of the considered data series, and
-
four-points algorithm (Figure 11b)
Z = P i + 1 P i Y = P i + 2 P i + 1 X = P i + 3 P i + 2
wherein i = 1, 2, 3, … M − 3, and M is the length of the considered data series.
These algorithms and their variations are based on the current analysis of the load history and the extraction of full cycles from it, where the full cycle manifests as a closed hysteresis loop.
Modification of the four-points criterion is possible by adding n points to the calculation algorithm, achieving a higher processing speed of the load history, but also requiring a larger memory for data storage. The newly defined algorithm (Figure 12) is a 4 + n points criterion, which is implemented if the 4 points criterion is not met.
No matter which the rainflow method algorithm is applied, the residuum of the uncounted half cycles inevitably appears. Methods of processing these residues are presented in [24].
The rainflow method is a nonlinear numerical algorithm, not a mathematical function [25]. The use of rainflow algorithms is limited to post-processing the results. Since the algorithms of this method require a series of data (points) and not only instantaneous measurements, this method is not suitable for real-time data processing and is more often used as a post-processing tool.
The Preisach hysteresis operator, as a continuous rate-independent operator, has found its application in the assessment of material damage depending on the number of load cycles. Its purely mathematical form allows its application in real-time, unlike the rainflow method which is limited to subsequent data processing. The results comparison of these two approaches is shown in [26,27], confirming the compatibility between the rainflow method and the Preisach hysteresis operator. The results compatibility of these two different approaches in determining the number of cycles in a random load history is shown by Equation (36) towards estimating fatigue damage.
D a c s = c s α , β N α , β = V a r P α , β
where c(s)(α,β) represents the cycle number determined by the rainflow counting method, N(α,β) is the number of cycles up to failure and Var(P(α,β)) is the variation of the Preisach function up to failure under cycles N(α,β).
The calculation of energy losses S and trapped energy would be greatly simplified by extracting hysteresis loops of the same dimensions from random load histories. The rainflow counting method in principle identifies the ranges, which correspond to closed hysteresis loops.

5. Experimental Results and Model Verification

The experiments presented in the paper were conducted on samples of structural mild steel class S275 at room temperature.
All experimental samples are cylindrical, with a total length of 190 mm and a circular cross-section of the measuring part, with a diameter of d = 10 mm. The dimensions of the test specimens are shown in Figure 13.
The samples were prepared according to the standard [28]. The influence of surface roughness was eliminated by the surface finishing of the sample.
Experimental equipment consisted of the SHIMADZU Servo Pulser machine allowing tension and compression fatigue loading and the extensometer SHOWA-SOKI TCK-1-IF, which was used for strain measurements. The gauge length was 25 mm (Figure 14). All tests were performed under constant velocity of 0.1 Hz and were rate independent.
A total of 4 specimens made of S275 grade structural steel were subjected to symmetrical cyclic stress and symmetrical load of the same range ε = ±1.5, but with different previous histories.
Two test specimens were exposed to the same load regime. A total of two different load histories were applied to the specimens during the tests. Each load history consisted of blocks of 5 symmetric, full cycles (Figure 15). The load application rate was constant in all tests and was 0.1 Hz. The test regime with the number of cycles for all samples subjected to cyclic load is shown in the Table 1.
A graphical representation of two different load histories in the form of a change in dilatation ε as an input as a function of time t is shown in Figure 15.
The given load history B represents a “shortening” of the stress history A, i.e., this load history is formed based on the previous one (A → B) by omitting the first ranges. However, the last, observed ranges in the cyclic load histories are the same. A comparison of hysteresis loops of the same amplitudes, formed at different load histories, is shown in Figure 16.
Experiments have shown remarkable matching of hysteresis loops with same dilatations ± ε(t), independent of the stress history (A or B). A comparison of the numerical model with the results of the experiments is shown in Figure 17. The comparison of results was performed outside the damage zone.
Comparison of the values obtained using the Preisach model defined by Equation (11) and the results of the experiments showed negligible deviations in the zone of nonlinear response. As the deviations are small, it can be concluded that the proposed model provides satisfactory solutions for defining the behavior of structural mild steel exposed to cyclic loading.
During the experiment the surface temperature at the middle of specimen was recorded via a Flir E50bx thermal imaging infrared camera. Temperature increase due to plastic deformation is evident and it is shown in Figure 18. Results obtained in Section 3 and specifically with Equation (35) are within decimal agreement with the recorded one shown in Figure 18.

6. Conclusions

In this paper, the Preisach model of elasto-plastic cyclic loading of a mild construction steel S275 specimen was considered. It is shown that the Lüder dilatation disappears after first cycle. Comparison of theoretical and experimental results are evident. In addition using Preisach model the locked-in energy and heat energy are calculated. After that the temperature field within the sample is obtained. The measurements of temperature during cyclic loading were performed using a thermal imaging infrared camera Flir 50bx. The calculated and measured temperatures are almost the same.
The paper clearly shows the capability of the Presisach model to analyze all aspects and phenomena associated with cyclic loading of structures in the elastoplastic region. The results obtained in this study clearly show that the model is simple enough to solve complicated structural analysis problems of cyclically loaded parts of structures. The results and findings of this paper are also background for scientists dealing with embedded sensors such are FBG, necessary for structural health monitoring of so-called smart structures.
In addition, the experimentalists now have a theoretical basis for measurements of elastic and plastic strains under severe loadings using regenerated fiber Bragg gratings sensors, as well as using a nested split-ring resonator probe.

Author Contributions

Methodology, D.S. and M.C.; validation, D.S., P.K. and C.D.; formal analysis, D.S.; resources, M.C. and D.S.; writing—original draft preparation, D.S.; writing—review and editing, P.K., C.D. and M.C. All authors read, discussed, and commented on the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Key R&D Program of China (No. 2018YFF0214705), the Key R&D Project of Jiangxi Science and Technology Department (No.20192BBH80021), and the Changzhou Science and Technology Program (No. CZ20200003).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

References

  1. Preisach, F. Über die magnetische Nachwirkung. Z. Phys. 1935, 94, 277–302. [Google Scholar] [CrossRef]
  2. Everett, D.H.; Whitton, W.I. A General Approach to Hysteresis. Trans. Faraday Soc. 1952, 48, 749–757. [Google Scholar] [CrossRef]
  3. Lubarda, A.V.; Sumarac, D.; Krajcinovic, D. Hysteretic response of ductile materials subjected to cyclic loads. Recent Adv. Damage Mech. Plast. ASME Publ. AMD 1992, 132, 145–157. [Google Scholar]
  4. Lubarda, A.V.; Sumarac, D.; Krajcinovic, D. Preisach Model and Hysteretic Behavior of Ductile Materials. Eur. J. Mech. A Solids 1993, 12, 445–470. [Google Scholar]
  5. Sumarac, D.; Stosic, S. The Preisach model for the cyclic bending of elasto-plastic beams. Eur. J. Mech. A Solids 1996, 15, 155–172. [Google Scholar]
  6. Šumarac, D.; Petrašković, Z. Hysteretic behavior of rectangular tube (box) sections based on Preisach model. Arch. Appl. Mech. 2012, 82, 1663–1673. [Google Scholar] [CrossRef]
  7. Knežević, P.; Šumarac, D.; Perović, Z.; Dolićanin, Ć.; Burzić, Z. A Preisach Model for Monotonic Tension Response of Structural Mild Steel with Damage. Period. Polytech. Civ. Eng. 2020, 64, 296–303. [Google Scholar] [CrossRef] [Green Version]
  8. Linder, M.; Stadler, A.; Hamann, G.; Fischer, B.; Jakobi, M.; Heilmeier, F.; Bauer, C.; Volk, W.; Koch, A.W.; Roths, J. Fiber Bragg Sensors Embedded in Cast Aluminum Parts: Axial Strain and Temperature Response. Sensors 2021, 21, 1680. [Google Scholar] [CrossRef] [PubMed]
  9. Ozbey, B.; Erturk, V.; Demir, V.H.; Altintas, A.; Kurc, O. A Wireless Passive Sensing System for Displacement/Strain Measurements in Reinforced Concrete Members. Sensors 2016, 16, 496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Zhang, Y.T.; Ao, T.; Jiao, W.; Cui, Y.H. Prediction of the Lüders band in fine grained steel strips under uniaxial tension. Comput. Mater. Sci. 2008, 41, 547–552. [Google Scholar] [CrossRef]
  11. Yoshida, F. A constitutive model of cyclic plasticity. Int. J. Plast. 2000, 16, 359–380. [Google Scholar] [CrossRef]
  12. Raghu, N.; Srinivasan, N.; Venkatraman, B. Study on the Deformation Band Characteristics in Mild Steel Using Digital Image Correlation. J. Multidiscip. Eng. Sci. Technol. JMEST 2014, 1, 400–403. [Google Scholar]
  13. Richard, W.; Richard, P.V.; Hertzberg, J.L. Deformation and Fracture Mechanics of Engineering Materials, 4th ed.; John Wiley & Sons Inc.: Hoboken, NJ, USA, 1996. [Google Scholar]
  14. Shi, Y.; Wang, M.; Wang, Y. Experimental and constitutive model study of structural steel under cyclic loading. J. Constr. Steel Res. 2011, 67, 1185–1197. [Google Scholar] [CrossRef]
  15. Iwan, W.D. On a Class of Models for the Yielding Behavior of Continuous and Composite Systems. J. Appl. Mech. 1967, 34, 612–617. [Google Scholar] [CrossRef]
  16. Mayergoyz, I.D. Mathematical Models of Hysteresis; Springer: Berlin/Heidelberg, Germany, 1991. [Google Scholar]
  17. Steinmetz, C.P. On the law of hysteresis. Proc. IEEE 1984, 72, 197–221. [Google Scholar] [CrossRef]
  18. Mills, A.F. Heat Transfer, 2nd ed.; Prentice-Hall: Upper Saddle River, NJ, USA, 1999. [Google Scholar]
  19. Holman, J.P. Heat Transfer, 10th ed.; McGraw-Hill: New York, NY, USA, 2010. [Google Scholar]
  20. Mercer, C.A.; Livesey, J. Statistical counting methods as a means of analysing the load histories of light bridges. J. Sound Vib. 1973, 27, 399–410. [Google Scholar] [CrossRef]
  21. Schijve, J. The Analysis of Random Load-time Histories with Relation to Fatigue Tests and Life Calculations: This Report Has Been Presented as a Paper to the 2nd. ICAF-AGARD Symposium 1961, Paris; Nationaal Lucht-en Ruimtevaartlaboratorium: Amsterdam, The Netherlands, 1960. [Google Scholar]
  22. Dowling, N.E. Fatigue Failure Predictions for Complicated Stress-Strain Histories. J. Mater. JMLS 1972, 7, 71–87. [Google Scholar] [CrossRef] [Green Version]
  23. Matsuishi, M.; Endo, T. Fatigue of metals subjected to varying stress. Jpn. Soc. Mech. Eng. 1968, 2, 37–40. [Google Scholar]
  24. Marsh, G.; Wignall, C.; Thies, P.R.; Barltrop, N.; Incecik, A.; Venugopal, V.; Johanning, L. Review and application of Rainflow residue processing techniques for accurate fatigue damage estimation. Int. J. Fatigue 2016, 82, 757–765. [Google Scholar] [CrossRef] [Green Version]
  25. Berglind, J.J.B.; Wisniewski, R.; Soltani, M. Fatigue load modeling and control for wind turbines based on hysteresis operators. In Proceedings of the American Control Conference (ACC) 2015, Chicago, IL, USA, 1 July 2015; pp. 3721–3727. [Google Scholar]
  26. Brokate, M.; Dreßler, K.; Krejčí, P. Rainflow counting and energy dissipation for hysteresis models in elastoplasticity. Eur. J. Mech. Solid 1996, 15, 705–737. [Google Scholar]
  27. Brokate, M.; Sprekels, J. Phase Transitions and Hysteresis. In Hysteresis and Phase Transitions; Springer: New York, NY, USA, 1996; pp. 150–174. [Google Scholar]
  28. ASTM. Standard Test Methods for Tension Testing of Metallic Materials; ASTM E8/E8M-21; ASTM International: West Conshohocken, PA, USA, 2021. [Google Scholar]
Figure 1. Elementary hysteresis operator.
Figure 1. Elementary hysteresis operator.
Sensors 21 03546 g001
Figure 2. The Lüders band phenomenon occurring under monotonous loading of mild steel.
Figure 2. The Lüders band phenomenon occurring under monotonous loading of mild steel.
Sensors 21 03546 g002
Figure 3. (a) Quintuple linear working diagram and (b) mechanical model.
Figure 3. (a) Quintuple linear working diagram and (b) mechanical model.
Sensors 21 03546 g003
Figure 4. (a) The stress–strain diagram of structural mild single crystal steel under cyclic loading and (b) mechanical single crystal model.
Figure 4. (a) The stress–strain diagram of structural mild single crystal steel under cyclic loading and (b) mechanical single crystal model.
Sensors 21 03546 g004
Figure 5. Model of polycrystalline material: parallel connection of an infinite number of single elements with yield limits Yic,minYiYic,max.
Figure 5. Model of polycrystalline material: parallel connection of an infinite number of single elements with yield limits Yic,minYiYic,max.
Sensors 21 03546 g005
Figure 6. The Preisach triangle for the material model defined by Equation (11).
Figure 6. The Preisach triangle for the material model defined by Equation (11).
Sensors 21 03546 g006
Figure 7. The wiping-out property of the Preisach model: (a) Load history u1–u9; (b) Load history u1–u′9.
Figure 7. The wiping-out property of the Preisach model: (a) Load history u1–u9; (b) Load history u1–u′9.
Sensors 21 03546 g007
Figure 8. The congruency property of the Preisach model: (a) αβ diagrams for thr input u1; (b) αβ diagrams for thr input u2; (c) The coincidence for loops for input u1 and u2.
Figure 8. The congruency property of the Preisach model: (a) αβ diagrams for thr input u1; (b) αβ diagrams for thr input u2; (c) The coincidence for loops for input u1 and u2.
Sensors 21 03546 g008
Figure 9. Geometry sketch of a cylinder.
Figure 9. Geometry sketch of a cylinder.
Sensors 21 03546 g009
Figure 10. Distribution of flux.
Figure 10. Distribution of flux.
Sensors 21 03546 g010
Figure 11. Rainflow cycle counting method: (a) three-points algorithm and (b) four-points algorithm.
Figure 11. Rainflow cycle counting method: (a) three-points algorithm and (b) four-points algorithm.
Sensors 21 03546 g011
Figure 12. Rainflow method algorithm with the n-points criterion.
Figure 12. Rainflow method algorithm with the n-points criterion.
Sensors 21 03546 g012
Figure 13. Test specimens dimensions.
Figure 13. Test specimens dimensions.
Sensors 21 03546 g013
Figure 14. Test equipment configuration.
Figure 14. Test equipment configuration.
Sensors 21 03546 g014
Figure 15. (a) Load history A; (b) Load history B.
Figure 15. (a) Load history A; (b) Load history B.
Sensors 21 03546 g015
Figure 16. A comparison of hysteresis loops of the same amplitudes ε = ±1.5, formed at different load histories.
Figure 16. A comparison of hysteresis loops of the same amplitudes ε = ±1.5, formed at different load histories.
Sensors 21 03546 g016
Figure 17. Comparison of numerical and experimental results.
Figure 17. Comparison of numerical and experimental results.
Sensors 21 03546 g017
Figure 18. Thermographic recording of the temperature of the specimen.
Figure 18. Thermographic recording of the temperature of the specimen.
Sensors 21 03546 g018
Table 1. Sample testing regimes.
Table 1. Sample testing regimes.
ε [%]±0.5±1±1.5
Number of cyclesIII-3555A.
III-4555
III-5/55B.
III-6/55
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sumarac, D.; Knezevic, P.; Dolicanin, C.; Cao, M. Preisach Elasto-Plastic Model for Mild Steel Hysteretic Behavior-Experimental and Theoretical Considerations. Sensors 2021, 21, 3546. https://doi.org/10.3390/s21103546

AMA Style

Sumarac D, Knezevic P, Dolicanin C, Cao M. Preisach Elasto-Plastic Model for Mild Steel Hysteretic Behavior-Experimental and Theoretical Considerations. Sensors. 2021; 21(10):3546. https://doi.org/10.3390/s21103546

Chicago/Turabian Style

Sumarac, Dragoslav, Petar Knezevic, Cemal Dolicanin, and Maosen Cao. 2021. "Preisach Elasto-Plastic Model for Mild Steel Hysteretic Behavior-Experimental and Theoretical Considerations" Sensors 21, no. 10: 3546. https://doi.org/10.3390/s21103546

APA Style

Sumarac, D., Knezevic, P., Dolicanin, C., & Cao, M. (2021). Preisach Elasto-Plastic Model for Mild Steel Hysteretic Behavior-Experimental and Theoretical Considerations. Sensors, 21(10), 3546. https://doi.org/10.3390/s21103546

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