gp.system

Syntax

l := gp.system(gp)

gp.system(gp) := l

Get/set the local coordinate system system used by a grid point during equations of motion. This is given as a matrix (list of vectors). The first entry is the local x direction, and so on. If an orthonormal coordinate system is not specified, FLAC will automatically adjust the axes to make them orthonormal.

Returns:

l - A list of vectors describing the local coordinate system of the grid point.

Accepts:

l - A list of vectors describing the local coordinate system of the grid point.

Arguments:

gp - a zone gridpoint pointer

Usage Example

The following example illustrates how gp.system can be used

model new
model large-strain off

zone create brick size 1 1 1
zone face skin
fish define localAxis ; switch x and z axis in local system
    local matA = matrix(3,3)
    matA(1,1) = 0
    matA(1,2) = 0
    matA(1,3) = 1
    matA(2,2) = 0
    matA(2,1) = 1
    matA(2,3) = 0
    matA(3,3) = 1
    matA(3,2) = 0
    matA(3,1) = 0
    loop foreach gp1 gp.list
        gp.system(gp1) = matA ; switch x and z axis only at bottom face
    endloop
end
[localAxis]
zone gp fix vel-x 1 range group 'bottom' ; velocity is applied in the z-direction