Custom Variable Mass 3DOF (Body Axes)

Implement three-degrees-of-freedom equations of motion of custom variable mass with respect to body axes

• 

Libraries:
Aerospace Blockset / Equations of Motion / 3DOF

Description

The Custom Variable Mass 3DOF (Body Axes) block implements three-degrees-of-freedom equations of motion of custom variable mass with respect to body axes. It considers the rotation in the vertical plane of a body-fixed coordinate frame about a flat Earth reference frame. For more information about the rotation and equations of motion, see Algorithms.

Ports

Input

F_x — Applied force along x-axis

scalar

Applied force along the body x-axis, specified as a scalar, in the units selected in Units.

Data Types: double

F_z — Applied force along z-axis

scalar

Applied force along the body z-axis, specified as a scalar, in the units selected in Units.

Data Types: double

M — Applied pitching moment

scalar

Applied pitching moment, specified as a scalar, in the units selected in Units.

Data Types: double

dm/dt — Rate of change of mass

scalar

Rate of change of mass (positive if accreted, negative if ablated), specified as a scalar.

Dependencies

To enable this port, select Include mass flow relative velocity.

Data Types: double

m — Mass

scalar

Mass, specified as a scalar, in the units selected in Units.

Data Types: double

dI/dt — Rate of change of inertia tensor

scalar

Rate of change of inertia tensor, I_yy, specified as scalar, in the units selected in Units.

Dependencies

To enable this port, set Mass type to Custom Variable.

Data Types: double

I — Inertia tensor

scalar

Inertia tensor, specified as a scalar, in the units selected in Units.

Dependencies

To enable this port, set Mass type to Custom Variable.

Data Types: double

g — Gravity

scalar

Gravity, specified as a scalar.

Dependencies

To enable this port, set Gravity source to External.

Data Types: double

V_re — Relative velocity

two-element vector

Relative velocity at which mass is accreted to or ablated from the body in body-fixed axes, specified as a two-element vector.

Dependencies

To enable this port, select Include mass flow relative velocity.

Data Types: double

Output

θ — Pitch attitude

scalar

Pitch attitude, within ±pi, returned as a scalar, in radians (θ).

Data Types: double

q — Pitch angular rate

scalar

Pitch angular rate, returned as a scalar, in radians per second.

Data Types: double

X_eZ_e — Location of body

two-element vector

Location of the body in the flat Earth reference frame, (Xe, Ze), returned as a two-element vector.

Data Types: double

U w — Velocity of body

two-element vector

Velocity of the body resolved into the body-fixed coordinate frame, (u, w), returned as a two-element vector.

Data Types: double

g — Gravity

scalar

Gravity, returned as a scalar.

Dependencies

To enable this port, select the Output mass properties for acceleration computation parameter.

Data Types: double

Parameters

Main

Units — Input and output units

Metric (MKS) (default) | English (Velocity in ft/s) | English (Velocity in kts)

Input and output units, specified as Metric (MKS), English (Velocity in ft/s), or English (Velocity in kts).

Units	Forces	Moment	Acceleration	Velocity	Position	Mass	Inertia
Metric (MKS)	Newton	Newton-meter	Meters per second squared	Meters per second	Meters	Kilogram	Kilogram meter squared
English (Velocity in ft/s)	Pound	Foot-pound	Feet per second squared	Feet per second	Feet	Slug	Slug foot squared
English (Velocity in kts)	Pound	Foot-pound	Feet per second squared	Knots	Feet	Slug	Slug foot squared

Programmatic Use

Block Parameter: units

Type: character vector

Values: Metric (MKS) | English (Velocity in ft/s) | English (Velocity in kts)

Default: Metric (MKS)

Axes — Body or wind axes

Body (default) | Wind

Body or wind axes, specified as Wind or Body.

Programmatic Use

Block Parameter: axes

Type: character vector

Values: Wind | Body

Default: Body

Mass type — Mass type

Custom Variable (default) | Fixed | Simple Variable

Mass type, specified according to the following table.

Mass Type	Description	Default for
Fixed	Mass is constant throughout the simulation.	3DOF (Body Axes) 3DOF (Wind Axes)
Simple Variable	Mass and inertia vary linearly as a function of mass rate.	Simple Variable Mass 3DOF (Body Axes) Simple Variable Mass 3DOF (Wind Axes)
Custom Variable	Mass and inertia variations are customizable.	Custom Variable Mass 3DOF (Body Axes) Custom Variable Mass 3DOF (Wind Axes)

The Custom Variable selection conforms to the previously described equations of motion.

Programmatic Use

Block Parameter: mtype

Type: character vector

Values: Fixed | Simple Variable | Custom Variable

Default: 'Custom Variable'

Initial velocity — Initial velocity of body

100 (default) | scalar

Initial velocity of the body, (V₀), specified as a scalar.

Programmatic Use

Block Parameter: v_ini

Type: character vector

Values: '100' | scalar

Default: '100'

Initial body attitude — Initial pitch altitude

0 (default) | scalar

Initial pitch attitude of the body, (θ₀), specified as a scalar.

Programmatic Use

Block Parameter: theta_ini

Type: character vector

Values: '0' | scalar

Default: '0'

Initial body rotation rate — Initial pitch rotation rate

0 (default) | scalar

Initial pitch rotation rate, (q₀), specified as a scalar.

Programmatic Use

Block Parameter: q_ini

Type: character vector

Values: '0' | scalar

Default: '0'

Initial incidence — Initial angle

0 (default) | scalar

Initial angle between the velocity vector and the body, (α₀), specified as a scalar.

Programmatic Use

Block Parameter: alpha_ini

