CN114428119B - 一种用于反演各向异性特征的复合材料弹性常数的方法 - Google Patents
一种用于反演各向异性特征的复合材料弹性常数的方法 Download PDFInfo
- Publication number
- CN114428119B CN114428119B CN202210064248.9A CN202210064248A CN114428119B CN 114428119 B CN114428119 B CN 114428119B CN 202210064248 A CN202210064248 A CN 202210064248A CN 114428119 B CN114428119 B CN 114428119B
- Authority
- CN
- China
- Prior art keywords
- mode
- guided wave
- elastic constant
- wave
- composite material
- 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
Links
- 239000002131 composite material Substances 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 33
- 239000006185 dispersion Substances 0.000 claims abstract description 34
- 238000012360 testing method Methods 0.000 claims abstract description 31
- 238000010586 diagram Methods 0.000 claims abstract description 16
- 238000001514 detection method Methods 0.000 claims abstract description 15
- 230000005284 excitation Effects 0.000 claims abstract description 11
- 235000019687 Lamb Nutrition 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 5
- 239000000463 material Substances 0.000 description 16
- 239000000523 sample Substances 0.000 description 12
- 238000005070 sampling Methods 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000006073 displacement reaction Methods 0.000 description 5
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000010206 sensitivity analysis Methods 0.000 description 3
- 229920000049 Carbon (fiber) Polymers 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 239000004917 carbon fiber Substances 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 239000007822 coupling agent Substances 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000011152 fibreglass Substances 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000009864 tensile test Methods 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- -1 chen Jiangyi Substances 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000004033 plastic Substances 0.000 description 1
- 229920003023 plastic Polymers 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- 238000009941 weaving Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/22—Details, e.g. general constructional or apparatus details
- G01N29/24—Probes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/011—Velocity or travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种用于反演各向异性特征的复合材料弹性常数的方法,它包括以下步骤:1、采用数值方法求解正交各向异性板中对称模态与反对称模态的频散方程,列出相应模态的频散曲线;2、通过分析弹性常数对频散曲线的影响,选择超声导波测试的激励频率;建立弹性常数与导波传播速度的对应关系图;3、利用超声导波检测方法测量对称模态与反对称模态在复合材料结构中的传播速度;4、利用弹性常数与导波传播速度的关系图反演得到对应的弹性常数。本发明的技术效果是:能对各向异性的复合材料进行无损测试,准确、快速地获得复合材料构件各主轴方向的弹性模量和剪切模量。
Description
技术领域
本发明属于复合材料力学性能测试领域,具体涉及一种利用超声导波检测确定复合材料弹性常数的方法。
背景技术
复合材料具有高比强度、高比模量和耐腐蚀、耐疲劳等优异性能,在航空航天、汽车船舶以及日常生活中得到了广泛的应用。在航空航天领域,复合材料构件尺寸较大,存在对这类大型构件进行弹性模量有关测试方法的特定需要。现有的超声测试方法采用体波或表面波进行测量,但对于复合材料而言,其弹性性能是各向异性的,无法通过体波法在大型构件上得到各主轴方向可靠的结果。
复合材料弹性模量的测量主要采用直条形结构件,在测试的长宽方向粘贴应变片,使用拉伸试验机进行计算确定。然而大型复合材料结构编织工艺复杂,且将其切割成直条形试件不仅成本昂贵,而且会造成很大浪费。
超声导波检测技术是无损检测方法之一,被广泛应用于金属、塑料、复合材料等材料的无损伤检测与物化特性分析。中国专利文献CN201611051177.X中公开了一种用于确定材料的弹性模量的空气耦合超声检测方法,包括以下步骤:1、根据对称模式和反对称模式的导波频散关系式,求解得到导波的第一理论频散曲线;2、利用第一理论频散曲线实验测得在给定的频厚积下待测板中的导波波长,得出导波的第一实测频散曲线;3、计算得出导波相速度随弹性模量的变化率和导波相速度随泊松比的变化率曲线;4、根据导波相速度变化率曲线筛选出多对计算点对;5、针对计算点对中的每一点,以其对应的弹性模量和泊松比重复步骤一、二以得到第二实测频散曲线;6、根据第二实测频散曲线利用相速度偏差的最小平方和优化算法求解待测材料的弹性模量。
该方法发明通过求解导波的理论频散曲线并与实测频散曲线进行比较分析得到待测材料的弹性模量,但该方法测量与计算过程复杂,且该方法基于各向同性材料(如:金属、树脂等均质材料)中的导波理论,能够测出其弹性模量与泊松比,但对于复合材料这类具有各向异性特征(至少有5个独立弹性常数)的材料与结构并不适用。
发明内容
针对现有技术存在的问题,本发明所要解决的技术问题就是提供一种用于反演各向异性特征的复合材料弹性常数的方法,它能对各向异性的复合材料进行无损测试,准确、快速地获得复合材料构件各主轴方向的弹性模量和面外的剪切模量,为航空航天领域大型复合材料结构设计和强度核算提供基础参数。
本发明所要解决的技术问题是通过这样的技术方案实现的,它包括以下步骤:
步骤1、采用数值方法求解正交各向异性板中对称模态与反对称模态的频散方程,列出相应模态的频散曲线;
步骤2、通过分析弹性常数对频散曲线的影响,选择超声导波测试的激励频率;建立弹性常数与导波传播速度的对应关系图;
步骤3、利用超声导波检测方法测量对称模态与反对称模态在复合材料结构中的传播速度;
步骤4、利用弹性常数与导波传播速度的对应关系图反演得到对应的弹性常数。
由于本方法基于各向异性材料中的导波理论,通过测试最低阶的对称模态(S0)与反对称模态(A0)波速,利用S0的波速在弹性常数与导波传播速度对应的关系图中对比查找反演得到主轴方向弹性模量(E11、E22);利用A0的波速在弹性常数与导波传播速度对应的关系图中对比查找得到面外剪切模量(G13、G23)。与背景技术中CN201611051177.X相比,本方法发明有更广泛的适用性,能测试各向异性的复合材料。
与现有弹性常数测试方法相比,本发明具有以下优点:
本发明在不对结构造成任何破坏的情况下,能测试各向异性材料的弹性常数;
本发明可以用测试板结构、梁结构、大尺寸薄壁结构以及变厚度结构等复杂成型件的等效弹性常数;
本发明操作简便,对测试环境要求不高,可用于工作现场测试,测试效率高,测试结果准确。
附图说明
本发明的附图说明如下:
图1为实施例的相速度频散曲线;
图2为实施例的群速度频散曲线;
图3为实施例的对称模态传播速度与弹性模量关系图;
图4为实施例的反对称模态传播速度与剪切模量关系图;
图5为实施例超声导波方法测试波速的实验系统示意图;
图6为复合材料板激励与采样点布置图;
图7为实施例一个接收点所采集的时域信号图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明:
本方法发明包括以下步骤:
步骤1、采用数值方法求解正交各向异性板中对称模态与反对称模态的频散方程,列出相应模态的频散曲线
根据文献“正交各向异性板中的兰姆波频散曲线”,蔡应龙,铁瑛,陈江义等,玻璃钢/复合材料,2018(10):4记载,正交各向异性材料中声波控制方程为:
式(1)、(2)中,ρ为密度,u1、u3是板中质点沿x轴与z轴的位移矢量分量,u1,11、u3,13、u3,33、u1,13表示位移矢量对坐标的微分,其中下标“,”代表微分符号,c11、c13、c33、c55是板的弹性常数,t为时间。
假设其位移解为:
u1=Aei(kx+kλz-ωt) (3)
u3=Bei(kx+kλz-ωt) (4)
式(3)、(4)中,k为波在x方向的波数,λ为待求系数,A和B为位移幅值,ω为波的圆频率,i是虚数单位。
将式(3)、(4)代入式(1)、(2)可以得到特征方程:
要使式(5)的解有意义,则系数矩阵的行列式应为0,即:
整理得,
式(7)中,λp,q的下标“p”与“q”分别对应下公式右端的“+”和“-”,
上式中,c=ω/k表示波的相速度。
通过整理并化简可得到正交各向异性板的频散方程为:
其中
式(10)中,h表示板厚,幂指数项“+1”和“-1”分表表示Lamb波的对称模态与反对称模态。
通过频散方程(10)可以发现:频散方程包含两个基本参数,波数k与相速度c,使用MATLAB确定其中一个参数值进行数值求解得到另一个参数,可计算得到对应材料的对称模态与反对称模态的频散曲线。
步骤2、通过分析弹性常数对频散曲线的影响,选择超声导波测试的激励频率;建立弹性常数与导波传播速度的对应关系图;
1、弹性常数敏感性分析
对于正交各向异性材料中的九个独立的弹性常数(E11、E22、E33、G12、G13、G23、ν12、ν13、ν23),使用控制变量法,每次只改变其中一个参数,其余八个参数保持不变,分析得到九个弹性常数与对称模态与反对称模态的相关性。通过分析发现:S0在低频段的传播速度仅与对应主轴方向弹性模量(E11、E22)主要相关,而A0在低频段的传播速度与对应的剪切模量(G13、G23)主要相关。
然后依据分析结果选取测量弹性常数所需的模态与激励频率范围。
2、建立导波传播速度与弹性常数间的对应关系图
用数值求解方法得到对称模态、反对称模态在材料中的传播速度与弹性模量、剪切模量一一对应关系图。
步骤3、利用超声导波检测方法测量对称模态与反对称模态在复合材料结构中的传播速度;
根据被测复合材料的基本参数与步骤2的弹性常数敏感性分析结果,使用超声检测方法测波速,选取合适的超声信号发射探头在结构中激励Lamb波信号,并在沿复合材料主轴方向选取等间距的接收点拾取信号,进而根据测得的时域信号进行信号处理分析计算得到所测结构中Lamb波传播的相速度。
对于碳纤维增强复合材料,当测得的信号(频厚积小于1.5MHz.mm时)传播速度小于3500m/s时为A0模态的速度,当测得的波速大于4000m/s时,接收到的信号则为S0模态。
与常规超声检测相比,Lamb波检测具有快速高效的特点,适合于板形结构的物化特性测试分析,因此本步骤采用Lamb波测试复合材料的弹性常数。
步骤4、利用弹性常数与导波传播速度之间的对应关系图反演得到对应的弹性常数
根据步骤3测得的对称模态或反对称模态传播速度,在步骤2得到的波速与弹性常数对应关系图中进行查找,即可反演得到被测材料的弹性常数。
实施例
一种0、90°铺层的复合材料:铺层顺序[0/90/0/90/0/90],厚度4.4mm;固化成型后的CFRP板密度为1400kg/m3。该复合材料单层板参数(厂家提供)为:C11=119GPa,C12=6.2GPa,C13=6.2GPa,C22=13.8GPa,C23=4.5GPa,C33=13.8GPa,C44=(C22-C23)/2,C55=5.7GPa,C66=5.7Gpa。
本实施例的反演过程如下:
步骤1、采用数值方法求解正交各向异性板中对称模态与反对称模态的频散方程,列出相应模态的频散曲线
文献“正交各向异性板中的兰姆波频散曲线”,蔡应龙,铁瑛,陈江义等,玻璃钢/复合材料,2018(10):4记载的二分数值方法计算得到相速度频散曲线参见图1、群速度频散曲线参见图2。图中的A0是Lamb波最低阶的反对称模态,而A1、A2、A3则是对应的高阶反对称模态;图中S0是Lamb波最低阶对称模态,S1、S2、S3是对应的高阶对称模态。
步骤2、通过分析弹性常数改变对频散曲线的影响,选择超声导波测试的激励频率;建立弹性常数与导波传播速度的对应关系图
由图2可知,除对称模态S0和反对称模态A0外,其他导波模态都存在截止频率(如A1的截止频率约为450kHz,S1的截止频率约为720kHz等),因此为减少多模态对检测结果的影响,因此本实施例选择小于截止频率的频率段进行波速测试。
对碳纤维增强复合材料(CFRP)中导波传播敏感性分析发现:对称模态S0在低频范围对Lamb波传播方向对弹性模量E敏感,对其他弹性常数不敏感;反对称模态A0在低频范围对Lamb波传播方向对剪切模量G敏感,而对其他弹性常数不敏感。
S0模态的传播速度在低频段仅与传播方向主轴的弹性模量相关,反演复合材料主轴方向弹性模量E,因此可以通过测量S0模态沿主轴方向的传播速度,可反演得到复合材料的两个主要弹性常数E11与E22。
A0模态在低频区域的传播速度仅与沿传播方向对应的剪切模量相关,而与其他弹性常数无关,因此可以通过测量A0模态在低频段的传播速度反演得到两个主要弹性常数G13与G23。
根据前述数值计算方法得到的结果建立导波传播速度与弹性常数间的关系图:
图3为S0模态传播速度与弹性模量E之间的关系图,
图4为A0模态传播速度与剪切模量G之间的关系图。
本实施例仅对弹性模量进行测试与验证。
步骤3、利用超声导波检测方法测量对称模态和反对称模态在复合材料结构中的传播速度;
用图5所示的超声检测装置测试复合材料中S0模态传播速度,该超声检测装置包括:信号发生器、功率放大器、超声探头、AE传感器、示波器和计算机,其中,超声探头和AE传感器接在复合材料上,用于时域信号的激励与接收。
测试过程:
实验器材:普源DG4062信号发生器,用于发射超声导波信号;
安泰ATA-300功率放大器,用于激励信号的放大;
RITEC有源放大器,用于接收信号的放大;
泰克DPO 3014示波器,用于采集超声导波信号;
计算机,用于操控LabVIEW进行信号采集与存储;
超声探头,中心频率100kHz,可用于发射或接收超声导波信号,本实施例中用于发射超声导波信号;
AE传感器,最佳工作频率100kHz~1MHz,可用于接收超声导波信号,主要用于接收结构中的离面位移信号;
横波耦合剂与纵波耦合剂,用于探头与被测物体之间进行声耦合;
测试方案:首先进行设备连接,信号发生器(DG4062)激励信号,经后置功率放大器(ATA300)进行信号放大,然后传输至超声探头,将探头放置在图6所示激励源位置,为减小探头与被测材料间的声阻抗差,在探头接触面涂抹均匀的耦合剂,并使用定制夹具进行固定,探头发射超声波并传入被测材料中,与结构耦合形成超声导波信号;然后在图6所示的采样点1~10进行信号采集,将AE传感器接触面均匀涂抹耦合剂,放置在采样点位置,并使用定制夹具进行固定,AE传感器采集的信号经前置放大器进行滤波与放大后将数据传输至示波器(DPO3014),将示波器通过网线与计算机连接,在计算机上进行信号采集与存储。测试过程的具体步骤如下:
根据前述参数敏感性分析,即对称模态在低频范围对Lamb波传播方向对弹性模量E敏感,本实施例的工作频率选择80kHz;
然后在信号发生器上发射中心频率为80kHz,周期数为5的Hanning窗调制正弦信号至奥林巴斯探头;
将AE传感器放置在采样点1位置进行信号拾取;
在计算机端操作示波器配套LabVIEW进行激励信号与接收信号的采集,采样点设置为100k,并将信号进行编号存储;
数据存储完成后,将AE传感器分别移至2~10号采样点位置,并重复上述步骤进行信号采集与存储,直至采样完成。
图7所示为某一采样点经MATLAB降噪处理后得到的时域信号;
接下来计算导波传播速度,
通过算法提取第一个波包到达时间ti,如图7所示;
通过波速计算公式Vi=Δdi/Δti,并对Vi取平均值得到超声导波S0模态在CFRP板中1方向(如图6所示方向)的传播速度V1=6955m/s,2方向的传播速度为V2=6498m/s。
步骤4、利用弹性常数与导波传播速度的关系图反演得到对应的弹性常数
通过图3所示传播速度与弹性常数关系图反演得到所测波速VS0所对应的弹性模量为E11=67.7GPa,E22=59.4GPa;
为验证测试结果的准确性,在速度测试完成后,将CFRP板切割成一组标准拉伸试验样件进行弹性模量测试。
测得1方向弹性模量为70.8GPa,2方向弹性模量为63.0GPa,与本发明方法测试结果误差小于6.0%,说明本发明测试准确有效。
Claims (1)
1.一种用于反演各向异性特征的复合材料弹性常数的方法,其特征是:包括以下步骤:
步骤1、采用数值方法求解正交各向异性板中对称模态与反对称模态的频散方程,列出相应模态的频散曲线;
所述正交各向异性板中对称模态与反对称模态的频散方程为:
式中,
h表示板厚,幂指数项“+1”和“-1”分表表示Lamb波的对称模态与反对称模态;
λp,q的下标“p”与“q”分别对应下公式右端的“+”和“-”;
其中,ρ为密度,c=ω/k表示波的相速度;k为波在x方向的波数,ω为波的圆频率,c11、c13、c33、c55是板的弹性常数;
步骤2、通过分析弹性常数对频散曲线的影响,得到超声导波传播速度与弹性模量关系为:对称模态在低频段的传播速度与对应主轴方向弹性模量相关,非对称模态在低频段的传播速度与对应的剪切模量相关;
选择超声导波测试的激励频率,建立弹性常数与导波传播速度的对应关系图;所述超声导波测试的激励频率选用小于导波高阶模态的截止频率;
步骤3、利用超声导波检测方法测量对称模态与反对称模态在复合材料结构中的传播速度,所述超声导波为Lamb波;
步骤4、利用弹性常数与导波传播速度对应关系图反演得到对应的弹性常数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210064248.9A CN114428119B (zh) | 2022-01-20 | 2022-01-20 | 一种用于反演各向异性特征的复合材料弹性常数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210064248.9A CN114428119B (zh) | 2022-01-20 | 2022-01-20 | 一种用于反演各向异性特征的复合材料弹性常数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114428119A CN114428119A (zh) | 2022-05-03 |
CN114428119B true CN114428119B (zh) | 2023-11-21 |
Family
ID=81312711
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210064248.9A Active CN114428119B (zh) | 2022-01-20 | 2022-01-20 | 一种用于反演各向异性特征的复合材料弹性常数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114428119B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116577417B (zh) * | 2023-07-13 | 2023-09-29 | 浙江大学 | 一种用于复合材料的自适应超声全聚焦缺陷成像方法 |
CN117214303B (zh) * | 2023-11-07 | 2024-01-30 | 华东交通大学 | 一种固体材料弹性模量的导波无损检测方法 |
Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1054763A1 (ru) * | 1982-06-16 | 1983-11-15 | Ленинградский Ордена Ленина Электротехнический Институт Им.В.И.Ульянова (Ленина) | Способ контрол упругих посто нных ограниченных пр моугольных пластин |
JPH05133861A (ja) * | 1991-11-12 | 1993-05-28 | Nkk Corp | 被検体の弾性定数測定法 |
JPH08193984A (ja) * | 1995-01-18 | 1996-07-30 | Toyota Autom Loom Works Ltd | 素材の異方性評価方法及びその装置 |
JP2004266878A (ja) * | 2004-06-28 | 2004-09-24 | Toyo Commun Equip Co Ltd | 三重モード圧電フィルタの周波数調整方法 |
CN104330300A (zh) * | 2014-11-24 | 2015-02-04 | 重庆大学 | 超高温陶瓷材料热-损伤耦合强度间接测量方法 |
CN104807570A (zh) * | 2015-05-04 | 2015-07-29 | 北京化工大学 | 一种超声兰姆波测量塑料薄板制品内应力的装置和方法 |
CN105808792A (zh) * | 2014-12-29 | 2016-07-27 | 北京强度环境研究所 | 一种贮箱晃动质量的数值计算方法 |
CN106055853A (zh) * | 2016-07-29 | 2016-10-26 | 四川大学 | 无初值依赖和角度依赖的正交各向异性弹性系数反演方法 |
CN106680375A (zh) * | 2016-11-25 | 2017-05-17 | 中国商用飞机有限责任公司 | 用于确定材料的弹性模量的空气耦合超声检测方法 |
CN106770648A (zh) * | 2016-11-21 | 2017-05-31 | 南京航空航天大学 | 一种复合材料刚度系数的无损测量方法 |
CN108426948A (zh) * | 2018-05-14 | 2018-08-21 | 南京航空航天大学 | 一种激发单一模态Lamb波的电磁超声换能器及其工作方法 |
CN110261485A (zh) * | 2019-07-26 | 2019-09-20 | 北京领示科技有限公司 | 一种超声波测量材料内部各处弹性模量及泊松比的方法 |
DE102018003311A1 (de) * | 2018-04-24 | 2019-10-24 | Diehl Metering Gmbh | Verfahren und Messeinrichtung zur Ermittlung einer Messinformation |
US10458958B1 (en) * | 2017-01-24 | 2019-10-29 | United States Of America As Represented By The Secretary Of The Air Force | Ultrasonic through-thickness modulus evaluation of materials |
WO2021010836A1 (en) * | 2019-07-18 | 2021-01-21 | Technische Universiteit Delft | Method and system for using wave analysis for speed of sound measurement |
CN112362745A (zh) * | 2020-11-02 | 2021-02-12 | 大连理工大学 | 一种均匀基底涂覆层弹性常数的超声测量与反演方法 |
CN112444563A (zh) * | 2020-11-25 | 2021-03-05 | 大连理工大学 | 一种基于超声背反射的横观各向同性材料损伤评价方法 |
CN112557220A (zh) * | 2020-11-12 | 2021-03-26 | 西北工业大学 | 高速冲击异形件正交异性弹性常数虚场同步表征法 |
CN113109447A (zh) * | 2021-04-20 | 2021-07-13 | 哈尔滨工业大学(深圳) | 一种用于复合材料的动态均质化超声全聚焦缺陷成像方法及系统 |
CN113640883A (zh) * | 2021-09-16 | 2021-11-12 | 中国石油大学(北京) | 测井速度预测、基于井震融合的孔隙度确定方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5629481B2 (ja) * | 2010-03-16 | 2014-11-19 | 富士重工業株式会社 | 損傷診断システム |
AR113330A1 (es) * | 2017-10-31 | 2020-04-08 | Univ La Republica | Equipo y método para determinar la elasticidad de sólidos blandos |
-
2022
- 2022-01-20 CN CN202210064248.9A patent/CN114428119B/zh active Active
Patent Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1054763A1 (ru) * | 1982-06-16 | 1983-11-15 | Ленинградский Ордена Ленина Электротехнический Институт Им.В.И.Ульянова (Ленина) | Способ контрол упругих посто нных ограниченных пр моугольных пластин |
JPH05133861A (ja) * | 1991-11-12 | 1993-05-28 | Nkk Corp | 被検体の弾性定数測定法 |
JPH08193984A (ja) * | 1995-01-18 | 1996-07-30 | Toyota Autom Loom Works Ltd | 素材の異方性評価方法及びその装置 |
JP2004266878A (ja) * | 2004-06-28 | 2004-09-24 | Toyo Commun Equip Co Ltd | 三重モード圧電フィルタの周波数調整方法 |
CN104330300A (zh) * | 2014-11-24 | 2015-02-04 | 重庆大学 | 超高温陶瓷材料热-损伤耦合强度间接测量方法 |
CN105808792A (zh) * | 2014-12-29 | 2016-07-27 | 北京强度环境研究所 | 一种贮箱晃动质量的数值计算方法 |
CN104807570A (zh) * | 2015-05-04 | 2015-07-29 | 北京化工大学 | 一种超声兰姆波测量塑料薄板制品内应力的装置和方法 |
CN106055853A (zh) * | 2016-07-29 | 2016-10-26 | 四川大学 | 无初值依赖和角度依赖的正交各向异性弹性系数反演方法 |
CN106770648A (zh) * | 2016-11-21 | 2017-05-31 | 南京航空航天大学 | 一种复合材料刚度系数的无损测量方法 |
CN106680375A (zh) * | 2016-11-25 | 2017-05-17 | 中国商用飞机有限责任公司 | 用于确定材料的弹性模量的空气耦合超声检测方法 |
US10458958B1 (en) * | 2017-01-24 | 2019-10-29 | United States Of America As Represented By The Secretary Of The Air Force | Ultrasonic through-thickness modulus evaluation of materials |
DE102018003311A1 (de) * | 2018-04-24 | 2019-10-24 | Diehl Metering Gmbh | Verfahren und Messeinrichtung zur Ermittlung einer Messinformation |
CN108426948A (zh) * | 2018-05-14 | 2018-08-21 | 南京航空航天大学 | 一种激发单一模态Lamb波的电磁超声换能器及其工作方法 |
WO2021010836A1 (en) * | 2019-07-18 | 2021-01-21 | Technische Universiteit Delft | Method and system for using wave analysis for speed of sound measurement |
CN110261485A (zh) * | 2019-07-26 | 2019-09-20 | 北京领示科技有限公司 | 一种超声波测量材料内部各处弹性模量及泊松比的方法 |
CN112362745A (zh) * | 2020-11-02 | 2021-02-12 | 大连理工大学 | 一种均匀基底涂覆层弹性常数的超声测量与反演方法 |
CN112557220A (zh) * | 2020-11-12 | 2021-03-26 | 西北工业大学 | 高速冲击异形件正交异性弹性常数虚场同步表征法 |
CN112444563A (zh) * | 2020-11-25 | 2021-03-05 | 大连理工大学 | 一种基于超声背反射的横观各向同性材料损伤评价方法 |
CN113109447A (zh) * | 2021-04-20 | 2021-07-13 | 哈尔滨工业大学(深圳) | 一种用于复合材料的动态均质化超声全聚焦缺陷成像方法及系统 |
CN113640883A (zh) * | 2021-09-16 | 2021-11-12 | 中国石油大学(北京) | 测井速度预测、基于井震融合的孔隙度确定方法及装置 |
Non-Patent Citations (2)
Title |
---|
正交各向异性板中的兰姆波频散曲线;蔡应龙,铁瑛,陈江义;《复合材料科学与工程》(第10期);第47-50页 * |
薄板结构中非线性低频超声Lamb波传播机理研究;丁湘燕,赵友选,胡 宁;《声 学 技 术》;第36卷(第5期);第455-466页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114428119A (zh) | 2022-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114428119B (zh) | 一种用于反演各向异性特征的复合材料弹性常数的方法 | |
US5408882A (en) | Ultrasonic device and method for non-destructive evaluation of polymer composites | |
CN114235241B (zh) | 基于莫尔圆应力原理的非线性超声应力检测方法及装置、复合材料平面应力检测方法 | |
CN106770648B (zh) | 一种复合材料刚度系数的无损测量方法 | |
CN113325075A (zh) | 一种用于金属薄板高周疲劳损伤的非线性波检测方法 | |
CN103926315B (zh) | 一种基于单纯形法的各向同性薄板材料弹性性质获取方法 | |
CN115683369A (zh) | 一种基于Lamb波双模态声时比的复合材料空耦超声平面应力检测方法 | |
CN105424810A (zh) | 纤维增韧陶瓷基复合材料均匀性评价方法 | |
CN114964580A (zh) | 一种基于空气耦合Lamb波的正交各向异性复合材料平面应力检测方法及可读存储介质 | |
CN112326786A (zh) | 基于电磁超声Lamb波S1模态群速度的金属板应力检测方法 | |
CN114994177B (zh) | 复合板材超声缺陷检测方法及装置和复合板材 | |
CN115166051A (zh) | 一种纤维增强树脂基复合材料弹性常数的测量方法 | |
CN109490417B (zh) | 一种金属材料平面各向异性超声检测方法 | |
CN106680375B (zh) | 用于确定材料的弹性模量的空气耦合超声检测方法 | |
Lobkis et al. | In-plane elastic property characterization in composite plates | |
Grondel et al. | The propagation of Lamb waves in multilayered plates: phase-velocity measurement | |
CN115753991B (zh) | 一种空耦超声高空间分辨率应力场测量装置的使用方法 | |
Hafizi et al. | An investigation of acoustic emission signal attenuation for monitoring of progressive failure in fiberglass reinforced composite laminates | |
Cha et al. | Compression wave velocity of cylindrical rock specimens: engineering modulus interpretation | |
Spalvier et al. | The stress-induced surface wave velocity variations in concrete | |
Valdes et al. | Delamination detection in composite laminates from variations of their modal characteristics | |
CN112945771A (zh) | 一种超声疲劳试验试样应力应变校准方法 | |
CN117214303B (zh) | 一种固体材料弹性模量的导波无损检测方法 | |
Ratnam et al. | Structural damage detection using combined finite-element and model Lamb wave propagation parameters | |
CN116358747B (zh) | 一种基于双晶空耦换能器的高空间分辨率平面应力测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |