Next Article in Journal
Damage Evolution Prediction during 2D Scale-Model Tests of a Rubble-Mound Breakwater: A Case Study of Ericeira’s Breakwater
Next Article in Special Issue
Hybrid Finite-Discrete Element Modeling of the Mode I Tensile Response of an Alumina Ceramic
Previous Article in Journal
Efficient Hydrodynamic Modelling of Urban Stormwater Systems for Real-Time Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Empirical Modeling of Transverse Displacements of Single-Sided Transversely Cracked Prismatic Tension Beams

Faculty of Civil Engineering, Transportation Engineering and Architecture, University of Maribor, Smetanova ulica 17, 2000 Maribor, Slovenia
Modelling 2022, 3(4), 481-498; https://doi.org/10.3390/modelling3040031
Submission received: 9 November 2022 / Revised: 8 December 2022 / Accepted: 12 December 2022 / Published: 16 December 2022
(This article belongs to the Special Issue Modeling Dynamic Fracture of Materials)

Abstract

:
While the effects of axial compression on beams have long been known, the effect of tensile axial loads on one-sided transversely cracked beams is less known. The crack namely shifts the position of the resultant of the axial normal stresses deeper into the uncracked part of the cross-section, and the crack tends to open, causing a transverse displacement. Therefore, this paper focuses on empirical modeling of the considered phenomenon for slender prismatic beams in order to establish a suitable 1D computational model based on detailed 3D FE mesh results. This goal can be achieved through the already established simplified model, where the crack is represented by an internal hinge endowed with a rotational spring. Several analyses of various beams differing in geometry, crack locations, and boundary conditions were executed by implementing 3D FE meshes to establish the appropriate model’s bending governing differential equation. After that, the corresponding parameter definitions were calibrated from the database of 3D FE models. By redefining the model’s input parameters, a suitable solution is achieved, offering a good balance between the results’ accuracy and the required computational effort. The functionality of the newly obtained solutions was verified through some comparative case studies that supplement the derivations.

1. Introduction

Beams represent a very widespread type of structural component, although they are actually quite basic structural elements. Although a beam is generally a long structural element of arbitrary cross-section, beams with constant cross-sections are the most common. Consequently, mainstream research is mostly limited to them.
Cracks can have a significant negative impact on structures, and their timely detection is very important to maintain the integrity of structures. Consequently, they have been intensively studied through recent decades. Crack investigations can be roughly divided into the following areas: studies of the effect of cracks on the structure’s response, detection and identification of crack parameters, and development of appropriate computational models of cracked elements.
Measured natural frequencies can be effectively and relatively easily used to detect the presence of cracks. Therefore, many initial studies focused on dynamic analyses to analyze the cracked structures’ dynamic parameters ([1,2,3,4,5]).
Many articles are devoted to the detection of cracks, the identification of crack parameters, and the further monitoring of their growth. They can be divided according to the structural elements (precast concrete members [6]) or structures they consider (dams [7], concrete structures [8] and bridges [9], pavements [10], as well as airplanes [11] and aircraft engine components [12]), or by the techniques and approaches used. Crack detection and identification were initially based on manual visual inspection of the structure. This subjective and scarcely efficient method was further replaced by approaches based on measurements of the natural frequencies of the cracked (and uncracked) structure. Recently, several other approaches have been studied intensively: computer vision techniques [13], CMOS line scan cameras [14], image processing and machine learning [15,16,17]; convolutional neural networks [18,19,20,21,22], multilayer ELM-based feature learning and classification, [23]; and soft elastomeric capacitors (SECs), [24].
Simultaneously with other studies, the development of appropriate computational models of cracked elements that support effective detection and identification was carried out. Although there are some analytical approaches [25,26], research has mainly focused on numerical models, where the finite element method has taken precedence. Different approaches to the derivation of the stiffness matrix and the corresponding load vector due to transverse loads can be found in many references, both for non-cracked uniform straight beams ([27,28]) as well as for cracked multi-stepped beams and beams with linearly-varying heights [29,30].
Axial displacements of straight isotropic beams due to axial loads and transverse displacements due to transverse loads are usually considered completely separately. A usual exception is situations where the influence of axial loads on transverse displacements is considered, which are mostly limited to compressive axial forces, mainly due to problems with buckling of the beams, which requires second-order analysis. However, transverse displacements occur when beams with a one-sided transverse crack (or cracks) are subjected to axial tensile forces [31]. A comparison of the transverse displacement lines for different types of beams showed that they differ in mathematical descriptions, as well as in the presence and distribution of shear force and bending moment. The simplest state occurs for a cantilever and a simply supported beam, as no bending moment occurs and, accordingly, both parts of the deflection line to the left and right of the crack are simple straight lines. However, for other boundary conditions where reactions occur the transverse displacement line is parabolic.
The observed phenomenon was mathematically and empirically modeled with a simple computational model in which a discrete internal moment, MN, was introduced at the crack location [31]. The appropriate internal moment lever ρ was determined empirically and the corresponding stiffness of the rotary spring, KrN, was obtained accordingly. Figure 1 shows the considered cracked beam with a rectangular cross-section of the length L, acting forces, the location L1 of the transverse crack, and the corresponding computational model with studied transverse displacements. Although the model was originally derived for prismatic beams, it was later applied also to beams with a variation linear of width [32].
In this paper, based on response values obtained from a series of analyses of more accurate 3D finite element models, the determinations of the two governing parameters of the simplified model are enhanced, leading to a noticeably better agreement of the results.
The new definitions can be directly incorporated into an existing beam finite element model of a cracked prismatic beam, producing the same results as the solutions of the simplified computational model governing differential equations.

2. Problem Description and Preliminary Studies

2.1. Governing Parameters of the Applied Mathematical Model

