Processing math: 100%

3D打印贝壳仿生复合材料的拉伸力学行为

侯祥龙 雷建银 李世强 王志华 刘志芳

边江, 曹学文, 杨文, 于洪喜, 尹鹏博. 入口压力对天然气超声速液化特性的影响[J]. 高压物理学报, 2018, 32(3): 031101. doi: 10.11858/gywlxb.20170639
引用本文: 侯祥龙, 雷建银, 李世强, 王志华, 刘志芳. 3D打印贝壳仿生复合材料的拉伸力学行为[J]. 高压物理学报, 2020, 34(1): 014102. doi: 10.11858/gywlxb.20190768
BIAN Jiang, CAO Xuewen, YANG Wen, YU Hongxi, YIN Pengbo. Influence of Inlet Pressure on Supersonic Liquefaction of Natural Gas Mixtures[J]. Chinese Journal of High Pressure Physics, 2018, 32(3): 031101. doi: 10.11858/gywlxb.20170639
Citation: HOU Xianglong, LEI Jianyin, LI Shiqiang, WANG Zhihua, LIU Zhifang. Tension Mechanical Behavior of 3D Printed Composite Materials Inspired by Nacre[J]. Chinese Journal of High Pressure Physics, 2020, 34(1): 014102. doi: 10.11858/gywlxb.20190768

3D打印贝壳仿生复合材料的拉伸力学行为

doi: 10.11858/gywlxb.20190768
基金项目: 国家自然科学基金(11602161,11772216);山西省1331重点创新团队项目
详细信息
    作者简介:

    侯祥龙(1988-),男,硕士,主要从事材料冲击动力学行为研究. E-mail: 592285804@qq.com

    通讯作者:

    李世强(1986-),男,博士,讲师,主要从事材料及结构冲击动力学行为研究.E-mail: lishiqiang@tyut.edu.cn

  • 中图分类号: O341

