摘要

采用Euler四元数表示的Kirchhoff方程来研究受力挤压作用下的弹性细杆的拓扑构形,进一步研究弹性细杆的力学性质;将得到的微分方程与约束条件组成微分代数方程后再转化为微分方程规范形式以便求解;为满足边界条件,应用数值打靶法求解边值条件,并将弹性细杆在力作用下的拉压过程用Matlab仿真出来.同时对由于误差导致的违约现象进行处理,并针对欧拉参数的特征,选取合适的修正系数以保持方程的稳定性.