The preliminary studies examined different types of beams with a one-sided transverse crack subjected to axial tensile forces [31,32], where a member which experiences simultaneous tension and bending is called a tension beam [33]. The beams considered differed from each other in terms of boundary conditions and geometry (although the height of beams in civil engineering practice is frequently greater than the width, examples can be found in the references where the width is much greater than the height [34]). The database was obtained through analyses of beams modeled with 3D finite element computational models, which allowed accurate crack modeling with a discrete approach. The analyses showed two different mathematical forms of the transverse displacement function. The calculated displacements of simply supported beams and cantilevers apparently followed a perfect linear distribution, while this distribution was parabolic for propped cantilevers.
The cause of transverse displacements due to tensile axial force was attributed to the shift of the resultant of the axial normal stresses in the cracked parts deeper into the non-cracked part of the cross-section at the location of the crack. A pair of equal parallel longitudinal forces Nx acting in opposite directions, but not through the same point, thus produces a local couple effect that tends to open the crack. This effect of opening the crack simultaneously causes transverse displacements along the element, and their magnitudes increase with the crack depth. Mathematically, this crack-opening phenomenon was modeled using the prolonged use of Okamura’s crack model [35]. This widely applied model in which each crack is represented by an internal hinge equipped with a rotational spring that takes into account the residual cross-sectional stiffness has proven to be efficient and reasonably accurate. To be able to use the genuine model to analyze also transverse displacements due to tensile axial forces, a so-called “internal” local bending moment MN was introduced. This moment, which depends both on the applied axial force Nx and on the relative depth δ of the crack, thus represents the first governing parameter of the upgraded model. The type (and consequently the sign) of this moment depends on the position of the crack. A crack on the bottom surface of the beam causes a sagging (or positive) bending moment, while the crack on the top surface causes a hogging (or negative) bending moment. The actual lever (or moment) arm of force ρ·d depends on the actual distribution of normal stresses in the longitudinal direction, which can be obtained adequately only by applying fracture mechanics theory. For simplicity, this arm was initially taken as half the crack depth d (which equals to δ·h, where h is the height of the cross-section), allowing the “internal” sagging moment MN to be simply written as:
M N = 1 2 · d · N x
The stiffness of the rotational spring, denoted by Kr, represents the second main governing parameter of the model. Several definitions for its evaluation can be found in the literature, which were derived for “classical” analyses of transverse displacements due to transverse loads. However, in the initial studies, none of the available definitions gave acceptable results when implemented within the upgraded model due to simplifying the definition of the internal bending moment MN. As a consequence, it was necessary to derive a new appropriate definition of stiffness, denoted as KrN.
Since Okamura’s model with its associated definition has already proven effective for standard bending analyses, the form of the new KrN function was assumed to be similar to the genuine definition. As a result, it was necessary to evaluate only the appropriate function fN(δ) in the denominator:
K r N = E I h · f N ( δ )
where EI and h represent the bending stiffness and height of the beam, respectively.
In order to include both model’s parameters in the analysis of transverse displacement functions v1(x) and v2(x), Figure 1, it is necessary to solve a system of two coupled differential equations of the 4. order, similar to the ordinary analysis of lateral displacements due to transverse loads. The solutions of these equations contain eight unknown integration constants.
Four of them are determined based on the actual boundary conditions of each considered beam case. Three further constants are determined from the crack continuity conditions at the crack location being the same as for the non-cracked condition (simple continuity conditions for displacement, bending moment and shear force).
The remaining constant is obtained from the slope discontinuity at the crack location L1, which now (compared to the authentic Okamura model) also includes the internal moment MN at the crack location in the additional term on the right side of the condition that connects the slopes to the right and left of the crack (φ2(L1) and φ1(L1), respectively):
φ 2 ( L 1 ) φ 1 ( L 1 ) = E I K r N · d 2 v 1 ( L 1 ) d x 2 + ξ ( δ )
where L1 represents crack location, i.e., its distance from the left end of the beam, and the ratio of the governing parameters of the model:
ξ ( δ ) = M N K r N
represents an additional discrete local change in slope.

2.2. Analysis of Transverse Displacements with the Simplified Computational Model

For further comparison with more accurate 3D models in order to identify the governing parameters of the model, GDE solutions for different types of beams were prepared. The following beam types were analyzed: a cantilever, a simply supported beam and a propped cantilever.
The analytical solutions obtained for the transverse displacement function for the cantilever clamped at the left end were:
v 1 ( x ) = 0                       0 x L 1
v 2 ( x ) = ( x L 1 ) · ξ ( δ )                       L 1 x L
The following linear functions followed for a simply supported beam:
v 1 ( x ) = ( L 1 L 1 ) · x · ξ ( δ )                       0 x L 1
v 2 ( x ) = ( x L 1 ) · L 1 · ξ ( δ )                       L 1 x L
Finally, the solutions derived for the propped cantilever were:
v 1 ( x ) = ( L L 1 ) · M N · x 2 · ( x 3 · L ) 2 · K r N · L 3 + 6 · E I · ( L L 1 ) 2                     0 x L 1
v 2 ( x ) = M N · ( L x ) · ( 2 · L 2 · ( L 1 x ) L 1 · x 2 + L · x · ( 2 · L 1 + x ) ) 2 · K r N · L 3 + 6 · E I · ( L L 1 ) 2                     L 1 x L
It can be seen from Equations (5)–(8) that the governing parameters MN and KrN appear indirectly (in a coupled form of a fraction), which does not allow their individual identification, but only the identification of their ratio ξ(δ). However, there are two noteworthy differences in the propped cantilever solutions. Namely, not only do both governing parameters appear in the uncoupled form, but reactions at both ends of the beam also occur, which introduces new parameters into the identification process. The magnitude of the vertical reactions having the opposite directions is given as:
V A = V B = 3 · E I · ( L L 1 ) · M N 3 · E I · ( L L 1 ) 2 + K r N · L 3

3. Enhanced Simplified Computational Model’s Data Acquisition and Calibration

3.1. Database for Determining ξ(δ) Coefficients

