### 实现简谐力作用下弹簧内力计算 对于简谐力 \( F(t) = A \sin(\omega t + \phi) \)，其中 \(A\) 是振幅，\(\omega\) 是角频率，\(t\) 是时间，\(\phi\) 是初相位。假设弹簧遵循胡克定律，则弹簧的内部力可以表示为： \[ F_s(x) = kx \] 这里 \(k\) 表示弹簧常数，而 \(x\) 则是相对于平衡位置的位移。 当考虑简谐外力施加于弹簧上时，系统的运动方程可以通过牛顿第二定律得出: \[ m\ddot{x} + c\dot{x} + kx = A \sin(\omega t + \phi)[^1] \] 为了简化问题并提供一个具体的代码例子，在此忽略阻尼项 (\(c=0\)) 并假定初始条件为零偏置和速度。下面是一个 Python 代码片段来模拟这种情况下的响应以及求解对应的弹簧内力： ```python import numpy as np from scipy.integrate import odeint import matplotlib.pyplot as plt def harmonic_oscillator(y, t, omega, phi, amplitude, k, mass): x, v = y dxdt = v dvdt = (amplitude * np.sin(omega * t + phi) - k * x) / mass return [dxdt, dvdt] # 参数设置 mass = 1.0 # 质量 kg stiffness = 4.0 # 弹簧刚度 N/m damping = 0 # 阻尼系数（本例设为0） frequency = 2.0 # 外部激励频率 Hz phase_shift = 0 # 初相位 radian force_amplitude = 1.0 # 力的幅度 N time_points = np.linspace(0, 10, 500) initial_conditions = [0, 0] # 开始时无位移也无速度 solution = odeint(harmonic_oscillator, initial_conditions, time_points, args=(2*np.pi*frequency, phase_shift, force_amplitude, stiffness, mass)) displacements = solution[:, 0] internal_forces = stiffness * displacements plt.figure() plt.plot(time_points, internal_forces) plt.title('Spring Internal Force Over Time') plt.xlabel('Time(s)') plt.ylabel('Force(N)') plt.grid(True) plt.show() ``` 该程序利用 `scipy` 库中的 ODE 解算器对上述微分方程进行了数值积分，并绘制了随时间变化的弹簧内力图象。