input-basic-units3

有限元计算程序

/*  =======================================================  */
/*  SI and US CONVERSIONS                                    */
/*                                                           */
/*  Tests variables and matrices in different Units Systems  */
/*  =======================================================  */

/* UNITS of LENGTH     */

print " LENGTH UNITS \n";

v =  1 km;
w =  1 mile;

print  "v    = ", v       , "\n";
print  "v    = ", v (in)  , "\n";
print  "v    = ", v (ft)  , "\n";
print  "v    = ", v (yard), "\n";
print  "v    = ", v (mile), "\n";

print "\n";

print  "w    = ", w       , "\n";
print  "w    = ", w (mm)  , "\n";
print  "w    = ", w (cm)  , "\n";
print  "w    = ", w (dm)  , "\n";
print  "w    = ", w (m )  , "\n";
print  "w    = ", w (km)  , "\n";

print "\n";
W = v + w;
V = v - w;
U = v * w;
X = v / w;
Y = sqrt(v*w);

print "v + w = ", W     , "\n";
print "v + w = ", W (mile), "\n";
print "v - w = ", V       , "\n";
print "v * w = ", U       , "\n";
print "v * w = ", U (ft^2), "\n";
print "v * w = ", U (m^2), "\n";
print "v / w = ", X       , "\n";
print "sqrt(v*w) = ", Y   , "\n";
print "sqrt(v*w) = ", Y(m), "\n";

print "\n";

/* ====================*/
/* UNITS of VOLUME     */
/* ====================*/
print "\n VOLUME UNITS \n";
x =  1       gallon;
y =  1       barrel;

print  "x = ", x     , "\n";
print  "x = ", x (m^3), "\n";
print  "x = ", x (cm^3), "\n";
print  "x = ", x (ft^3), "\n";
print  "x = ", x (in^3), "\n";
print  "x = ", x (barrel), "\n";

print "\n";

print  "y = ", y, "\n";
print  "y = ", y (m^3), "\n";
print  "y = ", y (cm^3), "\n";
print  "y = ", y (ft^3), "\n";
print  "y = ", y (in^3), "\n";
print  "y = ", y (gallon), "\n";

/* ====================*/
/* UNITS of MASS       */
/* ====================*/

print "\n";
print "\n MASS UNITS \n";

x = 1 kg;
y = 1 klb;

print  "x = ", x     , "\n";
print  "x = ", x (lb), "\n";
print  "x = ", x (klb), "\n";
print  "x = ", x (ton), "\n";

print "\n";
print  "y = ", y     , "\n";
print  "y = ", y (kg), "\n";

print "\n";
W = x + y;

print "x + y = ", W      , "\n";
print "x + y = ", W (klb), "\n";
 
print "\n";
/* ====================*/
/* UNITS of TIME       */
/* ====================*/
print " TIME UNITS \n";
x =  1  hr;

print  "x = ", x, "\n";
print  "x = ", x (min), "\n";
print  "x = ", x (sec), "\n";

/* =====================*/
/* UNITS of TEMPERATURE */
/* =====================*/

print "\n";
print "\n TEMPERATURE UNITS \n";

x = 1 deg_C;
y = 1 deg_F;
u = 10 mm/deg_C;   /* not properiate units */ 
v = 10 in/deg_F;   /* Worng */
U = 10 mm/DEG_C;   
V = 10 in/DEG_F;

print  "x = ", x        , "\n";
print  "x = ", x (deg_F), "\n";

print "y = ", y        , "\n";
print  "y = ", y (deg_C), "\n";
print "\n";

print  "u = ", u , "\n";
print  "v = ", v , "\n";
print  "U = ", U , "\n";
print  "V = ", V , "\n";

print "\n";

print "u*(1 deg_C) = ", u*(1 deg_C) (mm), "\n";
print "u*(1 deg_F) = ", u*(1 deg_F) (mm), " WRONG !!\n";
print "\n";
print "U*(1      DEG_C) = ", U*(1      DEG_C) (mm), "\n";
print "U*(1.8    DEG_F) = ", U*(1.8    DEG_F) (mm), "\n";
print "V*(0.5556 DEG_C) = ", V*(0.5556 DEG_C) (in), "\n";
print "V*(1      DEG_F) = ", V*(1      DEG_F) (in), "\n";

print "\n";
print "v*(1 deg_C) = ", v*(1 deg_C) (in), "\n";
print "v*(1 deg_F) = ", v*(1 deg_F) (in), " WRONG !!\n";
print "U*(1 DEG_C) = ", U*(1 DEG_C) (mm), "\n";
print "U*(1 DEG_F) = ", U*(1 DEG_F) (mm), "\n";
print "V*(1 DEG_C) = ", V*(1 DEG_C) (in), "\n";
print "V*(1 DEG_F) = ", V*(1 DEG_F) (in), "\n";

/* ===========================*/
/* UNITS of FREQUENCY & SPEED */
/* ===========================*/

print "\n FREQUENCY & SPEED UNITS \n";

y = 60 rpm;     /* rev. per min  */

print  "y = ", y, "\n";
print  "y = ", y (Hz), "\n";
print  "y = ", y (cps), "\n";

/* ===========================*/
/* UNITS of FORCE             */
/* ===========================*/

print "\n FORCE UNITS \n";
x =  1  kN;
y =  1 kips;

print  "x = ", x      , "\n";
print  "x = ", x (kgf), "\n";
print  "x = ", x (kips), "\n";
print  "x = ", x (lbf), "\n";
print "\n";
print  "y = ", y, "\n";
print  "y = ", y (kN), "\n";
print  "y = ", y (kgf), "\n";

/* ===========================*/
/* UNITS of PRESSURE          */
/* ===========================*/

print "\n PRESSURE UNITS \n";
x =  1  kPa;
y =  1  ksi;

print  "x = ", x, "\n";
print  "x = ", x (psi), "\n";

print  "y = ", y, "\n";
print  "y = ", y (kPa), "\n";

/* ===========================*/
/* UNITS of ENERGY            */
/* ===========================*/

print "\n ENERGY UNITS \n";
x =  1  kJ;
print  "x = ", x, "\n";
print  "x = ", x (Jou), "\n";

/* ===========================*/
/* UNITS of POWER             */
/* ===========================*/
print "\n POWER UNITS \n";
x =  1  kW;
print  "x = ", x, "\n";
print  "x = ", x (Watt), "\n";

/* ===========================*/
/* UNITS of PLANE ANGLE       */
/* ===========================*/
print "\n PLANE ANGLE UNITS \n";
x =  180 deg;
y =  PI;
print  "x = ", x, "\n";
print  "x = ", x (rad), "\n";
print  "y = ", y, "\n";
print  "y = ", y (deg), "\n";

U = [10 in/DEG_F, 2 in, 20 mm/DEG_C];
V = [2 DEG_F; 2; 1.8 DEG_F];
Z = U*V;
PrintMatrix(U, V, Z);