Modeling and Simulation of a Cracked Beam with Different Location Using FEM

Nowadays presence of crack in different engineering structures becomes a se-rious threat to the performance. Since most of the civil and mechanical structures may be damaged due to material fatigue, mechanical vibration, environmental attack and long-term service. Moreover, dynamical systems of a beam usually possess a non-linear character, which causes practical difficul-ties on the model-based damage detection techniques. This paper presents a novel approach to detect damage in a simply supported beam. In this study, a numerical simulation using the Finite Element Method (FEM) has been done to determine the frequencies to detect the crack in a concrete beam of length 0.12 m and width 0.015 m. A vibration-based model is employed to simulate the results by using COMSOL Multiphysics. At the tip, by performing the computational analysis it is found that the presence of cracks affects the natural frequencies of the concrete structure. It is observed that after applying load, the frequencies of the cracked beam have been changed with the varia-tion of the location of the crack for all the modes of vibration. It also found that maximum frequency reserved at the cracked point so it will also help us to detect different hidden defects in any structure. A comparison is also made with the experimental results. It is also found that the effect of crack is more near the fixed end than at the free end.


Introduction
Identification of structural crack location has gained increasing attentions from the scientific and engineering domains since the unpredicted structural failure may cause catastrophic, economic, and life loss [1]. So early detection of these defects seems pivotal for safety and economic reasons, as their detection can remarkably extend the existence of the structure, which increases its reliability at the same time [2]. The most important technique of structural health monitoring is to provide information on the anticipation of structures, simultaneously detects and locates structural fatigue. This needs ideas of the model of structures in great detail, which may not be possible all-time. As we know when a structure suffers from damage, its dynamic property (like as float, string or object) can change and it is observed that crack caused a stiffness reduction in modal natural frequencies. Consequently it results in the change within the dynamic response of the beam [3]. Within the recent years, with the advancement of science and technology, the study of detecting the crack has been increasing among the researchers. Being very commonly utilized in steel construction and machinery industries, health monitoring and thus the analysis of injury within the type of crack within the Beam structures pose a crucial role [4]. Since long efforts are on their to get a feasible solution for crack detection in beam structures, during this regard many approaches have so far being considered. Since there always exist constraints in ferroconcrete structures and buildings deteriorate overtime, cracking seems unavoidable and appears altogether kinds of structures, as an example, concrete wall, and beam, slab, and brick walls. The presence of the crack not only changes the regional stress and strain fields of the crack tip but also effects on structural dynamics [5]. When a structure suffers from damages, its dynamic properties can change. If these fatigues cannot be timely detected and repaired, subsequent fracture can bring catastrophic failure to the beam structures.
Many extensive researches on crack detection in different methods have been performed theoretically and experimentally till today. Kim and Zhao [6] studied a very distinctive crack detection technique employing a harmonic response, where the displacement and slope modes of a cracked cantilever beam are thought of 1st, that the approximate formula for displacement and slope response under single-point harmonic excitation comes. They conclude that the slope response incorporates a pointy amendment with the crack location and additionally the depth of the crack. A new idea has developed by Lee and Chung [7] that how to look out the lowest four natural frequencies of the cracked structure by FEM and additionally the approximate crack location is obtained by exploitation. Later, Owolabi et al. [8] have experimentally investigated of the possessions of cracks and damages on the dependability of structures. An analytical and experimental approach for the fault detection in cantilever beams like structure by vibration analysis has developed by Nahvi [9]. Presently numerous analytical, numerical and experimental techniques are in use for crack detection throughout a fiber-reinforced composite, laminated composites and non-composite structures for its vibration analysis. A numerous vibration based fatigue designation techniques bestowed by Sanjay K. Behera et al. [10] researches for cracked composite and non-composite structures. In the same case I. Goda [11], has applied numerical study exploitation finite parts is performed to research the free vibration response of laminated composite beams. They had a tendency to perform dynamic modeling of the laminated beams by associate eigenvalue analysis, exploitation associate eight-node bedded shell part to simulate the free vibrations. The major importance of their study was for the mechanical designer to conceive and optimize composite structures subjected to dynamic loadings [2].
Prasad et al. [12] discussed the effect of location of crack from the free end to the fixed end on crack growth rate along vibrating cantilever beam and a mathematical model was developed using dimensional analysis to find out the value of the crack growth rate along vibrating cantilever beam.
The above papers are studied the changes neighborhood in natural frequencies due to the presence of multiple cracks at intervals associate concrete beam. The target of this paper to look out a method for predicting crack parameters (crack depth associated its location) throughout the concrete beam from changes in natural frequencies. Constant studies are disbursed by using COMSOL Multiphysics software to judge modal parameters (natural frequencies and deflections) for various crack position parameters. The technique developed to detect crack by using the FEM software has become popular in recent years.

Mathematical Modeling
The vibration-based damage detection has become one of the commonly used tools for crack detection. This approach is mainly based on changes in dynamic characteristics, such as natural frequency and crack position parameter [2]. In this study we consider a Computational domain of a beam, where the dimensions of the beam and crack dimensions are taken from Yamuna and Sambasivarao [13].

Governing Equation
The equation of motion in matrix form for vibration of a beam under load is given by, (2) Computational Water, Energy, and Environmental Engineering In-plane, the load P(t) can be expressed in the form as shown below, where, P s = the static portion of P. P t = the amplitude of the dynamic portion of P and Ω = the frequency of excitation. Equation (2) represents an eigen value problem and the roots of the equation give rise to square of the natural frequency given by the equation,

