This work addresses the development of a new four-noded rectangular Mindlin plate bending element (MP4C) with a crack which consists of three degrees of freedom (DOF) at each corner node. The crack in the element is assumed to be not closed and nonpropagating. The crack affects the elastic strain energy and the flexibility matrix of the element, whereas the mass matrix remains unchanged. The complete element stiffness matrix is constructed as the inverse of the combined flexibility matrix of both noncracked and cracked elements. To evaluate the behavior of the proposed cracked Mindlin plate element, numerical examples are provided. They are based on developing user subroutines in ABAQUS. The finite element analysis results using the developed plate element are in excellent agreement with those reported in previous work. The cracked plate element developed in this paper provides a simple and robust approach to model the real service conditions in plate-like structures.

1. Introduction

Civil structures such as highway bridges are subject to a complex environment during their lifetime. In the past two decades, researchers at the University of Connecticut have been exploring different approaches to conduct long-term monitoring for application to bridges. These efforts include vibration monitoring on different types of concrete bridges, nondestructive load testing on steel girder bridges, environmental evaluation of a multispan continuous bridges, and the development of a variety of damage detection methods based on both experimental and analytical bridge models [1, 2]. This paper reports on the study of novel techniques that can be used for finite element analysis of cracked plate-like structures to consider their real service conditions.

A crack or local defect is known to affect the dynamic response of a structural member. The first attempts to quantify defects were made by Thomson [3] and Kirmsher [4]. In their investigations, the effect of a notch on the structure flexibility was simulated by a local bending moment or reduced section, with magnitudes measured by experiments. Since 1950s, as researchers started to realize that the local flexibility will reduce the stiffness of a structural member and subsequently reduce its natural frequency, the analysis of the local flexibility of a cracked region of structural element was largely considered for the structural damage identification.

In the following a few decades, the concept of local flexibility in crack regions has been extensively applied in developing techniques of finite elements with cracks. In those elements, the crack was represented by a localized compliance, and the stiffness matrix of the cracked element was expressed as a function of the crack depth and location. While the local flexibility approach has been developed for elements like truss, beam, and rod, few researches have dealt with plate elements with cracks. Guan-Liang et al. [5] discussed an approximate formulation of a cracked plate bending element by finding the strain energy under the occurrence of crack. They derived the element stiffness matrix from the integration of stress intensity factors (SIFs). Almost at the same time, Krawczuk [6] independently developed a similar cracked plate finite element with in-plane forces, and later extended his method to develop shell element with a through crack and plate element with an elastoplastic through crack [7, 8]. These existing plate elements, however, neglect the effects of transverse shear, and so they are insufficient for practical engineering applications when it is necessary to model plate bending behaviors. In engineering applications, many plate members must be modeled as moderately thick plate elements due to the relatively large thickness in comparison with other two dimensions.

Mindlin [9] presented a complete theory for the shear deformable plates in static and dynamics, referred to as the first-order shear deformation theory. This theory, which accounts for transverse shear deformations, is better suited for the solution of plate bending problems than the Kirchhoff plate theory. A wealth of further analytic, numerical, and experimental investigations done by other investigators with the Mindlin theory make it feasible for use as a rigorous cracked plate bending element. Thus this paper develops a new rectangular cracked Mindlin plate bending finite element using the local flexibility approach and explicitly considering three independent fracture modes [10]. This cracked element possesses four corner nodes with three engineering degrees of freedom (DOF) at each node. In the derivation of the stiffness matrix, selective numerical integration is used to overcome the shear-locking problem. This paper also explores a practical engineering analysis tool that can readily model the dynamic response of structures in the presence of discrete cracks, with reasonable engineering accuracy. This paper provides a framework for incorporating the developed elements into the commercial finite element software ABAQUS [11] through a user element subroutine. The behaviors of the developed noncracked and cracked elements are examined and compared with results presented by other researchers and elements available in ABAQUS. The examples presented in this paper address the performance of the proposed elements for use in problems with either thin or thick plates.

