how to create function for this problem..

i want to create a function for Double Spherical isolation bearing but still not able to do that.. can anyone create the function which can simulate this bearing..

This device is used to isolate the superstructre from substructure to reduce the earthquake damage....

thanks in advance

[Daniel kiracofe]

Suggest that you post the code that you have tried so far, and then maybe we can help to improve it.

https://www.mathworks.com/matlabcentral/answers/6200-tutorial-how-to-ask-a-question-on-answers-and-get-a-fast-answer

[Habib Muhammad]

function [ Nh,Kh,Beff] = isolation(W,D,h)

% Parameters for base isolator 
clear all
clc
global Nh W D H R1 R2 h
% Nh is the Horizontal component of the force 
W = 50e10;           % Vertical load supported by the bearing
D = 0.2;             % Lateral displacement between the upper and lower bearing plates
H = 0.3;             % Distance between the centers of the two spherical surfaces
R1= 20;              % Curvature radius of the sliding spherical surface
R2= 15;              % Curvature radius of the rotation spherical surface
h = 25;              % Depth of the slider
mu=0.01;             % Alpha Angle between the straight line of the spherical centers and 
% the vertical direction after lateral displacement appeared

Nh = W*D/H
H = R1 + R2 - h
Dy= diff(D)
F = Nh  + mu*W*(Dy)
% mu is the coefficient of friction pairs on the sliding spherical surface
% Dy is the horizontal velocity

Kh = W/H    % Lateral stiffness 
Beff=2/((D/mu*h)+1)*pi
end

[Habib Muhammad]

here is the function i want to develop.. i want to input the W D and h or load it from a txt file and in response i want to get Nh, Kh and Beff...