Journal of Applied Mathematics

Volume 2012 (2012), Article ID 734521, 17 pages

doi:10.1155/2012/734521

Research Article

Analysis of the Elastic Large Deflection Behavior for Metal Plates under Nonuniformly Distributed Lateral Pressure with In-Plane Loads

Jeom Kee Paik, Ju Hye Park, and Bong Ju Kim

The Ship and Offshore Research Institute, The Lloyd’s Register Educational Trust Research Centre of Excellence, Pusan National University, Busan 609-735, Republic of Korea

Received 21 February 2012; Accepted 22 August 2012

Academic Editor: Carl M. Larsen

Copyright © 2012 Jeom Kee Paik et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The Galerkin method is applied to analyze the elastic large deflection behavior of metal plates subject to a combination of in-plane loads such as biaxial loads, edge shear and biaxial inplane bending moments, and uniformly or nonuniformly distributed lateral pressure loads. The motive of the present study was initiated by the fact that metal plates of ships and ship-shaped offshore structures at sea are often subjected to non-uniformly distributed lateral pressure loads arising from cargo or water pressure, together with inplane axial loads or inplane bending moments, but the current practice of the maritime industry usually applies some simplified design methods assuming that the non-uniform pressure distribution in the plates can be replaced by an equivalence of uniform pressure distribution. Applied examples are presented, demonstrating that the current plate design methods of the maritime industry may be inappropriate when the non-uniformity of lateral pressure loads becomes more significant.

1. Introduction

Ships and ship-shaped offshore structures are composed of metal plate elements, and the accurate computation of nonlinear behavior of the plate elements in deck, bottom and side shells up to the ultimate limit state is a basic requirement for the structural safety assessment. The plate elements in ships and ship-shaped offshore structures are generally subjected to combined inplane and lateral pressure loads. Inplane loads include biaxial compression/tension, biaxial inplane bending and edge shear, as shown in Figure 1, which are mainly induced by overall hull girder bending and/or torsion of the vessel. Lateral pressure loads are due to water pressure and/or cargo. In rough weather, roll and/or pitch motions of vessels are typical as shown in Figure 2, and, subsequently, the distribution of lateral pressure loads often becomes non-uniform as shown in Figure 3. For simplicity, however, the maritime industry usually applies a simplified design method in which the non-uniform distribution of lateral pressure loads is replaced by an equivalence of uniformly distributed lateral pressure loads with an average magnitude of applied pressure loads.

Figure 1: In-plane load components applied in a plate element.

Figure 2: Roll motion of a vessel causing nonuniformity of lateral pressure load distribution in the hull cross-section.

Figure 3: Non-uniformly distributed lateral pressure loads in a plate element.

A large number of studies have been available in the literature, for example [1–6]. In the present paper, a mathematical algorithm is derived to analyze the elastic large deflection behavior, including buckling and postbuckling response, of plate elements under combined inplane and lateral pressure loads noted previously, with the emphasis on the non-uniformly distributed lateral pressure loads. The Galerkin method is applied to solve the nonlinear governing differential equations of elastic large deflection plate theory for plate elements. In the literature, the studies of metal plate buckling are also available, for example [7–13].

2. Elastic Large Deflection Analysis

The elastic large deflection behavior of a plate element with initial imperfections can be analyzed by solving two differential equations, one representing the equilibrium condition and the other representing the compatibility condition [5]:

Once the Airy stress function is defined, the membrane stress components at a certain position inside the plate may be expressed as follows:

To apply the Galerkin method for solving the nonlinear governing differential equations, the added deflection and initial deflection are assumed as follows:

where and are functions which satisfy the boundary conditions for the plate. and are unknown and known deflection coefficients, respectively.

Equations (2.4) and (2.5) are substituted into (2.2) to obtain the stress function . In this case, the particular solution of the stress function can be expressed as follows:

where the coefficients will be second-order functions with regard to the unknown deflection coefficients .

With the homogeneous solution of the stress function which satisfies the applied loading condition, the complete stress function may be given by

To compute the unknown coefficients , the Galerkin method is applied to the equilibrium equation (2.1), resulting in the following equation:

Equations (2.4), (2.5), and (2.7) into (2.8), and performing the integration over the whole volume of the plate, a set of third-order simultaneous equations with regard to the unknown coefficients will be obtained.

3. Application to the Elastic Large Deflection Analysis of Simply Supported Plates

