Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 548708, 13 pages

http://dx.doi.org/10.1155/2014/548708

Research Article

Geometric Nonlinear Meshless Analysis of Ribbed Rectangular Plates Based on the FSDT and the Moving Least-Squares Approximation

L. X. Peng,1,2 Yue-Ping Tao,³ Hong-Qiao Li,1,2 and Gui-Kai Mo¹

¹College of Civil Engineering and Architecture, Guangxi University, Nanning 530004, China

²Guangxi Key Laboratory of Disaster Prevention and Engineering Safety, Guangxi University, Nanning 530004, China

³Guangxi Polytechnic of Construction, Nanning 530003, China

Received 17 December 2013; Revised 14 January 2014; Accepted 21 January 2014; Published 13 March 2014

Academic Editor: Yumin Cheng

Copyright © 2014 L. X. Peng 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.

Based on the first-order shear deformation theory (FSDT) and the moving least-squares approximation, a new meshless model to study the geometric nonlinear problem of ribbed rectangular plates is presented. Considering the plate and the ribs separately, the displacement field, the stress, and strain of the plate and the ribs are obtained according to the moving least-squares approximation, the von Karman large deflection theory, and the FSDT. The ribs are attached to the plate by considering the displacement compatible condition along the connections between the ribs and the plate. The virtual strain energy formulation of the plate and the ribs is derived separately, and the nonlinear equilibrium equation of the entire ribbed plate is given by the virtual work principle. In the new meshless model for ribbed plates, there is no limitation to the rib position; for example, the ribs need not to be placed along the mesh lines of the plate as they need to be in FEM, and the change of rib positions will not lead to remeshing of the plate. The proposed model is compared with the FEM models from pieces of literature and ANSYS in several numerical examples, which proves the accuracy of the model.

1. Introduction

Ribbed plate has been widely used in engineering, such as bridges, ship hulls, and aviation, and it is a popular structure with obvious advantages. The ribs make the structure stiffer and allow it to achieve larger bearing capacity than flat plate with roughly the same weight. However, the ribs also bring difficulties to analysis, and the calculation of ribbed plate is more complicated than that of flat plates. Based on the fact that the ribs of many ribbed plates are attached to the plate with uniform spacing and close to one another, and that ribbed plates show different elastic characteristics in the two perpendicular directions, early researchers transformed the ribs to an addition layer to the plate and used the orthotropic model to approximate the ribbed plates [1]. Another early model was the grillage model [2]. The models were simple and fulfilled the demand of fast and easy computation in engineering. Therefore, they are still used in some design environments, where accurate analysis is not the first concern. However, because the models were introduced in the age of lacking computational tools and some approximations were adopted, they cannot give satisfying results in solving generalized ribbed plate problems. Due to the advances of computers and numerical methods in the past decades, a ribbed plate model which has more universality was introduced, regarding the ribbed plate as a composite structure of ribs and plate and analyzing them separately; combining them by imposing the displacement compatibility conditions between them, this is also the model which is accepted widely. Several methods have been developed, such as the Rayleigh-Ritz method [3–7] and the finite element methods (FEM) [8, 9].

