斜直圆筒内基于Euler-Rodrigues参数的细长杆螺旋屈曲研究 (1).pdfVIP

计 算 物 理 Vo1．29．No．4 第 29卷 第4期 2012年 7月 CHINESE JOURNAL OF COMPUTATIONAL PHYSICS July，2012 文章编号 ：1001．246X(2012)04-0549—08 斜直 圆筒内基于 Euler-Rodrigues参数 的 细长杆螺旋屈 曲研究 袁张贤， 王鑫伟 (南京航空航天大学 机械结构力学及控制国家重点实验室，江苏 南京 210016) 摘 要 ：基于 Euler．Rodrigues参数描述受约束细杆 的后屈 曲变形状态 ，建立相应 的能量泛函，导出非线性平衡微分 方程及接触力表达式 ．采用有限元法对能量泛函进行分析计算 ，通过与现有文献 中数值、理论结果 比较 ，验证 了公 式和求解过程 的正确性．结合算例 ，揭示 了正弦屈曲模态与螺旋屈 曲模态之间的转化过程 ，证实了在受约束细长杆 的后屈 曲响应 中存在不 同屈 曲模态之间的互相转变． 关键词 ：Euler—Rodrigues参数 ；斜直 圆筒；螺旋屈 曲；模态转变 ；有限元 中图分类号 ：0343．9；TB121 文献标识码 ：A 0 引言 受倾斜圆筒约束的细长杆的后屈曲问题具有极其重要的工程意义 ，吸引着大批 国内外学者．该 问题是一 个包含几何、接触 、曲率非线性且有多解的复杂问题 ，至今仍然没有被圆满地解决．自从 Lubinski 首次研究 了无重 圆杆 的螺旋屈 曲问题后 ，关于细长杆的正弦屈曲与螺旋屈 曲问题得到了广泛的关注 ，文献 [2—3]对 前人的工作进行了概括与总结 ，阐述 了该问题 的发展方向和研究近况．由于理论推导 的局限性 ，各种假设的 存在使得理论结果与工程应用还有一段距离．特别是考虑重力影 响时 ，问题 的复杂性大大提高，理论公式一 般仅适用于临界载荷的预测 -5]．Dawson 根据水平井 中计重钻杆屈 曲临界载荷公式得到 了斜直井 中计重 钻杆屈曲临界载荷公式．谈梅兰 基于最小位能原理由Rayleogh—Ritz法得到了计重弹性杆的非线性螺旋屈 曲临界载荷．更为复杂的后屈曲情况 ，广泛采用数值方法 进行研究分析．Mitchell 使用 Galerkin法并结 合最速下降法和 Newton—Raphson法研究了未考虑边界约束 的有重细杆 的正弦屈 曲和螺旋屈曲问题．刘峰 采用有限元法并结合增量加权迭代法研究 了考虑边界约束有重细杆 的屈曲问题．甘立 飞 采用微分求积单 元法研究了斜直井管柱的非线性屈曲问题．迄今为止 ，受约束计重细杆的正弦屈曲和螺旋屈曲问题 ，虽然都 分别得到了广泛的研究 以及较深的认识 ，但两者间的相互转变过程还不是十分 明了 ．采用有限元法和微 分求积法 仅成功地给 出了直井 内钻柱从正 弦屈 曲模态过渡 到 1／4螺旋屈 曲模态 的仿真结果． Munteanu 也仅在无重情况下模拟 出完整的螺旋屈曲模态．而恰是在过渡 阶段容易发生一些工程应用所关 心的问题，如细杆脱离圆筒 、屈曲模态发生转变等． 考虑到 Euler—Rodrigues参数在描述大位移特别是有限转动方面的优越性 ¨ ，本文采用 Euler—Rodrigues 参数描述细杆后屈 曲变形状态 ，由此得出计重细杆的能量泛函，并导出了相应的有限元列式进行求解．通过 与现有结果和理论分析结果的比较来验证公式推导、算法和 自编程序 的正确性．结合具体算例 ，揭示正弦屈 曲与螺旋屈曲之间的转变过程，同时给出后屈曲过程 中不同模态之问发生转变的数值结果． 1 理论分析 1．1 分析模型 考虑一倾斜刚性 圆筒，如图 1所示 ， 为圆筒中心线与垂线的夹角．圆筒 中有一各 向同性圆截面细长杆 ， 收稿 日期 ：2011—09—09；修 回日期 ：2011—11—30 基金项 目：国家 自然科学基金 (No资助项 目 作者简介 ：袁张贤(1988一)，男 ，硕士 ，主要从事工程问题的力学计算建模的研究 通讯作者 ：王鑫伟(1948一)，男 ，教授 ，主要从事计算力学研究 ，E-mail：wangx@nuaa．edu．cn 550