Hello everyone, still on vacation but I'm back in the Los Angeles area.

I received an email from retired engineer Djajadi Susastra, who has been using Hewlett Packard calculators since 1989. Susastra obtained a BSc in Mechanical Engineering that year. He recently purchased an HP Prime calculator and has asked me to share his programs on mechanical engineering (the last three deal with business and finance), which he learned how to program the Prime from this blog.

Djajadi, I am very thankful for the compliments and with pleasure, here are the programs he has sent:

I received an email from retired engineer Djajadi Susastra, who has been using Hewlett Packard calculators since 1989. Susastra obtained a BSc in Mechanical Engineering that year. He recently purchased an HP Prime calculator and has asked me to share his programs on mechanical engineering (the last three deal with business and finance), which he learned how to program the Prime from this blog.

Djajadi, I am very thankful for the compliments and with pleasure, here are the programs he has sent:

**CALCULATING CROSS SECTION AREA OF CYLINDER.**

EXPORT AREA()

BEGIN

LOCAL OD, ID, AA;

INPUT (OD,
"OD=" );

INPUT
(ID,"ID=" );

AA:=(OD^2-ID^2)/4*PI;

MSGBOX
("XSECTION AREA=" +AA);

END;

**CALCULATING MAXIMUM BENDING MOMENT ON A CYLINDER AT YIELD.**

EXPORT BENDING()

BEGIN

LOCAL OD,
ID, YI, I, BEN;

INPUT
(OD, "OD (IN)=" );

INPUT
(ID, "ID= (IN)=" );

INPUT
(YI, "YIELD (PSI)=");

I:=(OD^4-ID^4)/64*PI;

BEN:=YI*I/(12*(OD/2));

MSGBOX
("MAX BENDING ALLOWED (FT-LB)=" +BEN);

END;

**CALCULATING COLLAPSE PRESSURE ON 110 KSI YIELD PIPE (PER API 5C3).**

EXPORT
COLLAPSE110()

BEGIN

LOCAL OD,
ID, YI, TEN, DT, COL;

INPUT
(OD, "OD (IN)=" );

INPUT
(ID, "ID= (IN)=" );

YI:=110000;

DT:=OD/((OD-ID)/2);

IF
DT≤12.44 THEN

COL:=2*YI*(DT-1)/DT^2;

MSGBOX
("COLLAPSE PRESS (PSI)=" +COL);

END;

IF
12.44<DT AND DT≤20.41 THEN

COL:=YI*(3.181/DT-.0819)-2852;

MSGBOX
("COLLPASE PRESS (PSI)=" +COL);

END;

IF
20.41<DT AND DT≤26.22 THEN

COL:=YI*(2.066/DT-.0532);

MSGBOX
("COLLAPSE PRESS (PSI)=" +COL);

END;

IF
26.22<DT AND DT≤1000 THEN

COL:=46950000/(DT*(DT-1)^2);

MSGBOX
("COLLAPSE PRESS (PSI)=" +COL);

END;

IF
DT≥1001 THEN

MSGBOX
("NO CASE");

END;

END;

**CALCULATING COLLAPSE PRESSURE ON 80 KSI YIELD PIPE (PER API 5C3).**

EXPORT
COLLAPSE80()

BEGIN

LOCAL OD,
ID, YI, TEN, DT, COL;

INPUT
(OD, "OD (IN)=" );

INPUT
(ID, "ID= (IN)=" );

YI:=80000;

DT:=OD/((OD-ID)/2);

IF
DT≤13.38 THEN

COL:=2*YI*(DT-1)/DT^2;

MSGBOX
("COLLAPSE PRESS (PSI)=" +COL);

END;

IF
13.38<DT AND DT≤22.47 THEN

COL:=YI*(3.071/DT-.0667)-1955;

MSGBOX
("COLLPASE PRESS (PSI)=" +COL);

END;

IF 22.47<DT
AND DT≤31.02 THEN

COL:=YI*(1.998/DT-.0434);

MSGBOX
("COLLAPSE PRESS (PSI)=" +COL);

END;

IF
31.02<DT AND DT≤1000 THEN

COL:=46950000/(DT*(DT-1)^2);

MSGBOX
("COLLAPSE PRESS (PSI)=" +COL);

END;

IF
DT≥1001 THEN

MSGBOX
("NO CASE");

END;

END;

**CALCULATING COLLAPSE PRESSURE ON 75KSI YIELD PIPE (PER API 5C3).**

EXPORT
COLLAPSE75()

BEGIN

LOCAL OD,
ID, YI, TEN, DT, COL;

INPUT
(OD, "OD (IN)=" );

INPUT
(ID, "ID= (IN)=" );

