From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 3250/CH3/EX3.6/Ex3_6.sce | 23 +++++++++++++++++++++++ 1 file changed, 23 insertions(+) create mode 100755 3250/CH3/EX3.6/Ex3_6.sce (limited to '3250/CH3/EX3.6/Ex3_6.sce') diff --git a/3250/CH3/EX3.6/Ex3_6.sce b/3250/CH3/EX3.6/Ex3_6.sce new file mode 100755 index 000000000..7469536ab --- /dev/null +++ b/3250/CH3/EX3.6/Ex3_6.sce @@ -0,0 +1,23 @@ +clc +// Given that +r = 150 // Radius of the circular disc of lead in mm +Ti = 50 // Initial thickness of the disc in mm +Tf = 25 // Reduced thickness of the disc in mm +mu = 0.25// Cofficient of friction between the job and the dies +K = 4 // Avg shear yield stress of the lead in N/mm^2 +// Sample Problem 6 on page no. 122 +printf("\n # PROBLEM 3.6 # \n") +R = r*sqrt(2) +rs = (R - ((Tf/(2*mu)) * log(1/(mu*sqrt(3))))) +funcprot(0) +function p1 = f(x), p1 = (((sqrt(3))*K)*exp((2*mu/Tf)*(R-x)))*x, +endfunction +I = intg(rs,R,f) +funcprot(0) +function p2 = f(y), p2 = ((2*K/Tf)*(R-y) + ((K/mu)*(1+log(mu*sqrt(3)))))*y, +endfunction +I_ = intg(0,rs,f) +F = 2*%pi*(I+I_) +printf("\n The maximum forging force = %e N",F) + + -- cgit