单层厚壁圆筒Green函数及其动应力分布.docVIP

单层厚壁圆筒Green函数及其动应力分布 　　摘要通过建立轴对称冲击载荷作用下单层厚壁圆筒的Green函数，得到了问题的位移分布的解析解和相应的动应力分布表达式，此解法不仅避免了将动力学的一般解分解位拟静态解和自由振动解叠加的过程，而且解的形式适用于各种初始条件，便于数值计算和应力分布规律的分析.通过一个实例的计算，表明本文求解方法的正确性，并以此结果具体分析了单层厚壁圆筒的动应力分布规律. 　　关键词轴对称；圆柱壳；Green函数；应力分布 　　中图分类号O347.1，O174.6 文献标识码A 文章编2013 　　实心柱体或单层厚壁圆筒在平面轴对称作用下的动应力历程和分布规律是一个典型的弹性力学问题，对工程应用具有重要意义.圆柱壳除了在军事、爆炸加工等方面有广泛的应用外，在石油、化工、航空、原子能、海洋开发等工程领域也被广泛应用，它们在各种载荷作用下的动应力响应问题方面已有很多的研究和结果[17].文献[8]利用有限Hankel 变换的方法从理论上解决了单层厚壁圆筒受动载荷作用的响应问题，文献[9]利用有限Hankel变换和Laplace变换的方法得到了动载荷作用下单层后壁圆筒的弹性动应力，文献[10]推导出了圆筒在轴对称平面应变条件下的弹性动响应解.这些方法都有各自的优点，但也存在一些缺陷，例如解的形式不能适用于各种载荷，不便于进行数值计算和结果分析，文献[11]研究了单层厚壁圆筒弹性动应力的简化计算方法以利于工程直接应用. 　　本文则通过建立问题Green函数的方法，利用冲量原理和变量分离的方法对Green函数求解[12]，得到了问题位移的解析解和相应的动应力场表达式，此方法不仅避免了复杂的积分变换，而且解的形式使用于各种有效载荷，便于数值计算和轴对称冲击载荷作用下厚壁圆筒的动应力响应历程和分布规律的分析. 　　将问题的位移分布以其Green函数表示，问题的关键在于计算本征函数的模方，在使用Matlab进行数值计算时，需要求出特征方程的根.文献[13～14]提供了Bessel方程本征值问题特征方程根的Matlab算法，为本文的数值计算和结果分析提供了方法基础.本文的最后，通过一个实例的计算，表明了本文方法的正确性，同时作者还就两种不同的情形具体分析了平面轴对称冲击载荷作用下厚壁圆筒的动应力响应历程和分布规律并得到了一些有意义的结论. 　　4结论 　　1） 将轴对称冲击载荷作用下厚壁圆筒的位移分布用其Green函数表示并相应地得到问题的动应力分布解析解，此方法不仅避免了将问题分解为拟静态动力学解和动态解的叠加的过程，也避免了复杂的积分变换，动应力表达式的形式便于进行数值计算动应力分布规律的分析.这种方法可推广应用于其他问题的讨论. 　　2） 图1表示圆筒外壁受瞬时冲击时径应力随时间的变化规律，从图中可以看出，当应力波波头随时间由冲击壁向另一壁传播而尚未到达某点时，该点的动应力为零.当波头到达该点时，该点动应力值呈强间断跳跃.而当波头远离该点时，该点的动应力值平缓地趋于静载时的静应力值，随着内壁的反射作用，应力波呈现拉伸现象并迅速地趋于静载时的静应力值，直到冲击壁的下一个应力到来，再次呈现强跳跃.这一现象与[7]相同，从而验证了本文解的正确性. 　　3） 图2表示圆筒外壁受瞬时冲击时在同一时刻沿径向各点的径应力变化规律.从图中可以看出，在内外壁附近的点，径向动应力值均出现显著变化，并且最大动应力值出现在内壁附近. 　　4） 图3和图4分别表示圆筒外内径比为s=4时外壁受瞬时冲击时径应力随时间的变化规律和在同一时刻沿径向各点的径应力变化规律.随着圆筒厚度的增加，径向动应力先是出现缓慢的拉伸，然后随着应力波头的到来而出现跳跃，内壁反射波到达的时间加长，趋于静载时的静应力值的时间也相应加长了.最大压缩应力出现在冲击壁附近，而在另一壁附近则出现最大拉伸应力. 　　5） 图5到图8表示圆筒内壁受瞬时冲击时径应力随时间的变化规律和在同一时刻沿径向各点的径应力变化规律，所得到的结果与外壁瞬时冲击时的结论形似，不同的是，当圆筒较薄时，最大压缩动应力均出现在内壁附近. 　　参考文献： 　　[1]朱文辉，薛鸿陆，韩钧万，等.爆炸容器动力学研究进展评述[J].力学进展， 1996，26（1）：6878. 　　[2]吴洪飞.容中柱形壳和球壳的耐外压实验研究[J].哈尔滨工业大学学报， 2002，34（3）：367369. 　　[3]黄承义，郑际嘉，刘土光.两端简支圆柱壳在径向冲击力下的动力响应分析[J].华中理工大学学报， 1997，25（1）：5154. 　　[4]郑津洋，陈勇军，邓贵德，等.强动载荷下离散多层圆筒的弹性动力响应分析[J].浙江大学学报：工学版， 2005