rel="nofollow" href="#ulink_9f2627f1-42cf-56fb-a6c5-8883a50b2dbc">[1.8] Factorization of ODE [1.4] is (Vakhromeev and Kornev 1972)
The general solution of ODE [1.9] is given by
where functions W0 and W1,2 are the general solutions of the following equations:
[1.10]
For large frequencies (aω ≫1), the following estimates for the derivatives of functions W0 and W1,2 are obtained:
The behavior of these solution components is different: W1 is the rapidly oscillating function, and W2 is the sum of exponentials rapidly decreasing from the edges of the beam.
Therefore, the situation under consideration is fundamentally different from the case when the characteristic equation has small and large modulo roots, which is typical for boundary layer theory. In our case, we are talking about the separation of solutions, one of which oscillates at the same rate as the EE decays (i.e. the characteristic equation has large real and imaginary roots with moduli of the same order). The self-adjoint eigenvalue problem [1.1], [1.2] can be referred to as the boundary value problem with quasi-separable variables (Bolotin 1960a, 1960b, 1961a, 1961b, 1961c; Bolotin et al. 1950, 1961).
We proceed to the construction of the EE described by equation [1.11]. Taking into account the expression for the natural frequency [1.8], we obtain the following relations for EEs localized in the vicinity of the edges x = 0 and x = L, respectively:
We assume that the beam is so long such that EEs do not affect each other, i.e. exp(−πλ-1L) ≪ 1.
Now, it remains to find the quantities x0, λ and constants C1, C2 from the boundary conditions to determine the eigenmodes and eigenfrequencies:
Note that in original works by Bolotin, a slightly different matching procedure is used. Namely, the matching conditions [1.14], [1.15] are not given at the domain boundaries. They are set at some points, which are then determined from the solution of the system of transcendental equations along with other unknown constants.
Substituting expressions [1.7] and [1.12] into conditions [1.14], and expressions [1.7] and [1.13] into [1.15], we obtain
[1.16]
[1.17]
For λ and x0, we have the following expressions:
Finally, the formula for the natural oscillation frequencies of the clamped beam is written as follows:
The same formula follows from the analysis of the transcendental equation (Weaver et al. 1990). Its difference from the numerical solution is less than 1% even for the first natural oscillation frequency.
Using DEEM to solve various problems proved its high accuracy even when calculating the first eigenfrequencies (Bolotin 1960a, 1960b, 1961a, 1961b, 1961c, 1961d; Gavrilov 1961a, 1961b; Bolotin 1963, 1970; Elishakoff and Wiener 1976; Elishakoff and Steinberg 1979; Emmerling 1979; Bolotin 1984; Elishakoff et al. 1993, 1994). It confirms “Crighton’s principle” (Crighton 1994):
All experience suggests that asymptotic solutions are useful numerically far beyond their nominal range of validity, and can often be used directly, at least at a preliminary product design stage, for example, saving the need for accurate computation until the final design stage where many variables have been restricted to narrow ranges.
The problem of natural oscillations was considered above. In the case of forced oscillations, the response is given as a series expansion in terms of normal modes of natural oscillations. It is important that one can use only the rapidly oscillating part of solution W0, neglecting the dynamic edge effect (DEE) (Bolotin 1961b, 1961c, 1970, 1984).
1.3. Linear problems solved
DEEM in its original form can be effectively used for calculating natural, free and forced oscillations of rectangular plates, shells of revolution and shallow shells with a rectangular base.
In particular, DEEM is widely used for the analysis of rectangular isotropic (Bolotin et al. 1958; Bolotin 1961d, 1963; Kudryavtsev
