液体火箭是典型的时变系统，也是复杂的动力学耦合系统，其动力学问题对于火箭的设计与研制极为重要。本文采用多体动力学方法研究液体火箭的动力学建模与仿真问题，主要开展三方面工作：首先，提出了变质量液体火箭的多体动力学建模方法。先基于有限段方法建立未加注火箭的柔性多体系统模型，再提出变质量质点系的达朗贝尔-拉格朗日原理，并结合拉氏乘子法建立柔性变质量充液贮箱的多体建模单元，然后提出变质量贮箱时变晃动的等效力学模型和变质量质点的建模方法。最后，采用状态空间方程方法实现了时变控制在多体系统中的建模；采用约束力方程方法实现了火箭分离过程的多体建模；采用谐波叠加法生成随机样本，实现了火箭跨音速抖振的多体动力学建模；考虑地球外形及其自转，给出了各项时变外载荷的多体建模方法。液体火箭的动力学方程为微分代数方程组，可采用BDF方法和Newton-Raphson迭代法进行数值求解。针对所提出的建模方法给出若干验证算例。通过多体动力学仿真结果与文献结果和解析解的对比，验证了多体建模方法的正确性。其次，根据本文提出的建模方法建立大型捆绑式液体运载火箭的多体系统模型，开展火箭“大气飞行-助推分离”连续过程的精确动力学仿真，实现了姿控、弹道、推进剂消耗与晃动以及箭体振动的耦合分析。根据多体动力学仿真结果分析了飞行箭体的动力学行为，预示了飞行箭体的姿态运动、各贮箱的推进剂晃动、箭体飞行轨迹以及助推器的分离过程；根据多体仿真结果还分析了箭体结构振动和动载荷，给出了姿控喷管摇摆幅度、箭体关键部位动间隙以及飞行全程全箭的动载荷包络。跨音速脉动压力使得箭体姿态、喷管摆角和箭体动载荷在跨音速时段之内剧烈振荡，但对飞行轨迹无显著影响，而其对晃动的扰动影响将持续到跨音速时段之后。最后，采用多体动力学建模方法实现了大型捆绑式液体运载火箭“结构-姿控-Pogo”耦合大回路的稳定性分析，研究了火箭姿控回路和Pogo回路之间的相互耦合影响。基于火箭的多体系统模型分别开展姿控回路、Pogo回路以及耦合大回路三种情况下的稳定性分析，通过结果对比阐明了大型捆绑火箭的姿控回路可影响Pogo稳定性，空间耦合模态对姿控回路和Pogo回路间的耦合程度有关键影响，从而明确了耦合大回路稳定性分析的必要性。