summaryrefslogtreecommitdiff
path: root/3705/CH13
diff options
context:
space:
mode:
Diffstat (limited to '3705/CH13')
-rw-r--r--3705/CH13/EX13.1/Ex13_1.sce26
-rw-r--r--3705/CH13/EX13.2/Ex13_2.sce36
2 files changed, 62 insertions, 0 deletions
diff --git a/3705/CH13/EX13.1/Ex13_1.sce b/3705/CH13/EX13.1/Ex13_1.sce
new file mode 100644
index 000000000..caa62df4c
--- /dev/null
+++ b/3705/CH13/EX13.1/Ex13_1.sce
@@ -0,0 +1,26 @@
+
+clear//
+
+//Variable Declaration
+d=150 //Depth of the web in mm
+wf=100 //Width of the flange in mm
+df=20 //Depth of the flange in mm
+t=20 //Thickness of the web in mm
+
+//Calculations
+y_bar=10**-3*(((wf*df*(d+df*0.5))+(d*t*d*0.5))/(wf*df+d*t)) //Distance of Neutral Axis in m
+//Simplfying the computation
+a=wf*df**3*12**-1
+b=wf*df*((d+df*0.5)-y_bar*10**3)**2
+c=t*d**3*12**-1
+f=t*d*((d*0.5)-y_bar*10**3)**2
+I=(a+b+c+f)*10**-12 //Moment of inertia in mm^3
+
+//Limit Moment
+yp=(wf*df+d*t)/(2*t) //Plastic Neutral Axis in mm
+Myp=I/y_bar //Yielding will start at moment without the stress term to ease computation
+mom=10**-9*((t*yp**2*0.5)+(wf*df*(d-yp+10))+(t*25**2*0.5)) //Sum of 1st moments
+Ml_Myp=mom*Myp**-1 //Ratio
+
+//Result
+printf("\n The ratio ML/Myp= %0.3f ",Ml_Myp)
diff --git a/3705/CH13/EX13.2/Ex13_2.sce b/3705/CH13/EX13.2/Ex13_2.sce
new file mode 100644
index 000000000..051a9aaa2
--- /dev/null
+++ b/3705/CH13/EX13.2/Ex13_2.sce
@@ -0,0 +1,36 @@
+
+clear//
+//
+
+//Variable Declaration
+E_st=200 //Youngs Modulus of Steel in GPa
+sigma_st_yp=290 //Yielding Stress in MPa
+E_al=70 //Youngs Modulus of Aluminium in GPa
+sigma_al_yp=330 //Yielding Stresss of Aluminium in MPa
+A_st=900 //Area of steel rod in mm^2
+A_al=600 //Area of Aluminium rod in mm^2
+L_st=350 //Length of the steel rod in mm
+L_al=250 //Length of the aluminium rod in mm
+
+//Calculations
+//Limit Load
+P_st=sigma_st_yp*A_st*10**-3 //Load in limiting condition in kN
+P_al=sigma_al_yp*A_al*10**-3 //Load in limiting condition in kN
+P_L=P_st+2*P_al //Total Loading in kN
+
+//Elastic Unloading
+//Solving for Pst and Pal using matri approach
+A=([[1,2;L_st*(E_st*A_st)**-1,-L_al*(E_al*A_al)**-1]])
+B=([P_L;0])
+C=linsolve(A,B) //Loading in kN
+
+//Residual Stresses
+P_res_st=-C(1)-P_st //Residual Load in kN
+P_res_al=-C(2)-P_al //Residual Load in kN
+sigma_st=P_res_st/A_st //residual Stress in Steel in MPa
+sigma_al=P_res_al/A_al //residual Stress in Aluminium in MPa
+
+
+//Result
+printf("\n The Residual stresses are as follows")
+printf("\n Sigma_st= %0.1f MPa and sigma_al= %0.1f MPa",sigma_st*10**3,sigma_al*10**3)