一种确定冲击载荷等效静载的方法技术

本发明专利技术属于结构、设备抗冲击设计技术领域，涉及一种可以确定冲击载荷的等效静载荷的方法。步骤一，根据以下公式计算动态缩放系数λ：

Method for determining equivalent static load of impact load

The invention belongs to the technical field of structure and equipment shock resistance design, and relates to a method for determining equivalent static load of impact load. The dynamic scaling factor lambda is calculated according to the following formula:

【技术实现步骤摘要】
一种确定冲击载荷等效静载的方法
本专利技术属于结构、设备抗冲击设计
，涉及一种可以确定冲击载荷的等效静载荷的方法。
技术介绍
所述等效静载荷指，针对不便于甚至无法直接进行冲击仿真或试验的设备、结构，将其冲击载荷等效为静载荷，为冲击的等效静态仿真或试验提供输入。文献《冲击载荷作用下的结构设计方法》(导弹与航天运载技术，2007年第3期：33～37)给出了一种确定冲击载荷等效静载的公式：式中：TD为冲击载荷持续时间；ωn为系统固有圆频率；λ为动态缩放系数，定义为等效静载荷Ps与冲击载荷峰值P之比：当已知系统的固有圆频率ωn和冲击载荷持续时间TD之后，先求得动态缩放系数λ，那么冲击载荷的等效静载荷为冲击载荷峰值与动态缩放系数的乘积：λP。该文献给出了一种计算结构冲击强度的方法，可以方便地给出冲击载荷的等效静载荷。不足之处在于：该方法缺乏严谨的理论依据，动态缩放系数的求解公式采用近似值，并且均大于1。因此利用该方法得到的等效静载荷总大于冲击载荷峰值，使得安全系数过大，导致设计保守。
技术实现思路
本专利技术的目的：提供一种具有高真实度特征的冲击载荷等效静载确定方法。本专利技术的技术方案：一种确定冲击载荷等效静载的方法，其特征在于，所述的方法包括如下步骤：步骤一，根据以下公式计算动态缩放系数λ：其中，TD为冲击载荷的持续时间；fn为待冲击结构的基频；步骤二，根据以下公式计算冲击载荷的等效静载Ps：Ps＝λP其中，P为冲击载荷的峰值。优选地，通过振动模态试验或有限元仿真的方法获取冲击结构的基频fn。优选地，根据GJB150A.18-2009确定冲击载荷的持续时间TD。本专利技术的有益效果：本技术方案计算得出的λ值，根据TDfn的变化能够出现大于和小于1的情况，真实地等效了冲击载荷，提高了等效静载的真实度。另外，该方法思路简单、明确，理论依据充分，结果精度较高，应用范围广，适合多种类型的冲击载荷，例如：着陆、坠撞、水上迫降、离散元撞击、空投空降、弹射和拦阻等。利用冲击载荷的等效静载，方便对结构、设备进行考核，有助于解决工程相关问题，节约工时和费用，具备广阔的应用前景。附图说明图1受外力作用的单自由度动响应系统；图2受基础加速度作用的单自由度动响应系统；图3λ随TDfn的变化曲线；图4确定冲击载荷等效静载的流程图；图5确定TE的示意图；图6梁式结构示意图；图7冲击力时间历程；图8机身某处的压力时间历程。具体实施方式步骤一，建立动态缩放系数的理论计算公式。基于冲击载荷为后峰锯齿波形式、冲击主要诱发基频响应、忽略阻尼对冲击响应的影响这3种假设，根据冲击响应位移最大值与静位移一致的原则，推导出动态缩放系数的理论计算公式。具体推导过程如下：1.冲击载荷的表达式采用后峰锯齿波描述冲击载荷L(t)(外力或基础加速度)波形，表达式如下：式中：P为冲击载荷峰值，TD为冲击载荷持续时间。2.单自由度系统的动响应图1和图2显示了两种常见的单自由度动响应系统：一是受外力作用；二是受基础加速度作用。下面分别推导这两种系统的动态响应。2.1受外力作用的单自由度动响应系统零初始条件下，无阻尼单自由度无阻尼系统在任意力f(t)作用下，其响应位移x(t)可写成：式中：m为系统质量；ωn为系统固有圆频率，表达式如下：式中：fn为系统固有频率。将式(1)所示的外力代入式(2)可得：2.2受基础加速度作用的单自由度动响应系统零初始条件下，无阻尼单自由度无阻尼系统在任意基础加速度a(t)作用下，其相对响应求解方程可写成：其响应相对位移z(t)表达式为：同理，将式(1)所示的基础加速度代入式(6)可得：3.动态缩放系数定义动态缩放系数λ为等效静载荷Ps与冲击载荷峰值P之比：在静力fs作用下，外力系统的静态位移为：在静态基础加速度as作用下，基础加速度系统的静态相对位移为：根据位移等效原则即动态位移最大值等于静态位移，可建立动态缩放系数λ的表达式。对于受外力作用的单自由度动响应系统，令式(4)的最大值等于式(9)，可得：根据式(8)的定义，可得动态缩放系数为：对于受基础加速度作用的单自由度动响应系统，令式(7)的最大值等于式(10)，可得：根据式(8)的定义，可得动态缩放系数为：比较式(12)和式(14)可知，受外力作用系统和受基础加速度作用系统，两者的动态缩放系数是一致的。将式(12)中的λ(t)对时间t求一次导数可得：由式(15)可知，当0≤t≤TD时，λ’(t)≥0，λ(t)为单调递增函数；当t>TD时，λ’(t)先大于零后谐频震荡，即λ(t)先增加后震荡减小。因此，λ(t)最大值发生在t>TD的第一个峰值时刻。求解λ(t)在t>TD时段的最大值，即可得动态缩放系数：令a＝TDωn，进一步将式λ(t)的下式展开整理得：式中：为相位，表达式如下：综合式(16)和式(17)可得，动态缩放系数为：根据式(3)，将式(19)化为工程上常用的固有频率fn的形式：可见，λ是与冲击力作用时间和结构固有频率乘积TDfn的函数。图3显示了根据式(20)得到的动态缩放系数λ随TDfn的变化曲线。可见：a)λ先随TDfn增大而增大，然后震荡衰减并逐渐趋近于1；b)TDfn在0～0.65范围内，λ一直随之增大，最大值约为1.26；c)当TDfn≤0.37时，λ值小于1。步骤二，建立确定冲击载荷等效静态的流程图4显示了确定冲击载荷等效静载的实施流程。对于不是严格后峰锯齿波形式的冲击载荷，优选地，采用GJB150A.18-2009附录C“冲击有效持续时间”规定的方法确定TD。具体方法如下：1.确定TETE定义为在冲击波形中，包括所有幅值超过最大峰值的1/3数据的最小时间长度。确定TE的示意图见图5。2.按下式计算TD：TD＝2.5TE可以采用振动模态试验或有限元仿真方法获取待冲击结构的固有频率，基频fn是指结构的第一阶固有频率。对于实际结构，优选地，采用振动模态试验确定fn。实施例一：以某梁式结构为例，说明基于本方法确定冲击载荷等效静载的实施过程。1.建立梁式结构的有限元模型图6显示了某梁式结构示意图。梁长2m，与OY轴负向呈45o夹角，横截面为等直圆环，圆环外径为0.05m，内径为0.04m。材料为合金钢，材料弹性模量为200GPa，泊松比为0.3，密度为7800kg/m3。在梁靠近原点O一端1/10处布置一线性弹簧，仅提供Y向刚度，用于调节梁绕OZ轴的弯曲频率。采用MSC.PATRAN建立该结构的有限元模型。2.计算有限元模型的固有频率分别设弹簧刚度值为1×105N/m，1×106N/m，1×107N/m，1×108N/m和1×109N/m。采用MSC.PATRAN中的SOL103实模态分析模块计算不同弹簧刚度下下的固有频率。表1显示这5个刚度值下，梁绕OZ轴的第一阶固有频率fn。表1不同刚度下的频率比较3.计算有限元模型的响应位移在梁自由端施加沿X轴正向的集中冲击力时间历程，力曲线如图7所示。可见，满足后峰锯齿波形式，冲击力持续时间TD＝0.05s，最大值P＝100N。利用MSC.PATRAN中的SOL112模态法计算上述5种不同刚度下自由端的位移时域响应，提取最大响应幅值，列于表2。表2不同刚度下的位移幅值比较4.基于有限元模型计算动态缩放系数利用MSC.本文档来自技高网...

【技术保护点】
一种确定冲击载荷等效静载的方法，其特征在于，所述的方法包括如下步骤：步骤一，根据以下公式计算动态缩放系数λ：

【技术特征摘要】
1.一种确定冲击载荷等效静载的方法，其特征在于，所述的方法包括如下步骤：步骤一，根据以下公式计算动态缩放系数λ：其中，TD为冲击载荷的持续时间；fn为待冲击结构的...

【专利技术属性】
技术研发人员：张玉杰，冯震宙，黄超广，
申请(专利权)人：中国航空工业集团公司西安飞机设计研究所，
类型：发明
国别省市：陕西,61