CN101446528B - 纤维增强复合材料二维随机孔隙模型的建立方法 - Google Patents
纤维增强复合材料二维随机孔隙模型的建立方法 Download PDFInfo
- Publication number
- CN101446528B CN101446528B CN2008102303592A CN200810230359A CN101446528B CN 101446528 B CN101446528 B CN 101446528B CN 2008102303592 A CN2008102303592 A CN 2008102303592A CN 200810230359 A CN200810230359 A CN 200810230359A CN 101446528 B CN101446528 B CN 101446528B
- Authority
- CN
- China
- Prior art keywords
- dimensional random
- model
- porosity
- reinforced composite
- fiber reinforced
- 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
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
一种纤维增强复合材料二维随机孔隙模型的建立方法,属于复合材料无损检测与评价领域。这种纤维增强复合材料二维随机孔隙模型的建立方法首先对纤维增强复合材料标准样品进行解剖,通过金相法统计孔隙率及孔隙分布特征,根据统计结果建立二维随机介质模型M(x,z);然后利用极值搜索法改造二维随机介质模型M(x,z),建立二维随机孔隙模型,统计二维随机孔隙模型的孔隙率及孔隙分布特征;最后对二组孔隙率及孔隙分布特征进行数据比较,若二者差异小于2%,则得到纤维增强复合材料二维随机孔隙模型,该方法可精确建立纤维增强复合材料二维随机孔隙模型,利用该模型可开展纤维增强复合材料孔隙的弹性性能分析,为纤维增强复合材料孔隙率超声无损检测提供分析依据。
Description
技术领域
本发明涉及一种纤维增强复合材料二维随机孔隙模型的建立方法,其属于复合材料无损检测与评价技术领域。
背景技术
迄今为止,用于建立含孔隙的纤维增强复合材料二维几何模型的方法,均基于以下类似假设:所有孔隙为球形、半径相等,均匀分布于介质中;以孔隙率1.5%为界假定孔隙为球状或圆盘状且尺寸符合不同分布等。然而,实际纤维增强复合材料中的孔隙尺寸变化范围很大,形状各异,分布也并不均匀。已有几何模型不能全面、准确地反映孔隙形貌及其随机特性,导致孔隙率超声检测实验和理论分析之间存在较大偏差。
发明内容
为了克服上述现有技术中存在的问题,本发明提供一种纤维增强复合材料二维随机孔隙模型的建立方法。该方法可精确建立纤维增强复合材料二维随机孔隙模型,利用该模型可开展纤维增强复合材料孔隙的弹性性能分析,为纤维增强复合材料孔隙率超声无损检测提供分析依据。
本发明的技术方案是:一种纤维增强复合材料二维随机孔隙模型的建立方法如下具体步骤:
(1)对纤维增强复合材料标准样品进行解剖,通过金相法统计由孔隙长度、宽度和宽长比表示的孔隙率及孔隙分布特征;
(2)根据统计结果初步选定自相关函数、自相关长度参数建立二维随机介质模型M(x,z);
(3)选定孔隙率、分区边长参数,利用极值搜索法改造所得的二维随机介质模型M(x,z),建立二维随机孔隙模型;
(4)统计所得的二维随机孔隙模型的孔隙率及孔隙分布特征;
(5)比较步骤(1)与(4)的数据,若二者差异小于2%,则得到纤维增强复合材料二维随机孔隙模型,否则改变建立参数,重复步骤(2)~(5)。
所述极值搜索法包括如下具体步骤:
(1)选定孔隙率P及分区边长D;
(2)将所得到的二维随机介质模型M(x,z),以分区边长D划分为若干区域,并在每个区域内进行以下操作:
①获取该区域所有节点数量L;
②将区域内各节点按数值由大到小排列,放入一维数组H;
③获取H中前L*P个节点在二维随机介质模型M(x,z)中的坐标,放入一维数组;
④将M(x,z)中符合一维数组中所记录坐标的节点设为孔隙,其余设为纤维增强复合材料。
本发明的有益效果是:这种纤维增强复合材料二维随机孔隙模型的建立方法首先对纤维增强复合材料标准样品进行解剖,通过金相法统计孔隙率及孔隙分布特征,根据统计结果建立二维随机介质模型M(x,z);然后利用极值搜索法改造二维随机介质模型M(x,z),建立二维随机孔隙模型,统计二维随机孔隙模型的孔隙率及孔隙分布特征;最后对二组孔隙率及孔隙分布特征进行数据比较,若二者差异小于2%,则得到纤维增强复合材料二维随机孔隙模型,该方法可精确建立纤维增强复合材料二维随机孔隙模型,利用该模型可开展纤维增强复合材料孔隙的弹性性能分析,为纤维增强复合材料孔隙率超声无损检测提供分析依据。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1为一种纤维增强复合材料二维随机孔隙模型的建立方法的流程图。
图2为碳纤维单向增强复合材料二维随机孔隙模型,建立参数为:高斯型自相关函数;a=7,b=10;P=1.2%;D=4。
具体实施方式
图1为一种纤维增强复合材料二维随机孔隙模型的建立方法的流程图。下面以碳纤维单向增强复合材料为例,介绍随机孔隙模型建立的具体实施方式:
(1)将碳纤维单向增强复合材料标准试样沿与纤维垂直方向解剖,使用金相显微镜观察并统计截面孔隙率和每个孔隙的长度,宽度,长宽比。重复上述步骤20次以上,以得到标准试样的体积孔隙率及孔隙特征的样本空间;
(2)选取高斯型自相关函数:其中a,b为自相关长度,初步取值为:a=7,b=10;将进行二维快速傅里叶变换得到生成二维随机场θ∈[0,2π];构造随机功率谱对随机功率谱做傅里叶逆变换,得到随机扰动δM(x,z);利用公式M(x,z)=M0+δM(x,z)得到随机介质模型。其中,M可为密度、拉梅参数或声速等参量,M0为M在大尺度不均匀性上的平均值,以纵波声速为例取值为M0=2.6km/s;
(3)选定孔隙率P=1.2%及分区边长D=4;将M(x,z)若干区域,并在每个区域内进行以下操作:
①获取该区域所有节点数量L;
②将区域内各节点按数值由大到小排列,放入一维数组H;
③获取H中前L*P个节点在二维随机介质模型M(x,z)中的坐标,放入一维数组;
④将M(x,z)中符合一维数组中所记录坐标的节点设为孔隙,其余设为复合材料。
(4)统计步骤(3)所得模型的孔隙率及孔隙分布特征,同步骤(1);
(5)优化建立参数。比较步骤(1)与(4)的数据,若二者差异小于2%则得到纤维增强复合材料二维随机孔隙模型,否则改变建立参数,重复步骤(2)~(5)。
图2示出了碳纤维单向增强复合材料二维随机孔隙模型,建立参数为:高斯型自相关函数;a=7,b=10;P=1.2%;D=4。
Claims (1)
1.一种纤维增强复合材料二维随机孔隙模型的建立方法,其特征在于:包括如下具体步骤:
(1)对纤维增强复合材料标准样品进行解剖,通过金相法统计由孔隙长度、宽度和宽长比表示的孔隙率及孔隙分布特征;
(2)根据统计结果初步选定自相关函数、自相关长度参数建立二维随机介质模型M(x,z);
(3)利用极值搜索法改造所得的二维随机介质模型M(x,z),建立二维随机孔隙模型,极值搜索法包括如下具体步骤:
选定孔隙率P及分区边长D,将所得到的二维随机介质模型M(x,z),以分区边长D划分为若干区域,并在每个区域内进行以下操作:
①获取该区域所有节点数量L,
②将区域内各节点按数值由大到小排列,放入一维数组H,
③获取H中前L*P个节点在二维随机介质模型M(x,z)中的坐标,放入一维数组Position,
④将M(x,z)中符合Position中所记录坐标的节点设为孔隙,其余设为纤维增强复合材料;
(4)统计所得的二维随机孔隙模型的孔隙率及孔隙分布特征;
(5)比较步骤(1)与(4)的数据,若二者差异小于2%,则得到纤维增强复合材料二维随机孔隙模型,否则改变建立参数,重复步骤(2)-(5)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102303592A CN101446528B (zh) | 2008-12-25 | 2008-12-25 | 纤维增强复合材料二维随机孔隙模型的建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102303592A CN101446528B (zh) | 2008-12-25 | 2008-12-25 | 纤维增强复合材料二维随机孔隙模型的建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101446528A CN101446528A (zh) | 2009-06-03 |
CN101446528B true CN101446528B (zh) | 2011-02-02 |
Family
ID=40742311
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008102303592A Active CN101446528B (zh) | 2008-12-25 | 2008-12-25 | 纤维增强复合材料二维随机孔隙模型的建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101446528B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102157015A (zh) * | 2011-04-24 | 2011-08-17 | 大连理工大学 | 纤维增强复合材料三维随机孔隙模型的建立方法 |
CN102507411B (zh) * | 2011-10-27 | 2013-11-06 | 哈尔滨工业大学 | 基于金相照片和图像分析技术分析铝基多孔复合材料中微孔分布的方法 |
CN102542587A (zh) * | 2012-01-17 | 2012-07-04 | 大连理工大学 | 一种孔隙尺寸离散度大的纤维增强复合材料二维随机孔隙模型建立方法 |
CN104318039B (zh) * | 2014-11-18 | 2017-12-15 | 张光亮 | 硬质合金自适应参数化建模的方法 |
CN105044146B (zh) * | 2015-07-02 | 2017-12-22 | 西安电子科技大学 | 一种两相复合材料的随机热均化分析方法 |
CN108152373B (zh) * | 2017-11-13 | 2021-07-16 | 大连理工大学 | 基于孔隙尺寸识别的含孔隙复合材料力学性能测量系统及方法 |
CN115859402B (zh) * | 2022-12-01 | 2023-09-12 | 哈尔滨师范大学 | 一种基于改进四参数随机生长法的植物纤维地膜结构模拟方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1546979A (zh) * | 2003-12-17 | 2004-11-17 | 东华大学 | 纤维压缩弯曲性能测量的方法及装置 |
CN1888860A (zh) * | 2006-07-24 | 2007-01-03 | 中国林业科学研究院木材工业研究所 | 植物纤维材料中化学成分含量模型的建立和含量测定方法 |
-
2008
- 2008-12-25 CN CN2008102303592A patent/CN101446528B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1546979A (zh) * | 2003-12-17 | 2004-11-17 | 东华大学 | 纤维压缩弯曲性能测量的方法及装置 |
CN1888860A (zh) * | 2006-07-24 | 2007-01-03 | 中国林业科学研究院木材工业研究所 | 植物纤维材料中化学成分含量模型的建立和含量测定方法 |
Non-Patent Citations (1)
Title |
---|
刘继忠等.碳纤维复合材料孔隙率的脉冲反射法超声衰减测试模型.《浙江大学学报》.2006,第40卷(第11期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN101446528A (zh) | 2009-06-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101446528B (zh) | 纤维增强复合材料二维随机孔隙模型的建立方法 | |
Van Doormaal et al. | Determination of permeability in fibrous porous media using the lattice Boltzmann method with application to PEM fuel cells | |
Park et al. | Application of lattice Boltzmann method to a micro-scale flow simulation in the porous electrode of a PEM fuel cell | |
CN110147561B (zh) | 一种含天然裂缝致密油气储层体积压裂缝网预测方法 | |
CN102542587A (zh) | 一种孔隙尺寸离散度大的纤维增强复合材料二维随机孔隙模型建立方法 | |
CN110133035B (zh) | 核磁共振耦合恒速压汞定量表征砂岩储层孔喉结构的方法 | |
CN103713316A (zh) | 一种基于岩石孔隙数字表征的速度预测方法及装置 | |
Wang | Structure effects on electro-osmosis in microporous media | |
CN102157015A (zh) | 纤维增强复合材料三维随机孔隙模型的建立方法 | |
Liao et al. | Effect of structural parameters on mass transfer characteristics in the gas diffusion layer of proton exchange membrane fuel cells using the lattice boltzmann method | |
CN108303360A (zh) | 一种煤岩三维孔隙网络结构参数表征方法 | |
CN108921945A (zh) | 结合居中轴线与实体模型的孔隙网络模型构建方法 | |
Zhu et al. | A fractal electrical conductivity model for water-saturated tree-like branching network | |
Delisée et al. | 3D morphology and permeability of highly porous cellulosic fibrous material | |
Yu et al. | Numerical study on the effect of microporous layer crack changes on water management in gas diffusion layer of proton exchange membrane fuel cell | |
CN113281149B (zh) | 一种节理岩体的表征单元体积尺度综合取值方法 | |
CN105334293A (zh) | 多重介质中致密气的流态识别方法及其装置 | |
CN110632677B (zh) | 一种碳酸盐岩地热储层分类方法 | |
CN103235879A (zh) | 基于小生境主从式并行遗传算法的双相介质参数反演方法 | |
Liu et al. | Production and construction quality management system of prefabricated buildings based on BIM technology | |
Zhao et al. | Percolation laws of a fractal fracture-pore double medium | |
Afrookhteh et al. | Stochastic reconstruction of carbon fiber paper gas diffusion layers of PEFCs: A comparative study | |
CN114460639B (zh) | 页岩油储层渗透率的预测方法及装置 | |
CN113758791A (zh) | 岩石的力学参数确定方法 | |
Hossain et al. | Non-Darcy free convective flow along a vertical cylinder embedded in a porous medium with surface mass flux |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |