有限元法在金刚石合成中的应用进展

王健康 李尚升 宋艳玲 李露 于昆鹏 韩飞 宿太超 胡美华 吴玉敏

潘强, 张继春, 肖清华, 邹新宽, 石洪超. 动能弹对混凝土靶侵彻深度的PSO-SVM预测[J]. 高压物理学报, 2018, 32(2): 025102. doi: 10.11858/gywlxb.20170577
引用本文: 王健康, 李尚升, 宋艳玲, 李露, 于昆鹏, 韩飞, 宿太超, 胡美华, 吴玉敏. 有限元法在金刚石合成中的应用进展[J]. 高压物理学报, 2019, 33(1): 013101. doi: 10.11858/gywlxb.20180550
PAN Qiang, ZHANG Jichun, XIAO Qinghua, ZOU Xinkuan, SHI Hongchao. Prediction of Penetration Depth of Projectiles into Concrete Targets Based on PSO-SVM[J]. Chinese Journal of High Pressure Physics, 2018, 32(2): 025102. doi: 10.11858/gywlxb.20170577
Citation: WANG Jiankang, LI Shangsheng, SONG Yanling, LI Lu, YU Kunpeng, HAN Fei, SU Taichao, HU Meihua, WU Yumin. Progress in the Application of Finite Element Method in Synthetic Diamonds[J]. Chinese Journal of High Pressure Physics, 2019, 33(1): 013101. doi: 10.11858/gywlxb.20180550

有限元法在金刚石合成中的应用进展

doi: 10.11858/gywlxb.20180550
基金项目: 河南省自然科学基金(182300410279);河南省科技攻关项目(182102210311,172102210283);河南理工大学材料工程专业学位研究生专业实践示范基地(2016YJD03)
详细信息
    作者简介:

    王健康(1992-),男,硕士研究生,主要从事金刚石大单晶合成研究. E-mail:527719003@qq.com

    通讯作者:

    李尚升(1966-),男,博士,副教授,主要从事单晶金刚石研究. E-mail: lishsh@hpu.edu.cn

  • 中图分类号: O521.1; TB34