2. Fundamental Relations in the Mindlin Plate Theory

Kirchhoff elements were introduced earlier by the influence of the classical plate theory. Mindlin elements appeared later with the assumption that a straight line originally normal to the plate midsurface remains straight but not normal after deformation. A shear factor was introduced to take account of the fact that the shear strain through the plate thickness is not uniformly distributed. Figure 1 shows the typical nodal DOFs in the Mindlin plate element. The inclusion of transverse shear effects in the classical plate theory indicates that the cross-sectional rotations 𝜓𝑥 and 𝜓𝑦 are no longer expressed solely in terms of the deflection 𝑤. Thus, there exist three basic reference quantities, namely, 𝑤, 𝜓𝑥 and 𝜓𝑦, at the midsurface.

2.1. Displacement Fields

Considering the DOF in the Mindlin plate element as shown in Figure 1, the displacement fields in the Cartesian coordinate system are given as [12]𝑢=𝑢(𝑥,𝑦,𝑧)=−𝑧𝜓𝑥(𝑥,𝑦),𝑣=𝑣(𝑥,𝑦,𝑧)=−𝑧𝜓𝑦(𝑥,𝑦),𝑤=𝑤(𝑥,𝑦,𝑧)=𝑤(𝑥,𝑦),(1) where 𝑤, 𝜓𝑥, and 𝜓𝑦 are the lateral (𝑧-direction) deflection and the two normal rotations of the midsurface, respectively.

2.2. Strain-Curvature Relations

For small displacements and rotations, strains are given by [12]:𝜀𝑥=−𝑧𝜓𝑥,𝑥,𝜀𝑦=−𝑧𝜓𝑦,𝑦,𝛾𝑥𝑦𝜓=−𝑧𝑥,𝑦+𝜓𝑦,𝑥,𝛾𝑦𝑧=𝑤,𝑦−𝜓𝑦,𝛾𝑧𝑥=𝑤,𝑥−𝜓𝑥,(2) where a comma denotes differentiation with respect to the following subscript. Equation (2) is the basis of the Mindlin plate theory, which allows nonzero transverse shear strains. The bending curvatures are given by𝜅𝑥=𝜓𝑥,𝑥,𝜅𝑦=𝜓𝑦,𝑦,𝜅𝑥𝑦=𝜓𝑥,𝑦+𝜓𝑦,𝑥.(3) Both curvatures and transverse shear strains appeared in the vector {𝜅} as follows:𝜅{𝜅}=𝑥,𝜅𝑦,𝜅𝑥𝑦,,−𝛾𝑧𝑥,−𝛾𝑦𝑧𝑇.(4)

2.3. Constitutive Relations

For homogeneous, isotropic, and linearly elastic material, the moment-curvature relations for the Mindlin plate are𝑀𝑐=[𝐷]{𝜅},(5) in which {𝑀𝑐} is the moment vector including {𝑀𝑥,𝑀𝑦,𝑀𝑥𝑦,𝑄𝑥,𝑄𝑦}𝑇 and [𝐷] is the rigidity matrix, analogous to flexural stiffness 𝐸𝐼 of a beam. The forces and associated stresses on the plate are described in later sections.

3. Mindlin Plate Element (MP4)

Figure 2 shows a four-node Mindlin plate element with 3 DOF at each node. The lateral displacement and rotation fields of the MP4 element with 𝑛 nodes can be given byâŽ§âŽªâŽ¨âŽªâŽ©ğ‘¤ğœ“ğ‘¥ğœ“ğ‘¦âŽ«âŽªâŽ¬âŽªâŽ­=𝑛𝑖=1âŽ›âŽœâŽœâŽœâŽğ‘ğ‘–000𝑁𝑖000ğ‘ğ‘–âŽžâŽŸâŽŸâŽŸâŽ âŽ§âŽªâŽ¨âŽªâŽ©ğ‘¤ğ‘–ğœ“ğ‘¥ğ‘–ğœ“ğ‘¦ğ‘–âŽ«âŽªâŽ¬âŽªâŽ­[𝑁]or{𝑢}={𝑑},(6) where the same shape functions are used for each displacement component.

