Question 1.3: Direct stresses of 160 N/mm^2 (tension) and 120 N/mm^2 (comp...

Direct stresses of 160 N / mm ^{2} (tension) and 120 N / mm ^{2} (compression) are applied at a particular point in an elastic material on two mutually perpendicular planes. The principal stress in the material is limited to 200 N / mm ^{2} (tension). Calculate the allowable value of shear stress at the point on the given planes. Determine also the value of the other principal stress and the maximum value of shear stress at the point. Verify your answer using Mohr’s circle.

Repeat the derivations presented using the Symbolic Math Toolbox in MATLAB®. Do not recreate Mohr’s circle.

The Blue Check Mark means that this solution has been answered and checked by an expert. This guarantees that the final answer is accurate.
Learn more on how we answer questions.

The stress system at the point in the material may be represented as shown in Fig. 1.13 by considering the stresses to act uniformly over the sides of a triangular element ABC of unit thickness. Suppose that the direct stress on the principal plane AB is \sigma . For horizontal equilibrium of the element,

\sigma AB \cos \theta=\sigma_{x} BC +\tau_{x y} AC

 

which simplifies to

\tau_{x y} \tan \theta=\sigma-\sigma_{x}  (i)

Considering vertical equilibrium gives

\sigma AB \sin \theta=\sigma_{y} AC +\tau_{x y} BC

 

or

\tau_{x y} \cot \theta=\sigma-\sigma_{y}  (ii)

Hence, from the product of Eqs. (i) and (ii),

\tau_{x y}^{2}=\left(\sigma-\sigma_{x}\right)\left(\sigma-\sigma_{y}\right)

 

Now, substituting the values \sigma_{x}=160 N / mm ^{2}, \sigma_{y}=-120 N / mm ^{2}, \text { and } \sigma=\sigma_{1}=200 N / mm ^{2}, we have

\tau_{x y}=\pm 113 N / mm ^{2}

 

Replacing cot\theta in Eq. (ii) with 1//tan\theta from Eq. (i) yields a quadratic equation in \sigma:

\sigma^{2}-\sigma\left(\sigma_{x}-\sigma_{y}\right)+\sigma_{x} \sigma_{y}-\tau_{x y}^{2}=0  (iii)

The numerical solutions of Eq. (iii) corresponding to the given values of \sigma_{x}, \sigma_{y} and \tau_{x y} are the principal stresses at the point, namely,

\sigma_{1}=200 N / mm ^{2}

 

given

\sigma_{ II }=-160 N / mm ^{2}

 

Having obtained the principal stresses, we now use Eq. (1.15) to find the maximum shear stress, thus

\tau_{\max }=\frac{\sigma_{ I }-\sigma_{ II }}{2}  (1.15)

\tau_{\max }=\frac{200+160}{2}=180 N / mm ^{2}

 

The solution is rapidly verified from Mohr’s circle of stress (Fig. 1.14). From the arbitrary origin O , OP _{1}, and OP _{2} are drawn to represent \sigma_{x}=160 N / mm ^{2} and \sigma_{y}=-120 N / mm ^{2}. The mid-point C of P _{1} P _{2} is then located. Next, OB =\sigma_{1}=200 N / mm ^{2} is marked out and the radius of the circle is then CB. OA is the required principal stress. Perpendiculars P _{1} Q _{1} and P _{2} Q _{2} to the circumference of the circle are equal to \pm \tau_{x y} (to scale), and the radius of the circle is the maximum shear stress.

Using the element shown in Fig. 1.13, derivations of the principal stresses and maximum shear stress are obtained through the following MATLAB file:

Capture
1
The 'Blue Check Mark' means that either the MATLAB code/script/answer provided in the answer section has been tested by our team of experts; or the answer in general has be fact checked.

Learn more on how do we answer questions.

% Declare any needed symbolic variables
syms sig tau_xy sig_x sig_y theta AB BC AC
% Define known stress values
sig_x = sym(160);
sig_y = sym(-120);
sig_val = sym(200);

% Define relationships between AB, BC, and AC
BC = AB*cos(theta);
AC = AB*sin(theta);

% For horizontal equalibrium of the element eqI ¼ sig*AB*cos(theta)-sig_x*BC-tau_xy*AC;

% For vertical equalibrium of the element eqII ¼ sig*AB*sin(theta)-sig_y*AC-tau_xy*BC;

% Solve eqI and eqII for tau_xy
tau_xyI = solve(eqI,tau_xy);
tau_xyII = solve(eqII,tau_xy);

% Take the square-root of tau_xyI times tau_xyII to get tau_xy
tau_xy_val = sqrt(tau_xyI*tau_xyII);

% Substitite the given value of sig into tau_xy
tau_xy_val = subs(tau_xy_val,sig,sig_val);

% Solve eqI for theta and substitute into eqII
eqI = simplify(eqI/cos(theta));
theta_I = solve(eqI,theta);
eqIII = subs(eqII,theta,theta_I);

% Substitute the value of tau_xy into eqIII and solve for the principal stresses (sig_p)
sig_p = solve(subs(eqIII,tau_xy,tau_xy_val),sig);
sig_I = max(double(sig_p));
sig_II = min(double(sig_p));

% Calculate the maximum shear stress using Eq. (1.15)
tau_max = (sig_I-sig_II)/2;

% Output tau_xy, the principal stresses, and tau_max to the Command Window
disp([‘tau_xy = +/-’ num2str(double(tau_xy_val)) ‘N/mm^2’])
disp([‘sig_I =’ num2str(sig_I) ‘N/mm^2’])
disp([‘sig_II =’ num2str(sig_II) ‘N/mm^2’])
disp([‘tau_max =’ num2str(tau_max) ‘N/mm^2’])

\tau_{\max }=\frac{\sigma_{ I }-\sigma_{ II }}{2}  (1.15)

 

The Command Window outputs resulting from this MATLAB file are as follows:

tau_xy = +/- 113.1371 N/mm^2
sig_I = 200 N/mm^2
sig_II = -160 N/mm^2
tau_max = 180 N/mm^2

Related Answered Questions