Progress in the Application of Finite Element Method in Synthetic Diamonds

  • 摘要: 金刚石以其优异的性能广泛应用于国防工程、机械加工、电子科技等领域,其需求量也日益增大。有限元法适用于复杂几何结构和物理问题的模拟分析,由此开辟了有限元法应用于金刚石合成和相应设备优化的新途径。阐述了有限元方法在六面顶压机及金刚石合成腔体工艺方面的应用进展。首先,考虑静力、应力强度、应力分布和形变等影响因素,对铰链梁和工作缸进行模拟分析,运用有限元法对顶锤的作用、破坏机理及新型顶锤设计进行探讨;其次,总结有限元法在金刚石腔体内的温度场、压力场、电学场等研究中的应用进展;最后,对有限元法在金刚石合成中的应用前景进行展望。

     

  • 超高分子量聚乙烯(ultra-high molecular weight polyethylene,UHMWPE)纤维复合材料层合板具有高强度、低质量密度等特点,在装甲防护领域得到了广泛应用[1-4]。UHMWPE作为防弹头盔、防弹衣等防护装备常用材料,主要通过分层、屈曲、纤维断裂、纤维拉伸耗散弹体的动能,达到抵抗子弹、破片侵彻的目的[5-6]。然而,目前尚未清晰阐明UHMWPE纤维复合材料在弹道冲击过程中的动态压缩机制[7]。因此,研究UHMWPE纤维复合材料层合板的动态压缩力学性能对防护材料抗弹性能优化及防护装备设计具有重要意义。

    许多学者尝试将UHMWPE层合板的力学性能与其抗冲击性能联系起来,以揭示层合板潜在的失效机制[8-11]。Lassig等[12]通过一系列准静态拉伸和剪切实验表征了UHMWPE材料的强度属性和冲击行为,推导了材料的本构参数并将其应用于仿真代码中,但该模型对于不同冲击场景及尺寸效应有待进一步研究。Zhu等[13]研究了UHMWPE材料在高应变率下的压缩行为,结果显示,材料厚度并不影响面外压缩方向的弹性模量。O’masta等[14]研究了铺层角度对UHMWPE板压缩性能的影响,发现样品的抗压强度随着夹角的减小而降低。孙颖等[15]研究了UHMWPE/乙烯基酯2.5D角联锁机织复合材料在不同应变率下的应力-应变曲线、压缩强度和能量吸收,结果表明,屈服应力、压缩模量随着应变率的增加而增大。Liu等[16]探究了UHMWPE纤维复合材料在面内压缩载荷下的动态变形机理,数值模拟显示,层间分层在中等应变水平下对复合材料峰值强度和压缩响应的影响较小。Sharifpour等[17]研究了UHMWPE层合板局部失效机制和早期裂纹微观力学评估。虽然高度各向异性纤维的纵向压缩性能已经得到了广泛的研究,但其横向力学性能有待进一步深入研究[18]。张羲黄等[19]通过压缩实验发现环氧树脂玻璃钢材料具有明显的应变率效应,且面内加载下材料的失效原因为分层损伤。潘月秀等[20]研究发现T800碳纤维/环氧树脂单向复合材料在动态压缩下发生脆性断裂,且具有明显的应变率强化效应和界面效应。叶卓然等[21]对比了UHMWPE纤维、芳纶纤维和碳纤维3种高性能纤维的力学性能,指出UHMWPE纤维在比强度、比模量等方面均优于其他纤维材料。目前,针对UHMWPE纤维复合材料层合板的研究主要通过实验和数值模拟分析来研究其抗弹性能[22-25]和失效模式[26-28],关于UHMWPE纤维复合材料层合板在高应变率下的压缩力学性能及微观失效机制研究鲜有报道。

    综上所述,研究UHMWPE纤维复合材料层合板的压缩性能并分析其失效模式对开发新型抗弹材料、建立材料在高应变率下的本构模型及失效准则具有重要意义。为此,本研究拟开展UHMWPE纤维复合材料层合板在静、动态压缩载荷下的力学性能实验,分析材料在不同应变率下的抗压强度和微观失效模式,以期为开发UHMWPE纤维复合材料本构模型和多尺度数值模拟提供依据。

    实验材料为湖南中泰特种装备有限公司生产的UHMWPE纤维复合材料层合板,层合板尺寸为400 mm×400 mm×5 mm,板材的面密度约为3.8 kg/m2。材料层间粘接强度较低,切割试样过大会使材料分层。为减小实验过程中样品尺寸的影响,静、动态实验样品均采用相同尺寸,样品形状对其力学响应的影响较小[29],结合霍普金森压杆(split Hopkinson pressure bar,SHPB)杆件形状,选用直径为10 mm、高度为5 mm的圆柱样品,样品结构及制备过程见图1图1(a)为样品的结构及加载方向,通过水射流切割成如图1(c)所示形状,再用小刀切除连接部分,得到实验所用样品,如图1(d)所示。

    图  1  样品的加工过程:(a) 样本结构,(b) 切割样品,(c) 加工后的样品, (d) 实验样品
    Figure  1.  Sample processing: (a) sample structure, (b) cutting sample, (c) post-processing samples, (d) experimental samples

    将样品进行分组后开展静态和动态压缩实验,每组实验重复3次,以保证结果的可重复性。实验装置如图2所示。通过华龙万能材料实验机(WDW-100C)进行静态压缩实验,加载速度设为2、5、10、20 mm/min,观察其压缩性能和失效状态,实验布置如图2(a)所示。

    图  2  实验装置:(a) 静态压缩实验装置,(b) SHPB实验装置
    Figure  2.  Experimental setup: (a) static compression experimental device, (b) SHPB experimental device

    应变率 ˙ε 通过加载速度v与样品的高度h计算得到

    ˙ε=vh
    (1)

    工程应力σe通过面外载荷F和原始样品的压缩面积A0计算得到,工程应变εe通过压缩位移Δh和样品的高度h计算得到,即

    {σe=F/A0εe=Δh/h
    (2)

    真实压缩应力σ和真实压缩应变ε[30]

    {σ=σe(1εe)ε=ln(1εe)
    (3)

    使用SHPB进行动态实验,杆的材料为35CrMnSi,撞击杆尺寸为16 mm×210 mm,撞击速度最高达80 m/s,入射杆、透射杆、缓冲杆的尺寸均为16 mm×1200 mm。为研究材料在不同动态应变率下的抗压性能,设定充气压力为0.25、0.45、0.65 MPa,具体实验布局如图2(b)所示。

    根据SHPB实验的一维应力波假定,样品两端面受力平衡,则

    εi+εr=εt
    (4)

    式中:εiεrεt分别为杆的入射、反射和透射应变。

    工程应力和工程应变可由反射波和透射波计算得到,将所得结果代入式(3)即可得到真实应力和真实应变,此时工程应力、工程应变、应变率可分别表示为[30]

    {σe=AEA0εtεe=2c0hh0εrdτ˙ε=2c0hεr
    (5)

    式中:A为当前样品当前的压缩面积,E为杆的弹性模量,c0为压杆中应力波的传播速度。

    图3所示,采用日立高新技术公司的SU3500型扫描电子显微镜(scanning electron microscopy,SEM)对实验后的样品进行微观形貌观察以及失效模式分析。由于样品材料不导电,因此在材料厚度方向进行喷金处理。

    图  3  扫描电子显微镜(SU3500)
    Figure  3.  Scanning electron microscope (SU3500)

    每种加载速度开展3次重复实验,并记录样品在加载过程中的变形情况。4种静态压缩速率(6.7×10−3、1.7×10−2、3.3×10−2和6.7×10−2 s−1)下材料的真实应力-真实应变曲线如图4所示,其中曲线为3次重复实验的平均值,误差带为各应力的标准差。

    图  4  准静态实验结果
    Figure  4.  Quasi-static experimental results

    图4可知,样品在压缩过程中呈现先凹后凸的变化趋势。当应变率在6.7×10−3~6.7×10−2 s−1区间时,UHMWPE纤维复合材料层合板未表现出应变率效应;当应变率为6.7×10−3 s−1时,样品出现应力滞后现象,可能是较低的加载速率导致样品的失效模式发生了变化。根据应力-应变曲线的变化趋势,可将其分为两个阶段,通过曲线求导近似,以紫色线作为分界线。在图4紫色线下方的第Ⅰ阶段,由于材料在制备过程中不可避免地存在缝隙和空腔[31],在样品压缩过程中这些缝隙和空腔闭合,导致应力-应变曲线缓慢上升,呈凹形上升现象;在紫色线上方第Ⅱ阶段,此时样品被压实,压缩作用下样品内部产生间接张力,最终单个铺层在纵向发生位移,出现片状挤出现象,曲线呈现凸形平台。

    将实验视频进行图像截取,以观察样品失效形貌。由图4图5可知,当应力小于150 MPa时,样品处于逐渐压紧状态,应力-应变曲线增长得较快,4种应变率下均满足该变化规律,侧面反映出实验设备和所制样品稳定可靠。继续加载样品时,单个铺层会从纵向挤出,应力增长速度逐渐变缓,随着应变率增加,曲线斜率并无明显变化,应变率效应并不显著,与文献[32]的结论一致。

    图  5  静态压缩下样品的变形
    Figure  5.  Deformation of samples under static compression

    图5可知,当应变率为6.7×10−3 s−1时,样品逐渐压紧,其表面从光滑变为粗糙,反光性能变差,纤维层逐渐从样品中挤出,即材料在压缩载荷下达到图4曲线的极限点,该点的应力值为材料在面外压缩方向的压缩强度。在面外压缩载荷下,UHMWPE纤维复合材料层合板因其正交结构出现了面内挤出现象,该现象是材料纤维与基体之间的层间粘接强度较低导致的,通过层间位移达到材料泄压。何业茂等[33]通过实验发现方形样品在面外压缩下出现的分层破坏现象是层合板在压缩载荷下纤维层间的剪切力大于其剪切强度导致的。本实验中,圆柱形样品也出现了分层破坏,且纤维层被大量挤出,其原因是面外压缩导致纤维轴向面内拉伸,进而使层状材料被压出。

    动态压缩实验设定的加载气压为0.25、0.45、0.65 MPa,对应的应变率分别为2.05×103、4.08×103和5.27×103 s−1。从图6(a)~图6(c)中可知,3种高应变率下材料的力学响应变化趋势一致,均具有线弹性阶段,随后出现微屈服段,最后压力上升直至断裂失效。曲线中的峰值应力为材料在相应应变率下的压缩强度。图6(d)表明其压缩强度呈现明显的应变率相关性。3种工况下,材料曲线均出现二次屈曲现象,第一次材料屈曲为材料内部间隙压实,第二次屈曲为材料面外压缩诱导的面内位错分层。当应变率为5.27×103 s−1时,材料曲线呈现3段屈曲,并且材料屈曲与其分层相关。动态压缩下的力学响应情况与静态压缩下的材料失效明显不同,高应变率下材料失效主要为分层和纤维断裂。由图6(e)材料不同方向的外部形貌可知,随着应变率的增大,样品周围的纤维碎屑逐渐增多且分层现象显著。当应变率为2.05×103 s−1时,样品的主要破坏形态为铺层的局部滑移和纤维扭结,并未出现纵向分层,所以其失效应变较小。随着应变率的增加,图6(e)中材料的分层现象更加明显。

    图  6  动态实验结果: (a) 0.25 MPa,(b) 0.45 MPa,(c) 0.65 MPa,(d)真实应力-应变曲线,(e) 动态压缩下样品的变形
    Figure  6.  Dynamic experimental results: (a) 0.25 MPa, (b) 0.45 MPa, (c) 0.65 MPa, (d) true stress-strain curves, (e) deformation of sample under dynamic compression

    表1列出了不同应变率下的实验结果,参考文献[34]的方法,将6.7×10−3 s−1应变率下的压缩强度作为静态压缩基准。动态压缩载荷下,随着应变率的增加,动态抗压强度和动态增强因子(dynamic increase factor,DIF)均呈现明显的动态增强效应。由于动态冲击时间短,材料缓冲位移小,因此材料受压后纤维铺层产生间接张力,进而形成分层来耗散能量。动态载荷下3种工况均出现了分层开裂和滑移现象,并且应变率越高,材料分层效应越明显,纤维崩裂的碎屑越多,材料的可变形位移越大,动态峰值应变也越大。

    表  1  动态压缩实验结果
    Table  1.  Dynamic compression experimental results
    Static compressive
    strength/MPa
    Loading air
    pressure/MPa
    Strain rate/s−1Dynamic compressive
    strength/MPa
    DIFDynamic peak
    strain
    424.520.252.05×103334.260.78100
    424.520.454.08×103499.971.18210
    424.520.655.27×103649.771.53220
    Note: (1) Static compressive strength is the maximum stress at a loading speed of 2 mm/min;
    (2) Dynamic compressive strength is the peak stress achieved by the sample under dynamic compressive loading;
    (3) DIF is the dynamic compressive strength divided by the static compressive strength;
    (4) Dynamic peak strain is the strain value corresponding to the dynamic compressive strength.
    下载: 导出CSV 
    | 显示表格

    图7显示了静、动态压缩(7种应变率)下材料的压缩强度和峰值应变。由图7(a)可知:4种静态应变率加载下,材料的压缩强度和峰值应变相近,无明显的应变率效应;但应变率为6.7×10−3 s−1与应变率为1.7×10−2~6.7×10−2 s−1相比时,实验曲线出现应力滞后现象,可能是加载速率变化过大导致应力波在材料内部的传递形式发生变化,从而导致失效模式不同;当应变率为6.7×10−3 s−1时,峰值应变为0.46,比1.7×10−2~6.7×10−2 s−1应变率下样品的峰值应变高4%~7%。因此,在应变相同的情况下,应变率为6.7×10−3 s−1时的压缩强度更低。动态压缩下样品的压缩强度和峰值应变均具有应变率相关性。在2.05×103、4.08×103、5.27×103 s−1 3种动态速率压缩下,材料的峰值应变显著低于4种静态加载工况,这是由于高应变率压缩下材料没有足够多的时间传递载荷和变形,导致材料局部发生分层失效。静态和动态加载下,材料压缩强度均逐渐升高。但是当应变率为2.05×103 s−1时,压缩强度显著低于静态加载工况,原因可能是0.25 MPa气压不足以使样品产生滑移分层,仅有局部的压实段和少量的纤维挤出,与准静态压缩样品相比,应变相差4倍以上,故其压缩强度较低。高应变率加载下,样品的峰值应变逐步增大,表明该材料具有应变率强化效应。

    图  7  不同应变率下的压缩强度和峰值应变:(a) 压缩强度,(b) 峰值应变
    Figure  7.  Compressive strength and peak strain at different strain rates: (a) compressive strength, (b) peak strain

    图8显示了静态压缩实验后样品的微观形貌,其中:图8(a)~图8(d)为样品在4种应变率下的整体失效形貌,图8(e)~图8(f)为典型的局部失效形貌。由图8可知,材料的整体压缩响应表现为铺层挤出,根据应变速率的不同,铺层挤出程度也不尽相同。当应变率为6.7×10−3 s−1时,样品压缩端和固定端的铺层挤出程度相对均匀,在持续加载下材料内部逐渐压实,最终出现铺层的大幅挤出,导致图4所示曲线出现屈服。其原因是在加载过程中压力和变形有足够多的时间在整个样品中传递,样品厚度方向均发生纤维层纵向位移,因此该工况下材料的峰值应变较大,压缩强度最小。当应变率为1.7×10−2、3.3×10−2、6.7×10−2 s−1时,样品单侧(压缩端)铺层挤出层数多,且3种工况下微观尺度的铺层变形形貌相同,因此图4应力-应变曲线变化趋势基本一致。由于加载速率提高,厚度方向的加载使样品非均匀受力,加载端纤维层的纵向位移较大,致使峰值应变较小。在微观层面,4种静态应变率下样品的纤维失效模式相同,均为面外加载下纤维铺层发生纵向位移。由于加载速率的差异,铺层位移存在差异。当应变率为6.7×10−3 s−1时,由于其加载速率较其余3种工况低,样品厚度方向铺层的纵向位移较均匀,铺层纤维未发生明显的拉伸。在1.7×10−2、3.3×10−2和6.7×10−2 s−1 3种工况下,纵向位移的铺层均发生较严重的纤维拉伸,因此其压缩强度比应变率为6.7×10−3 s−1时大。从图8(e)和图8(f)中的纤维尺度来看,横向的压缩作用使铺层在纤维方向上产生间接拉伸作用,层间剪切力大于层压板的剪切强度,纤维被逐渐拉长直至断裂。由此可见,UHMWPE纤维复合材料层合板在横向压缩过程中会诱导材料在纤维方向产生拉伸载荷。从图8(f)中纤维逐渐变细甚至断裂可知,在面外压缩载荷下材料主要通过纤维断裂承载冲击。Attwood等[35]也通过理论和数值模拟解释了这种面外压缩产生的材料纤维轴向位移形成的剪切滞后区。

    图  8  静态加载后UHMWPE纤维复合材料层合板样品的微观形貌
    Figure  8.  Microstructure of UHMWPE fiber composite laminate samples after static loading

    图9显示了动态加载后样品的微观形貌,其中:图9(a)~图9(c)显示了样品在3种应变率下发生的分层、位错滑移形貌。如图9(a)所示,当应变率为2.05×103 s−1时,样品在其高度方向上只有一层发生了明显的铺层滑移现象,且为局部失效,因此其压缩强度最低。随着加载气压的升高,如图9(b)~图9(c)所示,分层现象更严重,主要表现为材料的位错滑移,端部纤维形貌与静态载荷下的拉伸状态不同,发生了缠绕和扭结,主要原因是动态载荷下除了压缩诱导的层间剪切外,反射波加剧了材料分层破坏,进一步揭示了材料在动态载荷下的强度增强效应。与静态加载不同,在瞬间冲击作用下材料变形将由等温变形转变为绝热变形,产生的温度效应将会影响材料的失效形式。当应变率为2.05×103 s−1时,样品的峰值应变显著小于其余2种工况,因此其失效形貌(见图9(a))也与其他2种工况(见图9(b)和图9(c))存在较大差异,较小的峰值应变也导致其压缩强度较低。当应变率为5.27×103 s−1时,如图9(c)所示,样品中不仅出现较多的材料分层,还发生明显的纤维缠绕,因此其压缩强度最高。从图9(d)和图9(f)的局部视图中可以发现,铺层发生脱粘和撕裂现象,与静态压缩实验的铺层挤出现象不同,铺层未发生明显卷曲。随着应变率的增加,材料的分层数逐渐增多,纵向分层位移也逐步增大。当应变率为2.05×103 s−1时,只有少量纤维挤出,对应的图6(a)曲线第Ⅱ阶段较为平滑;当应变率为4.08×103 s−1时,材料出现2处分层,对应的图6(b)曲线双段屈曲效应明显;当应变率为5.27×103 s−1时,应力-应变曲线近似呈3段屈曲现象,如图6(c)所示。通过材料分层数可知,应力-应变曲线的多次屈曲是由滑移分层导致的,可通过改进纤维和基体属性调控其力学性能。

    图  9  动态加载后UHMWPE纤维复合材料层合板样品的微观形貌
    Figure  9.  Microstructure of UHMWPE fiber composite laminate samples after dynamic loading

    (1) UHMWPE纤维复合材料层合板在静态压缩载荷下无应变率效应,真实应力-应变曲线呈现先凹后凸的变化趋势。当应变率为6.7×10−3 s−1时,层合板有足够长的响应时间,样品两端均呈现纤维铺层挤出,应力-应变曲线出现应力滞后现象;当应变率为1.7×10−2、3.3×10−2、6.7×10−2 s−1时,样品均为单侧纤维挤出,压缩强度和峰值应变接近,无明显的应变率效应。

    (2) UHMWPE纤维复合材料层合板动态压缩载荷存在应变率效应,其应力-应变曲线呈现“多屈曲”效应,屈曲状态与材料的纵向位错分层相联系,压缩强度随应变率的增加而升高。

    (3) UHMWPE纤维复合材料层合板横向压缩载荷在静、动态下呈现不同的失效模式:静态压缩载荷下材料的主要力学响应为纤维拉伸、断裂和铺层滑移;动态载荷下除了压缩诱导的层间剪切外,反射波也会导致层裂破坏,材料失效主要表现为纤维缠绕、扭结、铺层脱粘以及材料的纵向位错滑移,宏观表现为大开口的分层破坏。

  • 图  铰链梁设计变量[15]

    Figure  1.  Design parameters of hinge beam[15]

    图  碳化钨顶锤工程图[27]

    Figure  2.  Engineering drawing for the WC anvil[27]

    图  锤面温度分别为150、220和300 ℃的热平衡状态下的顶锤温度分布[32]

    Figure  3.  Temperature distributions on the anvil under thermal equilibrium state when the initial temperatures of the anvil are (a) 150, (b) 220 and (c) 300 ℃, respectively[32]

  • [1] 罗中平, 姚亚林, 夏志华, 等. 磁性对人造金刚石单晶机械性能的影响 [J]. 高压物理学报, 2009, 23(3): 236–240 doi: 10.3969/j.issn.1000-5773.2009.03.012

    LUO Z P, YAO Y L, XIA Z H, et al. Study on the influence of magnetism of synthetic diamond single crystal on their mechanical performances [J]. Chinese Journal of High Pressure Physics, 2009, 23(3): 236–240 doi: 10.3969/j.issn.1000-5773.2009.03.012
    [2] 张文凯, 彭放, 郭振堂, 等. 高压烧结镀Cr、Ti膜金刚石/铜复合材料热导率研究 [J]. 高压物理学报, 2012, 26(3): 306–312 doi: 10.11858/gywlxb.2012.03.010

    ZHANG W K, PENG F, GUO Z T, et al. Research on thermal conductivity of diamond with Cr, Ti coating/copper composite materials by sintering under high pressure [J]. Chinese Journal of High Pressure Physics, 2012, 26(3): 306–312 doi: 10.11858/gywlxb.2012.03.010
    [3] 曹海要, 战再吉. 铜/金刚石复合材料电磁轨道烧蚀特性的实验研究 [J]. 高压物理学报, 2016, 30(4): 317–322 doi: 10.11858/gywlxb.2016.04.008

    CAO H Y, ZHAN Z J. Experimental study of Cu/diamond composite electromagnetic rail ablation characteristics [J]. Chinese Journal of High Pressure Physics, 2016, 30(4): 317–322 doi: 10.11858/gywlxb.2016.04.008
    [4] 陆振勤, 胡毅胜, 郭滇生, 等. 金刚石合成设备的技术进展 [J]. 超硬材料工程, 2011, 23(3): 23–25 doi: 10.3969/j.issn.1673-1433.2011.03.006

    LU Z Q, HU Y S, GUO D S, et al. Technical progress of diamond synthesis equipment [J]. Superhard Material Engineering, 2011, 23(3): 23–25 doi: 10.3969/j.issn.1673-1433.2011.03.006
    [5] 张聪. 有限单元法在六面顶压机结构分析及优化设计中的应用 [D]. 长春: 吉林大学, 2010: 18-23.
    [6] 王勖成. 有限单元法基本原理和数值方法 [M]. 2版. 北京: 清华大学出版社, 1997: 1-2.
    [7] 黄志新. ANSYS Workbench 14.0超级学习手册 [M]. 北京: 人民邮电出版社, 2013: 1-3.
    [8] 韩奇钢, 马红安, 李瑞, 等. 基于有限元法确立碳化钨顶锤的破裂判据 [J]. 高压物理学报, 2010, 24(1): 1–5 doi: 10.3969/j.issn.1000-5773.2010.01.001

    HAN Q G, MA H A, LI R, et al. The criterion for crack of tungsten carbide anvil based on finite element method [J]. Chinese Journal of High Pressure Physics, 2010, 24(1): 1–5 doi: 10.3969/j.issn.1000-5773.2010.01.001
    [9] 韩奇钢, 马红安, 肖宏宇, 等. 基于有限元法分析宝石级金刚石的合成腔体温度场 [J]. 物理学报, 2010, 59(3): 1923–1927

    HAN Q G, MA H A, XIAO H Y, et al. Finite element method study on the temperature distribution in the cell of large single crystal diamond [J]. Acta Physica Sinica, 2010, 59(3): 1923–1927
    [10] 韩奇钢. 人造金刚石的制备方法及其超高压技术 [J]. 高压物理学报, 2015, 29(4): 313–320 doi: 10.11858/gywlxb.2015.04.012

    HAN Q G. Preparation methods and ultra-high pressure technologies of synthetic diamonds [J]. Chinese Journal of High Pressure Physics, 2015, 29(4): 313–320 doi: 10.11858/gywlxb.2015.04.012
    [11] 王海阔, 贺端威, 许超, 等. 基于国产铰链式六面顶压机的大腔体静高压技术研究进展 [J]. 高压物理学报, 2013, 27(5): 633–661 doi: 10.11858/gywlxb.2013.05.001

    WANG H K, HE D W, XU C, et al. Development of large volume-high static pressure techniques based on the hinge-type cubic presses [J]. Chinese Journal of High Pressure Physics, 2013, 27(5): 633–661 doi: 10.11858/gywlxb.2013.05.001
    [12] 王良文, 刘建设, 汪曙光, 等. UDS-Ⅲ型金刚石压机结构的有限元分析 [J]. 机械设计, 2005, 22(3): 38–41 doi: 10.3969/j.issn.1001-2354.2005.03.014

    WANG L W, LIU J S, WANG S G, et al. Finite element analysis on the structures of UDS-Ⅲ typed diamond press [J]. Journal of Machine Design, 2005, 22(3): 38–41 doi: 10.3969/j.issn.1001-2354.2005.03.014
    [13] 刘竹丽, 梁颖, 秦东晨. 基于有限元法的铰链梁疲劳强度分析与设计研究 [J]. 金刚石与磨料磨具工程, 2006(4): 38–42 doi: 10.3969/j.issn.1006-852X.2006.04.012

    LIU Z L, LIANG Y, QIN D C. Investigation on fatigue strength and design of hinge sleeve based on FEM [J]. Diamond & Abrasives Engineering, 2006(4): 38–42 doi: 10.3969/j.issn.1006-852X.2006.04.012
    [14] 张聪, 马红安, 韩奇钢, 等. 高压下国产六面顶压机铰链梁和工作缸的应力分析 [J]. 高压物理学报, 2010, 24(5): 321–325 doi: 10.11858/gywlxb.2010.05.001

    ZHANG C, MA H A, HAN Q G, et al. Stress analysis on hinge sleeve and cylinder of china-type cubic-anvil high pressure apparatus under high pressure [J]. Chinese Journal of High Pressure Physics, 2010, 24(5): 321–325 doi: 10.11858/gywlxb.2010.05.001
    [15] 汪曙光. 基于Ansys Workbench六面顶压机优化分析 [J]. 金刚石与磨料磨具工程, 2015(2): 69–73

    WANG S G. Optimization analysis of cubic diamond press based on Ansys Workbench [J]. Diamond & Abrasives Engineering, 2015(2): 69–73
    [16] 牛炜霖. 聚晶金刚石取放机械手设计及压机顶锤应力分析 [D]. 郑州: 郑州大学, 2014: 15-17.
    [17] 韩奇钢. 有限元法在六面顶顶锤优化设计中的应用 [D]. 长春: 吉林大学, 2008: 22-24.
    [18] HAN Q G, LI M Z, JIA X P, et al. Finite element analysis and design of cubic high-pressure anvils based on the principle of lateral support [J]. Review of Scientific Instruments, 2010, 81(12): 123901. doi: 10.1063/1.3514091
    [19] 韩奇钢, 马红安, 李瑞, 等. 六面顶压机硬质合金顶锤应力与破裂机理的有限元分析 [J]. 重型机械, 2007(3): 27–31 doi: 10.3969/j.issn.1001-196X.2007.03.008

    HAN Q G, MA H A, LI R, et al. The stress and fracture mechanism analysis on cemented carbide anvils of cubic apparatus based on FEM [J]. Heavy Machinery, 2007(3): 27–31 doi: 10.3969/j.issn.1001-196X.2007.03.008
    [20] 韩奇钢, 贾晓鹏, 马红安, 等. 基于三维有限元法模拟分析六面顶顶锤的热应力 [J]. 物理学报, 2009, 58(7): 4812–4816 doi: 10.3321/j.issn:1000-3290.2009.07.068

    HAN Q G, JIA X P, MA H A, et al. Finite element simulations of thermal-stress on cemented tungsten carbide anvil used in cubic high pressure apparatus [J]. Acta Physica Sinica, 2009, 58(7): 4812–4816 doi: 10.3321/j.issn:1000-3290.2009.07.068
    [21] HAN Q G, ZHANG Q, LI M Z, et al. An effective solution for the best set of beveling parameters of the cubic high-pressure tungsten carbide anvil [J]. Chinese Physics Letters, 2012, 29(11): 116201. doi: 10.1088/0256-307X/29/11/116201
    [22] ADAMS D M, SHAW A C. A computer-aided design study of the behaviour of diamond anvils under stress [J]. Journal of Physics D, 1982, 15(9): 1609–1635. doi: 10.1088/0022-3727/15/9/006
    [23] 于歌, 韩奇钢, 李明哲, 等. 新型圆角式高压碳化钨硬质合金顶锤的有限元分析 [J]. 物理学报, 2012, 61(4): 63–67

    YU G, HAN Q G, LI M Z, et al. Finite element analysis of the high-pressure tungsten carbide radius-anvil [J]. Acta Physica Sinica, 2012, 61(4): 63–67
    [24] HAN Q G, YANG W K, ZHU P W, et al. Design and performance of tapered cubic anvil used for achieving higher pressure and larger sample cell [J]. Review of Scientific Instruments, 2013, 84(7): 073902. doi: 10.1063/1.4812351
    [25] MOSS W C, HALLQUIST J O, REICHLIN R, et al. Finite element analysis of the diamond anvil cell: Achieving 4.6 Mbar [J]. Applied Physics Letters, 1986, 48(19): 1258–1260. doi: 10.1063/1.96996
    [26] HAN Q G, JIA X P, QIN J M, et al. FEM study on a double-beveled anvil and its application to synthetic diamonds [J]. High Pressure Research, 2009, 29(3): 449–456. doi: 10.1080/08957950903109835
    [27] HAN Q G, LI M Z, JIA X P, et al. Modeling of effective design of high pressure anvils used for large scale commercial production of gem quality large single crystal diamond [J]. Diamond and Related Materials, 2011, 20(7): 969–973. doi: 10.1016/j.diamond.2011.05.017
    [28] HAN Q G, MA H A, HUANG G F, et al. Hybrid-anvil: a suitable anvil for large volume cubic high pressure apparatus [J]. Review of Scientific Instruments, 2009, 80(9): 096107. doi: 10.1063/1.3227239
    [29] KONDRAT’YEV A I, VOHRA Y K. Finite-element modeling of stresses and strains in a diamond anvil cell device: case of a diamond-coated rhenium gasket [J]. High Pressure Research, 2007, 27(3): 321–331. doi: 10.1080/08957950701557573
    [30] KONDRAT Y A, MURPHY M J, WEIR S T, et al. Comparison of finite element modeling and experimental pressure distribution in a diamond anvil cell [C]//APS Southeastern Section Meeting, 2002.
    [31] 刘芳, 刘永刚, 谢鸿森. 六面顶压机高压腔体温度场的数值模拟 [J]. 高压物理学报, 2012, 26(2): 135–140 doi: 10.11858/gywlxb.2012.02.003

    LIU F, LIU Y G, XIE H S. Numerical simulation of temperature field in sample assembly of cubic press [J]. Chinese Journal of High Pressure Physics, 2012, 26(2): 135–140 doi: 10.11858/gywlxb.2012.02.003
    [32] 李瑞, 马红安, 韩奇钢, 等. 基于ANSYS的金刚石合成腔内的温度场分析 [J]. 吉林大学学报(工学版), 2008, 38(3): 535–538

    LI R, MA H A, HAN Q G, et al. Analysis of temperature field in diamond synthetic chamber based on ANSYS [J]. Journal of Jilin (University Engineering and Technology Edition), 2008, 38(3): 535–538
    [33] LI Z C, JIA X P, HUANG G F, et al. FEM simulations and experimental studies of the temperature field in a large diamond crystal growth cell [J]. Chinese Physics B, 2013, 22(1): 363–367.
    [34] ZHANG Q, LI R, GU X, et al. Thermal analysis of the growth process of synthetic diamond in the large volume cubic press apparatus with large deformation of high pressure cell [J]. Journal of Crystal Growth, 2015, 420: 80–83. doi: 10.1016/j.jcrysgro.2015.03.036
    [35] 卢金斌, 穆云超, 王秦生, 等. 基于ø69 mm腔体温度场数值模拟的高品级金刚石合成工艺优化 [J]. 金刚石与磨料磨具工程, 2010, 30(6): 50–53 doi: 10.3969/j.issn.1006-852X.2010.06.011

    LU J B, MU Y C, WANG Q S, et al. Optimization of diamond synthesis process based on numerical simulation on temperature field in ø69 chamber [J]. Diamond & Abrasives Engineering, 2010, 30(6): 50–53 doi: 10.3969/j.issn.1006-852X.2010.06.011
    [36] 刘衍聪, 郭欣, 伊鹏, 等. 六面顶压机合成腔中心温度轨迹控制方法研究 [J]. 人工晶体学报, 2011, 40(1): 27–32 doi: 10.3969/j.issn.1000-985X.2011.01.006

    LIU Y C, GUO X, YI P, et al. Study on temperature tracking control approach of cubic press center [J]. Journal of Synthetic Crystals, 2011, 40(1): 27–32 doi: 10.3969/j.issn.1000-985X.2011.01.006
    [37] LI R, MA H A, HAN Q G, et al. Simulation of pressure distribution in a pyrophyllite high-pressure cell by finite-element analysis [J]. High Pressure Research, 2007, 27(2): 249–257. doi: 10.1080/08957950701385819
    [38] LI R, XU B, ZHANG Q, et al. Finite-element analysis on pressure transfer mechanism in large-volume cubic press [J]. High Pressure Research, 2016, 36(4): 575–584. doi: 10.1080/08957959.2016.1238915
    [39] LI R, MA H A, YIN B H, et al. Finite element analysis on pressure-transmitting behavior of pyrophyllite block based on ANSYS/LS-DYNA [J]. Journal of Jilin University, 2008, 38(2): 292–297.
    [40] 李瑞, 马红安, 臧传义, 等. 有限元法计算模拟金刚石合成腔内的电学场 [J]. 金刚石与磨料磨具工程, 2006(3): 4–7 doi: 10.3969/j.issn.1006-852X.2006.03.002

    LI R, MA H A, ZANG C Y, et al. Computer simulation of the electric field in diamond synthetic chamber using FEM [J]. Diamond & Abrasives Engineering, 2006(3): 4–7 doi: 10.3969/j.issn.1006-852X.2006.03.002
    [41] 肖宏宇, 臧传义, 马红安, 等. 碳素扩散场对塔状晶体合成的影响 [J]. 超硬材料工程, 2008, 20(6): 1–4 doi: 10.3969/j.issn.1673-1433.2008.06.001

    XIAO H Y, ZANG C Y, MA H A, et al. Effects of carbon diffusion field on the synthesis of tower shape diamond crystals [J]. Superhard Material Engineering, 2008, 20(6): 1–4 doi: 10.3969/j.issn.1673-1433.2008.06.001
  • 加载中
图(3)
计量
  • 文章访问数:  8250
  • HTML全文浏览量:  3565
  • PDF下载量:  41
出版历程
  • 收稿日期:  2018-05-02
  • 修回日期:  2018-05-28

目录

/

返回文章
返回