Obtaining wave kinematic data in an Abaqus/Aqua analysis

Utility routines GETWAVE, GETWAVEVEL, GETWINDVEL, and GETCURRVEL are provided to access the fluid kinematic data for an Abaqus/Aqua analysis.

These routines can be used only from within user subroutine UEL.

The following topics are discussed:

GETWAVE (get wave kinematics)
GETWAVEVEL, GETWINDVEL, and GETCURRVEL (get wave, wind, and current velocities)

GETWAVE (get wave kinematics)

Utility routine interface

       PARAMETER(MWCOMP=number of wave components)
       DIMENSION WAMP(MWCOMP),WPERD(MWCOMP),WXLAMB(MWCOMP),
     1 WPHI(MWCOMP),WOFF(3),WANG(2,MWCOMP)
       ...
       CALL GETWAVE(MWCOMP,NWCOMP,WAMP,WPERD,WXLAMB,WPHI,WOFF,WANG,
     1 ELEVB,ELEVS,JWTYPE,JRCD)

Variables returned from the utility routine

NWCOMP

Number of wave components (always 1 for Stokes wave theory).

WAMP

Array containing the amplitude of the wave components.

WPERD

Array containing the period of the wave components.

WXLAMB

Array containing the wavelength of the wave components.

WPHI

Array containing the phase angle of the wave components.

WOFF

Used only for gridded wave data (JWTYPE=2), when WOFF gives the position of the origin of the gridded coordinate system with respect to the global system.

WANG(2,*)

For Stokes fifth-order wave theory WANG(1,1) and WANG(2,1) are the direction cosines of wave travel. For Airy wave theory WANG(1,K1) and WANG(2,K1) are the direction cosines of the direction of travel of the K1th wave. For gridded wave data WANG(1,1) and WANG(2,1) are the direction cosines of the wave data grid. In all cases these direction cosines are with respect to the global coordinate system.

ELEVB

User-defined elevation of the seabed.

ELEVS

User-defined elevation of the still water surface.

JWTYPE

Integer flag indicating the wave type, as follows:

JWTYPE=0	Airy wave theory
JWTYPE=1	Stokes fifth-order wave theory
JWTYPE=2	Wave data obtained from gridded values

JRCD

The error flag JRCD is returned from GETWAVE as 0 if all the wave kinematic data are read correctly and as −1 if an error occurred (for instance, NWCOMP is greater that MWCOMP).

GETWAVEVEL, GETWINDVEL, and GETCURRVEL (get wave, wind, and current velocities)

Utility routine interface

      CALL GETWAVEVEL (NDIM, X, V, A, LERROR, NOEL, XINTERMED)
      CALL GETWINDVEL (NDIM, X, V, NOEL, XINTERMED)
      CALL GETCURRVEL (NDIM, X, V, NOEL, XINTERMED)

Variables to be provided to the utility routine

NDIM: Dimensionality of the element. It should be set to 2 for two-dimensional cases (for example, beams in a plane) and 3 for three-dimensional cases (for example, beams in space).
X(1..NDIM): Global coordinates of the point.

Variables returned from the utility routine

V(1..NDIM): Velocity components in the global coordinate system.
A(1..NDIM): Wave acceleration components in the global coordinate system. This variable is returned by GETWAVEVEL only.
LERROR: For gridded wave data LERROR is returned as 0 if the current point is within the grid or above the crest; it is returned as 1 if the point is outside the bounds of the grid. For Airy and Stokes waves LERROR is always returned as 0. If LERROR is returned as 1, the global coordinates of the nearest grid point are returned in X. LERROR is returned by GETWAVEVEL only.
NOEL: Element number.
XINTERMED(NDIM): An array containing the intermediate configuration coordinates of the load integration point. For nonstochastic analysis this array is not used. In a stochastic analysis the wave field is based upon this configuration. Additional details are found in UWAVE.