To begin with, three prismatic beams, B1, B2, and B3, which differ in their geometric properties, were calculated in detail for different boundary conditions. For each type of beam, analyzes were performed with three different boundary conditions (cantilever clamped at the left end (C), simply supported beam (SS) and propped cantilever clamped at the left end (PC)). Geometric data of these beams (length L, width b and height h), together with the analyzed boundary conditions, are given in Figure 2, Figure 3 and Figure 4, representing nine basic (without detailed information on crack locations) combinations. The Young’s modulus E and the axial load Nx were 30 GPa and 1 MN, respectively. Furthermore, in each beam, transverse cracks were introduced individually at a distance of 0.5 m along each type of beam. At each location, six cracks of different depths were modeled. Relative depths between 0.1 and 0.6 were applied in steps of 0.1. Each crack was modeled on the upper surface of the beams, and as a result, its impact on the beam was reflected by a negative internal moment MN.
Each 3D FE model consisted of 2000 elements in the axial direction, 30 in the vertical direction and 2 in the transverse horizontal direction. A total of 744 such models were analyzed (where almost 600,000 equations were solved in each analysis). The calculated vertical displacement at the crack point of each analyzed structure allowed the computation of the ratio ξ(δ) of the governing parameters of the model. The main results of each analysis were the vertical displacement at the crack location and the reactions in the supports (where they occurred). Ratio values for the same crack depths obtained for different locations were then compared to examine the effect of crack location on the identified ratio ξ(δ). The average results from all analyzed locations for each considered relative crack depth for beam type B1 are given in Table 1 for cantilever clamped at the left end (B1C), simply supported beam (B1SS) and propped cantilever clamped at the left end (B1PC). In the table, the results for the two additionally analyzed beams with modified lengths, a cantilever of 7 m length, B1C7, and a simply supported beam of length 6 m, B1SS6, are also included (they were analyzed only for initial testing of the hypothesis that beam length does not affect the identified values.).
The identified values of ξ(δ) were certainly not ideally the same for each relative depth δ, but their differences were clearly small enough to be easily attributed to the numerical modeling of the structure. Furthermore, since these values were identified at different crack locations L1, their good agreement indicates that they depend solely on the relative depth of the crack δ and not on the location of the crack itself. Therefore, the identification process was repeated for beam type B2, and the obtained average results for all three considered boundary conditions are presented in Table 2.
Although these values were not directly comparable to those from beam B1 due to the different cross-sectional dimensions, they also showed good agreement for each relative crack depth δ analyzed.
Consequently, the analyses were also performed for beam type B3, and the results for all three analyzed boundary conditions are presented in Table 3.
As with the previous beams, the average values of the three structures with different boundary conditions for each relative crack depth were in good agreement.
However, these values were again incomparable with the previously obtained ones due to the different cross-sectional dimensions of the analyzed beams.
Nevertheless, it can be reliably concluded from Table 1, Table 2 and Table 3 that the ratio ξ(δ) clearly depends on the relative depth of the crack δ, but at the same time it is independent of the location of the crack. Yet, in order to enable a direct comparison of the results for different cross-sections, it was necessary to convert the obtained ξ(δ) values to a common basis.

3.2. Conversion of Identified ξ(δ) Values into a Common Comparable Platform

In order to further compare the results for different cross-sections, Equation (4) was transformed into the following form in which the internal moment MN in the numerator was defined by an upgraded definition compared to the simple definition given in Equation (1):
ξ ( δ ) = ρ ( δ ) · δ · h · N x E I h · f N ( δ )
For a rectangular cross-section, the expression takes form:
ξ ( δ ) = 12 · ρ ( δ ) · f N ( δ ) · δ · N x E A
where EA represents the axial stiffness of the beam.
By introducing a new dimensionless variable representing the normalized local slope change:
ξ ¯ ( δ ) = ρ ( δ ) · f N ( δ ) · δ
the following dimensionless relation can be obtained:
ξ ¯ ( δ ) = E A 12 · δ · N x · ξ ( δ )
which now allows a direct comparison of the obtained values for different cross-sections. Therefore, the values from Table 1, Table 2 and Table 3 were transformed accordingly to the common platform and are presented in Table 4. Table 4a–c show results for beam types B1, B2 and B3, respectively.
It can be seen from all the tables that the calculated values for each relative depth δ now show a rational agreement.
Nevertheless, for better insight into the dispersion of the results from all obtained values for each relative crack depth δ, two extreme values (minimum value, ξ ¯ m i n , and maximum value,   ξ ¯ m a x ) were identified from the 127 values calculated for each relative depth. Further, the mean ξ ¯ m e a n and median ξ ¯ m e d i a n values as well as their ratios were also calculated for all relative crack depths, Table 5.
It can be seen from Table 5 that the mean and median values (columns 4 and 5) match almost perfectly, as their ratio (column 6) is almost equal to 1 for all relative crack depths. From the last two columns of Table 5, it is further evident that also the deviations of the extreme values from the mean values (Δn and Δp) are rather minor, as they do not exceed ±1.7%. Therefore, it is obvious that from an engineering point of view, the same ξ ¯ ( δ ) function can be applied to different cross-sections. Its polynomial approximation is thus given as:
ξ ¯ ( δ ) = 0.0595 + 1.3840 · δ 10.3141 · δ 2 + 48.3954 · δ 3 90.0572 · δ 4 + 75.7498 · δ 5       0 δ 0.6
Figure 5 shows the piecewise linear approximation of the ξ ¯ ( δ ) data (black dashed line) and the polynomial approximation derived by Equation (16) (red line).

3.3. Decoupling the Governing Parameters of the Model