Finite Element Analysis
In this analysis we consider a beam with two degrees of freedom (slope and deflection) per node [14] which is shown in Figure 1.
The governing equations of a deformation curved beams can be simplified by determining the force, moment, deflection and twist along the fifth metatarsal. This can be done subjected to both a point wise and a distributed load by using Young's Modulus and Moments of Inertia given by the following equation, where, E = Young's Modulus. I = Moments of Inertia.
The deflection of the beam due to applied loads on the domain are shown in Figure 2. Now, If we write these values in the matrix form of the above system then we find, i.e. The Element stiffness matrix due to bending,

[ ] [ ] [ ][ ]
The stiffness matrix K crack or K c of a cracked beam element: From equilibrium condition as in Figure 3. Where, Hence the stiffness matrix K crack or K c of a cracked beam element can be obtained, The cracked element stiffness matrix becomes,   21  22  21  2  22  21  22  21  12  11  22  12  21  11   22  22  21  22  21  11 22  12 21   12  12  11  12  11 1 e e e e e e e e e e C C L C

Boundary Conditions
Since two end sided of beam are fixed and no deformation will occur so q 1 = θ 1 = q 3 = θ 3 = 0. At that point there is no bending moment, shear force so the loading condition F 1 = F 3 = M 1 = M 3 = 0. Only the load applied in the middle portion, F 2 = 500 N and the bending moment M 2 and the cracked position of the beam will remain active. There are three faces present bounding the calculation domain which are thin Elastic Layer (Boundary 6) is shown in Figure 4.
For Symmetry thin Elastic Layer, 0 n u ⋅ = (16)  Table 1 and Table 2. As our domain with defect is complicated thus computer processor capacity becomes a significant issue for the computational study. The finer mesh is used along the whole computational model for numerical simulation. We had used 16 GB DDR3 RAM, Intel core i5 processor based computer for our simulation.

Numerical Results and Discussions
In this study, we have investigated the frequency of the concrete beam containing double crack using finite element method. For our simulation, we construct a solid concrete beam and have used different parameter values according to Table 3 and Table 4. In Figure 6 it is shown that a load of 500 N is applied vertically on the top phase of the concrete beam. Figure 7 shows that the deflection and the phase of the computational domain after applying load. After applying the load, it is found that the load affects the body and frequency variability is observed essentially in the affected area. Figure 8 shown that different magnitudes of stress applied in the domain and corresponding deflection due to load. It is observed that, after applying load on the beam, a vibration on the structural body is created. The natural frequencies become maximum at the crack position due to the presence of vibration on that point. The maximum portion of load is absorbed in the middle part on the body, so there will be maximum absorption of load on that point. It is also noted that maximum load creates much vibration, which is the major causes to creating crack on any structural body. Figure 9 shows that the slices of the Load distribution at different crack positions of the domain. We observed the load is maximum at the bottom of the domain. It is also found that the load was applied at the crown edge of the domain and distributed to bottom edge especially at the cracked position with high   vibration and differs the frequency. But there is a difference is formed for the (f) that, after the load the beam distorted and deflected. So due to double deflection in this situation, the applied loads distributed into two end sides. Finally we observed that, at the cracked position vibration and frequencies are increased as the load increased and the presence of crack affects the natural frequency of the structure.    position of crack. It is also seen that as long as the crack changes its position (goes to end points) its frequencies are also increased gradually. Finally from the graph it is clear to us that natural frequencies of the beam are directly affected by the location of the cracks.

Validation of the Study
In Table 5 we collect data and tried to compare with the results. According to Figure 11 it is observed that when the position of the crack moves from the fixed end towards the free end of the beam, the effect of the crack also decreases gradually. In order to check the accuracy of the present analysis, we tried to compare our experiment with the experiment done by Priyadarshini A [14] and Kisa et al. [15] the case is considered to validate the program. The presence and position of the crack are typically detected from the comparison of the basic modes of cracked beam. It should be mentioned here that they have used different models and different computational method for their experiments. By comparing with the study of Priyadarshini A [14] and Kisa et al. [15] our results shows an agreement with their results though initially different approach for first two cracks. The analysis was performed on the first six mode shapes. Results of this research on the effectiveness of the damage detection technique applied to higher vibration modes lead to the conclusions that, in practical the presence of crack affects the natural frequency of the structure. The changes of the natural frequency are directly influenced by the different location of the crack.

Conclusion
A numerical simulation is done for a cracked concrete beam with different locations to detect the cracks. The flexibility matrix method is used to calculate the stiffness of the cracked beam by using FEM based vibration model using COMSOL Multiphysics [16]. In this study it is observed that, the frequencies of the cracked beam have been changed with the change of the location of the crack for all the modes of vibration. The effect of crack is higher at near the fixed end than at the free end. From our simulation it is found that the frequencies increase proportional to load distribution and at cracked position the body gets much frequency. Finally we come to the conclusion that the presence of crack and different amounts of loads has a great influence on changing the mode shapes and frequencies of beam. The results confirm that detection of crack size, crack location in a concrete beam depends on natural frequencies and mode shapes. For the established curvature function, if the load is higher, the value of the curvature modulus is also higher, which indicates that higher modes are more sensitive to the presence of the defect. In future our model can be used for crack inclination and different depths.