如何根据以下等式绘制 h / L 和 h / L0 之间的曲线:
[(h / L){((Dk ^ 4 / g)+ 1)tanh(kh)}] ^(1/2)=(h / L0)^(1/2){1-(TU /小时)(H / L)}
我自己尝试了但是我无法获得正确的曲线。我的代码如下:
D=10^11;
g=9.8;
T=1;
U=5;
h=1;
a=T*U/h;
L= (-3:0.01:3);
L0=(-3:0.01:3);
y=h./L;
x=h./L0;
k=(2*pi/h).*y;
y= (x.*(1-a.*y).^2)./(((16*D*(pi^4)*(y.^4)/(g*h^4))+1).*tanh(2*pi.*y));
plot(x,y,'b');
xlabel('h/L0');
ylabel('h/L');
此处g
是由重力引起的加速度,k
是波数。 D
是刚性因子,等于 D = EI / d 其中 d 是密度,而 EI 是刚性和 k = 2pi / L 。