CN106055766B - 一种随机损伤对点阵夹层板动力学特性影响的分析方法 - Google Patents

一种随机损伤对点阵夹层板动力学特性影响的分析方法 Download PDF

Info

Publication number
CN106055766B
CN106055766B CN201610357186.5A CN201610357186A CN106055766B CN 106055766 B CN106055766 B CN 106055766B CN 201610357186 A CN201610357186 A CN 201610357186A CN 106055766 B CN106055766 B CN 106055766B
Authority
CN
China
Prior art keywords
dot matrix
sandwich plate
matrix sandwich
damage
plate model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610357186.5A
Other languages
English (en)
Other versions
CN106055766A (zh
Inventor
宋宏伟
路玲玲
黄晨光
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Mechanics of CAS
Original Assignee
Institute of Mechanics of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Mechanics of CAS filed Critical Institute of Mechanics of CAS
Priority to CN201610357186.5A priority Critical patent/CN106055766B/zh
Publication of CN106055766A publication Critical patent/CN106055766A/zh
Application granted granted Critical
Publication of CN106055766B publication Critical patent/CN106055766B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/24—Sheet material

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明提供了一种随机损伤对点阵夹层板动力学特性影响的分析方法，包括，首先获取该损伤程度下缺失胞元的数目；再利用软件随机产生缺失胞元的位置坐标，同时将点阵夹层板模型中对应位置坐标处的胞元删除，再对点阵夹层板模型进行模态分析，得到点阵夹层板模型的固有频率；再利用固有频率计算点阵夹层板模型在指定固有频率下的损伤评价结果；最后设置多组点阵夹层板模型的随机损伤工况，并重复上述过程以得到每一随机损伤工况的损伤评价结果，最终得到随机损伤对点阵夹层板模型指定固有频率的平均影响。本发明能够获得随机损伤在不同条件下对点阵夹层板动力学特性的影响规律，为后期点阵夹层板的损伤识别提供思路与依据。

Description

一种随机损伤对点阵夹层板动力学特性影响的分析方法
技术领域
本发明涉及材料领域，特别是涉及一种随机损伤对点阵夹层板动力学特性影响的分析方法。
背景技术
作为一种新型轻质多功能结构，点阵夹层板具有高强、隔热、屏蔽辐射等显著特点，被认为是最有发展潜力的防热形式之一，在热防护领域得到了越来越多的关注与应用。当应用于高速飞行器时，点阵夹层结构在极端的气动力/热载荷作用下，容易发生局部熔穿、结构屈曲等损伤。在试件制备过程中，点阵夹层板不可避免会出现各类制备缺陷，如虚焊、脱焊、芯材扭转、弯曲、断裂等。点阵夹层板损伤与缺陷的出现会引起夹层板的力学特性变化，如改变结构固有振动特性，容易引发结构共振；降低结构局部刚度与强度，产生局部薄弱部位；改变结构临界屈曲温度，增加结构屈曲失效风险等。因此，有必要开展损伤与缺陷对点阵夹层板动力学特性的影响分析，为点阵夹层板损伤识别提供思路与依据。
目前，已经开展了损伤与缺陷对点阵夹层板振动特性影响的理论、数值与实验研究，但在现有工作中，损伤与缺陷都为预先给定的。在实际中，由于点阵夹层板服役环境复杂多变以及各类突发事故，点阵夹层板的损伤往往也是随机分布，如结构遭受意外冲击载荷。同时，由于点阵夹层板工艺不成熟，随机缺陷也会在制备中出现，如虚焊、缺焊等。不同边界条件、损伤程度和位置下的随机损伤对点阵夹层板力学特性影响异常复杂。2001年，Wallach与Gibson研究了随机损伤对夹层板静力学特性的影响。他们发现模量与强度随着损伤增大而线性降低。但是，迄今为止，关于随机损伤对点阵夹层板动力学特性影响规律的研究未见报道。
发明内容
本发明的目的是要提供一种随机损伤对点阵夹层板动力学特性影响的分析方法。
特别地，本发明提供一种随机损伤对点阵夹层板动力学特性影响的分析方法，其特征在于，包括如下步骤：
步骤100，根据给定点阵夹层板模型的损伤程度和总胞元数目，得到该损伤程度下缺失胞元的数目；
步骤200，利用软件随机产生缺失胞元在所述点阵夹层板模型中的位置坐标，再将所述点阵夹层板模型中该位置坐标处的胞元删除，再对所述点阵夹层板模型进行模态分析，得到损伤后所述点阵夹层板模型的固有频率，并输出固有频率结果；
步骤300，通过固有频率结果，计算所述点阵夹层板模型在指定固有频率下的损伤评价结果；
步骤400，设置多组所述点阵夹层板模型的随机损伤工况，并重复上述步骤以得到每一随机损伤工况下所述点阵夹层板模型的损伤评价结果，最终得到多组随机损伤对所述点阵夹层板模型在指定固有频率下的平均影响。
进一步地，所述步骤100中，缺失胞元的数目Nmissing由下式得到：
Nmissing＝η×NTotal
η为给定的损伤程度，NTotal为模型总胞元数目。
进一步地，所述步骤300中所述点阵夹层板模型指定固有频率的损伤评价结果通过下式得到：
Δfi为固有频率的评价参数，分别为损伤与完好点阵夹层板模型的第i阶固有频率。
进一步地，所述步骤400中，多组随机损伤工况对所述点阵夹层板模型在指定固有频率下的平均影响通过下式得到：
为相同边界条件与损伤程度下第t组随机损伤工况引起的第i阶固有频率的相对变化，N为随机工况的总数。
进一步地，所述相同边界条件包括针对所述点阵夹层板模型进行一边固支三边自由、两边固支两边自由和四边固支三种情况。
进一步地，所述步骤200中，随机产生缺失胞元的位置坐标的软件是MATLAB。
进一步地，所述步骤200中，是利用ANSYS将所述点阵夹层板模型中由MATLAB随机产生位置坐标处的胞元删除，并对所述点阵夹层板模型进行模态分析的。
本发明能够获得随机损伤在不同条件下对点阵夹层板动力学特性的影响规律，为后期点阵夹层板的损伤识别提供思路与依据。
附图说明
图1是本发明一个实施例分析方法的流程示意图；
图2是本发明一个实施例中随机损伤对点阵夹层板前六阶中第一阶的模态影响示意图；
图3是本发明一个实施例中随机损伤对点阵夹层板前六阶中第二阶的模态影响示意图；
图4是本发明一个实施例中随机损伤对点阵夹层板前六阶中第三阶的模态影响示意图；
图5是本发明一个实施例中随机损伤对点阵夹层板前六阶中第四阶的模态影响示意图；
图6是本发明一个实施例中随机损伤对点阵夹层板前六阶中第五阶的模态影响示意图；
图7是本发明一个实施例中随机损伤对点阵夹层板前六阶中第六阶的模态影响示意图。
具体实施方式
如图1所示，本发明一个实施例的随机损伤对点阵夹层板动力学特性影响的分析方法一般性地包括如下步骤：
步骤100，根据给定点阵夹层板模型的损伤程度和总胞元数目，得到该损伤程度下缺失的胞元数目；
先选定一块点阵夹层板作为测试模型。缺失胞元的数目Nmissing由下式得到：
Nmissing＝η×NTotal
η为给定的损伤程度，NTotal为模型总胞元数目。
步骤200，利用软件随机产生缺失胞元在点阵夹层板模型中的位置坐标，再将点阵夹层板模型中该位置坐标处的胞元删除，再对点阵夹层板模型进行模态分析，得到损伤后点阵夹层板模型的固有频率，并输出固有频率结果；
通过随机的方式设定点阵夹层板损伤的位置，该损伤可以由相应位置处缺失的胞元来表现，随机产生缺失胞元的位置坐标的软件可以是MATLAB。将缺失胞元的位置坐标赋给点阵夹层板模型的相应位置坐标，对该点阵夹层板模型进行模态分析的软件可以是ANSYS。
通过模态分析即可得到该点阵夹层板模型在此随机损伤状态下的所有固有频率。
本步骤中虽然给出了MATLAB和ANSYS两个指定软件，但在其它的实施例中，也可以采用其它任意一个可以实现上述功能的软件。
步骤300，通过固有频率结果，计算点阵夹层板模型在指定固有频率下的损伤评价结果；
该点阵夹层板模型指定固有频率的损伤评价结果可以通过下式得到：
其中，Δfi为固有频率的评价参数，分别为损伤与完好点阵夹层板模型指定的第i阶固有频率。
步骤400，设置多组点阵夹层板模型的随机损伤工况，并重复上述步骤以得到每一随机损伤工况下点阵夹层板模型的损伤评价结果，最终得到多组随机损伤对点阵夹层板模型指定固有频率的平均影响。
多组随机损伤对点阵夹层板指定固有频率的平均影响通过下式评价：
为相同边界条件与损伤程度下第t组随机损伤工况引起的第i阶固有频率的相对变化，N为随机工况的总数。
这里的相同边界条件是指对点阵夹层板模型的四个边是否固定的情况，一般包括针对点阵夹层板模型进行一边固支三边自由、两边固支两边自由和四边固支这三种情况。本实施例中所采用的点阵夹层板模型可以是金字塔型金属点阵夹层板。
本发明能够获得随机损伤在不同条件下对点阵夹层板动力学特性的影响规律，为后期点阵夹层板的损伤识别提供思路与依据。
实施例1
选择金字塔型金属点阵夹层板模型。在算例中总胞元数目NTotal＝400，损伤程度η为0.01～0.1，则缺失胞元的数目Nmissing取值范围为4～40。随机损伤的Nmissing个胞元位置由MATLAB程序生成，随机分布于点阵夹层板结构中。
在计算中，共考虑三种边界条件：一边固支三边自由(CFFF)，两边固支两边自由(CCFF)，四边固支(CCCC)，得到前六阶固有频率在随机损伤下的相对变化，结果如图2-7所示。在这里，我们只涉及边界条件对点阵夹层板固有特性的影响。
通过附图2-7，可以看出，在相同损伤程度条件下，随着约束条件越严格，固有频率的变化也越大。这是因为，随着约束条件越严格，结构整体振动模态越复杂；同时，一些能够反映结构局部振动的模态也会出现。由于复杂的振动模态与局部振动模态都对结构损伤更加敏感。因此，在相同损伤程度下，点阵夹层板固有频率的变化与约束条件严格性成正比。同时，通过附图2-7可以看出，在相同边界条件下，随着损伤程度的增大，点阵夹层板结构的固有频率变化也增大。这是因为，缺失胞元的数目越多，点阵夹层板结构刚度变化越大，根据动力学方程，点阵夹层板固有频率变化也越大。
通过上述分析结果可以看出，本发明提出的分析方法的数值结果与分析结果一致，说明提出的这种随机损伤对点阵夹层板振动特性影响的分析方法是正确的。
至此，本领域技术人员应认识到，虽然本文已详尽示出和描述了本发明的多个示例性实施例，但是，在不脱离本发明精神和范围的情况下，仍可根据本发明公开的内容直接确定或推导出符合本发明原理的许多其它变型或修改。因此，本发明的范围应被理解和认定为覆盖了所有这些其它变型或修改。

Claims (6)

1.一种随机损伤对点阵夹层板动力学特性影响的分析方法，其特征在于，包括如下步骤：
步骤100，根据给定点阵夹层板模型的损伤程度和总胞元数目，得到该损伤程度下缺失胞元的数目；
步骤200，利用软件随机产生缺失胞元在所述点阵夹层板模型中的位置坐标，再将所述点阵夹层板模型中该位置坐标处的胞元删除，再对所述点阵夹层板模型进行模态分析，得到损伤后所述点阵夹层板模型的固有频率，并输出固有频率结果；
步骤300，通过固有频率结果，计算所述点阵夹层板模型在指定固有频率下的损伤评价结果；
所述点阵夹层板模型指定固有频率的损伤评价结果通过下式得到：
Δfi为固有频率的评价参数，fi D,fi I分别为损伤与完好点阵夹层板模型的第i阶固有频率；
步骤400，设置多组所述点阵夹层板模型的随机损伤工况，并重复上述步骤以得到每一随机损伤工况下所述点阵夹层板模型的损伤评价结果，最终得到多组随机损伤对所述点阵夹层板模型在指定固有频率下的平均影响。
2.根据权利要求1所述的分析方法，其特征在于，
所述步骤100中，缺失胞元的数目Nmissing由下式得到：
Nmissing＝η×NTotal
η为给定的损伤程度，NTotal为模型所有胞元数目。
3.根据权利要求1所述的分析方法，其特征在于，
所述步骤400中，多组随机损伤工况对所述点阵夹层板模型在指定固有频率下的平均影响通过下式得到：
Δfi t为相同边界条件与损伤程度下第t组随机损伤工况引起的第i阶固有频率的相对变化，N为随机工况的总数。
4.根据权利要求3所述的分析方法，其特征在于，
所述相同边界条件包括针对所述点阵夹层板模型进行一边固支三边自由、两边固支两边自由和四边固支三种情况。
5.根据权利要求1所述的分析方法，其特征在于，
所述步骤200中，随机产生缺失胞元的位置坐标的软件是MATLAB。
6.根据权利要求5所述的分析方法，其特征在于，
所述步骤200中，是利用ANSYS将所述点阵夹层板模型中由MATLAB随机产生位置坐标处的胞元删除，并对所述点阵夹层板模型进行模态分析的。
CN201610357186.5A 2016-05-26 2016-05-26 一种随机损伤对点阵夹层板动力学特性影响的分析方法 Active CN106055766B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610357186.5A CN106055766B (zh) 2016-05-26 2016-05-26 一种随机损伤对点阵夹层板动力学特性影响的分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610357186.5A CN106055766B (zh) 2016-05-26 2016-05-26 一种随机损伤对点阵夹层板动力学特性影响的分析方法

Publications (2)

Publication Number Publication Date
CN106055766A CN106055766A (zh) 2016-10-26
CN106055766B true CN106055766B (zh) 2019-05-10

Family

ID=57175169

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610357186.5A Active CN106055766B (zh) 2016-05-26 2016-05-26 一种随机损伤对点阵夹层板动力学特性影响的分析方法

Country Status (1)

Country Link
CN (1) CN106055766B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111931319A (zh) * 2020-07-13 2020-11-13 中国科学院力学研究所 一种非线性支承管束在横向流中的振动特性分析方法

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106503385B (zh) * 2016-11-07 2019-06-21 中国石油大学(华东) 一种点阵夹层材料等效弹性模量的计算方法
CN111537607B (zh) * 2020-05-07 2023-06-23 太原科技大学 基于庞加莱截面突变特性的超声导波缺陷检测方法
CN111964857B (zh) * 2020-07-27 2023-04-11 江苏大学 一种轻材点阵夹芯结构高频动态响应测试方法
CN112131718B (zh) * 2020-09-07 2022-06-10 中国科学院力学研究所 一种基于结构时域响应的点阵夹层板损伤识别方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103884776A (zh) * 2014-03-28 2014-06-25 大连理工大学 一种提高随机损伤定位向量法监测结果准确性的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6644535B2 (en) * 2001-05-18 2003-11-11 Massachusetts Institute Of Technology Truss core sandwich panels and methods for making same

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103884776A (zh) * 2014-03-28 2014-06-25 大连理工大学 一种提高随机损伤定位向量法监测结果准确性的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"Effect of local damages on the buckling behaviour of pyramidal truss core sandwich panels";Wu Yuan，et al.;《Composite Structures》;20160419;第271-277页 *
"基于振动特性的点阵桁架夹心板损伤检测数值模拟";田淑侠，等;《2009年远东无损检测论坛论文精选》;20091210;第31卷(第12期);第928-930页 *
"符合材料蜂窝夹层板低速冲击损伤模拟及实验研究";顾卫平 等;《机械科学与技术》;20130731;第32卷(第07期);第1017-1021页 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111931319A (zh) * 2020-07-13 2020-11-13 中国科学院力学研究所 一种非线性支承管束在横向流中的振动特性分析方法
CN111931319B (zh) * 2020-07-13 2021-05-07 中国科学院力学研究所 非线性支承管束在横向流中的振动特性分析方法

Also Published As

Publication number Publication date
CN106055766A (zh) 2016-10-26

Similar Documents

Publication Publication Date Title
CN106055766B (zh) 一种随机损伤对点阵夹层板动力学特性影响的分析方法
Donders et al. The effect of spot weld failure on dynamic vehicle performance
CN105426595A (zh) 一种铝合金热弹塑性变形模拟用本构模型的建立方法
Okita et al. A new concept of digital twin of artifact systems: synthesizing monitoring/inspections, physical/numerical models, and social system models
CN103049623A (zh) 一种激光焊接热源模型的建立方法
CN109684724B (zh) 一种客车多材料车身粘铆复合钢骨架结构分析方法
US20230038640A1 (en) Method for calculating service life of material under action of thermal shock load
P. Jirandehi et al. General quantification of fatigue damage with provision for microstructure: A review
Liu et al. Crack growth analysis and fatigue life estimation in the piston rod of a Kaplan hydro turbine
CN113987889A (zh) 一种焊接结构超高周疲劳主s-n曲线拟合及寿命预测方法
CN107727350B (zh) 微纳卫星矢量振动试验方法
CN104057210A (zh) 基于强度理论的焊接设计方法
CN106055760B (zh) 一种金属点阵夹层板的损伤识别方法
Ashwill et al. Concepts to facilitate very large blades
Mashayekhi Prediction of All-Steel CNG Cylinders Fracture in Impact by Using‎ Damage Mechanics Approach
CN115493954A (zh) 轨道车体复合材料疲劳耐久性及损伤容限强度验证方法
Youn et al. Structural durability design optimisation and its reliability assessment
Zade Finite Element Analysis and fatigue analysis of crane hook with different materials
Liu Acoustic Emission Signal Analysis and Damage Mode Identification of Composite Wind Turbine Blades
Amarir et al. Durability analysis of welded rectangular profiles under damped loads for automotive utilization
CN113408065B (zh) 带有指向性的减振装置用于随机振动仿真的等效建模方法
Chuanliang et al. Shape memory alloy hybrid composites for improving impact properties
CN106290011A (zh) 一种用于测试隔磁片受冲击过程力学响应的方法
Amarir et al. Finite element simulation of welded rectangular profiles for electric vehicles
Maican et al. Nonlinear impact resistance simulation of a newly designed roof from an agricultural tractor cabin.

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant