2020iypt，第13题摩擦振子

结果：

改变相关变量后的变化情况：

请注意：该模拟结果仅仅具有参考意义而并无实际的物理意义

改变初始位置我们可以看到振幅发生了变化。

改变动摩擦因数我们看到波动向左侧压缩。

改变两转轴间的间距可以看到波动向右侧拉长。

相关参数在代码中给出。欢迎大家讨论，如果这篇文章对你有帮助的话记得三连支持一下；

谢谢！

matlab代码：

%%初步的动力学基础以及初步求解

clear all

global x l p d h g 

x=0.20;%两轴的一半间距

l=0.4;%重物的长度

p=1.5*10^3;%重物的密度

d=0.05;%重物的宽度

h=0.01;%重物的高度

g=9.8;%重力加速度

m=p*h*l*d;%重物的质量

u=0.46;%动摩擦因数

%——————————————————————————

a2=zeros(10000,1);

v2=zeros(10000,1);

x2=zeros(10000,1);

x2(1,1)=0.01;%重心偏移原点的位置

t=0.001;%时间步长

G=m*g;%重力

for i=2:10000

G=m*g;%重力

N1=((x*G-G*x2(i-1,1))/2+x2(i-1,1)*G)/x;%受力分析的方程

N2=(x*G-G*x2(i-1,1))/(2*x);     %受力分析的方程

a2(i,1)=(N2*u-N1*u)/m;            %受力分析的方程

v2(i,1)=v2(i-1,1)+a2(i,1)*t;  %速度

x2(i,1)=x2(i-1,1)+v2(i,1)*t;  %位移

end

hold on

plot(x2,'Color','b')

%plot(v2,'Color','r')

%plot(a2,'Color','c')