Skip to content Skip to navigation
Summary: Problem Set for Introductory Programming
Write a script that will ask for pressure value in psi and display the equivalent pressure in kPa with a statement, such as "The converted pressure is: ..."
% This script converts pressures from psi to kPa
% User is prompted to enter pressure in psi
clc % Clear screen
disp('This script converts pressures from psi to kPa:')
disp(' ') % Display blank line
psi=input('What is the pressure value in psi? ');
kPa=psi*6.894757; % Calculating pressure in kPa
disp(' ') % Display blank line
str = ['The converted pressure is: ', num2str(kPa), ' kPa.'];
disp(str);
This script converts pressures from psi to kPa:
What is the pressure value in psi? 150
The converted pressure is: 1034.2135 kPa.
Write a script that generates a table of conversions from Fahrenheit to Celsius temperatures for a range and increment entered by the user, such as
Enter the beginning temperature in F:
Enter the ending temperature in F:
Enter the increment value:
Test your script with 20 the beginning Fahrenheit value, 200 the ending Fahrenheit value and 20 the increment.
% This script generates a table of conversions
% From Fahrenheit to Celsius temperatures
clc % Clear screen
disp('This script generates a table of conversions from Fahrenheit to Celsius')
disp(' ') % Display blank line
lowerF=input('Enter the beginning temperature in F: ');
upperF=input('Enter the ending temperature in F: ');
increment=input('Enter the increment value: ');
Fahrenheit=[lowerF:increment:upperF]; % Creating a row vector with F values
Celsius=5/9*(Fahrenheit-32); % Converting from F to C
disp(' ') % Display blank line
str = ['Fahrenheit Celsius '];% Displaying table header
disp(str);
% Tabulating results in two columns, ' is being used to transpose row to column
disp([Fahrenheit' Celsius'])
This script generates a table of conversions from Fahrenheit to Celsius
Enter the beginning temperature in F: 20
Enter the ending temperature in F: 200
Enter the increment value: 20
Fahrenheit Celsius
20.0000 -6.6667
40.0000 4.4444
60.0000 15.5556
80.0000 26.6667
100.0000 37.7778
120.0000 48.8889
140.0000 60.0000
160.0000 71.1111
180.0000 82.2222
200.0000 93.3333
Pascal's Law states that pressure is transmitted undiminished in all directions throughout a fluid at rest. (See the illustration below). An initial force of 150 N is transmitted from a piston of 25 mm^2 to a piston of 100 mm^2. This force is progressively increased up to 200 N. Write a script that computes the corresponding load carried by the larger piston and tabulate your results.
 |