As it became apparent that the same functions ρ(δ) and fN(δ) could be used for different cross-sections, their separate evaluation was carried out. The research focus was on data on propped cantilevers, as it was clear from the displacement function of cantilever and simply supported beam (Equations (5)–(8)) that only the ratios ξ(δ) (but not separate individual functions) could be obtained. At first glance, the four parameters (displacement at the crack location, two vertical reactions, and a bending moment reaction) should allow direct calculation of both model parameters independently for each crack location L1. However, basic equilibrium in the vertical direction requires that both vertical reactions be of equal magnitude (but in opposite directions), reducing the number of available parameters by one. Furthermore, the equilibrium of moments further implies that the reaction bending moment is equal to the vertical reaction multiplied by the length L of the beam, which further reduces the number of parameters by one. Therefore, the number of equations available for each crack location was reduced to two, which was still sufficient to determine the two discrete unknowns. The solution for the rotational spring stiffness KrN from which the function fN(δ) could be further found, was thus derived as:
K r N = 3 · E I · ( L L 1 ) 2 · V B 3 · E I · ξ · ( L L 1 ) L 3 · V B
and the corresponding value ρ(δ) could be further evaluated from Equation (4). Therefore, in the first attempt, the analysis for each crack location provided both discrete values using Equations (4) and (17) sequentially. The results for beam type B1PC and relative crack depth of δ = 0.6 are given in the second and third columns of Table 6 (where the value of ξ = 8.650680 × 10−3 was used). The fourth column of the table presents the discrete values of the lever arm function normalized to the value at the first crack.
The last column of the table shows that the mutual matching of the results is quite good (the differences between the results are less than 1%) for the first third of the locations (i.e., up to the first third of the beam length from the clamped end). However, the increase in the discrepancies of the results becomes quite apparent as the crack location L1 approaches the simply supported end. An even greater dispersion of results was evident for shallow cracks, where the zone with noticeable deviations extended further across the complete beam field.
Therefore, in the second attempt, only vertical reactions were used as input data for parameters identification. However, since only one vertical reaction cannot allow identification, vertical reactions from two arbitrary crack locations were used to obtain the two required values for each relative crack depth δ. Thus, each individual reaction was used multiple times, as it could be used in combination with all the remaining reactions for other crack locations of the same relative depth δ. In this way, combining the Nc crack data made it possible to obtain significantly more pairs of solutions than Nc.
Thus, the three considered propped cantilevers of beam types B1, B2 and B3 and their crack locations (35 in overall) gave a total of 188 solution pairs (MN, KrN) for each crack depth δ considered.
The results exclusively for beam type B1PC and a relative crack depth of 0.6 are given in Table 7, where all 55 identified KrN values are presented, ranked in ascending order.
The differences between the results are still noticeable. Nevertheless, a comparison of Table 7 and Table 8 shows that the dispersion of the results in Table 7 is evidently smaller as the deviations of the results from the mean value ranged from −0.902% to 0.778%. Similar, but even slightly smaller discrepancies of range −0.660% and 0.749% were detected for the values of the internal moment MN. Thus, this approach to determining the two governing parameters not only provided a greater variety of results for KrN and MN, but also offered a more stable precision with a mean value that was very close to the median value.
Then, for each type of beam as well as for each relative crack depth δ, the discrete values of the normalized lever arms ρ(δ) were evaluated and further averaged. These values are presented in Table 8.
The weighted mean value for each relative crack depth δ was further used in an interpolation procedure to obtain a polynomial approximation function for ρ(δ).
ρ ( δ ) = 4.23 37.48 · δ + 169 · δ 2 399.65 · δ 3 + 477 · δ 4 226 · δ 5 0 δ 0.6
Figure 6, which proves the usefulness of the derived lever arm function, shows the piecewise linear interpolation of the original identified discrete values in the black dashed line, while the polynomial approximation function for ρ(δ) is plotted in the red line.
With the help of the known functions ξ ¯ ( δ ) and ρ(δ), it was further possible to calculate the missing function fN(δ) from the definition of the stiffness KrN of the rotational spring. The required function follows from Equation (14) in the form:
f N ( δ ) = ξ ¯ ( δ ) ρ ( δ ) · δ
The original function fN(δ) from Equation (19), initially given in the fractional form, was further replaced by a simple uniform polynomial function fNp(δ) for practical reasons only. Figure 7 compares the graph of the original function fN(δ) (black line) with its simple polynomial replacement fNp(δ) (red line). From a comparison of the two graphs, it is evident that the match is generally more than decent, except in the area of very shallow cracks (i.e., δ < 0.1).
For the case without cracks, however, the rotational spring stiffness KrN(δ) must be infinite, which consequently requires the function in the definition denominator to be zero. Therefore, in the derivation of the final function, a zero initial condition was additionally included, which resulted in a slightly modified function fNo(δ):
f N 0 ( δ ) = 0.127 · δ + 11.414 · δ 2 22.026 · δ 3 + 95.961 · δ 4 171.221 · δ 5 + 159.754 · δ 6 0 δ 0.6

4. Numerical Analyses

Initially, a comparison of transverse displacements for different locations L1 and relative crack depths δ was performed. The results for a propped cantilever type B1PC are presented in Table 9. The first value (in regular) for each crack location/relative crack depth pair represents the displacement result (in mm) from the 3D model, and the second value (in italics) was obtained using a simplified model using the newly derived functions. Although it is evident that the directly comparable values are not identical, the matching of the results is completely fine from an engineering point of view, especially after considering the significant difference between the computational efforts of the two models.
In addition, Figure 8 shows in black lines the transverse displacement lines for beam B1PC calculated by Equations (9) and (10) for relative crack depths δ from 0.1 to 0.6 and the crack located at the distance of 2 m from the left, clamped end. Both governing parameters of the model were calculated by implementing the presented solutions.
Further, these solutions are compared against the displacements from the 3D FE models, which are plotted as red dots in Figure 8. It is evident that a very good agreement has been obtained.
In the last step, the reactions on the supports (bending moment MA and vertical reaction VA in the left, clamped support, and vertical reaction VB in the right, simple support) obtained from both models were compared. Accompanying values for beams whose displacement lines are presented in Figure 8 are given in Table 10. The first value (in regular) for each relative crack depth δ represents the displacement result (in mm) from the 3D model, and the second value (in italics) was obtained by implementing a simplified model using the newly derived functions, Equations (18) and (20).
Despite the not-perfect agreement of the values from both computational models, the engineering matching of the results is evident. The average percentage error values for MA, VA, and VB are 0.411%, 1.347%, and −0.437%, respectively. The slightly higher discrepancy value for the vertical reaction VA is mainly influenced by the value −43.92 N (from the 3D FE model) for the shallowest crack, which also shows a visible discrepancy with the value of 45.53 N for the reaction VB (from the same model). However, based on the equilibrium of moments, it can be concluded with certainty that VA is a less accurate value than VB.

