Zonal Harmonic Gravity Model

Calculate zonal harmonic representation of planetary gravity

• 

Libraries:
Aerospace Blockset / Environment / Gravity

Description

The Zonal Harmonic Gravity Model block calculates the zonal harmonic representation of planetary gravity at a specific location based on planetary gravitational potential. This block provides a convenient way to describe the gravitational field of a planet outside its surface.

By default, the block uses the fourth order zonal coefficient for Earth to calculate the zonal harmonic gravity. It also allows you to specify the second or third zonal coefficient.

For information on the planetary parameter values for each planet in the block implementation, see Algorithms.

Examples

Developing the Apollo Lunar Module Digital Autopilot

Develop Apollo Lunar Module digital autopilot using Simulink^® and Aerospace Blockset™.

• Open Live Script

Limitations

The block excludes the centrifugal effects of planetary rotation and the effects of a precessing reference frame.

Ports

Input

X_ecef — Planet-centered planet-fixed coordinates

m-by-3 matrix

Planet-centered planet-fixed coordinates, specified as an m-by-3 matrix, from the center of the planet in the selected length units. If Planet model has a value of Earth, this matrix contains Earth-centered Earth-fixed (ECEF) coordinates.

Data Types: double

Output

g_ecef — Gravity values

m-by-3 array

Gravity values, returned as an m-by-3 array, in the x-axis, y-axis and z-axis of the planet-centered planet-fixed coordinates, in the selected length units per second squared.

Data Types: double

Parameters

Units — Input units

Metric (MKS) (default) | English

Input units, specified as:

Units	Position	Equatorial Radius	Gravitational Parameter
Metric (MKS)	Meters	Meters	Meters cubed per second squared
English	Feet	Feet	Feet cubed per second squared

Programmatic Use

Block Parameter: units

Type: character vector

Values: 'Metric (MKS)' | 'English'

Default: 'Metric (MKS)'

Degree — Degree of harmonic model

4 (default) | 2 | 3

Degree of harmonic model, specified as.

• 2 — Second degree, J2. Most significant or largest spherical harmonic term, which accounts for the oblateness of a planet.

• 3 — Third degree, J3.

• 4 — Fourth degree, J4 (default).

Programmatic Use

Block Parameter: degree

Type: character vector

Values: '2' | '3' | '4'

Default: '4'

Action for out-of-range input — Out-of-range input behavior

Warning (default) | ErrorNone

Out-of-range input behavior, specified as:

Value	Description
None	No action.
Warning	Warning in the Diagnostic Viewer, model simulation continues.
Error	MATLAB® returns an exception, model simulation stops.

Programmatic Use

Block Parameter: action

Type: character vector

Values: 'None' | 'Warning' | 'Error'

Default: 'Warning'

Planet model — Planetary model

Mercury (default) | Venus | Earth | Moon | Mars | Jupiter | Saturn | Uranus | Neptune | Custom

Planetary model, specified as Mercury, Venus, Earth, Moon, Mars, Jupiter, Saturn, Uranus, Neptune, or Custom.

Selecting Custom enables you to specify your own planetary model.

• Selecting Mercury, Venus, Moon, Uranus, or Neptune limits the degree to 2.

• Selecting Mars limits the degree to 3.

Dependencies

Selecting Custom enables the Equatorial radius, Gravitational parameter and J values rate parameters.

Programmatic Use

Block Parameter: ptype

Type: character vector

Values: 'Mercury' | 'Venus' | 'Earth' | 'Moon' | 'Mars' | 'Jupiter' | 'Saturn' | 'Uranus' | 'Neptune' | 'Custom'

Default: 'Earth'

Equatorial radius — Planetary equatorial radius

6378136.3 (default) | scalar

Planetary equatorial radius, specified as a scalar, in the length units that the Units parameter defines.

Dependencies

To enable this parameter, set Planet model to Custom.

Programmatic Use

Block Parameter: R

Type: character vector

Values: scalar

Default: '6378136.3'

Gravitational parameter — Planetary gravitational parameter

398600441500000 (default) | scalar

Planetary gravitational parameter, specified as a scalar, in the length units cubed per second squared that the Units parameter defines.

Dependencies

To enable this parameter, set Planet model to Custom.

Programmatic Use

Block Parameter: GM

Type: character vector

Values: scalar

Default: '398600441500000'

J values — Zonal harmonic coefficients

[1.0826269e-03 -2.5323000e-06 -1.6204000e-06] (default) | 3-element array

Zonal harmonic coefficient, specified as a 3-element array.

Dependencies

To enable this parameter, set Planet model to Custom.

Programmatic Use

Block Parameter: jvalue

Type: character vector

Values: scalar

Default: '[1.0826269e-03 -2.5323000e-06 -1.6204000e-06]'

Algorithms

This block is implemented using the following planetary parameter values for each planet:

Planet	Equatorial Radius (Re) in Meters	Gravitational Parameter (GM) in m 3 /s 2	Zonal Harmonic Coefficients (J Values)
Earth	6378.1363e3	3.986004415e14	[ 0.0010826269 -0.0000025323 -0.0000016204 ]
Jupiter	71492e3	1.268e17	[0.01475 0 -0.00058]
Mars	3397.2e3	4.305e13	[ 0.001964 0.000036 ]
Mercury	2439.0e3	2.2032e13	0.00006
Moon	1738.0e3	4902.799e9	0.0002027
Neptune	24764e3	6.809e15	0.004
Saturn	60268e3	3.794e16	[0.01645 0 -0.001]
Uranus	25559e3	5.794e15	0.012
Venus	6052.0e3	3.257e14	0.000027

References

[1] Vallado, David, Fundamentals of Astrodynamics and Applications. New York: McGraw-Hill, 1997.

[2] Fortescue, P., J. Stark, G. Swinerd, eds.. Spacecraft Systems Engineering, 3d ed. West Sussex: Wiley & Sons, 2003.

[3] Tewari, A. Boston: Atmospheric and Space Flight Dynamics Modeling and Simulation with MATLAB and Simulink. Boston: Birkhäuser, 2007.

Extended Capabilities

C/C++ Code Generation
Generate C and C++ code using Simulink® Coder™.

Version History

Introduced in R2009b