Optimization of characteristics of multilayer spherical control joint-hinge stiffness
V. Gonca^{1} , S. Polukoshko^{2} , J. Shvab^{3} , A. Boiko^{4}
^{1, 3, 4}Riga Technical University, Institute of Mechanics, Ezermalas Str. 6, Riga, LV-1006, Latvia
^{2}Ventspils University College, Engineering Research Institute “VSRC”, Inzenieru Str. 101, LV-3601, Ventspils, Latvia
Journal of Vibroengineering, Vol. 15, Issue 1, 2013, p. 271-279.
Received 22 October 2012; accepted 4 March 2013; published 31 March 2013
In this work the analytical expression is derived on the basis of the variational method for the evaluation of angular stiffness characteristics of the spherical multilayered elastomeric package joint-hinge subjected to loading with hinge moment and hydrostatic pressure on the lateral surfaces. Joint-package consists of alternating metallic and thin elastomeric layers. Metallic plates-layers are assumed to be rigid. It is demonstrated that the hydrostatic pressure can both mitigate and "harden" the hinge angular stiffness characteristics depending on which side surface of the elastomeric layer it is applied. Obtained relationships allow solving problems of optimal design, operation and control program selection for compensating joints of this type.
Keywords: elastomers, packaged joint, optimization, hydrostatic pressure, Ritz’s method.
1. Introduction
Rubber and rubberlike materials (elastomers) are widely used in many sectors of industry [2, 4, 7]. Physical properties of elastomers, as polymeric materials, are qualitatively different from traditional construction materials because of their ability to maintain large elasticity deformations and small volume compressibility under deformation [3, 6]. Reinforced elastomeric structures (laminated elastomeric) consist of a large number of alternating thin layers of rubber and reinforcing layers of other, much more rigid than rubber, material. The connection of elastomer with reinforcing layer is usually done by means of vulcanization or gluing. Multilayer elastomeric structures have a special place: their axial compression stiffness is by several orders greater than the shear stiffness [6, 7, 10]. These structures are used in machine building, shipbuilding, civil engineering, aviation and aerospace due to its unique mechanical properties. Multilayer elastomeric structures successfully replace traditional technical systems, such as bearing, joints, compensating devices, shock-absorbers because of its important advantages: improving of machine dynamics, vibration and noise reduction, low shear and compression stiffness ratio. In practice, the packages of thin layers rubber-metal elements of different shapes are used: flat, cylindrical, conical, and others (Fig. 1). Number of layers may be different (at least three). In many applications of multilayer elastomeric structures their stiffness characteristics – the relationship between external forces imposed on the package and its displacement under load – are of great practical consequence.
Fig. 1. Examples of multilayer elastomeric structures
a) Rectangular
b) Ring
c) Cylindrical
d) Conical
e) Spherical
This paper considers a stiffening behavior of spherical multilayer package joint-hinge with perfectly rigid metallic reinforcing layers under influence of hinge moment and hydrostatic pressure exerted on the free side of the package. Strength and stability of elements are regarded as ensured. Flexible laminated spherical hinges are used in rocket vehicle to provide the control of direction of thrust vector by means of movable nozzle. Of all the mechanical deflection types, the movable nozzles are the most efficient to thrust vector control system [9]. Example of an elastomeric spherical hinge application in rocket technology is given in Fig. 2.
Reinforced elastomers have been in use for more than 50 years. There are a large number of publications concerning this structure. The theory of the elastomeric layer is developed [7]; methodology of determination of stress and deformations are proposed; compression, shear and tilting stiffness for round, square and infinite-strip form are determined theoretically [1, 2, 11]. For thin spherical elements there is no theoretical basis and method of stiffness calculation in the case of side follower load taking into account orientation of side surface changing [8, 10].
2. Problem definition, analytical model and method of analytical solution
A flexible hinge made up of alternate spherical segments of elastomeric and rigid metal shims is considered. The behavior of elastomeric layers is only taken into account, assuming that the metallic layers geometry makes them perfectly rigid at this scheme of loading. The joint-hinge is loaded with hinge moment ${M}_{h}$ (in $xoz$ plane) and with simultaneous action of hydrostatic pressure ${q}_{1}$ and ${q}_{2}$ on the free lateral surfaces of elastomeric layers, respectively, on the internal surface ($\phi ={\phi}_{1}$) and outside one ($\phi ={\phi}_{2}$). The purpose of this work is to obtain an analytical dependence for the angular stiffness characteristics of package-joint taking into consideration the influence of the pressure ${g}_{1}$ and ${q}_{2}$ on this characteristic. Loading scheme and the geometric parameters are shown in Fig. 3 in section $\theta =0$ of spherical coordinates $r$, $\phi $, $\theta $.
Fig. 2. Example of elastomeric hinge application
Fig. 3. Object geometry and loading scheme at $\theta =0$
Shear stiffness characteristics of joint-hinge are determined by means of Ritz's variation method using the linear theory of elasticity for an elastomeric material [6, 7]. At first, stiffness characteristic for one elastomeric layer is determined, i.e. joint with only one elastomeric layer is considered, then the package-joint response is studied.
2.1. Calculation of the stiffness characteristics of a single layer
The solution is carried out in two stages. Initially only hinge moment ${M}_{h}$ is imposed; loading and deformation scheme (in section $\theta =0$) is illustrated in Fig. 4. In the considered case the radial displacement $u$ and strains ${\epsilon}_{\phi}$ and ${\epsilon}_{\theta}$ are equal to zero. Geometrical boundary conditions for the displacement $w$ and $v$ (in the meridional and circumferential directions respectively) may be written as: if $r={R}_{1}$, $w=0$ and $v=0$, if $r={R}_{2}$, $w=\delta {R}_{2}\mathrm{c}\mathrm{o}\mathrm{s}\theta $ and $v=-\delta {R}_{2}\mathrm{c}\mathrm{o}\mathrm{s}\varphi \mathrm{s}\mathrm{i}\mathrm{n}\theta $.
Functions, satisfying these conditions and to the expected character of deformation are:
$\psi \left(r\right)=\frac{\left({r}^{3}-{R}_{1}^{3}\right){R}_{2}^{2}}{{r}^{2}({R}_{2}^{3}-{R}_{1}^{3})}.$
Sought dependence “moment – rotation angle” is defined by means of the Ritz method from the total potential energy functional minimum condition $\partial \Pi /\partial \delta =0$ [6]:
${\epsilon}_{r\varphi}=0.5\delta {R}_{2}r\frac{d}{dr}\left(\frac{\psi \left(r\right)}{r}\right)\mathrm{c}\mathrm{o}\mathrm{s}\theta ,\mathrm{}{\epsilon}_{r\theta}=-0.5\delta {R}_{2}r\frac{d}{dr}\left(\frac{\psi \left(r\right)}{r}\right)\mathrm{c}\mathrm{o}\mathrm{s}\varphi \mathrm{s}\mathrm{i}\mathrm{n}\theta .$
Obtained angular stiffness characteristic is:
where: $\alpha =\frac{{R}_{1}}{{R}_{2}}$.
Now the influence of pressure ${q}_{1}$ and ${q}_{2}$, acting respectively on the free sides of the elastomeric layer at $\phi ={\phi}_{1}$ and $\phi ={\phi}_{2}$ is_{}determined. First, the effect of the pressure ${q}_{2}$ is only taken into account with simultaneous action of hinge moment ${M}_{h}$. The lateral surface of the elastomeric layer $\phi ={\phi}_{2}$ in an arbitrary section $\theta $ is considered (Fig. 5).
Fig. 4. Loading and deformation scheme at $\theta =0$
Fig. 5. Deformation scheme in an arbitrary section
For thin elastomeric layers, even at small joint-hinge rotation angles $\delta $ reconfiguring of the elastomeric layers sides becomes significant. Since the pressure ${q}_{2}$ is follower load, which is always normal to the side surface of the elastomeric layer, significant change of its projections takes place relatively to the undeformed state of the elastomeric layer side. This leads to changing of the joint-hinge loading and must be taken into account in its analytical model. It is assumed that in view of the hinge geometry, its radial and flexural stiffness is significantly greater than the shear stiffness, i.e. radial and bending deformations, which can be induced by corresponding projections of pressure ${q}_{2}$_{}on the deformed surface, may be neglected. The elastomeric layer generator $AB$ under the joint-hinge moment moves meridionally (on $\phi $ - direction) at the angle $\beta =\angle ABD$ and in the circumferential direction $\theta $ at an angle $\gamma =\angle ABC$ and takes the position ${A}_{1}B$. Pressure ${q}_{2}$ remains normal to the deformed side of the elastomeric layer. Relative to the original surface $\phi ={\phi}_{2}$ vector ${q}_{2}$ of pressure can be decomposed into radial, circumferential and meridional direction components:
From Fig. 5: $AB={R}_{2}-{R}_{1},$$AD=\delta {R}_{2}\mathrm{c}\mathrm{o}\mathrm{s}\theta ,$$BC=\delta {R}_{2}\mathrm{c}\mathrm{o}\mathrm{s}{\varphi}_{2}\mathrm{s}\mathrm{i}\mathrm{n}\theta ,$$\mathrm{t}\mathrm{a}\mathrm{n}\left(\beta \right)=\frac{\delta \mathrm{c}\mathrm{o}\mathrm{s}\left(\theta \right)}{1-\alpha},$$\mathrm{t}\mathrm{a}\mathrm{n}\left(\gamma \right)=\frac{\delta \mathrm{c}\mathrm{o}\mathrm{s}\left({\varphi}_{2}\right)\mathrm{s}\mathrm{i}\mathrm{n}\left(\theta \right)}{1-\alpha},$ where $\alpha =\frac{{R}_{1}}{{R}_{2}}.$
Hinge shear takes place in $xoz$ plane. To determine the effect of ${q}_{2}$ pressure on the angular stiffness of swivel joint, it is necessary to calculate the moment $M\left({q}_{2}\right)$ in$xoz$ plane as the results of pressure ${q}_{2}$, acting on the lateral surface of the deformed elastomeric layer at $\phi ={\phi}_{2}$:
As soon as ${q}_{x}={q}_{u}\mathrm{c}\mathrm{o}\mathrm{s}\theta +{q}_{v}\mathrm{s}\mathrm{i}\mathrm{n}\theta \text{,}$${q}_{z}={q}_{u}\mathrm{s}\mathrm{i}\mathrm{n}\theta +{q}_{v}\mathrm{c}\mathrm{o}\mathrm{s}\theta \text{,}$$d{F}^{*}=\frac{r\mathrm{s}\mathrm{i}\mathrm{n}{\varphi}_{2}drd\theta}{\mathrm{c}\mathrm{o}\mathrm{s}\beta \mathrm{c}\mathrm{o}\mathrm{s}\gamma}\text{:}$
The integral (3) is not taken by elementary functions. Let us denote:
For the sufficiently thin elastomeric layers ${a}^{2}{\mathrm{s}\mathrm{i}\mathrm{n}}^{2}\theta <1$. Expanding the integrand in (3) into series relatively to ${a}^{2}{\mathrm{s}\mathrm{i}\mathrm{n}}^{2}\theta $ and leaving the first three members, after integration we obtain:
From the equations (3) and (4):
${k}_{2}\left(\delta \right)=\frac{1}{3}\left(1+\alpha +{\alpha}^{2}\right)\sqrt{1+\frac{{\delta}^{2}}{{\left(1-\alpha \right)}^{2}}}\left(1-\frac{3}{8}{a}^{2}-\frac{5}{64}{a}^{4}\right){\mathrm{c}\mathrm{o}\mathrm{s}}^{2}{\varphi}_{2}\mathrm{s}\mathrm{i}\mathrm{n}{\varphi}_{2}.$
Total angular hinge stiffness under the simultaneous action of hinge moment ${M}_{h}$ and pressure ${q}_{2}$ is defined from the condition of movable metallic layer equilibrium at $r={R}_{2}$ in the plane of the hinge bending $xoz$: ${M}_{h}-M\left({q}_{2}\right)=\delta \pi G{R}_{2}^{3}k.$ From this expression the value of hinge moment ${M}_{h}$, which provides a desired angle $\delta $ of hinge rotation at presence of the pressure ${q}_{2}$, is defined:
If the ${q}_{1}$ pressure on the side surface of an elastomeric layer $\phi ={\phi}_{1}$_{}is imposed, considering the scheme discussed above for ${q}_{2}$_{}(at $\phi ={\phi}_{2}$) and taking into consideration the opposite signs of the corresponding direction cosine to the surface $\phi ={\phi}_{1}$, we obtain the equation for the moment ${M}_{h}$, ensuring predetermined hinge rotation angle $\delta $:
${k}_{1}\left(\delta \right)=\frac{1}{3}\left(1+\alpha +{\alpha}^{2}\right)\sqrt{1+\frac{{\delta}^{2}}{{\left(1-\alpha \right)}^{2}}}\left(1-\frac{3}{8}{a}^{2}-\frac{5}{64}{a}^{4}\right){\mathrm{c}\mathrm{o}\mathrm{s}}^{2}{\varphi}_{1}\mathrm{s}\mathrm{i}\mathrm{n}{\varphi}_{1}.$
Consequently, the imposed pressure ${q}_{2}$_{}leads to increase of control hinge moment for a given rotation angle $\delta $ providing, i.e. to “mitigation” of angle stiffness; the imposed pressure ${q}_{1}$ at $\phi ={\phi}_{1}$ leads to decrease of control joint-hinge moment for a desired rotation angle providing, i.e. “toughen” the angular stiffness. Finally, under the simultaneous action of the pressure ${q}_{1}$ and ${q}_{2}$ on both sides of the elastomeric layer the control joint-hinge moment ${M}_{h}$ is equal to:
2.2. Calculation of stiffness characteristics of the package joint-hinge
For estimation of the stiffness characteristics of the packaged–joint consisting of N elastomeric layers, all of the above calculations are repeated. Deriving of the corresponding equations for each elastomeric layer one should keep in mind, that in a multilayer joint the moments caused by pressure ${q}_{1}$ and ${q}_{2}$ are transferred (and accumulated) from layer to layer – from mobile top layer of the packaged hinge to fixed lower layer. The introduced denotements for arbitrary $i$th layer are:
Then from expressions (7) and (8) for total packaged-joint we get a system of equations:
System (9) has $N+1$ equations and $N+1$ desired angular values $\delta $, ${\delta}_{1}$, ${\delta}_{2}$,…, ${\delta}_{i}$, ${\delta}_{N}$ and allows to find: angles of rotation ${\delta}_{i}$ separately for each layer and total angle of packaged joint-hinge $\delta $ for a given ${M}_{h}$, ${q}_{1}$ and ${q}_{2}$; hinge moment ${M}_{h}$, which provides the desired behavior of package – hinge for the given rotating angle $\delta $ and pressure ${q}_{1}$ and ${q}_{2}$ value; pressure value ${q}_{1}$ and ${q}_{2}$, providing the desired behavior mode when hinge moment ${M}_{h}$ and rotation angle $\delta $ of total hinge package are given.
Introducing the concept of the “middle” layer, from the system (9) for the sufficiently thin elastomeric layers the approximate, but sufficiently accurate analytical expression may be derived for the calculation of the angular packaged-joint stiffness for different combinations of values of ${M}_{h}$, ${q}_{1}$, ${q}_{2}$, $\delta $. For this purpose let us assume that all the elastomeric layers work as one “middle” layer, separated in the midsection of pack, the geometry of which is denoted as: ${R}_{1m}$, ${R}_{2m}$, ${\phi}_{1m}$, ${\phi}_{2m}$, ${k}_{m}$, ${k}_{1m}$, ${k}_{2m}$, $\delta m=\delta /N$ (where $N$ – the number of elastomeric layers in packaged joint). In this case, from the system (9) for the total angle of rotation of the packaged joint we have:
From the equations (8) - (10) the connection of control hinge moment ${M}_{h}$ with angle of rotation of all package hinge $\delta $ and pressure ${q}_{1}$ and ${q}_{2}$ is derived:
Here we consider small angles of joint rotation $\delta $, which provide hinge deformation (angles $\beta $), for which the linear relationships between the angular deformation and shear stresses in the material of the elastomeric layers are still valid. For shear and torsion it is valid for deformation up to 40-50 %. This condition imposes a limit on the number of layers of elastomeric and at a given angle $\delta $ and a given elastomeric layer geometry.
Let us consider separately the metallic and elastomeric layers. For small angles of rotation $\delta $ on the surface of the metal layer ${R}_{1m}<r<{R}_{2m}$, and $\phi ={\phi}_{2}$, we calculate the moment about axis $y$:
where $x=r\mathrm{s}\mathrm{i}\mathrm{n}\varphi \mathrm{c}\mathrm{o}\mathrm{s}\theta \text{,}$$z=r\mathrm{c}\mathrm{o}\mathrm{s}\varphi \text{,}$ the intensity of surface force on the metal layer for surface $\phi ={\phi}_{2}$ as the results of pressure ${q}_{2}$ is equal:
${Z}_{x}{|}_{\varphi ={\varphi}_{2}}={q}_{r}\mathrm{c}\mathrm{o}\mathrm{s}\varphi -{q}_{\varphi}\mathrm{s}\mathrm{i}\mathrm{n}\varphi ={q}_{2}\mathrm{s}\mathrm{i}\mathrm{n}{\varphi}_{2}\text{,}$
Then the equation (12) may be rewritten:
As it follows from (13), the presence of pressure does not influence joint shear provided that the metal layer is completely rigid (not deformable). It should be noted that for large angles of rotation $\delta $, which lead to an essentially nonlinear angular deformities in the elastomeric layer packet hinge, this influence is not zero. This is due to the fact that surface loads ${q}_{r}$, ${q}_{\varphi}$, ${q}_{\theta}$ on the metal layer at $\phi ={\phi}_{2}$ do not change since this metallic layer is not deformed, but arms vary: $x=r\mathrm{s}\mathrm{i}\mathrm{n}{\varphi}_{2}\mathrm{c}\mathrm{o}\mathrm{s}\theta +r\delta \mathrm{c}\mathrm{o}\mathrm{s}{\varphi}_{2}\mathrm{s}\mathrm{i}\mathrm{n}\theta $, $z=r\mathrm{c}\mathrm{o}\mathrm{s}\varphi +r\delta \mathrm{c}\mathrm{o}\mathrm{s}\theta $.
3. Numerical examples and experiment results
The example of theoretical calculation and results of experiment are given below for the multilayer joint-hinge with geometry: ${R}_{11}=$ 390 mm, ${R}_{2N}=$ 417 mm, ${h}_{elast}=$ 3 mm, ${h}_{met}=$ 3 mm, $N=$ 5, ${\phi}_{1}=$ 47°, ${\phi}_{2}=$ 57°, $\delta =$ 6° and shear modulus of the elastomeric layers $G=$ 460 kPa (Fig. 6).
Fig. 6. Testing packaged joint-hinge geometry
Theoretical calculations were performed in accordance with analytical expression (11), where the moment dependence on pressures ${q}_{1}$ and ${q}_{2}$ is linear due to small rotation angle and linear dependence between stress and strain. In Fig. 7 plots of joint-hinge moments ${M}_{h}$ dependence on hydrostatic side pressure ${q}_{1}$ and ${q}_{2}$_{}are shown for the simultaneous action of hinge moment ${M}_{h}$ and the pressure ${q}_{1}$ or ${q}_{2}$ (separetely (Fig. 7(a)) and together (Fig. 7(b)). In Fig. 8 plots of joint-hinge moments ${M}_{h}$ dependence on rotation angle $\delta $ and the hydrostatic side pressure ${q}_{1}$ and ${q}_{2}$ separately _{}are presented.
Fig. 7. Plots of joint - hinge moment ${M}_{h}$_{}(kNm) dependence on the side pressures ${q}_{1}$_{}and ${q}_{2}$ (MPa) for $\delta =$ 6°
a)${M}_{h}={M}_{h}\left({q}_{1}\right)$, ${M}_{h}={M}_{h}\left({q}_{2}\right)$
b)${M}_{h}={M}_{h}\left({q}_{1},{q}_{2}\right)$
Experimental investigation of the stiffness characteristics of this joint were performed at the Moscow Institute of Thermal Technology. Testing was fulfilled under static rotation of joint with fixation of angle (Fig. 9). Because of the complexity of the experiment, only the pressure ${q}_{1}$ was applied to the inner surface side and integral characteristics “moment-rotation” were measured. In Fig. 10 plots of experimental and theoretical joint stiffness characteristics $K$ (kNm/deg) dependence on the hydrostatic side pressure ${q}_{1}$ are shown. Due to the fact that shear modulus of rubber sample and joint elastomeric layer are different, results of experiment and calculations somewhat differ. In multilayer joint $G$ of elastomeric layer may reach 490 kPa.
Dependence between ${M}_{h}$ and ${q}_{1}$ is linear to implement the requirements of constant deflection angle $\delta $. Nonlinear dependence appears under simultaneous increasing of ${M}_{h}$ and ${q}_{1}$.
The calculations indicate and experimental data confirm that for sufficiently thin elastomeric layers pressure ${q}_{1}$_{}and ${q}_{2}$ acting on the free sides of elastomeric layers essentially affects the angular stiffness characteristics of the hinge - package that must be taken into account both during design and operation of such class of joint-hinges for the considered type of loading.
Fig. 8. Plots of joint - hinge moment ${M}_{h}$ (kNm) dependence on the side pressures ${q}_{1}$, ${q}_{2}$_{}(MPa) and rotation angle $\delta $ (rad)
a)${M}_{h}={M}_{h}\left({q}_{1},\delta \right)$
b)${M}_{h}={M}_{h}\left({q}_{2},\delta \right)$
Fig. 9. Scheme of experimental setup
Fig. 10. Plot of experimental and theoretical dependence of hinge rigidity on the side pressure ${q}_{1}$
4. Conclusions
Analytical expression is derived on the basis of the variational method for estimation of angular stiffness characteristics of the spherical multilayer elastomeric package-joint during its loading with hinge moment and hydrostatic pressure on the lateral surfaces. The theoretical calculations are verified experimentally, which indicate that at the shear of spherical packaged- joint with rather thin elastomeric layers under the hinge moment ${M}_{h}$ the presence of hydrostatic pressure ${q}_{1}$_{}and ${q}_{2}$, acting on the free sides of the elastomeric layers, influences the angular stiffness characteristics of the packaged hinge. The derived analytical expressions for the angular stiffness characteristics allow to take into account the effect of the above noted in the design and operation of this type of joints for the considered loading case. In the future research work it is planned to include the case of geometric nonlinearity.
Acknowledgements
This work has been supported by the ERDF’s Grant, within the project “SATTEH”, No. 2010/0189/2DP/2.1.1.2.0/10/APIA/VIAA/019, being implemented in Engineering Research Institute “Ventspils International Radio Astronomy Centre” of Ventspils University College (VIRAC).
References
- Amin A. F. M. S., Lion A., Sekita S., Okui Y. Nonlinear dependence of viscosity in modeling the rate-dependent response of natural and high damping rubbers in compression and shear: experimental identiﬁcation and numerical veriﬁcation. International Journal of Plasticity, Elsevier, Vol. 22, 2006, p. 1610 1657, www.elsevier.com/locate/ijplas [Search CrossRef]
- Bauman J. T. Fatigue, Stress and Strain of Rubber Components: Guide for Design Engineers. Munich, Carl Hanser Verlag, 2008. [Search CrossRef]
- Frolov N. N., Moldovanov S. J., Lozovoy S. B. Mechanics of Thin Rubber Elements. Kuban State Technological University, Krasnodar, Publishing House “Jug”, 2011. [Search CrossRef]
- Gent A. N. Engineering with Rubber: How to Design Rubber Components. 2-nd Ed., Munich, Carl Hanser Verlag, 2001. [Search CrossRef]
- Kelly J. M., Konstantinidas D. A. Mechanics of Rubber Bearings for Seismic and Vibration Isolation. John Wiley & Sons, UK, 2011. [Search CrossRef]
- Lavendel E. E. Design Methods of Workpieces from Highly Elastic Materials. Riga, Zinatne, 1990. [Search CrossRef]
- Malkov V. M. Mechanics of Multilayered Elastomeric Structures. St-Petersburg, St.-Petersburg University Press, 1998. [Search CrossRef]
- Mtenga P. V. Elastomeric Bearing Pads under Combined Loading. Final Report, Contract No: BC352_16, 2007, www.dot.state.fl.u..._STR/FDOT_BC352_16_rpt.pdf [Search CrossRef]
- Pisacane V. L. Fundamentals of Space System. 2-nd Edition, Oxford University Press, 2005. [Search CrossRef]
- Semenenko V. P. Some tasks of performance of rubber joints. State Space Agency of Ukraine, Scientific and Technical Collection Space Technology, Missile Armament, Issue 1, 2012, p. 149 161. [Search CrossRef]
- Tsai H.-C. Flexure analysis of circular elastic layers bonded between rigid plates. International Journal of Solids and Structures, No. 40, 2003, p. 2975 2987. [Search CrossRef]