5. Conclusions

In the article, the influence of the axial tensile force on the transverse displacements of slender prismatic beams with a single-sided crack was analyzed. A simplified model of a cracked beam was used to calculate transverse displacements and reactions. For this model, where the effect of the axial tensile force is modeled by an internal moment MN acting at the crack location, improved governing parameters have been redefined based on extensive studies.
In the first step, many different beams were analyzed with 3D finite element models, which allow for accurate modeling of cracks. The study covered three different beams differing in geometrical properties, and for each beam type, analyses were made for different boundary conditions. For each analyzed combination of beam type and boundary conditions, transverse cracks were individually modeled at distances of 0.5 m along the entire length of each beam, and 6 cracks of different depths were modeled separately for each location.
Then, for each combination of boundary conditions, an analytical solution of the transverse displacement differential equation for the simplified cracked beam model was obtained. From the analytical expression for the transverse displacement at the crack location, the values of the ratio of the model’s governing parameters were obtained for all analyzed beams. These ratios depended on the associated geometrical data of the beam as well as the depth of the crack. However, after the transformation of the values to a common platform, it became evident that the ratio is truly independent of the beam geometry.
A new universal definition of the lever arm ρ(δ) as a function of relative crack depth δ was then derived. This definition replaced the rather basic definition from the initial studies and consequently allowed more realistic modeling of the internal moment MN.
The rotational spring stiffness KrN definition has also been updated to complete the improved model’s parameters.
The performed numerical tests showed that the model with improved definitions of its parameters can be effectively implemented for structural analyses through solutions of standard bending GDE without any axial force directly included in the governing differential equation. The presented expressions namely lead to reliable engineering results based on values from more thorough (and computationally more demanding) 3D FE models. These models initially served as a database generator and then also for independent verification.
However, although the discussed model represents an effective alternative to 3D FE models, it should be clearly noted that the presented solutions are limited to cracks with a relative crack depth value of up to 0.6.

Funding

This research was partially funded by the Slovenian Research Agency (grant number P2-0129(A)).

Data Availability Statement

Not applicable.

Acknowledgments

The author acknowledges partial financial support from the Slovenian Research Agency (research core funding No. P2-0129 (A) “Development, modelling and optimization of structures and processes in civil engineering and traffic”).

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

SymbolDescription
bbeam width
B1beam No. 1
B2beam No. 2
B3beam No. 3
CCantilever
dcrack depth
fNfunction in the denominator of rotational spring stiffness definition
hbeam height
EYoung modulus
EAaxial stiffness
EIflexural rigidity
Krrotational spring stiffness for the transverse displacement due to pure bending
KrNrotational spring stiffness for the transverse displacement due to axial tensile force
Lbeam length
L1location of the crack i.e., distance from the left end of the beam
MNinternal local bending moment
Nxaxial force
PCpropped cantilever
SSsimply supported beam
v1function of transverse displacements to the left of the crack
v2function of transverse displacements to the right of the crack
xCoordinate
δrelative crack depth
ξlocal change in slope
ξ ¯ m i n identified minimum average value of local change in the slope
ξ ¯ m a x identified maximum average value of local change in the slope
ξ ¯ m e a n mean identified average value of local change in the slope
ξ ¯ m e d i a n median identified average value of local change in the slope
Δn deviation   of   ξ ¯ m i n   a gainst   the   mean   identified   average   value   ξ ¯ m e a n  
Δp deviation   of   ξ ¯ m a x   a gainst   the   mean   identified   average   value   ξ ¯ m e a n
ρnormalized lever arm function
φ1slope to the left of the crack
φ2slope to the right of the crack

