APPENDIX 1
FIGURE A1. Plot of superellipse exponents k (xaxis) against superellipse correction factors C_{m} (crosssectional areas, mass, and the contribution of a frustrum's point mass to I of an entire body) and C_{I} (for an individual frustrum's I) relative to an ellipse (yaxis). The polynomials enable area calculations with exponents of greater than one significant figure (Table A1).
TABLE A1. Superellipse exponents (k) and coefficients for determining mass (C_{mass}) and mass moment of inertia (C_{I}) relative to these quantities for an elliptical frustrum (k=2.0). C_{mass} is for conversion of both mass and the contribution of a frustrum's point mass to I of a body part from its center of rotation. C_{I} is for conversions of I_{x,y,z} for any individual frustrum. The inverse of a coefficient gives mass or I for an ellipse if the quantity is known for a shape with another exponent k.
Exponent k 
Superellipse: 
Superellipse: 
C_{mass} 
C_{I} 

1 
0.4292 
0.4244 
1.1 
0.5644 
0.5006 
1.2 
0.6641 
0.5731 
1.3 
0.7417 
0.6412 
1.4 
0.8032 
0.7048 
1.5 
0.8527 
0.7639 
1.6 
0.8931 
0.8188 
1.7 
0.9265 
0.8695 
1.8 
0.9551 
0.9164 
1.9 
0.9789 
0.9598 
2 
1.0000 
1.0000 
2.1 
1.0177 
1.0372 
2.2 
1.0338 
1.0716 
2.3 
1.0475 
1.1036 
2.4 
1.0599 
1.1333 
2.5 
1.0705 
1.1609 
2.6 
1.0805 
1.1866 
2.7 
1.0894 
1.2105 
2.8 
1.0972 
1.2328 
2.9 
1.1041 
1.2537 
3 
1.1108 
1.2732 
APPENDIX 2.
Derivation of mass moment of inertia I for superelliptical frustra
Solid Edge and other modeling software calculate inertial properties (mass, center of mass, and mass moment of inertia). The geometry for these properties will be especially accurate with NURBs ellipsoids deformed to drape over a skeleton (Mallison, 2010, 2011), for example, as in an animal's thoracic region. For a model with nondeformed ellipsoid cross sections, it is useful to estimate how differing crosssectional shapes would affect mass m and mass moments of inertia I about orthogonal axes.
Symmetrical, curved cross sections, called superellipses, have shapes described by equation (5).
5)
Here, x and y are variables, a and b are semimajor and semiminor axes, and exponent k governs the shape's curvature. In an ellipse, k=2.
For an object with superelliptical cross sections, the following equation gives mass moment of inertia about a vertical axis y, with n segments (each i), and radii from dorsal to ventral (DV) and lateral to medial (LM).
6)
A given segment's density is _{i }, its length is l_{i}, and and are respective averages of dorsoventral and mediolateral radii of the segment's anterior and posterior faces. The distance r_{i} is the radius from the center of rotation of a body to the center of mass m_{i} of segment i. Critically, C is a unique constant for a given superelliptical cross section of exponent k from equation 5 (as long as k is the same for front and back cross sections).
F_{max} for a superellipse, and hence for I and m for objects with superelliptical sections, can be calculated using exponent k. Mosen (2009) provides a general equation for F_{max} (although not for I or m), using the inverse of k. The equation becomes simpler when k is calculated and substituted algebraically.
7)
B is a beta function incorporating gamma functions defined as
8)
Substituting all terms and expanding equation 3, we obtain mass moment of inertia of an individual superellipsebased frustrum i about a vertical axis.
9)
Mass and I_{y} vary predicably with the exponent used to define each superellipse cross section, and can be calculated by multiplying the ellipse values by a correction factor. Table A1 presents correction factors for superelipsesection frustra of varying exponents k, including that of an ellipse (k=2). Figure A1 plots these coefficients, and fitted polynomials that enable calculation of areas with exponents of higher significant figures. Common crosssectional shapes for terrestrial vertebrates, with exponents of 2.3–2.5 (Motani 2001), result in masses 5–7% greater than for elliptical cross sections. This discrepancy supports Motani's (2001) findings that elliptical cross sections substantially underestimate body volume. I y of individual frustra increase more rapidly than mass as exponents exceed k=2, from 10–16% greater as k ranges from 2.3 to 2.5.
APPENDIX 3.
Congruent axes for joint markers in MSC Adams.
Joint markers on different bodies in Adams must be coincident. Aligning their positions, coordinate axes and axis orientations requires four steps.
1. Rightclick on the joint and obtain "Info". This dialog box identifies the i and j markers.
2. Copy the angles of the i marker's global axis orientation.
3. Rightclick the j marker, and select "Modify".
4. Paste the i marker's angles in place of the listed j marker's.