Chapter 6: Applications of Double Integration
Section 6.6: Second Moments
Example 6.6.10
Find the moments of inertial Ix and Iy, the total mass m, and the radii of gyration Rx and Ry of the lamina with density ρ=3+2 x+3 y and whose shape is that of R, the region bounded by the ellipse x2+4 y2=1.
See Example 6.5.10.
Solution
Mathematical Solution
With ρr,θ=3+r2 cosθ+3 sinθ and φθ=1/4−3 cos2θ, the calculations for the moments of inertia and the radii of gyration are given in Table 6.6.10(a), where the relevant integrations are implemented in polar coordinates. An expression for the ellipse is given by r=φθ.
m=∫02 π∫0φθr⋅ρ ⅆr ⅆθ = 3 π2
Ix=∫02 π∫0φθr⋅ρ⋅r sinθ2 ⅆr ⅆθ = 3 π32
Iy=∫02 π∫0φθr⋅ρ⋅r cosθ2 ⅆr ⅆθ = 3 π8
Rx=Iy/m=3 π/83 π/2=14=1/2
Ry=Ix/m=3 π/323 π/2=116=1/4
Table 6.6.10(a) Moments of inertia and radii of gyration
The polar-coordinate representation of the ellipse, namely, r=φθ, is obtained by converting the Cartesian representation to polar, as follows. Since x2+4 y2=1 ⇒ r2cos2θ+4 sin2θ=1,
r=1cos2θ+41−cos2θ=14−3 cos2θ
Maple Solution - Interactive
The total mass is calculated in Table 6.6.10(b) via a task template that integrates over an ellipse via the modified int command in the Student VectorCalculus package.
Tools≻Tasks≻Browse:
Calculus - Vector≻Integration≻Multiple Integration≻2-D≻Over an Ellipse
Integrate fx,y over an Ellipse
fx,y=
Equation of ellipse:
From θ= to θ=
Table 6.6.10(b) Total mass calculated with a task template that integrates over an ellipse
The ellipse is expressed in polar coordinates as
12⁢tan⁡θ2+114+tan⁡θ2
=sec2θ1+4 tan2θ
=1/cos2θ+4 sinθ2
=1/cos2θ+41−cos2θ
=1/4−3 cos2θ
Similarly, the moments of inertia could be obtained with this same task template by appropriately modifying the integrand fx,y. Alternatively, as shown in Table 6.6.10(c), the appropriate calculations can be implemented from first principles.
Define the density ρr,θ
Context Panel: Assign Name
ρ=3+r 2 cosθ+3 sinθ→assign
Compute the total mass m
Calculus palette: Iterated double-integral template
Context Panel: Evaluate and Display Inline
Context Panel: Assign to a Name≻m
∫02 π∫01/4−3 cos2θr⋅ρ ⅆr ⅆθ = 32⁢π→assign to a namem
Obtain the moments of inertia
Context Panel: Assign to a Name≻Ix (or Iy, as appropriate)
∫02 π∫01/4−3 cos2θr⋅ρ⋅r sinθ2 ⅆr ⅆθ = 332⁢π→assign to a nameIx
∫02 π∫01/4−3 cos2θr⋅ρ⋅r cosθ2 ⅆr ⅆθ = 38⁢π→assign to a nameIy
Obtain the radii of gyration
Rx=Iy/m
Rx=12
Ry=Ix/m
Ry=14
Table 6.6.10(c) Moments of inertial and radii of gyration from first principles
Maple Solution - Coded
Initialize
Tools≻Load Package: Student Multivariate Calculus
Loading Student:-MultivariateCalculus
Define the density function ρ.
ρ≔3+2 x+3 y:
Total mass
Use the MultiInt command in the Student MultivariateCalculus package.
MultiIntρ,x,y=Ellipsex2+4 y2=1,output=integral;m≔ MultiIntρ,x,y=Ellipsex2+4 y2=1
32⁢π
Moments of Inertia
Use the MultiInt command in the Student MultivariateCalculus package. b
MultiInty2⋅ρ,x,y=Ellipsex2+4 y2=1,output=integral;Ix≔ MultiInty2⋅ρ,x,y=Ellipsex2+4 y2=1
332⁢π
MultiIntx2⋅ρ,x,y=Ellipsex2+4 y2=1,output=integral;Iy≔ MultiIntx2⋅ρ,x,y=Ellipsex2+4 y2=1
38⁢π
Radii of gyration
Rx=sqrtIy/m
Ry=sqrtIx/m
<< Previous Example Section 6.6 Next Chapter >>
© Maplesoft, a division of Waterloo Maple Inc., 2024. All rights reserved. This product is protected by copyright and distributed under licenses restricting its use, copying, distribution, and decompilation.
For more information on Maplesoft products and services, visit www.maplesoft.com
Download Help Document