本项目开发一套基于Newmark-beta数值积分算法的结构动力方程求解器，旨在精确计算工程结构在地震、风载或机械振动等动态荷载作用下的时程响应。该程序的核心功能详细包括：1. 动力方程求解：针对单自由度（SDOF）及多自由度（MDOF）线性弹性体系，建立通用的二阶常微分动力平衡方程数值解法。2. 算法灵活配置：完整实现Newmark-beta逐步积分公式，允许用户通过调整参数（gamma和beta）来选择不同的积分方案，如无条件稳定的“平均加速度法”（gamma=0.5, beta=0.25）或条件稳定的“线性加速度法”（gamma=0.5, beta=1/6），以适应不同刚度特性的结构分析需求。3. 自动化计算流程：程序初始化阶段自动计算积分常数并组装有效刚度矩阵，在时间步进循环中递归计算每一时刻的有效荷载向量，进而求解当前步的位移，并根据运动学关系更新速度和加速度。4. 阻尼模型处理：内置瑞利阻尼（Rayleigh Damping）计算模块，可根据用户提供的模态阻尼比和参考频率自动构造系统的比例阻尼矩阵。5. 数据分析与可视化：能够完整记录结构关键节点在全时域内的动力响应数据，并绘制位移、速度、加速度的时程曲线，自动捕捉并输出最大响应值，为结构抗震性能评估和设计优化提供直接依据。
