diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /698/CH2/EX2.9/P9_stresses_at_section.sce | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '698/CH2/EX2.9/P9_stresses_at_section.sce')
| -rw-r--r-- | 698/CH2/EX2.9/P9_stresses_at_section.sce | 70 |
1 files changed, 70 insertions, 0 deletions
diff --git a/698/CH2/EX2.9/P9_stresses_at_section.sce b/698/CH2/EX2.9/P9_stresses_at_section.sce new file mode 100644 index 000000000..776d56ded --- /dev/null +++ b/698/CH2/EX2.9/P9_stresses_at_section.sce @@ -0,0 +1,70 @@ +clc
+//Example 2.9
+//Stresses at a section
+//------------------------------------------------------------------------------
+
+//Given Data:
+//Loads
+//Load acting vertically downward
+Py=-2000//N (compressive)
+//Load acting horizontally
+Px=1000//N
+
+//Dimensions
+//Diameter
+d=0.05//m
+//length of neck region
+a=0.125//m
+//distance between midpoint of section and line of action of vertical load
+b=0.2//m
+//distance between midpoint of section and line of horizontal load
+c=0.25//m
+//Moment of inertia
+I=(%pi/64)*(d^4)
+//Polar moment of inertia
+J=(%pi/32)*(d^4)
+//Area of cross section
+A=(%pi/4)*(d^2)
+//------------------------------------------------------------------------------
+
+//Printing steps and result file to .txt
+res9=mopen(TMPDIR+'9_stresses_at_section.txt','wt')
+//Bending moment due to horizontal load about midpoint
+Mx=Px*c
+mfprintf(res9,"Bending moment due to horizontal load about midpoint Mx= %0.2f Nm\n",Mx)
+//Bending moment due to vertical load about midpoint
+My=Py*b
+mfprintf(res9,"Bending moment due to vertical load about midpoint My= %0.2f Nm\n",My)
+//Total bending moment
+M=sqrt((Mx^2)+(My^2))
+mfprintf(res9,"The total bending moment:\n\tM=sqrt((Mx^2)+(My^2))\n")
+mfprintf(res9,"The total bending moment M = %0.2f Nm\n\n",M)
+//Normal stress in X direction
+Sx=Py/A-((M*(d/2))/I)
+mfprintf(res9,"The normal stress is given by\n\tSx=Py/A-((M*(d/2))/I)\n")
+mfprintf(res9,"Normal stress acting at the section is %0.2f MN/m^2\n\n",Sx*(10^-6))
+
+//Couple moment due to horizontal load about midpoint
+T=Px*b
+//Shear stress
+Txy=(T*(d/2))/J
+mfprintf(res9,"The shear stress is given by\n\tTxy=(T*(d/2))/J\n")
+mfprintf(res9,"Normal stress acting at the section is %0.2f MN/m^2\n\n",Txy*(10^-6))
+//------------------------------------------------------------------------------
+
+//Extreme streeses:
+//Maximum normal stress
+Sn_max=(Sx/2)+sqrt(((Sx/2)^2)+(Txy^2))
+mfprintf(res9,"Maximum normal stress is %0.2f MN/m^2\n",Sn_max*(10^-6))
+//Minimum normal stress
+Sn_min=(Sx/2)-sqrt(((Sx/2)^2)+(Txy^2))
+mfprintf(res9,"Minimum normal stress is %0.2f MN/m^2\n\n",Sn_min*(10^-6))
+S=max((abs(Sn_max)),(abs(Sn_min)))
+mfprintf(res9,"The maximum numerical normal stress is %0.2f MN/m^2",S*(10^-6))
+//Maximum shear stress
+T_max=sqrt(((Sx/2)^2)+(Txy^2))
+mfprintf(res9,"\nThe maximum shear stress is %f MN/m^2",T_max*(10^-6))
+mclose(res9)
+editor(TMPDIR+'9_stresses_at_section.txt')
+//------------------------------------------------------------------------------
+//-----------------------------End of program-----------------------------------
\ No newline at end of file |