Not many nonlinear analyses of ribbed plates can be found in pieces of literature, and most of them were based on the FEM [10–13], which benefits from the good adaptability and high accuracy of the method. However, no method is perfect, and FEMs also have disadvantages. The FEMs rely on the meshes that discretize problem domain to construct their approximated solution, but the large deformation of problem domain always leads to mesh disorder, and, therefore, time-consuming and accuracy-suffering remeshing is unavoidable, which brings difficulties to both programming and analysis. And for ribbed plate problems, most FEMs require that the ribs are placed along the mesh lines and any change in their positions will lead to the remeshing of the plate domain to accommodate the change. If the layout of ribs needs to be optimized, there may be hundreds of times of remeshing before obtaining the ideal result. And if the optimization is carried out under the consideration of the large deformation of a ribbed plate, the number of remeshing may become a dramatic figure combing the iteration from the nonlinear analysis. Meshless or meshfree method [14–19] is a numerical method which bases their approximated solution entirely on a set of nodes distributed in a problem domain. Without a mesh, the meshless methods overcome the aforementioned difficulties that FEM encountered with the meshes. The moving least-squares approximation originated in data fitting. Nayroles et al. [14] were the first to use a moving least-squares procedure to develop a meshless approximation. By introducing moving least-squares interpolants to construct the trial and test functions for the variation principle (weak form), Belytschko et al. [15] improved the method proposed by Nayroles et al. [14] and proposed the element-free Galerkin method (EFG). Nevertheless, due to the fact that the shape function of most meshless methods lacks Kronecker delta properties, and that the unknowns of the governing equation are nodal parameters other than nodal displacements, the displacement compatible conditions between the components of a composite structure cannot be implemented directly in meshless methods as they can in FEMs when the structure is analyzed, which limits the application of meshless methods in engineering. Recently, the analysis of plate and composite structure with meshless methods has made some progress. Lei et al. [20, 21] analyzed buckling and large deformation of functionally graded plate using the element-free kp-Ritz method. Zhang et al. [22] used a local Kriging meshless method to study the thermal buckling of functionally graded plates. The author Peng and his coworkers have proposed meshless methods to solve the linear bending, free vibration, and elastic buckling problems of ribbed plates with a derived transformation equation to address the nodal parameter issue [23, 24]. However, the equation did not consider all necessary displacement compatible conditions, which leads to failure in solving the large deformation problem of ribbed plates. The objective of this paper is to propose a meshless model to study the geometric nonlinear behaviors of ribbed plates from the perspective of composite structure. Based on the first-order shear deformation theory (FSDT), the moving least-squares approximation (MLS) and von Karman’s large deflection theory, the displacement field, nonlinear strains, and nonlinear equilibrium equations of the plate and ribs are derived. A new equation that transforms the nodal parameters of the ribs to those of the plate is introduced, and the equation allows the displacement compatible condition between the plate and the ribs to be implemented directly. Because of the meshless characteristics of the proposed model, the ribs need not to be placed along the mesh lines of the plate, and the change of rib positions will not lead to remeshing of the plate. Mesh disorder due to the large deformation of problem domain is avoided, as well. Some numerical examples are utilized to demonstrate the accuracy of the proposed model. The calculated results are compared with the results from ANSYS and pieces of literature. The proposed meshless model of ribbed plate can provide a substantial ground for future optimization of rib layout under the consideration of large deformation.

2. Moving Least-Squares Approximation

In MLS [15], a function defined in a domain Ω can be approximated by in the subdomain . is defined as

where are the monomial basis functions, is a factor that measures the domain of influence of the nodes, is the number of basis function, and are their coefficients. In this paper, the quadratic basis (, in 1D); (, in 2D) are used for the ribs and plates, respectively. The unknown coefficients can be determined by minimizing a weighted discrete norm

where or is the weight function, = 0 outside Ωx, is the number of nodes in that makes the weight function > 0, and are the nodal parameters. Minimizing Γ with respect to ,

we obtain

where

Therefore, (1) can be expressed in a standard form as

where are the shape functions.

3. Meshless Model of a Ribbed Plate

The meshless model of a ribbed plate, shown in Figure 1, is composed of a plate and an -rib. The plate and the rib are discretized by a set of nodes. The degree of freedom (DOF) of every node of plate is , where , , and are the nodal translations of the plate in , and directions, respectively. and are the rotation about the -axis and the -axis, respectively. The DOF of every node of -rib is . The rib is assumed to be made from the same material as the plate. The Young’s modulus is and the Poisson’s ratio is . If there is a -rib, we can derive similar equations for -rib as those for the -rib.

Figure 1: Meshless model of a ribbed plate.

3.1. Displacement Field Approximation of a Ribbed Plate

Based on the FSDT [25, 26], the displacement field of a plate is given as

According to the MLS approximation, the functions , and can be expressed in a discrete form

where are the nodal parameters of the th node of the plate, is the number of nodes of the plate, and are independent of . Similarly, the displacement field of the -rib is

where are the nodal parameters of the -rib and is the number of nodes of the -rib. The shape functions and are obtained from (6), and the cubic spline function is used as the weight function.

3.2. A New Transformation Equation of the Nodal Parameters

Along the axis of the -rib (Figure 1), we take a normal section parallel to -axis, as shown in Figure 2.

Figure 2: Section of plate and x-rib.

