您需要使用元素的操作。
clc;清晰的所有;
清晰的所有
关闭所有%赫兹接触计算器一个球体
clc;
格式长
F = 5000;
v1 = 0.3;
v2 = 0.3;
E1 = 700000
E2 = 700000;
d1 = 10;
d2 = 100;
b = (1 / d1) + (1 / d2)
m1 =(个v1 ^ 2) / E1
m2 = (1-v2 ^ 2) / E2
a0 = 0.76241947
一个= (0.182:0.1:2.286);
z = [];
U = [];
Q = [];
好吧= [];
Pmax = [];
qxy = [];
Pmax = (3 * F) /(2 *π* (^ 2));
z = 0.48 . *;
U = (2。* (1 + (z ^ 2 /。^ 2)))% qxy方程的简化——>不应该只提供一个输出!
Q = (1. / U);% qxy方程的简化
好吧= (1 - (abs (z / a)。*(每股(1. / (z / a)))))% qxy方程的简化——>不由于矩阵划分错误! !
% < - - - - - - -上述方程是改变
qxy = -Pmax。* (((OK)。* (1 + v1)) - q)
情节(Pmax)
%我试图让球体的赫兹接触应力的计算器
%我需要找到“qxy”值为每个迭代的“一个”,但我不明白
%为什么“u”变量仅创建一个输出。之后,我还需要画出“qxy”值。
%我会很感激你的解决我的问题和反馈
%的代码会很感激。提前谢谢你