精确Cosserat弹性杆动力学的分析力学方法.pdfVIP

物理学报 Acta Phys. Sin. Vol. 62, No. 4 (2013) 044601 精确Cosserat 弹性杆动力学的分析力学方法* 1† 12 2 薛纭 翁德玮 陈立群 1) ( 上海应用技术学院机械工程学院, 上海 201418 ) 2) ( 上海大学, 上海市应用数学和力学研究所, 上海 200072 ) ( 2012年7月11 日收到; 2012年8月29 日收到修改稿) 以脱氧核糖核酸和工程中的细长结构为背景, 大变形大范围运动的弹性杆动力学受到关注. 将分析力学方法运 用到精确Cosserat 弹性杆动力学, 旨在为前者拓展新的应用领域, 为后者提供新的研究方法. 基于平面截面假定, 在 弯扭基础上再计及拉压和剪切变形形成精确Cosserat 弹性杆模型. 用刚体运动的概念描述弹性杆的变形, 导出弹性 杆变形和运动的几何关系; 在定义截面虚位移及其变分法则的基础上, 建立用矢量表达的d’Alembert-Lagrange 原理, 在线性本构关系下化作分析力学形式, 并导出Lagrange 方程和Nielsen 方程, 定义正则变量后化作Hamilton 正则方 程; 对于只在端部受力的弹性杆静力学, 导出了将守恒量预先嵌入的Lagrange 方程, 并讨论了其首次积分. 从弹性杆 的d’Alembert-Lagrange 原理导出积分变分原理, 在线性本构关系下化作Hamilton 原理. 形成的分析力学方法使弹 性杆的全部动力学方程具有统一的形式, 为弹性杆动力学的对称性和守恒量的研究及其数值计算铺平道路. 关键词: 精确Cosserat 弹性杆, 分析动力学方法, 变分原理, Lagrange 方程 PACS: 46.70.Hg, 45.20.Jj, 45.10.Db DOI: 10.7498/aps.62.044601 方法已有所研究, 由于忽略了弹性杆截面因变形导 1 引言 致随弧坐标的平移, 即拉压和剪切变形, 使得变分 原理中与主矢有关的项以及关于主矢的3 个动力 弹性杆力学具有悠久的历史和广泛的应用背 学方程不具有分析力学的形式, 并且因主矢与主矩 景. 20 世纪70 年代, 弹性杆被作为脱氧核糖核酸 的耦合使得主矩方程的表达也不尽如人意8−12 , (DNA) 等生物大分子的力学模型研究其位形的拓 给应用带来了不便. 在弯扭基础上再计及轴向拉压 扑性质、平衡及其稳定性, 进一步推动了弹性杆力 和截面剪切变形的弹性杆模型是Kirchhoff 模型的 学的发展1−4 . 鉴于应用背景的复杂和多样性, 弹 改进1516 , 这被称之为Cosserat 弹性杆或精确模 性杆力学面临新的问题: 其一是极端细长导致小 型, Kirchhoff 弹性杆是其特例. 应变累积成超大位移, 形成异常复杂的几何形态; 本文引入了关于主矢的本构方程后形成的精 其二是全新的约束形式, 如自接触等, 对组蛋白的 确Cosserat 弹性杆的分析力学方法, 其变分原理和 缠绕包括了经典约束具有的单面、非定常、非完 动力学方程完全具有分析力学的形式, 这既有助于 25 整等特征 . 因分析力学的建模方法在约束处 充分利用分析力学的现有成果为弹性杆力学提供 67 理和数值计算上的优势 , 文献 [8—12] 建立了