In the present paper, the procedure described in [5] is applied. It is assumed that the plate element is simply supported along the four edges where support members such as longitudinal stiffeners and transverse frames are located. The simply supported edge conditions for the plates should satisfy the following conditions

To satisfy the boundary condition, the added deflection function and the initial deflection can be assumed in Fourier series:

where and are the unknown and the known coefficients, respectively. The condition of combined loads, namely biaxial loads, biaxial inplane bending and edge shear and lateral pressure loads are given as follows:

For simplicity in expressing the various functions, the following abbreviations are involving sine or cosine terms:

To find the stress function which should satisfy (2.2), (3.2) and (3.3) are substituted into (2.2) and the following equation yields:

A particular solution for the Airy stress function can be obtained as follows:

Equation (3.7) is substituted into (2.2), the coefficients , , and , as follows:

Equation (3.8) is substituted into (3.7), and a particular solution is obtained as follows:

By considering the condition of load application, the homogeneous solution for the stress function is given by

The Airy stress function is then expressed by the sum of the particular and homogeneous solution as follows:

To compute the unknown coefficients, the Galerkin method is applied to (2.1) as follows:

Substitution of (3.11) into (2.1) and then (2.1) to (3.12) after a derivation, a set of third order simultaneous equations for the unknown deflection coefficients is obtained, as follows:

where = lateral pressure loads, and the coefficients to and to are given in Appendix.

In (3.13), lateral pressure will take the following form in general.

When uniform pressure loads are applied, that is, with , is taken as constant value.

4. Applied Examples and Discussion

Some applied examples are now presented. A simply supported rectangular plate of under uniaxial compression and lateral pressure loads is considered, in which lateral pressure loads are applied first, followed by the monotonical application of axial compressive loads. The elastic modulus and Poisson’s ratio are = 205.8 GPa and = 0.3. The elastic buckling strength of this plate under uniaxial compression in the direction is determined as follows:

where = buckling coefficient which is taken as = 4.

Also, the lateral pressure loads are taken as and in the present examples. The initial deflection of the plate is selected as follows (see Figure 4):

where .

Figure 4: Configuration of the initial deflection in the plate.

The added deflection of the plate due to applied loads can be expressed as follows.

In the present study, the maximum number of amplitudes in the shorter direction of the plate is assumed as , but that in the longer direction of the plate is varied to investigate the best selection of in terms of the accuracy and efficiency of the computations.

Figures 5(a) to 5(d) present the relation between the axial compressive loads versus total deflection at the center of the plate, that is, at and , where the magnitude of lateral pressure loads and the maximum number of deflection amplitudes are varied. It is found that may be sufficient enough in terms of accuracy and efficiency of the resulting computations. Table 1 represents the deformed shapes of the plate with different magnitudes or shapes of lateral pressure loads when = 0 or 1.2 where is taken. It is observed that the deflection modes of the plate are changed due to the application of axial compressive loads. This is because the plate tends to deflect as per the original buckling pattern associated with axial compressive loads.

Table 1: Deformed shapes of the plate when .

(MPa)

(MPa)

0.1
0.1

0.1
0.15

0.1
0.2

0.1
0.3

0.1
0.4

(a)

(b)

(c)

(d)

Figure 5: (a) The axial compressive load versus total deflection curves at  MPa. (b) The axial compressive load versus total deflection curves at  MPa and  MPa. (c) The axial compressive load versus total deflection curves at  MPa and  MPa. (d) The axial compressive load versus total deflection curves at  MPa and  MPa.

Figure 6 presents a comparison of the elastic large deflection responses of the plate between non-uniform and uniform pressure loads, the latter being taken as the average value of the pressure loads as per the current practice of the maritime industry.

Figure 6: The axial compressive loads versus total deflection at and .

It is found that the current practice of the maritime industry with an average magnitude of applied pressure loads underestimates the lateral deflection of the plate compared to the real condition of the pressure load application, that is, with a non-uniform distribution. The difference between the two cases becomes larger and larger as the non-uniformity of lateral pressure loads becomes more significant.

The underestimation of the plate deflection due to external pressure loads gives rise to the overestimation of plate strength performance which may lead to unsafe design at optimistic side.

5. Concluding Remarks

The aim of the present paper has been to analyze the elastic large deflection behavior of metal plates subject to combined inplane and lateral pressure loads. When lateral pressure loads applied in plate elements are non-uniform, the current practice of the maritime industry applies some simplified design methods in which the non-uniform pressure distribution in the plates is replaced by an equivalence of uniform pressure distribution.

