60 HIGHER-ORDER ELEMENTS
torsional rigidity J obtaining
Γ4-560 for finite elements
J = 2
νάΩ =
Λ
Ω
[4-487 for boundary elements
The approximate solutions are about 10 % of the analytical values.
Notice that better agreement is obtained for the problem variable
than for the torsional rigidity, which is computed approximately from
the approximate values of the problem variables.
Table
3.1 VALUE OF TORSIONAL RIGIDITY AND V FUNCTION FOR INTERNAL POINTS
r
Finite elements Boundary elements
i
iii
x
act
/f
. /i.
(linear) (linear)
0
035
0
0-75
0
0-75
0
0-350
0-414
0-638
0-566
0-782
0-638
0-800
0-341
0-392
0-627
0-561
0-790
0-665
0-793
0-334
0-401
0-626
0-557
0-772
0-629
0-791
Rigidity 5026 4-560 4-487
Example 3.2
Figure 3.4 describes the aerofoil shape called NACA 0018 which can
be studied using boundary elements. The surface of the aerofoil was
* 100
Figure 3.4 NACA 0018 aerofoil
divided into a series of elements. Because of symmetry only half the
aerofoil needs to be considered. The problem can be expressed in
terms of a stream function u for a perfect fluid, such that the velocities
are