References

  1. Shifrin, E.I.; Ruotolo, R. Natural frequencies of a beam with an arbitrary number of cracks. J. Sound Vib. 1999, 222, 409–423. [Google Scholar] [CrossRef]
  2. Khiem, N.T.; Lien, T.V. A simplified method for natural frequency analysis of a multiple cracked beam. J. Sound Vib. 2001, 245, 737–751. [Google Scholar] [CrossRef]
  3. Li, Q.S. Free vibration analysis of non-uniform beams with an arbitrary number of cracks and concentrated masses. J. Sound Vib. 2002, 252, 509–525. [Google Scholar] [CrossRef]
  4. Fernández-Sáez, J.; Navarro, C. Fundamental frequency of cracked beams in bending vibrations: An analytical approach. J. Sound Vib. 2002, 256, 17–31. [Google Scholar] [CrossRef]
  5. Lin, H.P.; Chang, S.C. Forced responses of cracked cantilever beams subjected to a concentrated moving load. Int. J. Mech. Sci. 2006, 48, 1456–1463. [Google Scholar] [CrossRef]
  6. Lee, S.; Jeong, M.; Cho, C.-S.; Park, J.; Kwon, S. Deep Learning-Based PC Member Crack Detection and Quality Inspection Support Technology for the Precise Construction of OSC Projects. Appl. Sci. 2022, 12, 9810. [Google Scholar] [CrossRef]
  7. Feng, C.; Zhang, H.; Wang, H.; Wang, S.; Li, Y. Automatic Pixel-Level Crack Detection on Dam Surface Using Deep Convolutional Network. Sensors 2020, 20, 2069. [Google Scholar] [CrossRef] [Green Version]
  8. Ali, L.; Alnajjar, F.; Jassmi, H.A.; Gocho, M.; Khan, W.; Serhani, M.A. Performance Evaluation of Deep CNN-Based Crack Detection and Localization Techniques for Concrete Structures. Sensors 2021, 21, 1688. [Google Scholar] [CrossRef]
  9. Kim, I.-H.; Jeon, H.; Baek, S.-C.; Hong, W.-H.; Jung, H.-J. Application of Crack Identification Techniques for an Aging Concrete Bridge Inspection Using an Unmanned Aerial Vehicle. Sensors 2018, 18, 1881. [Google Scholar] [CrossRef] [Green Version]
  10. Ren, J.; Zhao, G.; Ma, Y.; Zhao, D.; Liu, T.; Yan, J. Automatic Pavement Crack Detection Fusing Attention Mechanism. Electronics 2022, 11, 3622. [Google Scholar] [CrossRef]
  11. Venugopal, A.; Mohammad, R.; Koslan, M.F.S.; Shafie, A.; Ali, A.b.; Eugene, O. Crack Growth Prediction on Critical Component for Structure Life Extension of Royal Malaysian Air Force (RMAF) Sukhoi Su-30MKM. Metals 2021, 11, 1453. [Google Scholar] [CrossRef]
  12. Rodella, J.; Dhondt, G.; Köster, P.; Sander, M.; Piorun, S. Determination of the Crack Propagation Direction in Mixed-Mode Missions due to Cyclic Loading. Appl. Sci. 2021, 11, 1673. [Google Scholar] [CrossRef]
  13. Hamishebahar, Y.; Guan, H.; So, S.; Jo, J. A Comprehensive Review of Deep Learning-Based Crack Detection Approaches. Appl. Sci. 2022, 12, 1374. [Google Scholar] [CrossRef]
  14. Zhang, W.; Zhang, Z.; Qi, D.; Liu, Y. Automatic Crack Detection and Classification Method for Subway Tunnel Safety Monitoring. Sensors 2014, 14, 19307–19328. [Google Scholar] [CrossRef] [PubMed]
  15. Munawar, H.S.; Hammad, A.W.A.; Haddad, A.; Soares, C.A.P.; Waller, S.T. Image-Based Crack Detection Methods: A Review. Infrastructures 2021, 6, 115. [Google Scholar] [CrossRef]
  16. Golding, V.P.; Gharineiat, Z.; Munawar, H.S.; Ullah, F. Crack Detection in Concrete Structures Using Deep Learning. Sustainability 2022, 14, 8117. [Google Scholar] [CrossRef]
  17. Chen, X.; Li, J.; Huang, S.; Cui, H.; Liu, P.; Sun, Q. An Automatic Concrete Crack-Detection Method Fusing Point Clouds and Images Based on Improved Otsu’s Algorithm. Sensors 2021, 21, 1581. [Google Scholar] [CrossRef]
  18. Yang, C.; Chen, J.; Li, Z.; Huang, Y. Structural Crack Detection and Recognition Based on Deep Learning. Appl. Sci. 2021, 11, 2868. [Google Scholar] [CrossRef]
  19. Hammouch, W.; Chouiekh, C.; Khaissidi, G.; Mrabti, M. Crack Detection and Classification in Moroccan Pavement Using Convolutional Neural Network. Infrastructures 2022, 7, 152. [Google Scholar] [CrossRef]
  20. Teng, S.; Liu, Z.; Chen, G.; Cheng, L. Concrete Crack Detection Based on Well-Known Feature Extractor Model and the YOLO_v2 Network. Appl. Sci. 2021, 11, 813. [Google Scholar] [CrossRef]
  21. Choi, D.; Bell, W.; Kim, D.; Kim, J. UAV-Driven Structural Crack Detection and Location Determination Using Convolutional Neural Networks. Sensors 2021, 21, 2650. [Google Scholar] [CrossRef] [PubMed]
  22. Xu, X.; Zhao, M.; Shi, P.; Ren, R.; He, X.; Wei, X.; Yang, H. Crack Detection and Comparison Study Based on Faster R-CNN and Mask R-CNN. Sensors 2022, 22, 1215. [Google Scholar] [CrossRef]
  23. Wang, B.; Li, Y.; Zhao, W.; Zhang, Z.; Zhang, Y.; Wang, Z. Effective Crack Damage Detection Using Multilayer Sparse Feature Representation and Incremental Extreme Learning Machine. Appl. Sci. 2019, 9, 614. [Google Scholar] [CrossRef] [Green Version]
  24. Yan, J.; Downey, A.; Cancelli, A.; Laflamme, S.; Chen, A.; Li, J.; Ubertini, F. Concrete Crack Detection and Monitoring Using a Capacitive Dense Sensor Array. Sensors 2019, 19, 1843. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Yankelevsky, D.Z.; Karinski, Y.S.; Feldgun, V.R. Analytical Modeling of Crack Widths and Cracking Loads in Structural RC Members. Infrastructures 2022, 7, 40. [Google Scholar] [CrossRef]
  26. Huang, X.; Shi, C.; Ruan, H.; Zhang, Y.; Zhao, W. Stable Crack Propagation Model of Rock Based on Crack Strain. Energies 2022, 15, 1885. [Google Scholar] [CrossRef]
  27. Smith, I.M.; Griffiths, D.V. Programming the Finite Element Method, 3rd ed.; John Wiley & Sons: Chichester, UK, 1997; pp. 21–27. [Google Scholar]
  28. Reddy, J.N. An Introduction to the Finite Element Method, 2nd ed.; McGraw-Hill: Boston, MA, USA, 1993; pp. 143–154. [Google Scholar]
  29. Skrinar, M. Elastic beam finite element with an arbitrary number of transverse cracks. Finite Elem. Anal. Des. 2009, 45, 181–189. [Google Scholar] [CrossRef]
  30. Skrinar, M. Computational analysis of multi-stepped beams and beams with linearly-varying heights implementing closed-form finite element formulation for multi-cracked beam elements. Int. J. Solids Struct. 2013, 50, 2527–2541. [Google Scholar] [CrossRef] [Green Version]
  31. Skrinar, M. Simple model for the vertical displacements computation of single cracked cantilever under tension loads. Z. Angew. Math. Mech. 2000, 80, 551–552. [Google Scholar] [CrossRef]
  32. Skrinar, M. Transverse Displacements of Transversely Cracked Beams with a Linear Variation of Width Due to Axial Tensile Forces. In Engineering Design Applications III, 1st ed.; Öchsner, A., Altenbach, H., Eds.; Springer International Publishing: Cham, Switzerland, 2020; Volume 1, pp. 183–203. [Google Scholar]
  33. Tauchert, T.R. Energy Principles in Structural Mechanics, 1st ed.; McGraw-Hill: Boston, MA, USA, 1974; p. 2. [Google Scholar]
  34. Bowles, J. Foundation Analysis and Design, 5th ed.; McGraw-Hill: Boston, MA, USA, 1996; pp. 518–519. [Google Scholar]
  35. Okamura, H.; Liu, H.W.; Chu, C.-S.; Liebowitz, H. A cracked column under compression. Eng. Fract. Mech. 1969, 1, 547–564. [Google Scholar] [CrossRef]
