Modeling of square wave

Hello everyone,

I am modeling the position varying gear mesh stiffness function with the following code inside ODE solver.

z(k) = mod(theta_p,2*pi/t_np);

if z(k)>=(2-c)*2*pi/t_np

k_m(k) = k_max;

else

k_m(k) = k_min;

end

where c, k_max, k_min, t_np are the constants:

k_min = 5e7/2; %minimum stiffness;

k_max = 5e7; %maximum stiffness;

c = 1.63; %engagement rate;

t_np = 25; %number of teeth on pinion;

The goal of this function is to make a square wave with t_np number of maximums over 2*pi interval.

This is how the time/varying gear mesh stiffness should look like:

When the gear have a fault, the stiffness of the faulty gear tooth is compromized, so when the broken tooth is in contact, the stiffness will decrease. I am struggling to figure out how to make the value of one of the minima less than others on each 2*pi intervals. Any suggestions highly appreciated.

Thank you!

[Steve Miller]

If you are interested in modeling gears with faulty behavior, you may be interested in the Simple Gear block in Simscape Driveline.

--Steve