基于结构动能的动态载荷等效静态转化方法.docxVIP

基于结构动能的动态载荷等效静态转化方法 在真实的物理世界中，大多数结构都是在动态负荷的作用下工作的。由于动态响应计算的复杂性和高耗时，结构动态响应的优化设计难以实施。韩国科学家g.j.公园总结并发展了基于有限模型的位移的动态载荷等效方法。致动函数。在基于等效场条件的情况下，动态载荷对应于静态载荷。等效原理是结构等效前后的位移场是相同的。提出基于等效补偿变换方法的应用，并将位移补偿的动态响应优化进行。 在G.J.Park文献中,虽然说明了静态载荷等效的关键时间点,即位移最大时间点作为等效载荷转化的时间点,然而结构在动态载荷的作用下,最大位移点与最大应力点有一个时间差;另外,结构在动态载荷的作用下,除了在动态力作用方向位移比较大之外,其他方向都有位移,如果只关心发生最大位移的方向而不考虑其他方向,从一定意义上来说不合理;如果只考虑应力也不合理. 总之,采用应力等效存在以下不足:①最大应力仅仅说明结构的局部受力状况,不能反映应力的分布情况;②由于动态分析采用有限元法计算,而结构最大应力与单元类型及单元尺寸有关,尤其在应力集中的局部位置,由于应力变化比较大,因此单元类型与尺寸对其影响非常大;③大量计算表明,一个结构在不同的单元类型与尺寸下获得的最大应力值差别非常明显.同理,采用位移等效方法也存在类似的问题.因此本文提出基于结构势能动态载荷等效方法. 1 应力刚度变化 结构在动态载荷作用下,若忽略结构阻尼,能量满足如下关系: 式中:为动能;为结构势能(其中为体积应变能,为应力刚度能,也称为外力势能;为外力功;S为应力刚度矩阵或应力刚化矩阵(stress stiffn ess matrix);U为位移向量;为速度向量;M为质量矩阵;K为刚度矩阵;F(t)为动态载荷;u(t)为结构在力方向的位移. 由于结构线性振动时,动能比较小,因此将式(1)近似为 2 动态载荷作用下的关键点识别 Courant1943年的工作就体现了有限元法的思想,1960年Clough第一次提出了有限元法的概念.有限元法为计算分析结构动力学问题提供了较为有效的方法,然而,结构动态分析过程非常复杂且耗时较大.因此,为有效解决动态响应优化问题,考虑用静态载荷等效代替动态载荷,关键问题是首先要确定在哪个时刻进行等效,即在动态载荷作用下结构关键点的识别. 关键点是指动态载荷作用下的结构最危险的时刻.一直以来,学术界普遍认为关键点即结构在动态载荷作用下,结构位移绝对值最大的时刻,或者部分研究中将动态载荷最大的时刻作为关键点.作者在前期的研究中发现,结构在动态载荷作用下位移达到最大值的时刻与应力达到最大值的时刻有一个时间差,另一方面,动态载荷最大的时刻,结构位移与应力都远远没有达到最大值,因此不论以位移绝对值最大或者动态载荷最大的时刻作为关键点都不合适.因此,作者提出基于动态应力解空间谱元离散的关键点识别方法,并在此基础上进行动态载荷等效静态转化. 3 结构的模态叠加法求解 结构在动态载荷作用下,基于结构势能进行动态载荷等效静态转化方法如下. 步骤1建立有限元模型. 运动微分方程为 式中:M=ΦmTΦ,C=ΦTcΦ,K=ΦTkΦ,F=ΦTf(t),q为模态位移向量;为模态加速度向量;为模态速度向量;f(t)为时间的任意函数;Φ为模态矩阵;m为质量系数;c为阻尼系数;k为刚度系数.在本文中假设c=0. 步骤2瞬态动力学分析. 采用模态叠加法.先求解式(3),获得q,然后得到结构的节点位移x=qΦ,再由公式σk(t)=DBkqkΦ获得第k阶模态应力,通过公式得到结构动态应力. 步骤3识别关键时间点. 获得结构动态应力解空间,将其在高斯—勒让德-罗巴托点谱元离散,构造时间点与其对应的动态应力解空间矩阵,应用Lagrange插值,得到高精度的近似函数,最后识别关键时间点.关键时间点识别的数学模型: 式中:m为解空间离散的单元数;k为每一个单元离散GLL点数;σi(tk)为第i个单元的第k个GLL点对应的动态应力;t0为仿真开始时间,通常为0;t1为仿真结束时间. 步骤4计算关键时间点对应的结构位移分布、结构势能. 根据步骤③中得到的关键时间点,通过公式(1)计算结构在这一点的结构势能及位移分布. 步骤5建立结构动态载荷等效静态转化数学模型 式中: 式中:为静态力作用下整体结构势能;为动态力作用下整体结构势能;为静态力作用下某单元结构势能;为动态力作用下某单元结构势能. 步骤6求解. 由于结构等效载荷是未知的,其特性极为复杂,一般为非线性多峰值函数,因此采用全局优化算法——区域细分算法(DIRECT)进行搜索结构等效静态载荷. 4 方法应用 4.11 密度和结构试验 本算例的动态载荷为常数,结构弹性模量为69 GPa,密度为2 710 kg/m3,截面积为1.94×10m-42,泊松比为