% This script computes the load carried by the larger piston in a hydraulic system
clc % Clear screen
disp('This script computes the load carried by the larger piston in a hydraulic system')
disp(' ') % Display blank line
initialF=150;
finalF=200;
increment=10;
area1=25;
area2=100;
F1=[initialF:increment:finalF]; % Creating a row vector with F1 values
F2=F1*area2/area1; % Calculating F2 values
disp(' ') % Display blank line
str = [' F1 F2 '];% Displaying table header
disp(str);
disp([F1' F2']) % Tabulating results in two columns, ' is being used to transpose row to column
This script computes the load carried by the larger piston in a hydraulic system
F1 F2
150 600
160 640
170 680
180 720
190 760
200 800
Modify your script in previous problem so that the user provides the following input:
Enter the initial force in N:
Enter the final force in N:
Enter the increment value:
Enter the area of small piston in mm^2:
Enter the area of big piston in mm^2:
Test your script with 150, 200, 10, 25 and 100 with respect to each input variable.
% This script computes the load carried by the larger piston in a hydraulic system
clc % Clear screen
disp('This script computes the load carried by the larger piston in a hydraulic system')
disp(' ') % Display blank line
initialF=input('Enter the initial force in N: ');
finalF=input('Enter the final force in N: ');
increment=input('Enter the increment value: ');
area1=input('Enter the area of small piston in mm^2: ');
area2=input('Enter the area of big piston in mm^2: ');
F1=[initialF:increment:finalF]; % Creating a row vector with F1 values
F2=F1*area2/area1; % Calculating F2 values
disp(' ') % Display blank line
str = [' F1 F2 '];% Displaying table header
disp(str);
disp([F1' F2']) % Tabulating results in two columns, ' is being used to transpose row to column
This script computes the load carried by the larger piston in a hydraulic system
Enter the initial force in N: 150
Enter the final force in N: 200
Enter the increment value: 10
Enter the area of small piston in mm^2: 25
Enter the area of big piston in mm^2: 100
F1 F2
150 600
160 640
170 680
180 720
190 760
200 800
Write a script to solve the Stress-Strain problem in the Problem Set
The m-file contains the following:
% This script uses readings from a Tensile test and
% Computes Strain and Stress values
clc % Clear screen
disp('This script uses readings from a Tensile test and')
disp('Computes Strain and Stress values')
disp(' ') % Display a blank line
Specimen_dia=12.7; % Specimen diameter in mm
% Load in kN
Load_kN=[0;4.89;9.779;14.67;19.56;24.45;...
27.62;29.39;32.68;33.95;34.58;35.22;...
35.72;40.54;48.39;59.03;65.87;69.42;...
69.67;68.15;60.81];
% Gage length in mm
Length_mm=[50.8;50.8102;50.8203;50.8305;...
50.8406;50.8508;50.8610;50.8711;...
50.9016;50.9270;50.9524;50.9778;...
51.0032;51.816;53.340;55.880;58.420;...
60.96;61.468;63.5;66.04];
% Calculate x-sectional area im m2
Cross_sectional_Area=pi/4*((Specimen_dia/1000)^2);
% Calculate change in length, initial lenght is 50.8 mm
Delta_L=Length_mm-50.8;
% Calculate Stress in MPa
Sigma=(Load_kN./Cross_sectional_Area)*10^(-3);
% Calculate Strain in mm/mm
Epsilon=Delta_L./50.8;
str = ['Specimen diameter is ', num2str(Specimen_dia), ' mm.'];
disp(str);
Results=[Load_kN Length_mm Delta_L Sigma Epsilon];
% Tabulated results
disp(' Load Length Delta L Stress Strain')
disp(Results)
This script uses readings from a Tensile test and
Computes Strain and Stress values
Specimen diameter is 12.7 mm.
Load Length Delta L Stress Strain
0 50.8000 0 0 0
4.8900 50.8102 0.0102 38.6022 0.0002
9.7790 50.8203 0.0203 77.1964 0.0004
14.6700 50.8305 0.0305 115.8065 0.0006
19.5600 50.8406 0.0406 154.4086 0.0008
24.4500 50.8508 0.0508 193.0108 0.0010
27.6200 50.8610 0.0610 218.0351 0.0012
29.3900 50.8711 0.0711 232.0076 0.0014
32.6800 50.9016 0.1016 257.9792 0.0020
33.9500 50.9270 0.1270 268.0047 0.0025
34.5800 50.9524 0.1524 272.9780 0.0030
35.2200 50.9778 0.1778 278.0302 0.0035
35.7200 51.0032 0.2032 281.9773 0.0040
40.5400 51.8160 1.0160 320.0269 0.0200
48.3900 53.3400 2.5400 381.9955 0.0500
59.0300 55.8800 5.0800 465.9888 0.1000
65.8700 58.4200 7.6200 519.9844 0.1500
69.4200 60.9600 10.1600 548.0085 0.2000
69.6700 61.4680 10.6680 549.9820 0.2100
68.1500 63.5000 12.7000 537.9830 0.2500
60.8100 66.0400 15.2400 480.0403 0.3000
Modify the script, you wrote above and plot an annotated Stress-Strain graph.
Edited script contains the plot commands:
% This script uses readings from a Tensile test and
% Computes Strain and Stress values
clc % Clear screen
disp('This script uses readings from a Tensile test and')
disp('Computes Strain and Stress values')
disp(' ') % Display a blank line
Specimen_dia=12.7; % Specimen diameter in mm
% Load in kN
Load_kN=[0;4.89;9.779;14.67;19.56;24.45;...
27.62;29.39;32.68;33.95;34.58;35.22;...
35.72;40.54;48.39;59.03;65.87;69.42;...
69.67;68.15;60.81];
% Gage length in mm
Length_mm=[50.8;50.8102;50.8203;50.8305;...
50.8406;50.8508;50.8610;50.8711;...
50.9016;50.9270;50.9524;50.9778;...
51.0032;51.816;53.340;55.880;58.420;...
60.96;61.468;63.5;66.04];
% Calculate x-sectional area im m2
Cross_sectional_Area=pi/4*((Specimen_dia/1000)^2);
% Calculate change in length, initial lenght is 50.8 mm
Delta_L=Length_mm-50.8;
% Calculate Stress in MPa
Sigma=(Load_kN./Cross_sectional_Area)*10^(-3);
% Calculate Strain in mm/mm
Epsilon=Delta_L./50.8;
str = ['Specimen diameter is ', num2str(Specimen_dia), ' mm.'];
disp(str);
Results=[Load_kN Length_mm Delta_L Sigma Epsilon];
% Tabulated results
disp(' Load Length Delta L Stress Strain')
disp(Results)
% Plot Stress versus Strain
plot(Epsilon,Sigma)
title('Stress versus Strain Curve')
xlabel('Strain [mm/mm]')
ylabel('Stress [mPa]')
grid
More about this module: Metadata | Downloads | Version History
This work is licensed by Serhat Beyenir under a Creative Commons Attribution License (CC-BY 3.0), and is an Open Educational Resource.
Last edited by Serhat Beyenir on Nov 17, 2011 1:16 pm -0600.