一种预估复材层合板鸟撞吸能强弱的方法
【技术领域】
[0001] 本发明涉及飞机鸟撞实验领域，特别是涉及一种预估复材层合板鸟撞吸能强弱的 方法。
【背景技术】
[0002] 吸能强弱是评价材料抗鸟撞优劣的重要指标之一。该指标通常由材料应力一应变 曲线所包围的面积，即应变能来表征。但应变能的计算需要事先测量应力一应变曲线，而不 同冲击应变率下的曲线数值并不一致，这对材料参数测试试验提出较高要求。
【发明内容】
[0003] 为了解决上述问题，本发明提供了一种预估复材层合板鸟撞吸能强弱的方法，该 方法基于应力-应变的线弹性假设，建立复材层合板应变能与等效弹性模量的关系。无 需测量层合板的应力-应变曲线，直接根据理论计算得到的弹性模量来预估其鸟撞吸能强 弱，主要包括如下步骤：
[0004]S1、获取预估复材层合板纵向弹性模量、横向弹性模量、泊松比、剪切模量；
[0005] S2、获取预估复材层合板的铺层角度及层数；
[0006]S3、计算预估复材层合板的弹性模量E。
[0007]S4、根据表征吸能强弱的应变能，当其处于线弹性变形段时，应变与应力满足：
[0008]
[0009] 其中，0为应力，e为应变，0 (e)为材料的应力一应变曲线，进而预估复材层合 板的鸟撞吸能强弱，即等效弹性模量越小，吸能越强。
[0010] 优选的是，所述表征吸能强弱的应变能W表达式如下：
[0011]
[0012] 其中，ea为材料的最大应变。
[0013] 在上述方案中优选的是，所述材料的应变能与弹性模量成反比，如下所示：
[0014]
[0015] 在上述方案中优选的是，所述弹性模量与单个铺层厚度成正比。
[0016] 该方法思路简单、明确，容易程序化，预估精度高，可以方便快捷地预估复合材料 层合板的鸟撞吸能强弱，对结构抗鸟撞设计具有指导作用。
【附图说明】
[0017] 图1是根据本发明一实施例的预估复材层合板鸟撞吸能强弱的方法流程图。
[0018] 图2是图1所示实施例的图鸟体动能随时间的变化曲线示意图。
【具体实施方式】
[0019] 为使本发明实施的目的、技术方案和优点更加清楚，下面将结合本发明实施例中 的附图，对本发明实施例中的技术方案进行更加详细的描述。在附图中，自始至终相同或类 似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本发明 一部分实施例，而不是全部的实施例。下面通过参考附图描述的实施例是示例性的，旨在用 于解释本发明，而不能理解为对本发明的限制。基于本发明中的实施例，本领域普通技术人 员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。下 面结合附图对本发明的实施例进行详细说明。
[0020] 在本发明的描述中，需要理解的是，术语"中心"、"纵向"、"横向"、"前"、"后"、"左"、 "右"、"竖直"、"水平"、"顶"、"底" "内"、"外"等指示的方位或位置关系为基于附图所示的方 位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元 件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明保护范围的 限制。
[0021] 如图1所示，为本发明预估复材层合板鸟撞吸能强弱的方法流程图，如图所示，本 发明涉及的方法包括：
[0022] S1、获取预估复材层合板纵向弹性模量、横向弹性模量、泊松比、剪切模量；
[0023]S2、获取预估复材层合板的铺层角度及层数；
[0024] S3、计算预估复材层合板的弹性模量E。
[0025]S4、根据表征吸能强弱的应变能，当其处于线弹性变形段时，应变与应力满足：
[0026]
[0027] 其中，〇为应力，e为应变，〇 (e)为材料的应力一应变曲线，进而预估复材层合 板的鸟撞吸能强弱，即等效弹性模量越小，吸能越强。
[0028] 需要指出的是，具体理论依据和评估方法如下：
[0029] 表征吸能强弱的应变能W表达式如下：
[0030]
[0031] 式中：〇为应力；e为应变；〇 (e)为材料的应力一应变曲线；ea为材料的最大 应变。
[0032] 在线弹性变形段，应变与应力满足：
[0033]
[0034] 式中：E为材料的弹性模量。
[0035] 这里认为线弹性变形段是影响材料吸能强弱的首要因素，因此将式（2)代入式 (1)，整理可得
[0036]
[0037] 由式⑶可知，当外应力一定时，材料的应变能与弹性模量成反比。对于复材层合 板这种各项异性材料，可以根据等效工程常数来评估其应变能大小。
[0038] 根据式（3)可以预估预估复材层合板的鸟撞吸能强弱，即等效弹性模量越小，吸 能越强。
[0039] 根据上述所示的吸能强弱预估理论，对3种碳纤维复材层合板的吸能性能进行分 析。其中，碳纤维工程常数见表1。
[0040] 表1碳纤维工程常数 [0041 ]
[0042] 表2显示了这3种层合板的铺层情况及其等效弹性模量计算结果。单个铺层厚度 为 0?125mm。
[0043] 表2碳纤维层合板等效弹性模量计算结果
[0044]
[0045] 根据表2的计算结果，可以初步判定：这3种碳纤维复材层合板的吸能强弱顺序 为，2#>1#>3#。
[0046] 以典型机翼前缘结构为仿真对象，对理论预估结果的准确性进行验证。该结构包 括蒙皮、隔板、角材和底板等。其中蒙皮为碳纤维复材层合板，分别选用表2所示的3种铺 层；其余结构为铝合金。材料的力学参数由试验测得，包含弹塑性变形段。结构模型均为壳 单元，鸟体模型采用SPH模型。
[0047] 采用PAM-CRASH2G-2006商用软件中的接触碰撞模块进行鸟撞仿真分析。仿真时 间为6ms。图2显示了鸟体动能随时间的变化曲线。可见，3种碳纤维复材层合板的吸能强 弱顺序为，2#>1#>3#。这和理论预估结果一致。
[0048] 作为本发明另一实施例，根据上述所示的吸能强弱预估理论，对3种GLARE(Glass ReinforcedAluminumLaminates)层合板的吸能性能进行分析。GLARE层板是由高强度错 合金薄板和高强度玻璃纤维增强树脂层交替层压构成的一种特殊复合材料。GLARE层合板 铺层由铝合金和玻璃纤维组成，相关工程常数见表3。
[0049] 表3铝合金和玻璃纤维工程常数
[0050]
[0051] 表4显示了这3种GLARE层合板的铺层情况及其等效弹性模量计算结果。
[0052] 表4GLARE层板弹性模量计算结果
[0053]
[0054] 根据表4的计算结果，可以初步判定：这3种GLARE层板的吸能强弱顺序为， 1#>2#>3#〇
[0055] 上述实施例中等效弹性模量的计算公式如下：
[0056] 复材二维单向层板通常包含4个独立的弹性常数：El、E2、vl2和G12。其二维刚 度矩阵中的参数可表示为
[0061] 其中：v21可根据下式求得：
[0062]
[0063] 对于含有多种铺层的复材层合板，其每层的刚度参数可表示为
[0070] 式中：9为铺层角。
[0071] 在得到每层的刚度参数之后，该层合板的刚度计算公式如下：
[0075] 式中：A。为拉伸刚度；B^为耦合刚度；D。为弯曲刚度；N为铺层总数；ZA单层法 向坐标。
[0076] 得到拉伸刚度Ax]之后，可根据下列式子求解单层结构的等效工程常数：
[0080] 式中：v为泊松比；E为弹性模量；G为剪切模量；t为层合板总厚度。
[0081] 最后需要指出的是：以上实施例仅用以说明本发明的技术方案，而非对其限制。 尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依 然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替 换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精 神和范围。
【主权项】
1. 一种预估复材层合板鸟撞吸能强弱的方法，其特征在于，包括： 51、 获取预估复材层合板纵向弹性模量、横向弹性模量、泊松比、剪切模量； 52、 获取预估复材层合板的铺层角度及层数； 53、 计算预估复材层合板的弹性模量E。 54、 根据表征吸能强弱的应变能，当其处于线弹性变形段时，应变与应力满足：其中，〇为应力，e为应变，〇 (e)为材料的应力一应变曲线，进而预估复材层合板的 鸟撞吸能强弱，即等效弹性模量越小，吸能越强。2. 如权利要求1所述的预估复材层合板鸟撞吸能强弱的方法，其特征在于：所述表征 吸能强弱的应变能W表达式如下：其中，ea为材料的最大应变。3. 如权利要求2所述的预估复材层合板鸟撞吸能强弱的方法，其特征在于：所述材料 的应变能与弹性模量成反比，如下所示：4. 如权利要求1所述的预估复材层合板鸟撞吸能强弱的方法，其特征在于：所述弹性 模量与单个铺层厚度成正比。
【专利摘要】本发明公开了一种预估复材层合板鸟撞吸能强弱的方法，属于飞机鸟撞实验领域。基于应力-应变的线弹性假设，建立复材层合板应变能与等效弹性模量的关系，用以预估复材层合板鸟撞吸能的强弱，本发明首先获取预估复材层合板纵向弹性模量、横向弹性模量、泊松比、剪切模量，并获取预估复材层合板的铺层角度及层数；进而计算预估复材层合板的弹性模量，从而预估复材层合板的鸟撞吸能强弱，即等效弹性模量越小，吸能越强。本发明无需测量层合板的应力-应变曲线，直接根据理论计算得到的弹性模量来预估其鸟撞吸能强弱。
【IPC分类】G06F17/50
【公开号】CN105022894
【申请号】CN201510496599
【发明人】张玉杰, 冯震宙, 黄超广
【申请人】中国航空工业集团公司西安飞机设计研究所
【公开日】2015年11月4日
【申请日】2015年8月13日