拉压不同模量问题非线性分析高性能算法研究

来源 :大连理工大学 | 被引量 : 6次 | 上传用户:sunto0724
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
实际工程中,诸如混凝土、玻璃纤维、石墨等材料都呈现出拉、压性能不同的力学行为;且随着温度的升高,材料的拉、压性能差异逐渐增大。近年来,随着材料制造技术的迅猛发展和工业生产规模的不断扩大,大量的复合材料被广泛地应用到航空、航天、土木工程等领域。一些复合材料不仅表现出各向异性的力学特性,在单个材料方向上同样也呈现出拉伸和压缩不同的力学行为。具有拉、压不同模量的材料也被称作双模量材料。对于这类材料而言,若仍采用经典的胡克定律对其进行力学分析,势必会引入很大的分析误差。这时,便需要采用拉、压不同模量弹性理论(也称作“双模量弹性理论”)。大量研究表明基于双模量弹性理论的材料力学行为数值分析通常存在收敛性困难。从计算方法的角度讲,发展针对拉、压不同模量材料力学分析的高性能算法显得十分必要。另一方面,由绳索和杆件组成的张拉整体结构,以及薄膜结构都具有轻质、可折叠等特性,可用作可展开的小卫星和太阳能帆板等结构的部件。绳索和薄膜均能承受较大的拉伸载荷;而当承受压缩载荷时,绳索会松弛,薄膜则会局部起褶。此类柔性结构的力学分析呈现出很强的非线性特征。若直接进行后屈曲分析,算法鲁棒性往往欠佳。因此,寻找合适的力学模型,并发展高效、稳定的求解方法对于实际工程分析具有重要意义。本文首先介绍了拉、压不同模量材料、tensegrity结构和薄膜褶皱分析的研究背景和现状,特别讨论了拉、压不同模量问题的研究意义和难点。分别建立了拉、压不同模量桁架、平面连续体小变形和大位移、小应变力学分析的参变量变分原理和数值方法,推导了适用于clustered tensegrity结构和经典tensegrity结构的统一切线刚度矩阵,并发展了高效稳定的计算方法,实现了对结构主动激励和展开过程的数值模拟。本文主要内容包含以下五个部分:第一,建立了适用于拉、压不同模量桁架小变形静、动力分析的参变量变分原理和保辛算法。推导了拉、压不同模量杆件在不同应力状态下的参变量统一本构方程,并证明了其与原本构方程之间的等价性。基于参数最小势能原理,推导出问题的有限元平衡方程和互补方程。新算法不需要更新刚度矩阵,表现出较高的计算效率。此外,基于Hamilton力学体系发展了拉、压不同模量杆件的动力参变量变分原理和保辛算法,将算法扩展到拉、压不同模量桁架动力响应分析。第二,建立了适用于拉、压不同模量桁架和tensegrity结构大位移、小应变分析的参变量共旋算法。利用共旋坐标法将杆件的大位移、小应变变形分解为整体坐标系下的刚体运动和局部坐标系下的小变形;绳索抗拉、不抗压的力学行为由局部坐标系下的参变量统一本构方程描述。推导了含有参数变量的杆件单元切线刚度矩阵。绳索和杆件单元的节点内力首先在局部坐标系下由参变量统一本构方程计算得到,而后由转换矩阵将其转换到整体坐标系下,集成结构节点内力向量。在非线性求解过程中,参变量共旋算法避免了对杆件和绳索拉、压状态的直接判断,表现出稳定的收敛性。算法也能准确预测tensegrity结构中松弛绳索的分布情况。第三,推导了针对clustered tensegrity结构非线性分析的切线刚度矩阵,并发展了高效的计算方法。由于clustered绳索的存在, clustered tensegrity结构比经典tensegrity结构拥有更多的内部机构运动数目,表现出更“柔”的力学特性,因而具有主动易控和可展开(折叠)的优良特性。本文基于共旋坐标法,推导了一个新的切线刚度矩阵。与经典tensegrity结构的切线刚度矩阵相比,该刚度矩阵包含了描述滑轮运动的附加刚度部分,且能够退化到经典结构的刚度矩阵,因而对于有限元程序实现和软件集成都十分方便。利用新的切线刚度矩阵,实现了对clustered tensegrity结构的主动激励和展开过程的数值模拟。本文方法表现出比力密度法和动力松弛法更好的收敛性。第四,针对拉、压不同模量二维连续体力学分析,建立了单元体在四种不同主应力状态下的参变量统一本构方程和参数势能表达式,并证明了统一本构方程的等价性。对于数值离散后的势能表达式,由参变量变分原理推导出问题的有限元平衡方程和互补方程,发展了相应的数值算法。算法执行过程中,避免了直接地判断单元的拉、压应力状态,因而不需要更新刚度矩阵。与传统算法相比,新算法克服了刚度矩阵随应力状态变化而容易发生奇异这一难点,展现出良好的收敛性。进一步考虑惯性力项,建立了拉、压不同模量平面材料的动力参变量变分原理和参数势能表达式,发展了新的动力积分算法。静、动力分析算法可对薄膜的褶皱分布区域及其变化过程进行定性预测。第五,建立了适用于拉、压不同模量二维连续体大位移、小应变分析的参变量共旋算法。在小变形分析的基础上,进一步考虑几何非线性因素;将参变量变分原理与共旋坐标法结合,推导了含有参数变量的单元切线刚度矩阵。由于在局部坐标系下采用了参变量统一本构方程,计算单元节点内力时无需通过判断单元当前所处的主应力状态来更新局部刚度矩阵,很大程度上改善了算法的收敛性。数值算例体现了参变量共旋算法相对于传统算法的收敛性优势。基于材料修正思想,提出了一个基于主应力的薄膜褶皱分析模型,并将参变量共旋算法成功地应用到薄膜起褶数值分析。两个薄膜褶皱分析的典型算例证明了新模型和新算法的可行性。
其他文献
中国国有企业委托代理关系的典型特征是委托人和代理人的政府化.其主要特征是委托人和代理人都由政府任命.这种政府化的委托人和代理人对企业治理结构的影响主要表现为:第一,
目的 :探讨血浆D-二聚体检测在诊断小儿肺炎支原体肺炎中的应用价值。方法 :将2012年1月至2017年1月期间常州市金坛区人民医院儿科收治的86例肺炎支原体肺炎患儿作为研究对象,
前不久,静海县团泊镇召开了电镀行业专项整治工作会议,县领导袁景生及县工委、环保局主要负责同志出席会议,全镇20家电镀企业厂长、经理参加了会议。
近年来,在多方面利好因素的影响下,使得我国城市建设工作得到了全面的落实,在这种发展形势下,使得人们对城镇室内燃气设施安全问题越发的关注。室内燃气设施的安全性往往会对
由于飞机要在高空中高速飞行,其制造过程中,对于生产工艺具有较为严格的要求,飞机的质量关系到乘客的生命财产安全。本文就数控加工技术在飞机制造中的应用进行探究,分析使用
东阳市目前正在积极开展电镀行业污染综合整治行动,争取在11月底前完成污染整治。
电沉积MnO2薄膜及其性能研究;Ti/SnO2/ZnO复合薄膜的制备及其光催化性能研究;硬铝合金表面微弧氧化着色工艺研究;三价铬硫酸盐电镀铬的发展现状;表面冶金涂层的磨损失效形式与对策
期刊
近些年来,水力发电企业的外部环境一直在发生着变化,导致企业经营者的考评管理和薪酬管理还存在一定的问题,如果不能及时处理这些问题,很有可能给企业的后续经营带来巨大的影
目的 :探讨用十二导联动态心电图检查诊断阵发性心房颤动的效果。方法 :随机选取2016年1月至2017年4月期间江苏省南通市肿瘤医院收治的80例阵发性心房颤动患者作为研究对象。根
以往认为心律失常的发生主要与心肌细胞兴奋性异常有关但近年来的研究表明细胞间电耦联障碍是心律失常的另一个重要原因之一,其所起的致心律失常作用甚至大于兴奋性异常。由心