With {𝑢}=[𝑤,𝜓𝑥,𝜓𝑦]𝑇, we obtain from (2) and  (4)âŽ§âŽªâŽªâŽªâŽ¨âŽªâŽªâŽªâŽ©ğœ“{𝑘}=𝑥,𝑥𝜓𝑦,𝑦𝜓𝑥,𝑦+𝜓𝑦,𝑥𝜓𝑥−𝑤,𝑥𝜓𝑦−𝑤,ğ‘¦âŽ«âŽªâŽªâŽªâŽ¬âŽªâŽªâŽªâŽ­=[𝜕][𝜕]=⎡⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣0𝜕{𝑢},where0𝜕𝜕𝑥000𝜕𝜕𝑦𝜕𝜕𝑦𝜕𝑥−𝜕𝜕𝑥10âˆ’ğœ•âŽ¤âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¥âŽ¦.𝜕𝑦01(7) From (6) and (7), we obtain for an element with 4 nodes[𝐵]{𝑘}=5×3𝑛[𝐵]=[𝜕𝑁]=⎡⎢⎢⎢⎢⎢⎢⎢⎣{𝑑},where][0𝑁1,𝑥0⋯000𝑁1,𝑦⋯𝑁𝑛,𝑦0𝑁1,𝑦𝑁1,𝑥⋯𝑁𝑛,𝑥−𝑁1,𝑥𝑁10⋯0−𝑁1,𝑦0𝑁1⎤⎥⎥⎥⎥⎥⎥⎥⎦.⋯0(8) The element stiffness matrix [𝐾] is[𝐾]=[𝐵]𝑇[𝐷𝐵]][𝑑𝐴.(9)

4. Cracked Mindlin Plate Element (MP4C)

Figure 3(a) illustrates a quadrilateral plate with a through crack of 2ğ‘Ž and subjected to twelve element forces 𝑆𝑛 (𝑛=1–12). These twelve forces are related by three equations of equilibrium, and consequently only nine of these forces are linearly independent. Figure 3(b) shows the independent force system (Przemieniecki [13]). The matrix [𝑇], given in Appendix A.1, was used to transform the nodal forces 𝑆1–𝑆12 to the independent ones 𝐹1–𝐹9. The stiffness matrix [𝐾𝑐] of the MP4C element can be now written as𝐾𝑐=[𝑇]𝑇𝐶𝑓−1[𝑇],(10) where [𝐶𝑓]−1 denotes the inverse matrix of [𝐶𝑓]=[𝐶𝑓]0+[𝐶𝑓]1 in which [𝐶𝑓]0 is the flexibility matrix of the noncracked element and [𝐶𝑓]1 is the flexibility matrix due to the presence of the crack. The details of the derivation of these two flexibility matrices are given in the next two sections.

4.1. Flexibility Matrix of the Noncracked Element

The elemental flexibility matrix is given by Przemieniecki [13]𝑓=𝑣[𝑒]𝑇[Φ𝑒]][𝑑𝑉,(11) in which matrix [𝑒] represents the approximate stress distribution due to unit force and [Φ] is the inverse matrix of flexural rigidity. Figure 4(a) depicts internal stresses on the plate cross-section. These stresses are distributed over the thickness of the plate and they produce bending moments, twisting moments, and vertical shear forces per unit of length in the 𝑥𝑦 plane, as shown in Figure 4(b). The in-plane stresses ğœŽğ‘¥, ğœŽğ‘¦, and 𝜏𝑥𝑦, vary linearly with 𝑧, whereas transverse shear stresses 𝜏𝑦𝑧 and 𝜏𝑧𝑥 usually vary parabolically through the thickness. The relationships between these stresses with the moments and forces are given byğœŽğ‘¥ğ‘¥=12𝑧𝑡3𝑀𝑥,ğœŽğ‘¦ğ‘¦=12𝑧𝑡3𝑀𝑦,𝜏𝑥𝑦=12𝑧𝑡3𝑀𝑥𝑦,𝜏𝑦𝑧=32𝑡1−2𝑧𝑡2𝑄𝑦,𝜏𝑧𝑥=32𝑡1−2𝑧𝑡2𝑄𝑥.(12)

Using (12) and assuming that the normal stresses vary linearly and the shearing stresses are constant along the plate edges, we can obtain the matrix [𝑒] as[𝑒]𝑇=12𝑧𝑡3⎡⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣02𝐿02(2−3𝜉)000𝐿2(2−3𝜉)(1−2𝜂)000𝐻2(−1+3𝜂)0000𝐻02(−1+3𝜂)(1−2𝜉)0000𝐿02(−1+3𝜉)000𝐿2(−1+3𝜉)(−1+2𝜂)000𝐿2(2−3𝜂)0000𝐿𝑡(2−3𝜂)(−1+2𝜉)000000128𝑧1−2𝑧𝑡2𝑡28𝑧1−2𝑧𝑡2⎤⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦.(13)in which 𝜉,𝜂 are the local coordinates for the element. Substituting (13) into (11), we can derive the flexibility matrix [𝐶𝑓]0 which is given in Appendix A.2.

4.2. Flexibility Matrix due to Crack

We have the following relationship due to crack in the element: 𝑐1𝑖,𝑗=𝜕2𝑈1𝜕𝐹𝑖𝜕𝐹𝑗,(14) where 𝐹𝑖 denotes independent nodal forces and 𝑈1 is the strain energy induced due to crack given in Ugural [14]𝑈1=1𝐸𝐴9𝑖=1𝐾2𝐼,𝑖+9𝑖=1𝐾2𝐼𝐼,𝑖+(1+𝑣)9𝑖=1𝐾2𝐼𝐼𝐼,𝑖𝑑𝐴1,(15) where 𝐾𝑗𝑖 are stress intensity factors for two modes of crack evolution, and 𝐴1 is the area of the crack. For a plate with uniform thickness 𝑡 containing a through crack of 2ğ‘Ž, and subjected to the 𝐹1–𝐹9 independent forces, the SIFs are given by Sih and Liebowitz [15]:𝐾∗𝐼=1âˆšî€œğœ‹ğ‘Žğ‘Žâˆ’ğ‘ŽğœŽğ‘¥ğ‘¥î‚™(0,𝑦)ğ‘Ž+ğ‘¦ğ¾ğ‘Žâˆ’ğ‘¦ğ‘‘ğ‘¦,∗𝐼𝐼=−1âˆšî€œğœ‹ğ‘Žğ‘Žâˆ’ğ‘Žğœğ‘¥ğ‘¦î‚™(0,𝑦)ğ‘Ž+ğ‘¦ğ¾ğ‘Žâˆ’ğ‘¦ğ‘‘ğ‘¦,∗𝐼𝐼𝐼=1âˆšî€œğœ‹ğ‘Žğ‘Žâˆ’ğ‘Žğœğ‘¦ğ‘§î‚™(0,𝑦)ğ‘Ž+ğ‘¦ğ‘Žâˆ’ğ‘¦ğ‘‘ğ‘¦,(16) where ğœŽğ‘¥ğ‘¥, 𝜏𝑥𝑦, and 𝜏𝑦𝑧 are normal, in-plane shear, and transverse shear stress acting on the crack face given byğœŽğ‘¥ğ‘¥(0,𝑦)=24𝑧𝐻𝑡312+3𝑦𝐻𝐹3+12−3𝑦𝐻𝐹7,𝜏𝑥𝑦(0,𝑦)=12𝑧𝑡3𝐹9,𝜏𝑦𝑧3(0,𝑦)=2𝑡1−2𝑧𝑡2𝐹9.(17) Substituting (17) into (16), the SIFs for an elastic nonpropagating crack are given by𝐾∗𝐼=24𝑧𝐻𝑡3𝐹312+3ğ‘Žî‚2𝐻+𝐹712−3ğ‘Žâˆš2ğ»î‚î‚„ğ¾ğœ‹ğ‘Ž,∗𝐼𝐼=−12𝑧𝑡3𝐹9âˆšğ¾ğœ‹ğ‘Ž,∗𝐼𝐼𝐼=32𝑡1−2𝑧𝑡2𝐹9âˆšğœ‹ğ‘Ž.(18)

The above SIFs are suitable for an infinite plate only. To take into account finite dimensions of the element, a geometric correction factor can be used [16]𝑌(𝑔)=1.0+0.018762ğ‘Žğ»î‚î‚€+0.18252ğ‘Žğ»î‚2+2.0242ğ‘Žğ»î‚3−2.43162ğ‘Žğ»î‚4.(19) Combining the geometric correction factor leads to𝐾𝐼=24𝑧𝐻𝑡3𝐹312+3ğ‘Žî‚2𝐻+𝐹712−3ğ‘Žâˆš2ğ»î‚î‚„ğ¾ğœ‹ğ‘Žğ‘Œ(𝑔),𝐼𝐼=−12𝑧𝑡3𝐹9âˆšğ¾ğœ‹ğ‘Žğ‘Œ(𝑔),𝐼𝐼𝐼=32𝑡1−2𝑧𝑡2𝐹9âˆšğœ‹ğ‘Žğ‘Œ(𝑔).(20)

The expression for the mode III SIF (𝐾𝐼𝐼𝐼) in (20) is first derived in this paper. We explicitly use the independent three fracture modes to evaluate the flexibility matrix of the cracked Mindlin plate element. For this approach, there is no need to incorporate fracture mode correction functions [5–8], used to transfer from the Kirchhoff theory to the Mindlin theory.

Substituting (15) and (20) into (14) leads to the following nonzero terms in the flexibility matrix of the MP4C element with crack:𝑐133=4𝑔𝑘−𝑔𝑘𝑔12+3𝑔42𝑌2,𝑐(𝑔)𝑑𝑔137=2𝑔𝑘−𝑔𝑘𝑔12+3𝑔412−3𝑔4𝑌2,𝑐(𝑔)𝑑𝑔177=4𝑔𝑘−𝑔𝑘𝑔12−3𝑔42𝑌2,𝑐(𝑔)𝑑𝑔199=𝑡1+(1+𝑣)2𝐻102𝑔𝑘−𝑔𝑘𝑔𝑌2.(𝑔)𝑑𝑔(21) The flexibility matrix [𝐶𝑓]1 derived in this procedure is given in Appendix A.3.

5. Numerical Examples

5.1. Example 1: Free Vibration of a Simply Supported Plate

In this example, the natural frequencies of free vibration of a simply supported square plate are calculated. The geometric data employed are the length of the edge of the plate 𝐿=0.1 m, the elastic modulus 𝐸=2.04×1011 N/m2, the Poisson ratio 𝜐=0.3, and the density 𝜌0=7860 kg/m3. Three values of plate thickness, 𝑡=0.001 m, 𝑡=0.005 m, and 𝑡=0.01 m, are considered representing the thin (𝑡/𝐿=1%), moderate (𝑡/𝐿=5%), and thick plate (𝑡/𝐿=10%), respectively.

The natural frequencies in three meshes (5×5, 7×7, and 10×10) have been computed and compared with the analytical solution and ABAQUS shell element (S4 and S4R) solutions. The analytical solution of the natural frequency 𝜔𝑚𝑛 for a simply supported rectangular plate is given [14]. The mode shapes of the first three modes are shown in Figure 5.

We use shell elements in ABAQUS’s library (2004) for comparison with MP4 element for the plate bending results. In ABAQUS, triangular shell element with 3 nodes (STRI3) is recommended to be used for Kirchhoff element [19]. Quadrilateral elements S4 and S4R, with a hypothesis of linear interpolation in the transverse direction, are commonly used as either Kirchhoff or Mindlin elements depending on the plate thickness. The S4 element is a fully integrated, isoparametric shell element with 6 DOFs at all four corner nodes. The S4R element has the same behavior as S4 except that it uses the reduced integration method, which means it contains only one center Gauss point. We used three different meshes (5×5, 7×7 and 10×10) for the plate when MP4, S4, and S4R elements were used. Free mesh was used for STRI3 elements. The three corresponding mesh patterns with STRI3 contain 20 elements and 17 nodes, 52 elements and 37 nodes, and 108 elements and 69, nodes respectively.

Tables 1, 2, and 3 compare the natural frequencies of the present MP4 plate element with other available reference solutions in a simply supported plate with thicknesses of 𝑡=0.001, 0.005, and 0.01 m, respectively. We found that the MP4 element, when used in vibration analysis, is expected to perform better than the STRI3, S4, and S4R elements, and that the best results for the natural frequencies can be obtained with 10×10 mesh of the plate. Table 4 compares the error of these four elements relative to thin plate theory and to thick plate theory for 10×10 mesh. It can be seen that the STRI3 element only works for a thin plate, and the S4R element works better than S4 for the thin plate and it is the reverse for the thick plate. We may see that the reduced integration is more accurate for the thin plate element whereas the full integration for thick plate element. The developed MP4 element performs excellent in vibration analysis for both thin and thick plate, especially for moderately thick plate. Note that when the MP4 elements are used for moderately thick plate (𝑡/𝐿=5%), the first and third mode of the plate are closer to the thick plate theory solution (𝛽=0.0029 and 0.0006), whereas the second mode is closer to the thin plate theory solution (𝛼=0.0032).

5.2. Example 2: Free Vibration of a Simply Supported Plate with a Through Crack

Figure 6 shows the plate with thickness of 0.001 m containing a crack in the middle of the plate. The crack is oriented parallel to the plate edges. A 5×5 mesh is used to model the plate, which includes 24 MP4 elements and one MP4C element shaded in the figure. Changes in the first four natural frequencies of the plate due to crack are investigated. It should be noted that in [8] the influence of the crack location in the plate on the changes of natural frequencies was analyzed. In this example, the most critical crack for a simply supported plate is located at the position of local maxima of bending moments. The ratio of natural frequency of cracked plate to one of the plates without crack, Δ𝜔𝑖=𝜔𝑐𝑖/𝜔𝑖, is used as the testing index, and 𝜔𝑖 is the 𝑖th natural frequency of noncracked plate, and 𝜔𝑐𝑖 is the 𝑖th natural frequency of the cracked plate. The relative crack length (2ğ‘Ž/𝐻) considered ranges from 5% to 20% of the plate width.

Table 5 compares the ratios (Δ𝜔𝑖) of natural frequencies of the developed MP4C element with other reference results for four modes. It is worth noting that the solutions in the references were obtained using the regular plate bending elements based on the Kirchhoff theory. The present results are in good agreement with the reference results. We observe that for the simply supported boundary condition, crack has the most significant effects on the basic bending mode of plate. Note that the values of Δ𝜔𝑖 decrease as the crack length increases in each single mode, and for higher modes, the effect of crack is not as obvious as that on the lower modes. The present results are close to those obtained using the Kirchhoff plate models since the transverse shear term only has slight influence in the flexibility matrix of the Mindlin plate. However, we observe that the present results are generally larger than those of other references, which indicates the crack has less effect in the Mindlin plate model than in the Kirchhoff plate model. The reason is that the Mindlin plate element is more flexible than the Kirchhoff plate element because the Kirchhoff plate element neglects the transverse deformation on its cross-section which actually adds stiffness. Further investigation of this phenomenon is given in next example.

5.3. Example 3: Free Vibration of a Cantilever Plate with a Through Crack

This example focuses on the effect of crack on the behavior of a cantilever square plate. The plate is made of aluminum with Young’s modulus 𝐸=6.9×1010 N/m2, Poisson ratio 𝑣=0.33, and density 𝜌0=2800 kg/m3. The length of an edge of the square plate was 0.24 m and its thickness was 0.00275 m. A through crack of length 2ğ‘Ž=0.0035 m was intentionally introduced at 𝑥=0.09 m, 𝑦=0.09 m, parallel to the clamped boundary (see Figure 7). In this case, a 4×4 mesh including 15 MP4 elements and one MP4C element is used.

Table 6 compares the ratios (Δ𝜔𝑖) of the first four natural frequencies of the present results with those of available results in Guan-liang et al. [5] and Krawczuk et al. [8]. The present results show very good agreement with experimental results in Guan-liang et al. [5] and other available numerical results. Similar to the previous example, the larger values of Δ𝜔𝑖 are observed in the present model. Table 7 compares the ratios (Δ𝜔𝑖) of the first four natural frequencies of the present results for four different thickness-to-length ratios: 𝑡/𝐿=0.1%, 1%, 5%, and 10%. We found that, with the same crack length, the thinnest plate (𝑡/𝐿=0.1%) is subjected to the least influence on natural frequency. Significant changes in natural frequencies are observed in the case of the moderately thick plate (𝑡/𝐿=5%) and the thick plate (𝑡/𝐿=10%) at the first and second modes, which shows that the crack effects on plates are modal dependent. The results in this example demonstrate that the plate thickness is a key factor when evaluating the effect of crack on plate bending behavior.

6. Conclusions

In this paper, we developed a new four-noded rectangular Mindlin plate finite element with a through crack. This element was used for vibration analyses based on plate bending. The performances of the MP4 and MP4C elements are compared with other available elements in the literature and also with ABAQUS shell elements (STRI3, S4, and S4R). Moreover, the developed elements are integrated into ABAQUS [1] as user elements.

From this study, we may conclude the following.(i)The formulation of the MP4 element, which is free of shear locking, is simple and straightforward. As evident through comparison with the thin and thick plate theories, the proposed MP4 element has shown better accuracy for natural frequencies than conventional shell elements available in ABAQUS. In particular, the proposed element yields better accuracy for moderately thick plates, which are more common types of plates in engineering practice.(ii)The local flexibility method is used to develop the cracked MP4C element based on the noncracked MP4 element. Mixed-mode SIFs have been independently introduced to fully represent the effect of crack on the flexibility matrix of the Mindlin plate element. We observe that dynamical behaviors of plate structures are influenced by the crack length, crack location, boundary conditions, and plate thickness.(iii)The present energy-based approach efficiently addresses the global vibration behavior of plate structures considering strain energy perturbation induced by crack and thus makes computational efficiency possible without the need for a highly dense mesh around crack tips. Moreover, the integration of MP4 and MP4C elements into FEM packages, such as ABAQUS, provides an efficient tool to model the static, dynamic, and stability behaviors of plate-like structures with multiple cracks.



The transformation matrix [𝑇] is[𝑇]=⎡⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣02𝐻200000−𝐿2−1110000000000000−1100−𝐻02𝐿200001−11000000000−1−100000000−𝐿02𝐻200−10000−1−1000001−10000000000−𝐻02𝐿1⎤⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦00001−1000000000110.(A.1)


The flexibility matrix of the analyzed element without the crack [𝐶𝑓]0 is derived in this paper𝐶𝑓0=12𝐸𝑡3⎡⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣4𝐻𝐿0𝑠𝑦𝑚4𝐻3𝐿−𝑣𝑣4𝐿𝐻−𝑣𝑣04𝐿−3𝐻2𝐻𝐿0−𝑣𝑣4𝐻𝐿02𝐻3𝐿−𝑣𝑣04𝐻3𝐿−𝑣𝑣−2ğ¿â„Ž0−𝑣𝑣4𝐿𝐻−𝑣𝑣02𝐿3𝐻−𝑣𝑣04𝐿𝑡3𝐻000000002(1+𝑣)𝐻𝐿1+25⎤⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦.(A.2)


The flexibility matrix of the analyzed element with crack [𝐶𝑓]1 is 𝐶𝑓1=6𝜋𝐸𝑡3⎡⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣0𝑠𝑦𝑚0000𝑐3300000000000000000𝑐37000𝑐770000000000000000𝑐99⎤⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦.(A.3)


ğ‘Ž:Half crack length
𝐴:Area of the plate element
𝐴1:Area of the crack
[𝐵]:Derivative of the shape function matrix
[𝐵𝑏]:Bending terms in [𝐵]
[𝐵𝑠]:Transverse shear terms in [𝐵]
𝑐1𝑖,𝑗:Terms in the flexibility matrix due to the crack
[𝐶]:Damping matrix in governing equation for structural dynamics
[𝐶𝑓]:Flexibility matrix of the cracked Mindlin plate element
[𝐶𝑓]0:Flexibility matrix of the Mindlin plate element without crack
[𝐶𝑓]1:Flexibility matrix due to the crack
{𝑑}:Nodal displacement vector in governing equation for structural dynamics
[𝐷]:Flexural rigidity matrix
ē:Stress distribution due to unit force in plate
𝐸:Elastic modulus of the plate
𝐹1–𝐹9:Independent nodal forces on the Mindlin plate element
𝐻:Width of the plate element
𝐾𝑖:Stress intensity factors, 𝑖=𝐼,𝐼𝐼,𝐼𝐼𝐼
{𝜅}:Curvature and strain column vector
[𝐾]:Stiffness matrix of the Mindlin plate element
[𝐾𝑐]:Stiffness matrix of the cracked Mindlin plate element
𝐿:Length of the plate element
{𝑀𝑐}:Moment and shear force column vector
𝑁𝑖(𝜉,𝜂):Shape function at node 𝑖 for the four-node Mindlin plate element
{𝑅ext}:External force vector in governing equation for structural dynamics
𝑆1–𝑆12:Nodal forces on the Mindlin plate element
𝑡:Thickness of the plate element
[𝑇]:Transformation matrix
𝑢:Displacement of the plate midsurface in positive 𝑥-direction
𝑈:Strain energy in the element
𝑈1:Additional strain energy due to the crack
𝜈:Displacement of the plate midsurface in positive 𝑦-direction
𝑤:Displacement of the plate midsurface in positive 𝑧-direction
𝑌(𝑔):Correction function for the infinite dimension of element to finite dimension
𝜀𝑥:Normal strain of the plate in 𝑥-direction
𝜀𝑦:Normal strain of the plate in 𝑦-direction
𝛾𝑥𝑦:Shear strain in 𝑥𝑦 plane of the plate
𝛾𝑦𝑧:Transverse shear strain in 𝑦𝑧 plane of the plate
𝛾𝑧𝑥:Transverse shear strain in 𝑧𝑥 plane of the plate
𝜌0:Mass density of the plate
𝜐:Poisson ratio of the plate
Δ𝜔𝑖:Ratio of frequencies of cracked plate to one of the plate without crack, Δ𝜔𝑖=𝜔𝑐𝑖/𝜔𝑖
𝜔𝑖:ith natural frequency of noncracked plate
𝜔𝑐𝑖:ith natural frequency of the cracked plate
ğœŽğ‘¥(0,𝑦):Normal stress acting on the crack face
𝜏𝑥𝑦(0,𝑦):Shear stress acting on the crack face
𝜓𝑥:Rotation of the plate midsurface pointing to negative 𝑦-direction
𝜓𝑦:Rotation of the plate midsurface pointing to positive 𝑥-direction
[Φ]:Inverse matrix of flexural rigidity.


This work has been carried out at the University of Connecticut as part of a research project on bridge monitoring directed by the middle author. The support of the Connecticut Department of Transportation and the U.S. Department of Transportation, Federal Highway Administration is gratefully acknowledged. The work was done with further support from the Connecticut Transportation Institute at the University of Connecticut. This work is an outcome of a course taught by the last author. The start-up support provided to him by the University of Connecticut is gratefully acknowledged.