Demonstrate DrawSailMeshVectors for different sail models.
This function calls SolarPressureForce on the input CAD model to
compute the solar force and resulting temperature.
Since version 7.
------------------------------------------------------------------------
See also Unit, DrawSailMeshVectors, CP1Props, SailMesh
------------------------------------------------------------------------
Contents
Define an offset sun vector
uSun = Unit([0.3*randn(2,1); 1]);
Draw vectors for the circular sail model
g = load('CircularSail');
DrawSailMeshVectors( g, uSun )
area =
79127
normal =
2.0108e-12
-2.4158e-13
78374
force =
0.0018796
0.0011384
-0.68529
Define an offset sun vector
uSun = Unit([1; 0.3*randn(2,1)]);
Draw vectors for the quadrant sail model
g = load('QuadSail_100');
DrawSailMeshVectors( g, uSun )
area =
10310
normal =
10000
9.5923e-14
4.3265e-13
force =
-0.058686
-0.0030077
-0.0028448
Draw vectors for the striped sail model
g = load('StripedSail');
DrawSailMeshVectors( g, uSun )
area =
491.85
normal =
450
-15.91
15.91
force =
-0.0023613
-6.2624e-06
-0.00012308
area =
491.85
normal =
450
15.91
15.91
force =
-0.0026457
-0.00049523
-0.00049152
area =
491.85
normal =
450
15.91
-15.91
force =
-0.0023684
-0.00014617
1.6603e-05
area =
491.85
normal =
450
-15.91
-15.91
force =
-0.0023904
-0.00033406
-0.00033057
Do the same thing by hand for a simplified circular mesh
theta = linspace(0,2*pi,8);
theta = theta(1:end-1);
rSail = 50;
x = rSail*cos(theta);
y = rSail*sin(theta);
dBillow = struct('b',0.8,'r',rSail);
[v,f] = SailMesh( x, y, 'HCircularBillow', dBillow, 1 );
uSun = Unit([0.3*randn(2,1); 1]);
[optical, infrared, thermal] = CP1Props;
DrawSailMeshVectors( struct('v',v,'f',f,'optical',optical,'thermal',thermal), uSun );
view(-10,15)