For a node of the -rib, there will be a corresponding point on the plate, and they have the same and coordinates. Their displacements follow

And necessarily, at the corresponding point of Node and Point in the contact surface between the plate and -rib, Point is

Remark 1. Point can be any point on the plate that corresponds to Node . For every node of the -rib, a corresponding point on the plate can be found. Similar to the process in [23, 24], the equations that transform the nodal parameters and of the -rib into the nodal parameters of the plate can be derived as

where

is an matrix, is an vector, is an square matrix, and is an vector. Every row of and matrices corresponds to a node of the rib, and, therefore, the matrices have rows. Equation (12) also gives equations as

or

where is the thickness of the plate and is the depth of the -rib. According to the FSDT, (16) can be written as

And because    (11),

The discrete form of (18) is

or

where , , and . Equation (20) leads to

For concentrically ribbed plates, just take . The combination of (13) and (21) gives a new equation that expresses the nodal parameters of the -rib in terms of the nodal parameters of the plate as follows:

where  , , , and is a matrix.

With this new transformation equation (22), the meshless model for ribbed plate is more applicable. No matter how the position of ribs changes, we only need to recalculate . Therefore, compared with a finite element model, this meshless model for ribbed plate is expected to have more advantages in future optimization analysis of rib position.

3.3. Strains and Stresses of the Plate and Rib

According to the FSDT and the von Karman theory, the strains of an isotropic plate are

For convenience, is rewritten as

where the linear component of the strain includes

where

The nonlinear component of the strain

where

The stress is

where

The strain of an isotropic -rib is

which can be rewritten as

where the linear component of the strain is

where

The nonlinear component of the strain is

where , .

The stress is

where

3.4. Nonlinear Formulations of the Ribbed Plate

The virtual work equation of the ribbed plate is

where is the sum of the internal and external force vectors, F ₁ is the sum of all of the loading vectors, is the virtual displacement and is the virtual strain of the plate, and is the virtual strain of the x-rib. We can write the relation between the strain and nodal parameters in an increment form as

Substituting (22) into (42), we have

The substitution of (41) and (43) into (40) gives us the nonlinear equilibrium equation of the entire ribbed plate as follows:

where

is the function of , where

is the function of , where

From (44), we have

Employing (29), (35), (41), and (43), we have

Equation (45) gives us

Substituting (49), (50), and (51) into (48), we obtain

where

where

We can obtain and for x-rib similarly. Therefore, (52) can be written as

where

From (44), we obtain

where

3.5. Solution to the Nonlinear Equilibrium Equations

This paper used the Newton-Raphson method to solve nonlinear equilibrium equation (57). The process is as follows.

(1)Take linear solution as the first approximated solution .(2)Substituting into (57) to calculate .(3)Employ (55) to determine .(4)Obtain the incremental displacements by

 and the improved solution

(5)Return to step and repeat steps to until is sufficiently small and is the final solution.

3.6. Enforcement of Essential Boundary Conditions

Due to a lack of Kronecker delta properties in the shape functions given in (6), it is difficult to impose the essential boundary conditions. The full transformation method that Chen et al. [17] introduced is adopted in this paper to enforce the essential boundary conditions.

4. Results and Discussion

4.1. Validation Studies

To show the convergence of the proposed model, and the influences of the support size of the nodes and the order of the basic functions, a clamped square plate subjected to a uniformly distributed pressure of 100 Pa is studied. The width of the plate is 1.8 m, and the thickness is 0.018 m. The Young’s modulus of the plate is  Pa and the Poisson’s ratio is .

The nonlinear deflection of the central point of the plate that is obtained by the proposed model under different support sizes (which are denoted by scaling factors, ) and different completeness orders of the basic functions is shown in Figure 3, and is compared with the solution that is given by the finite element software, ANSYS, using SHELL63 element.

Figure 3: Nonlinear deflection of the flat plate under different and N c.

In this paper, rectangular support is employed, and thus the scaling factors and are defined by

where , are the lengths of the rectangular support of nodes in the and directions, respectively, and , are the distances between the two neighboring nodes in the and directions, respectively. For convenience, we choose . From Figure 3, it can be observed that for a certain meshless scheme (in this case nodes) and all of the solutions for different completeness orders of the basic functions converge when the support size () is larger than 5. Higher completeness orders (N c) need a larger support size to make the solution converge.

