Governing equations and boundary conditions next up previous
Next: The fundamental solution Up: BEM FOR STRUCTURAL ANALYSIS Previous: BEM FOR STRUCTURAL ANALYSIS

Governing equations and boundary conditions

The theory of linear elastostatics may be found in numerous textbooks, like, for example, that of Timoshenko[7]. The basic equations are summarized here for the sake of completeness. The coordinates system ( tex2html_wrap_inline91 ), the displacements ( tex2html_wrap_inline93 ), the stresses ( tex2html_wrap_inline95 ), and the surface tractions ( tex2html_wrap_inline97 ) are shown in Figure 4.

  figure169
Figure 4:  The main coordinates, displacements, stresses, and tractions in three dimensions.

The equilibrium equations:

  eqnarray175

where tex2html_wrap_inline99 are the body force (per unit volume) components. The surface tractions are related to the stresses via:

equation179

where tex2html_wrap_inline101 are the components of the outward unit vector normal to the solid surface.

The strains are given in terms of the deformations as:

equation182

The stresses are related to the strains via the constitutive relationships:

equation188

where E is the modulus of elasticity, and tex2html_wrap_inline105 the Poisson's ratio; tex2html_wrap_inline107 is the Kronecher delta, and tex2html_wrap_inline109 .

For the sake of simplicity we will limit ourselves to two dimensions. The extension to three dimensions is straightforward. The next two sections summarize the fundamentals of the boundary element formulation in elastostatics, as described in more detail by Brebbia and Dominguez [8].

In two dimensions we may consider either the plain strain ( tex2html_wrap_inline111 ) or the plate stretching ( tex2html_wrap_inline113 ) situation. The latter can be considered equivalent to the former by taking as tex2html_wrap_inline115 and tex2html_wrap_inline117 . The following formulation corresponds to the plain strain condition. However, all the numerical examples that will follow correspond to the plate stretching condition.

  figure206
Figure 5:  The boundary conditions on the ``hydrofoil'' beam.

The ``hydrofoil'' beam shown in Figure 5 is considered. The pressures (tractions) acting on the sides of the beam are taken equal to those resulting from the hydrodynamic analysis of a cavitating 2-D hydrofoil with the same section. The pressures have been taken with reference to the pressure inside the cavity. They are also shown in Figure 5 together with the rest of the boundary conditions. The beam is taken to be clamped at the ``trailing edge'' of the hydrofoil. It should be noted that this 2-D beam problem has no physical equivalent in three dimensions, in which case a beam with a hydrofoil cross section is clamped to a body (like the hull of a ship or the hub of a propeller). However, we expect this 2-D problem to be representative for the investigation of the accuracy of the BEM with decreasing hydrofoil thickness to chord ratio. Notice that from the four quantities, tex2html_wrap_inline119 , only the two need to be specified on each boundary. The horizontal component of the pressure forces acting on the wetted side is ignored ( tex2html_wrap_inline121 ).


next up previous
Next: The fundamental solution Up: BEM FOR STRUCTURAL ANALYSIS Previous: BEM FOR STRUCTURAL ANALYSIS

Baris Gucun
Tue Mar 4 18:15:49 CST 1997