Figure 1. Transversely cracked tension beam and associated computational model.
Figure 1. Transversely cracked tension beam and associated computational model.
Modelling 03 00031 g001
Figure 2. Tension beam B1 with different boundary conditions: (a) cantilever (B1C), (b) simply supported beam (B1SS), (c) propped cantilever (B1PC).
Figure 2. Tension beam B1 with different boundary conditions: (a) cantilever (B1C), (b) simply supported beam (B1SS), (c) propped cantilever (B1PC).
Modelling 03 00031 g002
Figure 3. Tension beam B2 with different boundary conditions: (a) cantilever (B2C), (b) simply supported beam (B2SS), (c) propped cantilever (B2PC).
Figure 3. Tension beam B2 with different boundary conditions: (a) cantilever (B2C), (b) simply supported beam (B2SS), (c) propped cantilever (B2PC).
Modelling 03 00031 g003
Figure 4. Tension beam B3 with different boundary conditions: (a) cantilever (B3C), (b) simply supported beam (B3SS), (c) propped cantilever (B3PC).
Figure 4. Tension beam B3 with different boundary conditions: (a) cantilever (B3C), (b) simply supported beam (B3SS), (c) propped cantilever (B3PC).
Modelling 03 00031 g004
Figure 5. Comparison of piecewise linear and polynomial approximation of ξ ¯ ( δ ) function.
Figure 5. Comparison of piecewise linear and polynomial approximation of ξ ¯ ( δ ) function.
Modelling 03 00031 g005
Figure 6. Piecewise linear interpolation and polynomial approximation of the ρ(δ) function.
Figure 6. Piecewise linear interpolation and polynomial approximation of the ρ(δ) function.
Modelling 03 00031 g006
Figure 7. Comparison of functions fN(δ) and fNp(δ).
Figure 7. Comparison of functions fN(δ) and fNp(δ).
Modelling 03 00031 g007
Figure 8. Comparison of transverse displacements along the propped cantilever B1PC for various relative crack depths δ.
Figure 8. Comparison of transverse displacements along the propped cantilever B1PC for various relative crack depths δ.
Modelling 03 00031 g008
Table 1. Identified ξ(δ) ratios for the beam type B1 (10−3) at different boundary conditions.
Table 1. Identified ξ(δ) ratios for the beam type B1 (10−3) at different boundary conditions.
δB1CB1C7B1SSB1SS6B1PC
0.10.0790120.0787190.0793590.0791870.079344
0.20.3572660.3565430.3573560.3569920.358737
0.30.9351360.9336910.9334120.9328180.939049
0.42.0380212.0352382.0327072.0328762.046533
0.54.1674784.1619234.1801884.1738134.184832
0.68.5734378.5612198.5861398.5765248.609655
Table 2. Identified ξ(δ) ratios for the beam type B2 (10−3) at different boundary conditions.
Table 2. Identified ξ(δ) ratios for the beam type B2 (10−3) at different boundary conditions.
δB2CB2SSB2PC
0.10.0529960.0533800.053119
0.20.2409560.2414240.241507
0.30.6325490.6316890.634021
0.41.3799401.3759811.383200
0.52.8206352.8203902.827522
0.65.7933665.7901825.806879
Table 3. Identified ξ(δ) ratios for the beam type B3 (10−3) at different boundary conditions.
Table 3. Identified ξ(δ) ratios for the beam type B3 (10−3) at different boundary conditions.
δB3CB3SSB3PC
0.10.0319320.0323650.031999
0.20.1451270.1457760.145436
0.30.3810250.3807160.381835
0.40.8312480.8286090.832963
0.51.6989931.6946811.702461
0.63.4890853.4808943.495967
Table 4. (a). Identified ξ ¯ ( δ ) ratios for the beam type B1 at different boundary conditions. (b). Identified ξ ¯ ( δ ) ratios for the beam type B2 at different boundary conditions. (c). Identified ξ ¯ ( δ ) ratios for the beam type B3 at different boundary conditions.
Table 4. (a). Identified ξ ¯ ( δ ) ratios for the beam type B1 at different boundary conditions. (b). Identified ξ ¯ ( δ ) ratios for the beam type B2 at different boundary conditions. (c). Identified ξ ¯ ( δ ) ratios for the beam type B3 at different boundary conditions.
(a)
δB1CB1C7B1SSB1SS6B1PC
0.10.0158020.0157440.0158720.0158370.015869
0.20.0714530.0713090.0714710.0713980.071747
0.30.1870270.1867380.1866820.1865640.187810
0.40.4076040.4070480.4065410.4065750.409307
0.50.8334960.8323840.8360380.8347630.836966
0.61.7146871.7122441.7172281.7153051.721931
(b)
δB2CB2SSB2PC
0.10.0158990.0160140.015936
0.20.0722870.0724270.072452
0.30.1897650.1895070.190206
0.40.4139820.4127940.414960
0.50.8461910.8461170.848256
0.61.7380101.7370551.742064
(c)
δB3CB3SSB3PC
0.10.0159660.0161830.015999
0.20.0725630.0728880.072718
0.30.1905130.1903580.190917
0.40.4156240.4143040.416481
0.50.8494960.8473400.851231
0.61.7445421.7404471.747983
Table 5. Extreme (min and max), mean, and median values for different relative crack depths.
Table 5. Extreme (min and max), mean, and median values for different relative crack depths.
δ ξ ¯ m i n ξ ¯ m a x ξ ¯ m e a n ξ ¯ m e d i a n ξ ¯ m e a n m e d i a n Δn (%)Δp (%)
0.10.0157440.0161830.0159200.0158991.001346−1.1071.650
0.20.0713090.0728880.0720640.0722870.996929−1.0501.142
0.30.1865640.1909170.1887350.1895070.995928−1.1511.156
0.40.4065410.4164810.4113840.4127940.996583−1.1771.239
0.50.8323850.8512310.8420250.8461170.995164−1.1451.093
0.61.7122441.7479831.7301361.7370550.996017−1.0341.032
Table 6. Identified discrete governing parameters KrN and ρ for δ = 0.6 for propped cantilever B1PC.
Table 6. Identified discrete governing parameters KrN and ρ for δ = 0.6 for propped cantilever B1PC.
L1KrN [×106]ρ(δ) ρ ¯ (δ)
0.506.7831860.488993100.00%
1.006.7604880.48735799.665%
1.256.7557330.48701499.595%
1.506.7494930.48656499.503%
2.006.7151600.48408998.997%
2.506.6602060.48012898.187%
3.006.5575080.47272496.673%
3.506.3457790.45746193.552%
3.756.1568560.44384290.766%
4.005.8328510.42048485.990%
4.503.9926710.28782858.861%
Table 7. All identified values of KrN [×106] for beam type B1PC and a relative crack depth of 0.6.
Table 7. All identified values of KrN [×106] for beam type B1PC and a relative crack depth of 0.6.
6.8109856.8130966.8151886.8465706.853885
6.8563166.8584026.8586806.8629986.864356
6.8653866.8654916.8658546.8659376.866773
6.8698856.8703756.8707386.8716716.871884
6.8721176.8735376.8739806.8742766.874359
6.8747316.8756736.8759296.8763866.876517
6.8772356.8778826.8783926.8784336.878469
6.8791296.8794026.8802056.8804086.880548
6.8807336.8816666.8822246.8823276.882662
6.8829186.8848056.8852646.8852826.885317
6.8855616.8858446.9002516.9197946.926458
Table 8. Identified values of the discrete lever arm ρ(δ) values from the three considered beams.
Table 8. Identified values of the discrete lever arm ρ(δ) values from the three considered beams.
δBeam B1Beam B2Beam B3Mean Value
0.60.4931160.4927070.4926970.492831
0.50.5304530.5291200.5293970.529623
0.40.5973250.5962700.5968600.596784
0.30.7229600.7205360.7209340.721418
0.21.0042430.9849610.9826910.990277
0.11.8127291.7957921.8596381.817743
Table 9. Comparison of the transverse displacement vmax (mm) for beam type B1PC at the crack location from two different computational models.
Table 9. Comparison of the transverse displacement vmax (mm) for beam type B1PC at the crack location from two different computational models.
L1δ = 0.6δ = 0.5δ = 0.4δ = 0.3δ = 0.2δ = 0.1
0.50.363720.216180.110270.0545890.0225810.005462
0.359820.209850.114960.0567630.0227080.005150
1.01.327420.755900.397540.1936850.0775290.017843
1.336970.757970.406900.1981150.0785190.017712
1.52.713231.497590.780350.3751730.1482420.033651
2.737041.508350.794220.3817250.1500180.033677
2.04.271222.291041.176790.5591590.2189570.049336
4.308892.309831.194580.5674750.2213630.049484
2.55.699862.977631.507360.7089280.2756530.061804
5.749193.002731.527970.7184510.2784640.062025
3.06.663173.399991.698410.7919340.3062240.068420
6.719343.428251.720260.8019090.3091610.068659
3.56.821793.413311.685750.7806780.3006360.067024
6.877403.440461.706840.7902000.3033720.067219
4.05.878692.898701.418230.6535610.2510620.055937
5.923712.919371.436190.6615790.2532260.056015
4.53.629361.775050.860370.3959230.1524950.034046
3.653051.783550.872900.4008790.1531610.033847
Table 10. Reactions for the beam type B1PC.
Table 10. Reactions for the beam type B1PC.
δMA [Nm]VA [N]VB [N]
0.1−227.9−43.9245.53
−228.4−45.6845.68
0.2−1017.2−201.8203.4
−1021.7−204.3204.3
0.3−2607.3−517.4521.4
−2619.1−523.8523.8
0.4−5488.8−1095.61097.5
−5513.5−1102.71102.7
0.5−10,607.4−2120.02120.8
−10,660.8−2132.22132.2
0.6−19,811.4−3963.73963.0
−19,887.2−3977.43977.4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Skrinar, M. Empirical Modeling of Transverse Displacements of Single-Sided Transversely Cracked Prismatic Tension Beams. Modelling 2022, 3, 481-498. https://doi.org/10.3390/modelling3040031

AMA Style

Skrinar M. Empirical Modeling of Transverse Displacements of Single-Sided Transversely Cracked Prismatic Tension Beams. Modelling. 2022; 3(4):481-498. https://doi.org/10.3390/modelling3040031

Chicago/Turabian Style

Skrinar, Matjaž. 2022. "Empirical Modeling of Transverse Displacements of Single-Sided Transversely Cracked Prismatic Tension Beams" Modelling 3, no. 4: 481-498. https://doi.org/10.3390/modelling3040031

APA Style

Skrinar, M. (2022). Empirical Modeling of Transverse Displacements of Single-Sided Transversely Cracked Prismatic Tension Beams. Modelling, 3(4), 481-498. https://doi.org/10.3390/modelling3040031

Article Metrics

Back to TopTop