Secondly, we vary the meshless scheme and obtain the variations of the nonlinear central deflection under certain completeness order of the basic functions , which are shown in Figures 4, 5, and 6, respectively. The solution that is given by ANSYS is also in the figures for comparison. Figures 4 to 6 indicate that for certain , the solution converges when the number of nodes increases. For certain , the solutions for larger support sizes () converge before those for smaller support sizes do.

Figure 4: Variation of nonlinear deflection of the flat plate, .

Figure 5: Variation of linear deflection of the flat plate, .

Figure 6: Variation of linear deflection of the flat plate, .

From the studies, we find that when the order of basic functions and the support size for the plate, the solutions are precise enough with a relatively lower computational cost. A similar convergence study for ribs can be carried out, and the solution given under and the support size for the rib is found to be satisfying. Therefore, all of the following examples are calculated with , (for the plate), and (for the ribs).

4.2. Rectangular Plate with One Rib

A rectangular plate clamped at two opposite sides and with one rib (Figure 7) located at  mm is studied.

Figure 7: Rectangular plate with one rib.

The two other sides of the plate are free. Both the plate and the rib are made of the same material, with Young’s modulus  GPa and Poisson’s ratio . The plate is under a uniformly distributed load in the -direction, and the q ranges from 0.2 to 0.8 MPa. The central deflections of the plate obtained by the proposed model and given by Koko and Olson [10] under different loads are shown in Table 1. The present solution is based on the following discrete scheme: 13 × 13 for the plate and 13 for the rib.

Table 1: Central deflection of the rectangular plate with one rib under different loads.

Load (MPa)
Koko and Olson [10] (mm)
Present results (mm)
Relative errors

0.2
5.526
5.31654
3.8%

0.3
7.172
6.9253
3.4%

0.4
8.631
8.47243
1.8%

0.5
9.868
9.82917
0.4%

0.6
11.053
11.0406
0.1%

0.7
11.974
12.1373
−1.4%

0.8
12.961
13.1409
−1.4%

4.3. Square Plate with Two Cross Ribs

A simply supported square plate with two ribs is studied. The ribs are located at  m and  m, and they are made of the same material (Figure 8) as the plate.

Figure 8: Square plate with two cross ribs.

Young’s modulus and Poisson’s ratio are  GPa and , respectively. Under a uniformly distributed load in the -direction, the load-central deflection curves of the plate given by the proposed model and a large deflection analysis of FEM are shown in Figure 9. When the plate is fully clamped, the load-central stress curves of the top surface of the plate are shown in Figure 10. The solution of FEM is obtained using ANSYS software with a three-dimensional model (SOLID45 element is used). The number of discrete elements is 30416 (Figure 11), and two layers are used along the thickness. Figure 9 shows good agreement. In Figure 10, one can observe that the present results given by a discrete scheme of nodes have some differences from the results of ANSYS, but the results given by a discrete scheme of nodes agree well with the results of ANSYS. It illustrates that the proposed model need more discrete nodes to calculate stress than it need to calculate displacement.

Figure 9: Load-central deflection curve of the simply supported square plate with two ribs.

Figure 10: Load-stress curve of the clamped square plate with two ribs.

Figure 11: 3D finite element model of a square plate with two ribs.

If the thickness of the plate is increased to 0.1 m and the plate is clamped, the load-central deflection curves of the plate are shown in Figure 12.

Figure 12: Load-central deflection curve of the clamped square plate with two ribs.

4.4. Square Plate with One Rib and Different Discretization Schemes

A square plate with one rib located at  m (Figure 13) is studied. The rib is made of the same material as the plate. Consider  Pa and .

Figure 13: Square plate with one rib.

The plate is under a uniformly distributed load q that ranges from 0 to 0.1 MPa in the -direction. To test the new transformation equation (22) that was derived in Section 3.2, the central deflection of the plate is given by ANSYS and the proposed model (with different meshless schemes). In ANSYS, a three-dimensional model for the ribbed plate is used, and the number of discrete elements (SOLID45 element) is 1800 (Figure 14).

Figure 14: 3D finite element model of a square plate with one rib.

Firstly, the present results (Figure 17) are computed using 81 plate nodes (uniform and nonuniform distribution, as shown in Figures 15 and 16, resp.).

