CN105181480A - 一种多点弯曲测量复合材料面外刚度的方法 - Google Patents
一种多点弯曲测量复合材料面外刚度的方法 Download PDFInfo
- Publication number
- CN105181480A CN105181480A CN201510494369.7A CN201510494369A CN105181480A CN 105181480 A CN105181480 A CN 105181480A CN 201510494369 A CN201510494369 A CN 201510494369A CN 105181480 A CN105181480 A CN 105181480A
- Authority
- CN
- China
- Prior art keywords
- plane stiffness
- point bending
- testpieces
- compound substance
- point
- 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.)
- Granted
Links
Abstract
一种多点弯曲测量复合材料面外刚度的方法,属于航空结构力学领域。利用多点弯曲特别是五点弯曲这一简便的试验方法形成的既有存剪状态又有耦合状态的在一个试验件上形成的多种应力场分布,同时利用参数反演方法根据有限元模型结果进行优化反演计算可以一次得到多个面外刚度数据。本发明提供的多点弯曲测量复合材料面外刚度的方法测量结果准确度高,消除了各刚度之间的耦合效应。采用本发明提供的多点弯曲测量复合材料面外刚度的方法所得到的复合材料面外刚度数据结果用于有限元分析精度高,试验件制备简单,试验方法简便,提高了数据的精度,缩短了设计周期。
Description
技术领域
本发明涉及航空结构力学技术领域,具体而言,涉及一种多点弯曲测量复合材料面外刚度的方法。
背景技术
材料刚度是材料受外力作用时抵抗变形的能力,一般用弹性模量表示。复合材料是各向异性材料,不同材料方向上的刚度各不相同。根据复合材料力学理论,用工程弹性常数表示的柔度矩阵表示如下:
用9个独立弹性常数可以表征材料属性,即三个材料主方向上的弹性模量E1,E2,E3,三个泊松比γ12,γ13,γ23,三个平面内的剪切弹性模量G12,G13,G23。其中把3方向的弹性模量E3和G13、G23平面内的剪切弹性模量G13、G23称为面外刚度。在进行厚板和编织类的复合材料性能评估时需要使用准确的面外刚度数据进行基于三维的强度分析和失效分析。
现有的复合材料面外刚度(E3,G23,G13)测量方法仅仅针对单独项目,而且很难实现单独应力场,试验件制备复杂,要求高。具体存在的问题包括:面外拉压模量E3的测量目前还没有正式的标准,测量的方法和测量试验件尺寸尚在进一步研究中;面外剪切模量G23,G13虽然有ASTM正式标准,但测试方法不能实现纯剪力场,有其它载荷耦合的影响。试验件制备复杂,往往容易受试验件加工质量影响引起非正常破坏。目前五点弯曲方法有用来进行面内剪切模量G12的测量的专利,还没有进行复合材料面外刚度测量的研究和专利的公开。
现在亟需解决的技术问题是如何设计一种复合材料面外刚度的测量方法,该测量方法能够消除各刚度耦合效应,解决现有技术中在测量复合材料面外刚度中存在的缺陷。
发明内容
本发明的目的在于解决上述现有技术中的不足,提供一种能够消除各刚度耦合效应、一次性可得到多个复合材料面外刚度数据的多点弯曲测量复合材料面外刚度的方法。
本发明的目的通过如下技术方案实现:一种多点弯曲测量复合材料面外刚度的方法,包括如下步骤:
S1:实验准备,选用一般的万能材料试验机作为试验机,制作复合材料试验件,确定对试验件进行弯曲的点数并准备试验夹具;
S2:按照S1中选定的多点弯曲方法进行试验和测量并得出相应的数据;
S3:对S2中得到的数据进行参数反演,得出复合材料面外刚度。
上述方案中优选的是,S1中选取的试验机的最小的工作载荷为100KN。将试验机的工作载荷设定较大的数量级,保障测量条件能够满足多种测量需求,以适应较大范围数量级的测量。
上述任一方案中优选的是,S1中制作的复合材料试验件为规则的长方体结构,并以试验件的长、宽、高建立坐标系,定义:对试验件进行弯曲的点数为N。
上述任一方案中优选的是,定义:S1中试验件的长度取值为L,其中,L=45*(N-2)。
上述任一方案中优选的是,S2中在试验机上居中放置试验件,确定以试验件的高度方向为试验机加载方向,调整每组加载点和支持点之间的距离,对试验件施加压缩载荷,对试验件每个加载点截面上沿试验件高度方向设置测量点,测量点的数量大于等于三个。
上述任一方案中优选的是,通过数据处理软件进行参数反演计算。
上述任一方案中优选的是,通过载荷及对应测量的点的位移,计算得出对应的应变下面外刚度数据。
上述任一方案中优选的是,通过多个载荷和对应测量的的位移数据可得出应变-模量曲线。
本发明所提供的多点弯曲测量复合材料面外刚度的方法的有益效果在于,测量结果准确度高,消除了各刚度之间的耦合效应。采用本发明提供的多点弯曲测量复合材料面外刚度的方法所得到的复合材料面外刚度数据结果用于有限元分析精度高,试验件制备简单,试验方法简便。
附图说明
图1是按照本发明的多点弯曲测量复合材料面外刚度的方法的优选实施例流程示意图;
图2是按照本发明的多点弯曲测量复合材料面外刚度的方法的图1所示实施例的选用的试验件的示意图;
图3是按照本发明的多点弯曲测量复合材料面外刚度的方法的图1所示实施例S3的流程示意图;
图4是是按照本发明的多点弯曲测量复合材料面外刚度的方法的图1所示实施例迭代过程曲线图。
具体实施方式
为了更好地理解按照本发明方案的多点弯曲测量复合材料面外刚度的方法,下面结合附图对本发明的多点弯曲测量复合材料面外刚度的方法的一优选实施例作进一步阐述说明。
结合图1-图4,本发明提供的多点弯曲测量复合材料面外刚度的方法包括如下步骤:
S1:实验准备,选用一般的万能材料试验机作为试验机,制作复合材料试验件,确定对试验件进行弯曲的点数并准备试验夹具;
S2:按照S1中选定的多点弯曲方法进行试验和测量并得出相应的数据;
S3:对S2中得到的数据进行参数反演,得出复合材料面外刚度。
S1中选取的试验机的最小的工作载荷为100KN。S1中制作的复合材料试验件为规则的长方体结构,并以试验件的长、宽、高建立坐标系,定义:对试验件进行弯曲的点数为N。定义:S1中试验件的长度取值为L,其中,L=45*(N-2)。
S2中在试验机上居中放置试验件,确定以试验件的高度方向为试验机加载方向,调整每组加载点和支持点之间的距离,对试验件施加压缩载荷,对试验件每个加载点截面上沿试验件高度方向设置测量点,测量点的数量大于等于三个。
本发明提供的多点弯曲测量复合材料面外刚度的方法在S2步骤中通过试验获得实验数据,测得的数据通过数据处理软件进行参数反演计算。反演计算过程是通过特定编制的软件进行的,通过载荷及对应测量的点的位移,计算得出对应的应变下面外刚度数据。不仅如此,使用同样的特定编制的软件通过多个载荷和对应测量的的位移数据可得出应变-模量曲线。采用特定编制的软件进行反演计算的流程如图1所示。
本发明提供的多点弯曲测量复合材料面外刚度的方法可作为:复合材料面外刚度测量的试验方法;复合材料有限元三维计算所需刚度数据输入;随炉试验件性能测量的试验方法。
在具体的使用过程中,用五点弯曲方法测量复合材料(材料牌号:5228A/CCF300)的面外刚度E3,G13,G23,选用100KN万能材料试验机作为试验机。制备试验件,5228A/CCF300材料单层厚度0.125mm,试验件铺层[0]120,厚度15mm,宽度40mm,长度135mm。准备五点弯曲夹具。选择五点弯曲是因为五点弯曲可以在试验件上形成两个纯剪切面,该面弯矩为零,有利于试验数据的处理。
进行试验过程中,以试验件的长、宽、高建立坐标系,对应的坐标轴分别为X、Y、Z,试验载荷为0-10KN,加载过程中使用DIC测量设备测量如图2标注的A1、A2、A3、B1、B2、B3共计6个点的位移。选取任意载荷下的位移测量数据,如8KN时的位移如下表。
将试验数据输入数据处理软件,使用参数反演方法进行面外刚度计算,面外刚度初值可任意设置,本实施例中将面外刚度例设为[E3=15000MPa,G13=9000MPa,G23=6000MPa]作为初值。计算迭代过程见下表,所生成的迭代过程曲线如图3所示。
E3 | G13 | G23 | |
1 | 1.50E+04 | 9.00E+03 | 6.00E+03 |
2 | 8.75E+03 | 5.36E+03 | 3.92E+03 |
3 | 9.61E+03 | 5.96E+03 | 3.06E+03 |
4 | 9.651E+03 | 6.02E+03 | 3.098E+03 |
5 | 9.651E+03 | 6.02E+03 | 3.098E+03 |
最终测量并计算得到的面外刚度为E3=9651MPa,G13=6002MPa,G23=3098MPa。
以上结合本发明的多点弯曲测量复合材料面外刚度的方法具体实施例做了详细描述,但并非是对本发明的限制,凡是依据本发明的技术实质对以上实施例所做的任何简单修改均属于本发明的技术范围,还需要说明的是,按照本发明的多点弯曲测量复合材料面外刚度的方法技术方案的范畴包括上述各部分之间的任意组合。
Claims (8)
1.一种多点弯曲测量复合材料面外刚度的方法,其特征在于,包括如下步骤:
S1:实验准备,选用一般的万能材料试验机作为试验机,制作复合材料试验件,确定对试验件进行弯曲的点数并准备试验夹具;
S2:按照S1中选定的多点弯曲方法进行试验和测量并得出相应的数据;
S3:对S2中得到的数据进行参数反演,得出复合材料面外刚度。
2.根据权利要求1所述的一种多点弯曲测量复合材料面外刚度的方法,其特征在于,S1中选取的试验机的最小的工作载荷为100KN。
3.根据权利要求1所述的一种多点弯曲测量复合材料面外刚度的方法,其特征在于,S1中制作的复合材料试验件为规则的长方体结构,并以试验件的长、宽、高建立坐标系,定义:对试验件进行弯曲的点数为N。
4.根据权利要求3所述的一种多点弯曲测量复合材料面外刚度的方法,其特征在于,定义:S1中试验件的长度取值为L,其中,L=45*(N-2)。
5.根据权利要求1所述的一种多点弯曲测量复合材料面外刚度的方法,其特征在于,S2中在试验机上居中放置试验件,确定以试验件的高度方向为试验机加载方向,调整每组加载点和支持点之间的距离,对试验件施加压缩载荷,对试验件每个加载点截面上沿试验件高度方向设置测量点,测量点的数量大于等于三个。
6.根据权利要求1所述的一种多点弯曲测量复合材料面外刚度的方法,其特征在于,通过数据处理软件进行参数反演计算。
7.根据权利要求1所述的一种多点弯曲测量复合材料面外刚度的方法,其特征在于,通过载荷及对应测量的点的位移,计算得出对应的应变下面外刚度数据。
8.根据权利要求1所述的一种多点弯曲测量复合材料面外刚度的方法,其特征在于,通过多个载荷和对应测量的位移数据得出应变-模量曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510494369.7A CN105181480B (zh) | 2015-08-13 | 2015-08-13 | 一种多点弯曲测量复合材料面外刚度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510494369.7A CN105181480B (zh) | 2015-08-13 | 2015-08-13 | 一种多点弯曲测量复合材料面外刚度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105181480A true CN105181480A (zh) | 2015-12-23 |
CN105181480B CN105181480B (zh) | 2018-08-24 |
Family
ID=54903716
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510494369.7A Active CN105181480B (zh) | 2015-08-13 | 2015-08-13 | 一种多点弯曲测量复合材料面外刚度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105181480B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109724727A (zh) * | 2017-10-31 | 2019-05-07 | 中国航发商用航空发动机有限责任公司 | 燃气轮机曲面叶片残余应力测量方法及其装置 |
CN109991101A (zh) * | 2019-05-16 | 2019-07-09 | 曾庆雪 | 一种钢结构横梁弯曲度检测装置及检测方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100131256A1 (en) * | 2008-11-26 | 2010-05-27 | Livermore Software Technology Corporation | Spot Weld Failure Determination Method in a Finite Element Analysis |
CN102384875A (zh) * | 2011-11-09 | 2012-03-21 | 吉林大学 | 显微镜下拉压弯复合载荷模式材料力学性能测试装置 |
CN102721609A (zh) * | 2012-07-03 | 2012-10-10 | 济南科汇试验设备有限公司 | 万能试验机 |
CN102830020A (zh) * | 2012-07-18 | 2012-12-19 | 中国电力科学研究院 | 一种特高压瓷套管连接弯曲刚度测试装置及安装测试方法 |
CN103207118A (zh) * | 2012-01-16 | 2013-07-17 | 中国石油天然气集团公司 | 钢管弯曲变形中试验管大变形及应力实时测量装置及方法 |
CN103616299A (zh) * | 2013-12-12 | 2014-03-05 | 山东大学 | 用于双材料界面混合模态断裂韧性测试的试件及测试方法 |
CN103760037A (zh) * | 2014-01-23 | 2014-04-30 | 湖南大学 | 一种轴承式损伤生物力学三点弯曲试验装置和试验方法 |
CN204188470U (zh) * | 2014-11-11 | 2015-03-04 | 吉林大学 | 三点、四点弯曲作用下材料微观力学性能的原位测试仪器 |
-
2015
- 2015-08-13 CN CN201510494369.7A patent/CN105181480B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100131256A1 (en) * | 2008-11-26 | 2010-05-27 | Livermore Software Technology Corporation | Spot Weld Failure Determination Method in a Finite Element Analysis |
CN102384875A (zh) * | 2011-11-09 | 2012-03-21 | 吉林大学 | 显微镜下拉压弯复合载荷模式材料力学性能测试装置 |
CN103207118A (zh) * | 2012-01-16 | 2013-07-17 | 中国石油天然气集团公司 | 钢管弯曲变形中试验管大变形及应力实时测量装置及方法 |
CN102721609A (zh) * | 2012-07-03 | 2012-10-10 | 济南科汇试验设备有限公司 | 万能试验机 |
CN102830020A (zh) * | 2012-07-18 | 2012-12-19 | 中国电力科学研究院 | 一种特高压瓷套管连接弯曲刚度测试装置及安装测试方法 |
CN103616299A (zh) * | 2013-12-12 | 2014-03-05 | 山东大学 | 用于双材料界面混合模态断裂韧性测试的试件及测试方法 |
CN103760037A (zh) * | 2014-01-23 | 2014-04-30 | 湖南大学 | 一种轴承式损伤生物力学三点弯曲试验装置和试验方法 |
CN204188470U (zh) * | 2014-11-11 | 2015-03-04 | 吉林大学 | 三点、四点弯曲作用下材料微观力学性能的原位测试仪器 |
Non-Patent Citations (1)
Title |
---|
李家驹: "论复合材料弯曲实验", 《玻璃钢/复合材料》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109724727A (zh) * | 2017-10-31 | 2019-05-07 | 中国航发商用航空发动机有限责任公司 | 燃气轮机曲面叶片残余应力测量方法及其装置 |
CN109991101A (zh) * | 2019-05-16 | 2019-07-09 | 曾庆雪 | 一种钢结构横梁弯曲度检测装置及检测方法 |
CN109991101B (zh) * | 2019-05-16 | 2022-01-25 | 曾庆雪 | 一种钢结构横梁弯曲度检测装置及检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105181480B (zh) | 2018-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bauchau et al. | Euler-Bernoulli beam theory | |
Hanabusa et al. | Numerical verification of a biaxial tensile test method using a cruciform specimen | |
WO2019047529A1 (zh) | 一种纤维增强复合材料动态剪切本构模型的构建方法 | |
CN102564856B (zh) | 基于数字图像相关的塑性多缺陷材料m积分测量方法 | |
Liao et al. | 3-D FEM stress analysis and strength evaluation of single-lap adhesive joints subjected to impact tensile loads | |
He | Influence of boundary conditions on stress distributions in a single-lap adhesively bonded joint | |
CN104316388A (zh) | 一种对各向异性材料结构件进行疲劳寿命测定的方法 | |
CN109323791A (zh) | 基于增量切割法的复合材料板材的残余应力分布测量方法 | |
Gunwant et al. | Stress and displacement analysis of a rectangular plate with central elliptical hole | |
CN105181480A (zh) | 一种多点弯曲测量复合材料面外刚度的方法 | |
Kang et al. | Biaxial tensile failure of short carbon-fibre-reinforced PEEK composites | |
JP2020004357A (ja) | フロントカバーの最適化設計方法、及び装置 | |
Kavdir et al. | The experimental and numerical study on the mechanical behaviours of adhesively bonded joints | |
Dionísio et al. | Fracture mechanics approach to stress singularities in composite adhesive joints | |
Matthews et al. | Finite element analysis for large displacement J-integral test method for Mode I interlaminar fracture in composite materials | |
CN104636543A (zh) | 一种基于有限差分法的重型龙门铣床横梁重力变形预测方法 | |
CN104392032A (zh) | 基于有限元法的纱线材料参数识别方法 | |
CN111855057A (zh) | 基于盲孔法的复合材料板材的残余应力分布测量方法 | |
CN107024401A (zh) | 获取金属材料各向异性和拉压非对称性的方法及系统 | |
CN111339697B (zh) | 一种基于修正梁理论的复合材料i型分层桥联律确定方法 | |
Ramaswamy et al. | Software analysis on free edge stresses in the composite laminates using various approaches and its applications | |
JP2023065041A (ja) | 応力拡大係数推定方法及び疲労強度推定方法 | |
CN104765928B (zh) | 一种塑性成形摩擦特性测算方法 | |
CN108195532B (zh) | 一种梁结构裂纹等效刚度的测量方法 | |
Ping et al. | Analysis of generalized stress intensity factors of V-shaped notch problems by FEM |
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 |