Type: character vector

Values: '0' | scalar

Default: '0'

Initial position (x,z) — Initial location

[0 0] (default) | two-element vector

Initial location of the body in the flat Earth reference frame, specified as a two-element vector.

Programmatic Use

Block Parameter: pos_ini

Type: character vector

Values: '[0 0]' | two-element vector

Default: '[0 0]'

Gravity Source — Gravity source

Internal (default) | External

Gravity source, specified as:

External	Variable gravity input to block
Internal	Constant gravity specified in mask

Programmatic Use

Block Parameter: g_in

Type: character vector

Values: 'Internal' | 'External'

Default: 'Internal'

Acceleration due to gravity — Gravity source

9.81 (default) | scalar

Acceleration due to gravity, specified as a double scalar and used if internal gravity source is selected. If gravity is to be neglected in the simulation, this value can be set to 0.

Dependencies

• To enable this parameter, set Gravity Source to Internal.

Programmatic Use

Block Parameter: g

Type: character vector

Values: '9.81' | scalar

Default: '9.81'

Include mass flow relative velocity — Mass flow relative velocity port

off (default) | on

Select this check box to add a mass flow relative velocity port. This is the relative velocity at which the mass is accreted or ablated.

Programmatic Use

Block Parameter: vre_flag

Type: character vector

Values: off | on

Default: 'off'

Output mass properties for acceleration computation — Option to include mass properties for acceleration port

off (default) | on

Select this check box to enable ports for mass properties for acceleration. You can then use these ports as inputs for these blocks:

• 3DOF Acceleration — m output port

• 3DOF Angular Acceleration — I and dI/dt coeff ports.

Programmatic Use

Block Parameter: mass_flag

Type: character vector

Values: 'off' | 'on'

Default: off

State Attributes

Assign a unique name to each state. You can use state names instead of block paths during linearization.

• The number of names must match the number of states, as shown for each item, or be empty. Set all or none of the block states.

• To assign names to single-variable states, enter unique names between quotes, for example, 'q' or "q".

• To assign names to two-variable states, enter a comma-separated list surrounded by braces, for example, {'Xe','Ze'}.

• If a state parameter is empty (' '), no name is assigned.

• To assign state names with a variable in the MATLAB^® workspace, enter the variable without quotes. A variable can be a character vector, cell array of character vectors, or string.

Velocity: e.g., {'u, 'w'} — Velocity state name

'' (default) | comma-separated list surrounded by braces

Velocity state names, specified as a comma-separated list surrounded by braces.

Programmatic Use

Block Parameter: vel_statename

Type: character vector

Values: '' | comma-separated list surrounded by braces

Default: ''

Position: e.g., {'Xe', 'Ze'} — Position state name

'' (default) | comma-separated list surrounded by braces

Position state names, specified as a comma-separated list surrounded by braces.

Programmatic Use

Block Parameter: pos_statename

Type: character vector

Values: '' | comma-separated list surrounded by braces

Default: ''

Pitch angular rate e.g., 'q' — Pitch angular rate state name

'' (default)

Pitch angular rate state name, specified as a character vector or string.

Programmatic Use

Block Parameter: q_statename

Type: character vector | string

Values: '' | scalar

Default: ''

Pitch attitude: e.g., 'theta' — Pitch attitude state name

'' (default)

Pitch attitude state name, specified as a character vector or string.

Programmatic Use

Block Parameter: theta_statename

Type: character vector | string

Values: ''

Default: ''

Algorithms

The Custom Variable Mass 3DOF (Body Axes) block considers the rotation in the vertical plane of a body-fixed coordinate frame about a flat Earth reference frame.

The equations of motion are

$$\begin{array}{l}{A}_{xb}=\dot{u}=A_{xe}-qw\\ A_{zb}=\dot{w}=A_{ze}+qu\\ A_{xe}=\frac{\left(F_x-\dot{m}{u}_{re}\right)}{m}-g\sin \theta \\ A_{ze}=\frac{\left(F_z-\dot{m}{w}_{re}\right)}{m}+g\cos \theta \\ {\dot{X}}_e=u\cos \theta +w\sin \theta \\ {\dot{Z}}_e=-u\sin \theta +w\cos \theta \\ \dot{q}=\frac{M_y-{\dot{I}}_{yy}q}{I_{yy}}\\ \dot{\theta }=q\end{array}$$

where the applied forces are assumed to act at the center of gravity of the body. Input variables are F_x, F_z, M_y, $$\dot{m}$$ (dm/dt), m, $$\dot{I}$$ (dI_yy/dt), and I_yy. u_re, w_re, and g are optional input variables.

Extended Capabilities

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

Version History

Introduced in R2006a

R2026a: 3DOF equations of motion block acceleration ports removed

As a result of the removal of the acceleration ports, the Include inertial acceleration parameter has been replaced with Output mass properties for acceleration computation for Custom Variable Mass 3DOF (Body Axes).

Related to this change, the block output ports have changed:

Block	Output Ports Removed	Output Ports Added
Custom Variable Mass 3DOF (Body Axes)	dq/dt AxbAzb AxeAze	g

R2021b: Custom Variable Mass 3DOF (Body Axes) Block Changes

The 3DOF equations of motion have been updated. Existing models created prior to R2021b that contain 3DOF equations of motion blocks continue to run. If you replace a pre-R2021b version of a 3DOF equation of motion block with an R2021b or later version, your updated model might have a higher tendency for algebraic loops. For an example of how to remove algebraic loops using unit delays, see Remove Algebraic Loops. For further information about algebraic loops, see Identify Algebraic Loops in Your Model.