YI:=75000;

DT:=OD/((OD-ID)/2);

IF
DT≤13.6 THEN

COL:=2*YI*(DT-1)/DT^2;

MSGBOX ("COLLAPSE
PRESS (PSI)=" +COL);

END;

IF
13.6<DT AND DT≤22.91 THEN

COL:=YI*(3.054/DT-.0642)-1806;

MSGBOX
("COLLPASE PRESS (PSI)=" +COL);

END;

IF
22.91<DT AND DT≤32.05 THEN

COL:=YI*(1.99/DT-.0418);

MSGBOX
("COLLAPSE PRESS (PSI)=" +COL);

END;

IF
32.05<DT AND DT≤1000 THEN

COL:=46950000/(DT*(DT-1)^2);

MSGBOX
("COLLAPSE PRESS (PSI)=" +COL);

END;

IF
DT≥1001 THEN

MSGBOX
("NO CASE");

END;

END;

**CALCULATING VON MISSES STRESS (THICK WALL, SHIGLEY AND MITCHELL, API 6A)**

EXPORT
EQUI_STRESS()

BEGIN

LOCAL OD,
ID, P, S1, S2, S3, SE, T;

INPUT
(OD, "OD (IN)=" );

INPUT
(ID, "ID= (IN)=" );

INPUT (P,
"PRESSURE (PSI)=");

T:=(OD-ID)/2;

LOCAL RA;

RA:=OD/40;

IF
T>RA THEN

S1:=P*ID^2/(OD^2-ID^2);

S2:=P*(OD^2+ID^2)/(OD^2-ID^2);

S3:=-1*P;

SE:=(((S1-S2)^2+(S2-S3)^2+(S3-S1)^2)/2)^.5;

MSGBOX
("EQUIVALENT STRESS (PSI)=" +SE);

ELSE

MSGBOX
("THIN WALL");

END;

END;

**CALCULATING YIELD PRESSURE (API 5 C3).**

EXPORT
YIELDPRESS()

BEGIN

LOCAL OD,
ID, YI, PY;

INPUT
(OD, "OD (IN)=" );

INPUT
(ID, "ID= (IN)=" );

INPUT
(YI, "YIELD (PSI)=");

PY:=.875*(2*YI*(OD-ID)/2)/OD;

MSGBOX
("YIELD PRESSURE (PSI)=" +PY);

RETURN

END;

**CALCULATING SHEAR AREA AND BEARING AREA ON ACME THREADS.**

EXPORT
THREAD()

BEGIN

LOCAL D,
N, P, ES, KN, M, SA, BA;

INPUT (D,
"OD=" );

INPUT
(N,"NO OF THREAD PER INCH" );

P:=1/N;

ES:=D-.5*P-.008*D^.5-.006*D^.5-.03*(1/N)^.5;

KN:=D-P+.05*P;

M:=D-.05*P;

SA:=PI*KN*(.5+N*TAN(14.5)*(ES-KN));

BA:=(M^2-KN^2)/4*N*PI;

MSGBOX
("ACME 2G, SHEAR AREA=" +SA);

MSGBOX
("BEARING AREA=" +BA);

END;

**CALCULATING GROSS MARGIN (IN %) WHEN COST AND PRICE ARE KNOWN.**

EXPORT
COST()

BEGIN

LOCAL P,
C, GM;

INPUT (P,
"PRICE=" );

INPUT (C,
"COST=" );

GM:=100-(100*C/P);

MSGBOX
("GROSS MARGIN (%)=" +GM);

END;

**CALCULATING PRICE WHEN COST AND GROSS MARGIN (IN %) ARE KNOWN.**

EXPORT
GROSSMARGIN()

BEGIN

LOCAL G,
C, P;

INPUT (G,
"GM (%)" );

INPUT (C,
"COST=" );

P:=C/(1-G/100);

MSGBOX
("PRICE=" +P);

RETURN

END;

**CALCULATING FUTURE VALUE WHEN PRESENT VALUE, INTEREST PER YEAR, AND PERIOD IN YEAR ARE KNOWN.**

EXPORT
FUTUREVAL ()

BEGIN

LOCAL PV,
I, P, N, FV;

INPUT
(PV, "PRESENT VALUE=" );

INPUT
(I,"INTEREST (ANNUAL)=" );

INPUT (N,
"PERIOD (IN YEAR)=");

FV:=PV*(1+I/100)^N;

MSGBOX
("FUTURE VALUE=" +FV);

END

Eddie

This blog entry is property of Edward Shore. 2015.

Fantastic. These are PRIME examples for me to set up some welding engineering equations for myself

ReplyDeletegreat thank you

ReplyDelete