Figure 15: Nonuniform distribution of the 81 plate nodes.

Figure 16: Uniform distribution of the 9 × 9 plate nodes.

Figure 17: Central deflection of the clamped plate with one rib (81 nodes).

Secondly, the present results (Figure 20) are calculated by discretizing the plate with 121 nodes (uniform and nonuniform distribution, as shown in Figures 18 and 19, resp.).

Figure 18: Nonuniform distribution of the 121 plate nodes.

Figure 19: Uniform distribution of the 11 11 plate nodes.

Figure 20: Central deflection of the clamped plate with one rib (121 nodes).

Finally, the present results (Figure 23) are given by a discretization scheme of 169 plate nodes (uniform and nonuniform distribution, as shown in Figures 21 and 22).

Figure 21: Nonuniform distribution of the 169 plate nodes.

Figure 22: Uniform distribution of the 13 × 13 plate nodes.

Figure 23: Central deflection of the clamped plate with one rib (169 nodes).

In Figures 15, 16, 18, 19, 21, and 22, the small circles represent the nodes of the plate and the black dots denote the nodes of the stiffener. Figures 17, 20, and 23 display clearly that, in the cases of uniform distribution, Point is taken to be a node of the plate, while, in the cases of nonuniform distribution, Point is no longer a node of the plate (see Section 3.2). In Figure 17, the agreement is good but slight difference exists. When we increase the nodes of plate, the agreement becomes much better (Figures 20 and 23). The results in Figures 17, 20, and 23 demonstrate that the placement of the plate nodes nearly has no effects on the final solution and prove the accuracy of the new transformation equation (22) that was derived in Section 3.2.

5. Conclusions

This paper presents a meshless model, which is based on the FSDT and the MLS approximation to study the geometric nonlinear behaviors of ribbed plate structures. Considering a ribbed plate as a composite structure of plate and ribs, and starting from the large deflection theory of von Karman, the nonlinear behaviors of the plate and the ribs were studied, respectively. Then, employing the meshless advantages of the proposed model, the nonlinear governing equations of the plate and the ribs were superposed with a new equation derived for the nodal parameter transformation of the plate and the ribs, and the geometric nonlinear equilibrium equation of the entire structure is established. The advantages of the proposed model are that the ribs can be placed anywhere on a plate and any changes of their positions will not lead to the remeshing of the plate, which enhances computational efficiency in solving the optimization of rib layout under the consideration of nonlinear deformation. And the proposed model do not rely on mesh; therefore, mesh disorder due to the large deformation of problem domain is avoided. The present results are compared with those from three-dimensional FEM analysis and pieces of literature. Good agreement can be observed, which proves the accuracy of the proposed meshless model.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The work that is described in this paper has been supported by the grants awarded by the National Natural Science Foundation of China (Projects nos. 11102044, 51168003, and 51168005) and the Systematic Project of Guangxi Key Laboratory of Disaster Prevention and Structural Safety (Project no. 2012ZDX07).

References

S. Kendrick, “The analysis of a flat plated grillage,” European Shipbuilding, vol. 5, pp. 4–10, 1956.

H. A. Schade, “The orthogonally stiffened plate under uniform lateral load,” Journal of Applied Mechanics ASME, vol. 62, pp. 143–146, 1940.

C. L. Kirk, “Vibration of centrally stiffened rectangular plate,” Journal of the Royal Aeronautical Society, vol. 65, pp. 695–697, 1961.

K. M. Liew, Y. Xiang, and S. Kitipornchai, “Transverse vibration of thick rectangular plates—III. Effects of multiple eccentric internal ring supports,” Computers and Structures, vol. 49, no. 1, pp. 59–67, 1993.

K. M. Liew, Y. Xiang, S. Kitipornchai, and M. K. Lim, “Vibration of rectangular Mindlin plates with intermediate stiffeners,” Journal of Vibration and Acoustics, vol. 116, no. 4, pp. 529–535, 1994.

K. M. Liew, Y. Xiang, S. Kitipornchai, and J. L. Meek, “Formulation of Mindlin-Engesser model for stiffened plate vibration,” Computer Methods in Applied Mechanics and Engineering, vol. 120, no. 3-4, pp. 339–353, 1995.

