Making Animation with Scilab

A block shape in a general configuration as shown in Figure 1 can be defined by the position and the orientation of the reference frame attached to one of the vertices and the length of each side (Lx, Ly, Lz). If the position coordinate and the orientation matrix of the reference frame are * r* and

where **p**^{0}_{i} is the coordinate of the i-th vertex in the basic configuration. If the coordinates of the vertices are obtained, we can make the block shape by combining six polygons using `plot3d`

command in the same way as in the last section.

Figure 1: Block shape (in a general configuration)

A Scilab script to draw a block shape in a general configuration is shown in Code 1. The result of the drawing by this script is shown in Figure 2. Function `eulerXYZ()`

used in the script is a function that converts XYZ Euler angles to the corresponding rotation matrix and is defined as in Code 2. Also, `transver()`

is a function that transforms a vertices matrix according to the coordinate transformation formula shown above and is defined as in Code 3.

// make_block_general.sce clear; xdel(winsid()); exec('eulerXYZ.sci', -1); exec('transver.sci', -1); exec('genpat.sci', -1); // Reference position r = [1; 1; 1]; // Reference orientation R = eulerXYZ(-%pi/3, 0, %pi/6); // Side lengths Lx = 0.15; Ly = 0.05; Lz = 0.30; // Vertices vertices_0 = [ 0, 0, 0; // #1 Lx, 0, 0; // #2 0, Ly, 0; // #3 0, 0, Lz; // #4 Lx, Ly, 0; // #5 0, Ly, Lz; // #6 Lx, 0, Lz; // #7 Lx, Ly, Lz]; // #8 vertices = transver(vertices_0, r, R); // Faces faces = [ 1, 2, 5, 3; // #1 1, 3, 6, 4; // #2 1, 4, 7, 2; // #3 4, 7, 8, 6; // #4 2, 5, 8, 7; // #5 3, 6, 8, 5]; // #6 // Patches patches = genpat(vertices, faces); // Draw patches h_fig = figure; h_fig.background = 8; h_pat = plot3d(patches.x, patches.y, patches.z); h_pat.color_mode = 4; h_pat.foreground = 1; h_pat.hiddencolor = 4; // Axes settings xlabel("x"); ylabel("y"); zlabel("z"); h_axes = gca(); h_axes.isoview = "on"; h_axes.box = "off"; h_axes.rotation_angles = [63.5, -127]; h_axes.data_bounds = [0.8, 0.9, 0.8; 1.3, 1.4, 1.3]; xgrid;

Code 1: Scilab script to draw a block shape (in a general configuration)

function R = eulerXYZ(a1, a2, a3) // Convert XYZ Euler angles to rotation matrix R1 = [ 1, 0, 0; 0, cos(a1), -sin(a1); 0, sin(a1), cos(a1)]; R2 = [ cos(a2), 0, sin(a2); 0, 1, 0; -sin(a2), 0, cos(a2)]; R3 = [ cos(a3), -sin(a3), 0; sin(a3), cos(a3), 0; 0, 0, 1]; R = R1*R2*R3; endfunction

Code 2: Scilab function that converts XYZ Euler angles to the corresponding rotation matrix

function vertices = transver(vertices_0, r, R) // Transform vertices by translation vector r and rotation matrix R vertices = repmat(r', size(vertices_0, 1), 1) + vertices_0*R'; endfunction

Code 3: Scilab function that transforms coordinates of vertices

Figure 2: A block shape made of six polygons (in a general configuration)

Sponsor Link

© 2010-2020 Sky Engineering Laboratory

We provide consulting and engineering services on analysis and control of mechanical systems.

We provide consulting and engineering services on analysis and control of mechanical systems.