1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
|
#include <iostream>
#include <string>
#include <sstream>
#include <cmath>
#include <iomanip>
#include <fstream>
#include <string>
#include <cstdlib>
using namespace std;
#define PI 3.14159
int main()
{
int structure;
cout << "Enter your 1 for Ibeam ";
cin >> structure;
int num;
double b1,b2,b3,d1,d2,d3,c2,I1,I2,I3,I,w,x,E,LO,theta,delta;
string filename;
char *outputfilename = "output.txt";
ifstream inFile;
ofstream outputfile(outputfilename);
cout << "Please enter the name of the file you wish to open: ";
cin >> filename;
inFile.open(filename.c_str());
if (inFile.fail())
{
cout << "\nThe file named " << filename
<< " was not successfully opened"
<< "\n Please check that the file currently exists."
<<endl;
system("PAUSE");
exit(1);
}
cout << "\nThe file has been successfully opened for reading.\n";
if (!outputfile){
cout << "Error Opening Output File: " << outputfilename
<< " .\nCheck if file already exists. "<< endl;
return 0 ;
}
cout << "Output File Success fully Opened. " << endl;
inFile >> num;
inFile >> b1;
inFile >> b2;
inFile >> b3;
inFile >> d1;
inFile >> d2;
inFile >> d3;
inFile >> w;
inFile >> E;
inFile >> LO;
cout << "b1 from file: " << b1<<endl;
cout << "b2 from file: " <<b2<<endl;
cout <<"b3 from file: " <<b3<<endl;
cout <<"d1 from file: " <<d1<<endl;
cout <<"d2 from file: " <<d2<<endl;
cout <<"d3 from file: " <<d3<<endl;
cout <<"w from file: " <<w<<endl;
cout <<"E from file: " <<E<<endl;
cout <<"LO from file: " <<LO<<endl;
if (structure == 1)
{
c2= (((b1 * d1) * (d1 / 2)) + ((b2 * d2) * (d1 + (d2/2))) + ((b3 * d3) * (d1 + d2 + (d3/2)))) / ((b1 * d1) + (b2 * d2) + (b3 * d3));
cout << "c2: " << c2 << " m\n";
I1= ((b1 * (pow(d1,3))) / 12) + (b1 * d1) * (c2 - (d1 / 2)) * (c2 - (d1 / 2));
cout << "I1: " << I1 << " m\n";
I2= (((b2 * (pow(d2,3))) / 12) + ((b2 * d2) * ((d1 + (d2 / 2) - c2) * (d1 + (d2 / 2) - c2))));
cout << "I2: " << I2 << " m\n";
I3= ((b3 * (pow(d3,3))) / 12) + (b3 * d3) * (pow((d1 + d2 + (d3/2) - c2),2));
cout << "I3: " << I3 << " m\n";
cout << "Sum of I: " << I1 + I2 + I3 << " m\n";
for (double x = 0.0; x <= 1; x +=0.1)
{theta= ((((w * x) / (6 * E * (I1+I2+I3)))) * ((3 * (LO * LO)) - (3 * LO * x) + (x * x)));
cout << "theta: " << theta << " m\n";
delta= (((w * (x*x)) / (24 * E * (I1 +I2+I3)))) *(((6*(LO * LO))-(4*LO*(x))+((x)*(x))));
cout << "delta: " <<delta<< " m\n";
}
outputfile <<"b1= " <<b1<<endl;
outputfile <<"b2= " <<b2<<endl;
outputfile <<"b3= " <<b3<<endl;
outputfile <<"d1= " <<d1<<endl;
outputfile <<"d2= " <<d2<<endl;
outputfile <<"d3= " <<d3<<endl;
outputfile <<"c2= " <<c2<<endl;
outputfile <<"I1= " <<I1<<endl;
outputfile <<"I2= " <<I2<<endl;
outputfile <<"I3= " <<I3<<endl;
outputfile<<"Sum of I= " <<(I1+I2+I3)<<endl;
outputfile <<"w= "<<w<<endl;
outputfile <<"E= "<<E<<endl;
outputfile <<"LO= "<<LO<<endl;
outputfile << " Distance x(m) Deflection delta(m) Rotation theta "<<endl;
outputfile <<"\n " <<(.1*LO); outputfile <<" "<<((((w * ((.1*LO)*(.1*LO))) / (24 * E * (I1 +I2+I3)))) *(((6*(LO * LO))-(4*LO*(.1*LO))+((.1*LO)*(.1*LO))))); outputfile<<" "<<(((((w * (.1*LO)) / (6 * E * (I1+I2+I3)))) * ((3 * (LO * LO)) - (3 * LO * (.1*LO)) + ((.1*LO) * (.1*LO)))))<<endl;
outputfile <<"\n " <<(.2*LO); outputfile <<" "<<((((w * ((.2*LO)*(.2*LO))) / (24 * E * (I1 +I2+I3)))) *(((6*(LO * LO))-(4*LO*(.2*LO))+((.2*LO)*(.2*LO))))); outputfile<<" "<<(((((w * (.2*LO)) / (6 * E * (I1+I2+I3)))) * ((3 * (LO * LO)) - (3 * LO * (.2*LO)) + ((.2*LO) * (.2*LO)))))<<endl;
outputfile <<"\n " <<(.3*LO); outputfile <<" "<<((((w * ((.3*LO)*(.3*LO))) / (24 * E * (I1 +I2+I3)))) *(((6*(LO * LO))-(4*LO*(.3*LO))+((.3*LO)*(.3*LO))))); outputfile<<" "<<(((((w * (.3*LO)) / (6 * E * (I1+I2+I3)))) * ((3 * (LO * LO)) - (3 * LO * (.3*LO)) + ((.3*LO) * (.3*LO)))))<<endl;
outputfile <<"\n " <<(.4*LO); outputfile <<" "<<((((w * ((.4*LO)*(.4*LO))) / (24 * E * (I1 +I2+I3)))) *(((6*(LO * LO))-(4*LO*(.4*LO))+((.4*LO)*(.4*LO))))); outputfile<<" "<<(((((w * (.4*LO)) / (6 * E * (I1+I2+I3)))) * ((3 * (LO * LO)) - (3 * LO * (.4*LO)) + ((.4*LO) * (.4*LO)))))<<endl;
outputfile <<"\n " <<(.5*LO); outputfile <<" "<<((((w * ((.5*LO)*(.5*LO))) / (24 * E * (I1 +I2+I3)))) *(((6*(LO * LO))-(4*LO*(.5*LO))+((.5*LO)*(.5*LO))))); outputfile<<" "<<(((((w * (.5*LO)) / (6 * E * (I1+I2+I3)))) * ((3 * (LO * LO)) - (3 * LO * (.5*LO)) + ((.5*LO) * (.5*LO)))))<<endl;
outputfile <<"\n " <<(.6*LO); outputfile <<" "<<((((w * ((.6*LO)*(.6*LO))) / (24 * E * (I1 +I2+I3)))) *(((6*(LO * LO))-(4*LO*(.6*LO))+((.6*LO)*(.6*LO))))); outputfile<<" "<<(((((w * (.6*LO)) / (6 * E * (I1+I2+I3)))) * ((3 * (LO * LO)) - (3 * LO * (.6*LO)) + ((.6*LO) * (.6*LO)))))<<endl;
outputfile <<"\n " <<(.7*LO); outputfile <<" "<<((((w * ((.7*LO)*(.7*LO))) / (24 * E * (I1 +I2+I3)))) *(((6*(LO * LO))-(4*LO*(.7*LO))+((.7*LO)*(.7*LO))))); outputfile<<" "<<(((((w * (.7*LO)) / (6 * E * (I1+I2+I3)))) * ((3 * (LO * LO)) - (3 * LO * (.7*LO)) + ((.7*LO) * (.7*LO)))))<<endl;
outputfile <<"\n " <<(.8*LO); outputfile <<" "<<((((w * ((.8*LO)*(.8*LO))) / (24 * E * (I1 +I2+I3)))) *(((6*(LO * LO))-(4*LO*(.8*LO))+((.8*LO)*(.8*LO))))); outputfile<<" "<<(((((w * (.8*LO)) / (6 * E * (I1+I2+I3)))) * ((3 * (LO * LO)) - (3 * LO * (.8*LO)) + ((.8*LO) * (.8*LO)))))<<endl;
outputfile <<"\n " <<(.9*LO); outputfile <<" "<<((((w * ((.9*LO)*(.9*LO))) / (24 * E * (I1 +I2+I3)))) *(((6*(LO * LO))-(4*LO*(.9*LO))+((.9*LO)*(.9*LO))))); outputfile<<" "<<(((((w * (.9*LO)) / (6 * E * (I1+I2+I3)))) * ((3 * (LO * LO)) - (3 * LO * (.9*LO)) + ((.9*LO) * (.9*LO)))))<<endl;
outputfile <<"\n " <<(1*LO); outputfile <<" "<<((((w * ((1*LO)*(1*LO))) / (24 * E * (I1 +I2+I3)))) *(((6*(LO * LO))-(4*LO*(1*LO))+((1*LO)*(1*LO))))); outputfile<<" "<<(((((w * (1*LO)) / (6 * E * (I1+I2+I3)))) * ((3 * (LO * LO)) - (3 * LO * (1*LO)) + ((1*LO) * (1*LO)))))<<endl;
system("PAUSE");
return 0;
}
else if (structure>1 || structure<1)
{
cout <<"You did not type a valid number." <<endl;
system("PAUSE");
return 0;
}
system ("PAUSE");
return 0;
}
|