Tension Mechanical Behavior of 3D Printed Composite Materials Inspired by Nacre

  • 摘要: 采用硬质和软质双组分材料,通过调控两种基体材料的装配夹角,采用光固化3D打印技术制备了不同装配方式的仿贝壳珍珠层复合材料,开展了准静态拉伸实验,结合扫描电镜观察,分析了其拉伸力学性能、断裂及能量耗散机理。研究结果表明,保持胞元边长不变,随着面内装配角度增加,仿贝壳珍珠层复合材料的强度呈线性增加趋势,断裂应变呈线性减小的趋势;随着面外装配角度增大,断裂应变呈线性减小趋势,而强度在面外装配角小于45°时呈增强趋势,超过45°时趋于稳定;面外装配角度为45°时,材料的强度达到最大值。试样在断裂前主要通过硬质材料的拔出、软/硬相界面处微裂纹的生成及微裂纹在扩展过程中的合并和偏转等方式耗散能量。

     

  • 近年来随着液化天然气(LNG,Liquefied Natural Gas)产业在全球迅速发展,天然气的液化技术和设备也在不断发展完善、日渐成熟。我国LNG领域内的相关研究起步较晚,许多技术远远落后发达国家水平,在天然气液化工艺及装置的生产等方面缺乏自主产权。因此,开展天然气液化工艺及装置的研究,对于实现液化装置的国产化、高效化有十分重要的意义[1-3]

    超声速旋流分离技术是一种新兴的天然气加工处理技术,被较为广泛地用于天然气脱水、脱重烃、脱酸等方面,近年来开始逐渐应用于天然气液化方面[4-7]。天然气超声速液化的原理是:高压天然气混合物在Laval喷管内达到一定的温度、压力条件,开始凝结成核,最终凝结成液滴,在后续工艺中进行进一步气液分离。与传统的天然气液化技术相比,具有结构工艺简单、支持无人操作(适用于海底天然气处理)、对水合物抑制剂依赖性小、投资和运行成本低等优势[8-9]

    为了探究入口复杂多变的压力条件对天然气超声速液化特性的影响,对甲烷-乙烷气体混合物的超声速凝结流动特性进行研究,在凝结成核与生长理论的基础上建立了适用于甲烷-乙烷双可凝气体混合物的凝结流动数学模型,重点研究了入口压力对天然气混合物在Laval喷管内主要流动与凝结参数的影响规律。

    Laval喷管结构主要包括入口段、渐缩段、喉部及扩张段4部分[10-11]。各部分参数如表 1所示,L0为入口长度,r1为渐缩段入口半径,L1为渐缩段长度,rcr为喉部截面半径,L2为渐扩段长度,r2为渐扩段出口半径。为尽量减小流场涡流的影响,渐缩段采用双三次曲线设计,喉部采用一段平缓光滑的圆弧作为过渡曲线,渐扩段采用等膨胀率设计,膨胀率取为10 000 s-1。考虑到实验加工方便,保证曲面的精度,且能够更加直观地观察Laval喷管内部的流场分布情况,所设计Laval喷管截面采用矩形截面,三维结构如图 1所示。

    表  1  Laval喷管各部分参数
    Table  1.  Parameters of Laval nozzle
    L0/mm r1/mm L1/mm rcr/mm L2/mm r2/mm
    50.00 20.00 56.01 2.50 71.28 6.15
    下载: 导出CSV 
    | 显示表格
    图  1  Laval喷管三维结构
    Figure  1.  Structure of Laval nozzle

    采用欧拉双流体模型开展数值计算,控制方程主要包括气相流动方程和液相流动方程。在无滑移假设及欧拉双流体模型的前提下分别建立气相及液相流动控制方程组,液滴数目守恒方程及液滴半径、数目、湿度关系式分别添加到对应源相方程中[12-14]

    气相流动控制方程组

    ρvt+xj(ρvuj)=Sm (1)
    ρc2t+xj(ρc2uj)=Sm,c2 (2)
    t(ρvui)+xj(ρvujui)=pvxi+xj[μ(ujxi+uixj23δijujxj)]+xj(ρv¯uiuj)+Su (3)
    t(ρvE)+xj(ρvujE+ujpv)=xj(keffTxj+uiτeff)+Sh (4)

    液相流动控制方程组

    t(ρY)+xj(ρujY)=SY (5)
    ρNt+xj(ρNuj)=J (6)
    rd=33Y/(4πρLN) (7)

    式中:uiuj为时均速度分量,m/s;ρv为气相密度,kg/m3ρ为气液混合相密度,kg/m3pv为时均压力,Pa;μ为黏度,kg/(m·s);δij为Kronecker delta数;E为总能量,J/kg;keff为有效导热系数,W/(m·K);τeff为有效应力张量,无量纲;Y为液相质量分数,无量纲;rd为液滴半径,m;drd/dt为液滴生长速率,m/s;N为液滴数目,kg-1

    成核模型采用文献[15-16]中提出的双组分气体自发凝结成核模型修正方法。液滴生长过程采用Gyarmathy液滴生长模型,模型中液滴与气体间的传热系数[17-18]

    kr=λvrd11+28π1.5Pvγ1+γKn (8)

    依据传热、传质过程,可推导得到液滴生长速率计算模型

    drddt=λvρLhLV(TsT)(1rcrd)rd(1+28π1.5Pvγ1+γKn) (9)

    式中:λv为气体导热系数,W/(m·K);Pv为气体Prandtl数;γ为气体比热比;hLV为凝结潜热,J/kg;Ts为气体压力对应的饱和温度,K;Kn表示Kundsen数。由于双组分气体不存在压力对应的饱和温度Ts这一概念,将双组分相图中露点线类比于单组分中饱和曲线。

    针对气体状态方程的选择,由于低温气体已偏离理想气体假设,本研究采用了NIST真实气体模型进行计算。

    湍流发生时会导致流体之间相互交换动量、能量,也会造成浓度的改变。本研究建模时忽略相间速度的滑移,即液滴产生不影响湍流,因此只考虑气相的湍流方程。FLUENT中提供了以下几种湍流模型:S-A模型、标准k-ε模型、RNG k-ε模型、Realizable k-ε模型、k-ω模型以及雷诺应力模型。S-A模型主要应用流动分离区附近模拟,标准k-ε模型、RNG k-ε模型一般用于各向同性的均匀湍流,k-ω模型可用于带压梯度的流动模拟和跨声速激波模拟,雷诺应力模型主要用于龙卷风、燃烧室等强烈旋转流动的模拟。对Laval喷管内跨声速流动,采用k-ω模型可以获得较为理想的计算精度和计算速度,故本研究采用该模型进行数值计算。

    气体在Laval喷管中的流动属于高速可压缩流动,采用密度基进行求解,流动控制方程组、湍流动能方程、湍流耗散率方程均采用二阶迎风格式进行离散。

    根据双组分气体在Laval喷管内的高速可压缩的流动特性,入口和出口边界设置为压力入口边界和压力出口边界条件,对于气体在Laval喷管内的超声速流动,由于所有的流动参数都可从Laval喷管内部外推得到,故在出口处不进行相应设置,固体壁面边界设置为无滑移、无渗流、绝热边界条件。

    在数学模型中,由于气相方程添加了源相方程,液相方程定义了标量以及引入的真实气体方程,这些仅靠在FLUENT自带的模型和材料物性无法满足要求,需要编写相应的用户自定义函数(UDF)。本研究编写UDF时,分别定义DEFINE AJUST、DEFINE SOURCE和DEFINE PROPERTY 3个宏函数。DEFINE AJUST宏用来定义过饱和度、过冷度、成核速率、液滴生长率、液滴半径、液滴质量以及液滴表面张力等参数,DEFINE SOURCE宏用来定义控制方程中的质量、动量和能量源相,DEFINE PROPERTY用来定义数值计算中用到的真实气体的热力学参数如黏度系数、导热系数等。

    为验证所建立的双组分气体凝结数学模型及数值计算方法的准确性,采用本研究所设计的Laval喷管结构,在中国石油大学(华东)超声速气体凝结流动实验系统开展了水-乙醇双可凝组分气体凝结相变实验研究。实验条件为:Laval喷管入口压力0.586 MPa,入口温度288.05 K,气体湿度98.1%,水与乙醇摩尔体积比84:16,气体流量为323.78标方每小时,实验测得的Laval喷管沿程压力分布如图 2所示,可以看出,压力分布实验结果与数值计算结果吻合较好,说明本研究所建立的双组分气体超声速凝结流动特性数学模型及数值计算方法具有一定的准确性和可靠性。

    图  2  Laval喷管沿程压力分布数据对比
    Figure  2.  Comparison of pressure distribution data in Laval nozzle

    保持入口温度及组成(甲烷体积分数90%、乙烷体积分数10%)不变,研究不同的入口压力对Laval喷管内部甲烷-乙烷双组分气体凝结过程中压力、温度、成核率、液滴生长率、液滴半径、液相质量分数的影响。在数值计算中设定的入口温度为270 K,设定入口压力分别为5.5、6.0和6.5 MPa。Laval喷管内双组分气体凝结参数的变化趋势及对比如图 3~图 8所示。

    图  3  Laval喷管内气体压力分布
    Figure  3.  Gas pressure distribution in Laval nozzle
    图  4  Laval喷管内气体温度分布
    Figure  4.  Gas temperature distribution in Laval nozzle
    图  5  Laval喷管内成核率分布
    Figure  5.  Nucleation rate distribution in Laval nozzle
    图  6  Laval喷管内液滴半径分布
    Figure  6.  Droplet radius distribution in Laval nozzle
    图  7  Laval喷管内液滴生长率分布
    Figure  7.  Droplet growth rate distribution in Laval nozzle
    图  8  Laval喷管内液相质量分数分布
    Figure  8.  Liquid mass fraction distribution in Laval nozzle

    从压力与温度分布可以看出,气体进入Laval喷管后压力、温度不断降低,当达到一定过冷度时,气体发生凝结并释放潜热,但凝结突跃现象对压力造成的影响并不显著,压力在Laval喷管渐扩段减小到了一个比较稳定状态,温度在减小到最小值后又略微上升,这主要是由于液滴凝结释放潜热引起的。随着入口压力的增大,出口压力略微升高,温升位置有所提前,出口温度也越高,这是因为,随着入口压力的增大,液滴成核与凝结量也随之增多,释放的潜热也就越多。

    从成核率分布可以看出,保持其他条件一致,当压力发生变化时,成核速率的变化趋势几乎相同,在刚进入Laval喷管的一段距离内为零,在某一位置处开始,成核率从零开始突跃一直增大到峰值后迅速减小至零。随着入口压力从5.5 MPa增大到6.5 MPa,成核的发生位置(Wilson点)不断向前移动,逐渐向Laval喷管喉部靠拢,且成核率的最值逐渐增大。当压力为5.5 MPa时,成核发生位置为x=0.120 6 m,且在x=0.147 5 m处达到极限成核,为4.044×1020 m-3·s-1;当压力为6 MPa时,成核发生位置较5.5 MPa时向前移动,为x=0.119 9 m,极限成核位置也随压力的增大而前移,在x=0.139 3 m处达到8.062×1020 m-3·s-1;当压力继续增大为6.5 MPa时,成核发生位置较6 MPa时更加靠近喉部,为x=0.118 2 m,极限成核位置在x=0.132 8 m处,为9.015×1020 m-3·s-1

    从液滴半径分布可以看出,随着入口压力的增大,Laval喷管内平均液滴半径越大,出口液滴半径也随之增大,当压力为5.5、6.0、6.5 MPa时,对应最大液滴半径尺寸分别为415.86、447.88和477.44 nm。由此可知,压力的升高有利于液滴的生长。

    从液滴生长率分布可以看出,液滴生长率在气体刚进入Laval喷管时一直为零,当液滴开始发生成核凝结时液滴生长率开始突增,变化到最大值后又迅速减小,最终减小为零。综合图 6图 7还可以看出,随着入口压力的升高,在成核开始时液滴生长率较大,液滴半径增长速度较快,但一段距离后液滴生长率下降更大,液滴半径增长速度也明显放缓。

    从液相质量分数分布可以看出,伴随着混合气体的凝结成核,液相质量分数也不断增大,且随着入口压力的升高,Laval喷管出口处的湿度值随之增大,当压力为5.5 MPa时,湿度的最大值为3.989 2%,当压力增大到6.5 MPa时湿度最终增大到7.382 0%。

    (1) 建立了三维双组分天然气混合物超声速凝结流动数学模型,对Laval喷管内双组分混合物凝结流动进行了数值模拟,得出沿Laval喷管轴向的参数分布,通过开展双可凝组分气体凝结相变实验,对比发现数值模拟与实验结果基本一致,说明了所建立数学模型及计算方法的正确性。

    (2) 利用数值模型研究了入口参数对天然气混合物超声速液化特性的影响,结果表明,保持Laval喷管入口温度及组成不变,增大入口压力,混合气体成核位置前移,成核率、平均液滴半径、液相质量分数均随之增大,即入口压力越大,混合气体在Laval喷管内越易发生凝结。在实际生产中,可以通过调节入口压力来促进天然气的凝结,提高Laval喷管的液化效率。

  • 图  贝壳结构示意图[1]

    Figure  1.  Schematic diagram of shell structure[1]

    图  双相材料介观几何构型与组合方式以及结构优化示意图

    Figure  2.  Mesoscale structure and combination of biphase materials and structural-optimized diagrams

    图  硬质材料(a)与软质材料(b)在准静态拉伸下的应力-应变关系

    Figure  3.  Stress-strain curve of stiff (a) and soft (b) phases under quasi-static tensile

    图  (a) I-15~I-75及B-M试样准静态拉伸应力-应变曲线,(b) I-15~I-75及B-M试样断裂应变及拉伸强度

    Figure  4.  (a) Stress-strain curve of I-15–I-75 and B-M specimens under quasi-static tensile; (b) fracture strain and tensile strength of I-15–I-75 and B-M specimens

    图  (a) O-15~O-75及B-M试样准静态拉伸应力-应变曲线,(b) O-15~O-75及B-M试样断裂应变及拉伸强度

    Figure  5.  (a) Stress-strain curve of O-15–O-75 and B-M specimens under quasi-static tensile, and(b) fracture strain and tensile strength of O-15–O-75 and B-M specimens

    图  断裂前试样的能量吸收

    Figure  6.  Energy absorption of specimens before fracture

    图  试样的扫描电镜图像

    Figure  7.  SEM images of specimens

    表  1  模型及试件编号

    Table  1.   Number of models and specimens

    TypeUnit cellSpecimens
    B-M
    I-15
    I-30
    I-45
    I-60
    I-75
    O-15
    O-30
    O-45
    O-60
    O-75
    下载: 导出CSV
  • [1] 贾贤. 天然生物材料及其仿生工程材料[M]. 北京: 化学工业出版社, 2007.

    JIA X. Natural biomaterials and biomimetic engineering eaterials [M]. Beijing: Chemical Industry Press, 2007.
    [2] YARAGHI N A, KISAILUS D. Biomimetic structural materials: inspiration from design and assembly [J]. Annual Review of Physical Chemistry, 2018, 69: 23–57. doi: 10.1146/annurev-physchem-040215-112621
    [3] SONG F, SOH A K, BAI Y L. Structural and mechanical properties of the organic matrix layers of nacre [J]. Biomaterials, 2003, 24(20): 3623–3631. doi: 10.1016/S0142-9612(03)00215-1
    [4] GU G X, TAKAFFOLI M, BUEHLER M J. Hierarchically enhanced impact resistance of bioinspired composites [J]. Advanced Materials, 2017, 29(28): 1700060. doi: 10.1002/adma.201700060
    [5] 赵赫威, 郭林. 仿贝壳珍珠母层状复合材料的制备及应用 [J]. 科学通报, 2017, 62(6): 576–589.

    ZHAO H W, GUO L. Synthesis and applications of layered structural composites inspired by nacre [J]. Chinese Science Bulletin, 2017, 62(6): 576–589.
    [6] 赵建民, 麦康森, 张文兵, 等. 贝壳珍珠层及其仿生应用 [J]. 高技术通讯, 2003, 13(11): 94–98. doi: 10.3321/j.issn:1002-0470.2003.11.022

    ZHAO J M, MAI K S, ZHANG W B, et al. A review on characteristic of nacre and its potential for biomimetic applications [J]. High Technology Letters, 2003, 13(11): 94–98. doi: 10.3321/j.issn:1002-0470.2003.11.022
    [7] 李恒德, 冯庆玲, 崔福斋, 等. 贝壳珍珠层及仿生制备研究 [J]. 清华大学学报(自然科学版), 2001, 41(4/5): 41–47.

    LI H D, FENG Q L, CUI F Z, et al. Biomimetic research based on the study of nacre structure [J]. Journal of Tsinghua University (Science and Technology), 2001, 41(4/5): 41–47.
    [8] 梁艳, 赵杰, 王来, 等. 贝壳的力学性能和增韧机制 [J]. 机械强度, 2007, 29(3): 507–511. doi: 10.3321/j.issn:1001-9669.2007.03.031

    LIANG Y, ZHAO J, WANG L, et al. Mechanical properties and toughening mechanisms of mollusk shell [J]. Journal of Mechanical Strength, 2007, 29(3): 507–511. doi: 10.3321/j.issn:1001-9669.2007.03.031
    [9] MENIG R, MEYERS M H, MEYERS M A, et al. Quasi-static and dynamic mechanical response of Haliotis rufescens (abalone) shells [J]. Acta Materialia, 2000, 48(9): 2383–2398. doi: 10.1016/S1359-6454(99)00443-7
    [10] BARTHELAT F, LI C M, COMI C, et al. Mechanical properties of nacre constituents and their impact on mechanical performance [J]. Journal of Materials Research, 2006, 21(8): 1977–1986. doi: 10.1557/jmr.2006.0239
    [11] 张勋, 刘书海, 肖华平. 冷冻铸造技术制备仿贝壳层状结构陶瓷复合材料研究进展 [J]. 材料导报, 2017, 31(7): 99–112.

    ZHANG X, LIU S H, XIAO H P. Applying freeze-casting technique to the fabrication of nacre-like lamellar structured ceramic composites: a state-of-the-art review [J]. Materials Review, 2017, 31(7): 99–112.
    [12] 侯东芳, 周根树, 郑茂盛. 贝壳珍珠层断裂过程的原位观察及其增韧机制分析 [J]. 材料科学与工程学报, 2007, 25(3): 388–391. doi: 10.3969/j.issn.1673-2812.2007.03.016

    HOU D F, ZHOU G S, ZHENG M S. In situ SEM observation of crack propagation and analysis of the toughening mechanism in nacre [J]. Journal of Marerials Science & Engineering, 2007, 25(3): 388–391. doi: 10.3969/j.issn.1673-2812.2007.03.016
    [13] 马骁勇, 梁海弋, 王联凤. 三维打印贝壳仿生结构的力学性能 [J]. 科学通报, 2016, 61(7): 728–734.

    MA X Y, LIANG H Y, WANG L F. Multi-materials 3D printing application of shell biomimetic structure [J]. Chinese Science Bulletin, 2016, 61(7): 728–734.
    [14] BEGLEY M R, PHILIPS N R, COMPTON B G, et al. Micromechanical models to guide the development of synthetic ‘brick and mortar’ composites [J]. Journal of the Mechanics and Physics of Solids, 2012, 60(8): 1545–1560. doi: 10.1016/j.jmps.2012.03.002
  • 加载中
图(7) / 表(1)
计量
  • 文章访问数:  8678
  • HTML全文浏览量:  4082
  • PDF下载量:  65
出版历程
  • 收稿日期:  2019-04-29
  • 修回日期:  2019-05-11

目录

/

返回文章
返回