Y. Xiang, S. Kitipornchai, K. M. Liew, and M. K. Lim, “Vibration of stiffened skew Mindlin plates,” Acta Mechanica, vol. 112, no. 1–4, pp. 11–28, 1995.

J. R. O'Leary and I. Harari, “Finite element analysis of stiffened plates,” Computers and Structures, vol. 21, no. 5, pp. 973–985, 1985.

A. Deb and M. Booton, “Finite element models for stiffened plates under transverse loading,” Computers and Structures, vol. 28, no. 3, pp. 361–372, 1988.

T. S. Koko and M. D. Olson, “Non-linear analysis of stiffened plates using super elements,” International Journal for Numerical Methods in Engineering, vol. 31, no. 2, pp. 319–343, 1991.

D. V. Rao, A. H. Sheikh, and M. Mukhopadhyay, “A finite element large displacement analysis of stiffened plates,” Computers and Structures, vol. 47, no. 6, pp. 987–993, 1993.

B. Chattopadhyay, P. K. Sinha, and M. Mukhopadhyay, “Geometrically nonlinear analysis of composite stiffened plates using finite elements,” Composite Structures, vol. 31, no. 2, pp. 107–118, 1995.

A. H. Sheikh and M. Mukhopadhyay, “Geometric nonlinear analysis of stiffened plates by the spline finite strip method,” Computers and Structures, vol. 76, no. 6, pp. 765–785, 2000.

B. Nayroles, G. Touzot, and P. Villon, “Generalizing the finite element method: diffuse approximation and diffuse elements,” Computational Mechanics, vol. 10, no. 5, pp. 307–318, 1992.

T. Belytschko, Y. Y. Lu, and L. Gu, “Element-free Galerkin methods,” International Journal for Numerical Methods in Engineering, vol. 37, no. 2, pp. 229–256, 1994.

Y. Cheng and J. Li, “Complex variable meshless method for fracture problems,” Science in China G, vol. 49, no. 1, pp. 46–59, 2006.

J. S. Chen, C. Pan, C. T. Wu, and W. K. Liu, “Reproducing Kernel Particle Methods for large deformation analysis of non-linear structures,” Computer Methods in Applied Mechanics and Engineering, vol. 139, no. 1–4, pp. 195–227, 1996.

M. Peng, P. Liu, and Y. Cheng, “The complex variable element-free Galerkin (CVEFG) method for two-dimensional elasticity problems,” International Journal of Applied Mechanics, vol. 1, no. 2, pp. 367–385, 2009.

H. Gao and Y. Cheng, “A complex variable meshless manifold method for fracture problems,” International Journal of Computational Methods, vol. 7, no. 1, pp. 55–81, 2010.

Z. Lei, K. M. Liew, and J. L. Yu, “Buckling analysis of functionally graded carbon nanotube reinforced composite plates using the element-free kp-Ritz method,” Composite Structures, vol. 98, pp. 160–168, 2013.

Z. Lei, K. M. Liew, and J. L. Yu, “Large deflection analysis of functionally graded carbon nanotube-reinforced composite plates by the element-free kp-Ritz method,” Computer Methods in Applied Mechanics and Engineering, vol. 256, pp. 189–199, 2013.

L. W. Zhang, P. Zhu, and K. M. Liew, “Thermal buckling of functionally graded plates using a local Kriging meshless method,” Composite Structures, vol. 108, pp. 472–492, 2014.

L. X. Peng, K. M. Liew, and S. Kitipornchai, “Buckling and free vibration analyses of stiffened plates using the FSDT mesh-free method,” Journal of Sound and Vibration, vol. 289, no. 3, pp. 421–449, 2006.

L. X. Peng, K. M. Liew, and S. Kitipornchai, “Analysis of stiffened corrugated plates based on the FSDT via the mesh-free method,” International Journal of Mechanical Sciences, vol. 49, no. 3, pp. 364–378, 2007.

J. N. Reddy, Theory and Analysis of Elastic Plates, Taylor & Francis, London, UK, 1999.

K. M. Liew, Y. Xiang, and S. Kitipornchai, “Research on thick plate vibration: a literature survey,” Journal of Sound and Vibration, vol. 180, no. 1, pp. 163–176, 1995.