In the present paper, the Galerkin method was used to solve the nonlinear governing differential equations of plate elements under non-uniformly distributed lateral pressure loads in addition to inplane loads. Some applied examples are presented, demonstrating that the current practice of the maritime industry, that is, with an average magnitude of applied pressure loads as an equivalence, results in a great underestimation of lateral deflection calculations when the non-uniformity of lateral pressure loads becomes more significant. The underestimation of the plate deflection may lead to unsafe design of plates at optimistic side.

Thin plates buckle in elastic regime, while stocky plates may buckle in elastic-plastic or plastic regime. The present paper deals with elastic behavior only, and further studies are then recommended to take into account the effect of plasticity which is dominant in thick plates. Also, it is highly desirable for practical design purpose to develop a simpler method or design formulation.

Appendix

The coefficients to and to are given as follows:

Nomenclature

: | Length of the plate
: | Breadth of the plate
: | Plate bending rigidity  
: | Young’s modulus
: | Airy’s stress function
: | Maximum half-wave number of the assumed added deflection function in the and directions
: | In-plane bending moment in the direction ()
: | In-plane bending moment in the direction ()
: | Lateral (out-of-plane) pressure load on the surface area
: | Axial force in the direction
: | Axial force in the direction
: | Thickness of the plate
: | (Added) deflection of the plate due to the action of external loads
: | Initial deflection of the plate
: | Total deflection of the plate
: | Axis direction normal to the plane
: | Aspect ratio of the plate ( )
: | Poisson’s ratio
: | Bending stresses in the and directions
: | Normal stresses in the and directions
: | Mean stresses in the and directions
: | Shear stress.

Acknowledgments

This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (Grant no. K20903002030-11E0100-04610).

References

R. Szilard, Theory and Analysis of Plates: Classical and Numerical Methods, Prentice Hall, Upper Saddle River, NJ, USA, 1974.

S. P. Timoshenko and S. Woinowsky-Krieger, Theory of Plates and Shells, McGraw-Hill, London, UK, 1981.

J. K. Paik and A. K. Thayamballi, Ultimate Limit State Design of Steel-Plated Structures, John Wiley & Sons, Chichester, UK, 2003.

O. F. Hughes and J. K. Paik, Ship Structural Analysis and Design, The Society of Naval Architects and Marine Engineers, New Jersey, NJ, USA, 2010.

J. K. Paik, A. K. Thayamballi, S. K. Lee, and S. J. Kang, “Semi-analytical method for the elastic-plastic large deflection analysis of welded steel or aluminum plating under combined in-plane and lateral pressure loads,” Thin-Walled Structures, vol. 39, no. 2, pp. 125–152, 2001.

J. K. Paik and M. S. Lee, “A semi-analytical method for the elastic-plastic large deflection analysis of stiffened panels under combined biaxial compression/tension, biaxial in-plane bending, edge shear, and lateral pressure loads,” Thin-Walled Structures, vol. 43, no. 3, pp. 375–410, 2005.

G. Wang, H. Sun, H. Peng, and R. Uemori, “Buckling and ultimate strength of plates with openings,” Ships and Offshore Structures, vol. 4, no. 1, pp. 43–53, 2009.

U. N. Kim, I. H. Choe, and J. K. Paik, “Buckling and ultimate strength of perforated plate panels subject to axial compression: experimental and numerical investigations with design formulations,” Ships and Offshore Structures, vol. 4, no. 4, pp. 337–361, 2009.

P. P. Chaithanya, P. K. Das, A. Crow, and S. Hunt, “The effect of distortion on the buckling strength of stiffened panels,” Ships and Offshore Structures, vol. 5, no. 2, pp. 141–153, 2010.

K. M. El-Sawy and M. I. Martini, “Stability of biaxially loaded square plates with single central holes,” Ships and Offshore Structures, vol. 5, no. 4, pp. 283–293, 2010.

S. Benson, J. Downes, and R. S. Dow, “Ultimate strength characteristics of aluminium plates for high-speed vessels,” Ships and Offshore Structures, vol. 6, no. 1-2, pp. 67–80, 2011.

J. K. Paik, “Ultimate strength of steel plates with a single circular hole under axial compressive loading along short edges,” Ships and Offshore Structures, vol. 2, no. 4, pp. 355–360, 2007.

M. S. Kumar, P. Alagusundaramoorthy, and R. Sundaravadivelu, “Interaction curves for stiffened panel with circular opening under axial and lateral loads,” Ships and Offshore Structures, vol. 4, no. 2, pp. 133–143, 2009.