Scribd
Upload
270 views2 pages

MAtlab Code For Ball and Ball Contact Mechanics

This document contains code to calculate contact parameters for two steel balls in contact or a ball contacting a flat surface. The user inputs the radii of the balls, normal force, Young's moduli, and Poisson's ratios of the materials. The code then calculates the reduced radius of curvature, reduced Young's modulus, contact area, maximum and average contact pressures, maximum deflection, maximum shear stress, and depth of maximum shear stress.

Uploaded by

Shivshankar Mudaliar
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
270 views2 pages

MAtlab Code For Ball and Ball Contact Mechanics

This document contains code to calculate contact parameters for two steel balls in contact or a ball contacting a flat surface. The user inputs the radii of the balls, normal force, Young's moduli, and Poisson's ratios of the materials. The code then calculates the reduced radius of curvature, reduced Young's modulus, contact area, maximum and average contact pressures, maximum deflection, maximum shear stress, and depth of maximum shear stress.

Uploaded by

Shivshankar Mudaliar
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 2

disp('contact parameters for two steel balls')

prompt = 'enter the radius of the ball A in meters = ';


Ra = input(prompt);
prompt = 'enter the radius of the ball B in meters = ';
Rb = input(prompt);
prompt = 'enter the normal force in N = ';
W = input(prompt);
prompt = 'enter Youngs modulus for ball A in N/m^2= ';
Ea = input(prompt);
prompt = 'enter Youngs modulus for ball B in N/m^2 = ';
Eb = input(prompt);
prompt = 'enter Poisson ratio of ball A material = ';
va = input(prompt);
prompt = 'enter Poisson ratio of ball B material = ';
vb = input(prompt);
x= (1/Ra)+(1/Rb);
Rx=(1/x);
y=(1/Ra)+(1/Rb);
Ry=(1/x);
R_1= (1)/((1/Rx)+(1/Ry)); % reduced radius of
curvature (m)
xx=0.5*(((1-va^2)/(Ea))+((1-vb^2)/(Eb)));
E_1=(1/xx); % reduced Young's modulus
(N/m^2)
a=(3*W*R_1/E_1)^(1/3); % Contact area dimensions
(m)
Pmax=(3*W)/(2*pi*a*a); % Maximum contact pressure
(MPa)
Pavg=(W)/(pi*a*a); % Average contact pressure
(MPa)
def=(1.0397)*((W*W)/(E_1*E_1*R_1))^(1/3); % Maximum deflection (m)
shear_max=(Pmax/3); % Maximum shear stress (MPa)
z=0.638*a; % depth at which maximum
shear stress occurs(m)

disp('Contact between a ball and a flat Surface')


prompt = 'enter the radius of the ball in meters = ';
Ra = input(prompt);
prompt = 'enter the normal force in N = ';
W = input(prompt);
prompt = 'enter Youngs modulus for ball in N/m^2= ';
Ea = input(prompt);
prompt = 'enter Youngs modulus for flat surface in N/m^2 = ';
Eb = input(prompt);
prompt = 'enter Poisson ratio of ball material = ';
va = input(prompt);
prompt = 'enter Poisson ratio of flat surface material = ';
vb = input(prompt);
x=(1/Ra);
Rx=(1/x);
y=(1/Ra);
Ry=(1/x);
R_1= (1)/((1/Rx)+(1/Ry)); % reduced radius of
curvature (m)
xx=0.5*(((1-va^2)/(Ea))+((1-vb^2)/(Eb)));
E_1=(1/xx); % reduced Young's modulus
(N/m^2)
a=(3*W*R_1/E_1)^(1/3); % Contact area dimensions
(m)
Pmax=(3*W)/(2*pi*a*a); % Maximum contact pressure
(MPa)
Pavg=(W)/(pi*a*a); % Average contact pressure
(MPa)
def=(1.0397)*((W*W)/(E_1*E_1*R_1))^(1/3); % Maximum deflection (m)
shear_max=(Pmax/3); % Maximum shear stress (MPa)
z=0.638*a; % depth at which maximum
shear stress occurs(m)

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy