
Citation: | XIU Chengdong, WANG Changfeng, LI Bing, GUAN Renguo. Numerical Simulation Study on the Influence of Hard Phase Shape on the Fracture Behavior of Ti-Al3Ti Bionic Composites[J]. Chinese Journal of High Pressure Physics, 2023, 37(4): 044201. doi: 10.11858/gywlxb.20230629 |
复合材料板具有较高的强度质量比、良好的耐腐蚀性和优异的可设计性,被广泛应用于航空航天和工业制造等领域[1]。在实际使用中复合材料板经常受到不同形式的冲击荷载,从而产生振动和屈曲问题,因此冲击载荷作用下复合材料板的动力稳定性问题备受关注。
近年来关于复合材料板的研究越来越多,尤其是冲击荷载作用下复合材料板的动力稳定性问题研究[2],对工程部件结构设计和使用具有重要的意义。Sun等[3]研究了在加热环境中应力波对功能梯度圆柱壳轴向冲击屈曲的影响;毛柳伟等[4]对弹性直杆在应力波作用下的动力分叉屈曲进行了分析与探讨,提出了求解应力波作用下直杆动力屈曲的数值方法;Lepik[5]讨论了在应力波影响下轴向压缩的弹塑性梁的屈曲;Zhang等[6]分析了不确定初始几何缺陷对薄板屈曲的影响;Abdelaziz等[7]利用双曲线剪切变形理论,分析了在各种边界条件下复合材料板的弯曲变形和屈曲;Kouchakzadeh等[8]采用线性和旋转弹簧的均匀分布来模拟边界条件,对矩形层压复合板的屈曲进行了分析;Czapski等[9]通过数值和实验方法,研究了残余应力对压缩至破坏期间薄壁层压板屈曲性能的影响。
在实际工程中,复合材料板多应用于振动环境,其在动力响应下的动态特性和振动分析是必不可少的,因而对该类材料的振动屈曲研究至关重要。Kuo[10]研究了两种非均匀分布纤维复合材料板的振动屈曲问题,Villarreal等[11]对典型正交异性板的本征频率和振动屈曲进行了理论分析,Eftekhari等[12]提出了通过组合应用有限元方法和微分正交方法求解矩形板的振动屈曲问题,Rehman等[13]探讨了壳体结构的缺陷和损坏对结构振动屈曲的影响,Sayyad等[14]将三角剪切变形理论应用于复合板的变形和振动屈曲研究。
关于复合材料板的振动屈曲问题已开展了较多的研究,但大多未考虑应力波效应对振动屈曲的影响,而动力屈曲一般与应力波相联系且具有局部发生的特点,研究含初始缺陷的复合材料板能更好地揭示实际工程中复合材料板在不同工况下发生动力屈曲的机理。基于此,本研究利用Kirchhoff薄板理论和Hamilton原理,考虑应力波效应,建立含初始几何缺陷的四边简支复合材料板的振动控制方程,得到板的屈曲临界荷载表达式,在此基础上通过数值计算讨论初始几何缺陷、振型函数初相位、铺层角度、屈曲模态阶数和铺层层数对复合材料板振动屈曲临界荷载的影响,为工程实际提供理论依据。
复合材料板在x = La处为固定边界条件,其余3边为简支边界条件,在
根据Kirchhoff薄板理论及经典弹性理论,复合材料板的位移与应变、弹性曲面的曲率和扭率的表达式为
{w=w1+w0u=u0−z∂(w−w1)∂xv=v0−z∂(w−w1)∂y | (1) |
{ε0x=∂u0∂xε0y=∂v0∂yγ0xy=∂u0∂y+∂v0∂x | (2) |
{κx=−∂2(w−w1)∂x2κy=−∂2(w−w1)∂y2κxy=−2∂2(w−w1)∂x∂y | (3) |
{εx=ε0x+zκxεy=ε0y+zκyγxy=γ0xy+zκxy | (4) |
式中:u0、v0、w0分别为复合材料板在x、y、
板的内力(Nx、Ny、Nxy)与内力矩(Mx、My、Mxy)之间的关系为
[NxNyNxyMxMyMxy]=[A11A12A16B11B12B16A12A22A26B12B22B26A16A26A66B16B26B66B11B12B16D11D12D16B12B22B26D12D22D26B16B26B66D16D26D66][ε0xε0yγ0xyκxκyκxy] | (5) |
式中:
{Aij=n∑k=1¯Qkij(hk−hk−1)Bij=12n∑k=1¯Qkij(h2k−h2k−1)Dij=13n∑k=1¯Qkij(h3k−h3k−1) | (6) |
式中:
¯Q=P−1Q(P−1)T | (7) |
式中:P为坐标转换矩阵,Q为刚度矩阵。
P=[cos2θsin2θ2sinθcosθsin2θcos2θ−2sinθcosθ−sinθcosθsinθcosθcos2θ−sin2θ] | (8) |
Q=[Q11Q120Q12Q22000Q66] | (9) |
考虑材料为正交各向异性材料,设E1、E2、G12、μ12、μ21分别为板材料x和y方向的拉压弹性模量、剪切弹性模量、主泊松比和副泊松比,则有
Q11=E11−μ12μ21,Q22=E21−μ12μ21,Q12=μ12E21−μ12μ21=μ21E11−μ12μ21,Q66=G12 | (10) |
复合材料板在左端受面内冲击荷载N作用(见图1)时,应力波沿x方向在板内传播,其应力变化如图2所示。
当应力波传播至波阵面位置Lcr(临界长度)时,板发生振动屈曲,板的内力Nt和应力波波速c分别表示为
Nt={N0⩽x⩽Lcr0x>Lcr,c=√A11ρh | (11) |
板发生振动屈曲时的变形能可以表示为
U=12∫Lcr0∫Lb0(Nxε0x+Nyε0y+Nxyγ0xy+Mxκx+Myκy+Mxyκxy)dxdy=12∫Lcr0∫Lb0[Nx∂u0∂x+Ny∂v0∂y+Nxy(∂u0∂x+∂v0∂y)+Mx(−∂2w0∂x2)+My(−∂2w0∂y2)+Mxy(−2∂2w0∂x∂y)]dxdy | (12) |
发生振动屈曲时的动能(考虑转动惯量)可以表示为
T=12∫h/2−h/2∫Lcr0∫Lb0ρ(k)[(∂u∂t)2+(∂v∂t)2+(∂w0∂t)2]dxdydz=12∫Lcr0∫Lb0[I0(∂u0∂t)2+I0(∂v0∂t)2+I0(∂w0∂t)2−2I1∂u0∂t∂2w0∂x∂t−2I1∂v0∂t∂2w0∂y∂t+I2(∂2w0∂x∂t)2+I2(∂2w0∂y∂t)2]dxdy | (13) |
式中:
板发生振动屈曲时的外力功可以表示为
W=12∫Lcr0∫Lb0Nt(∂w∂x)2dxdy | (14) |
考虑Hamilton变分原理,即
δ∫t2t1(T−U+W)dt=0 | (15) |
将式(12)~式(14)代入式(15)并进行变分计算,可得
∂Nx∂x+∂Nxy∂x−I0∂2u0∂t2+I1∂3w0∂x∂t2=0 | (16) |
∂Ny∂y+∂Nxy∂y−I0∂2v0∂t2+I1∂3w0∂y∂t2=0 | (17) |
∂2Mx∂x2+∂2My∂y2−2∂2Mxy∂x∂y−Nt(∂2w0∂x2+∂2w1∂x2)=I0∂2w0∂t2+I1(∂3u0∂x∂t2+∂3v0∂x∂t2)−I2(∂4w0∂x2∂t2+∂4w0∂y2∂t2) | (18) |
对于正交各向异性正规对称正交铺设的复合材料板,其刚度矩阵满足[16]
D16=D26=0,A16=A26=0,Bij=0 | (19) |
根据 Kirchhoff 薄板理论及经典弹性理论,薄板中面在变形过程中没有伸长变形,将板的本构关系代入式(16)~式(18)中,略去含
D11∂4w0∂x4+(2D12+4D66)∂4w0∂x2∂y2+D22∂4w0∂y4+Nt(∂2w0∂x2+∂2w1∂x2)−I0∂2w0∂t2+I2(∂4w0∂x2∂t2+∂4w0∂y2∂t2)=0 | (20) |
设三边简支和应力波传播到Lcr(
w0(x,y,t)=Rij[siniπxLcr+ii+1sin(i+1)πxLcr]sinjπyLbsin(ωt+φ) | (21) |
由屈曲模态确定的缺陷分布形式是板结构最有可能发生的屈曲形式,能够很好地确定结构的缺陷敏感性[18]。对于复合材料板在制造过程中出现的初始几何缺陷,引入屈曲模态的
w1(x,y)=εRij[siniπxLcr+ii+1sin(i+1)πxLcr]sinjπyLb | (22) |
式中:i、j为屈曲模态阶数,
根据式(21)和式(22),利用棣莫弗公式对控制方程式(20)中的各项进行求导并化简,得到
D11(πLcr)4[i4+(i+1)4]sin(ωt+φ)+(2D12+4D66)(πLcr)2[i2+(i+1)2](jπLb)2sin(ωt+φ)+2D22(jπLb)4sin(ωt+φ)−Nt(πLcr)2[i2+(i+1)2][sin(ωt+φ)+ε]+2I0ω2sin(ωt+φ)+I2ω2{(πLcr)2[i2+(i+1)2]+2(jπLb)2}sin(ωt+φ)=0 | (23) |
根据式(23)可以得到N的表达式
N={D11(πLcr)4[i4+(i+1)4]+(2D12+4D66)(πLcr)2[i2+(i+1)2](jπLb)2+2D22(jπLb)4+2I0ω2+I2ω2{(πLcr)2[i2+(i+1)2]+2(jπLb)2}}sin(ωt+φ)(πLcr)2[i2+(i+1)2][sin(ωt+φ)+ε] | (24) |
板发生屈曲时,临界条件为
Ncr=D11(πLcr)2[i4+(i+1)4]sinφ[i2+(i+1)2](sinφ+ε)+(2D12+4D66)(jπLb)2sinφ(sinφ+ε)+2D22(j2πL2b)2L2crsinφ[i2+(i+1)2](sinφ+ε) | (25) |
利用MATLAB数值分析应力波未反射时初始几何缺陷、初相位、铺层角度、屈曲模态阶数、铺设厚度以及铺层层数对复合材料板振动屈曲临界荷载的影响,使用的材料参数见表1[21]。
图3显示了复合材料板的模态阶数i为1、2、3,j为1时的屈曲模态。当板的x方向模态增大时,x方向的屈曲模态第一峰值增大且波数增加,而y方向的屈曲模态呈正对称分布。模态阶数的增加使板振动的屈曲模态变得更复杂。
设置7组算例,分别以初始几何缺陷、初相位、铺层角度、x和y两个方向屈曲模态阶数、铺层层数及铺设厚度为变量进行算例分析,研究以上因素对板振动屈曲临界荷载的影响,算例参数见表2。
Group | Initial defect coefficient | Order of mode | Laying angle/(°) | Initial phase | Number of layers laid | Thickness of the plate/m | |
x direction | y direction | ||||||
A | Variable | i = 1 | j = 1 | [0, 0, 0, 0, 0] | π/2 | 5 | 0.01 |
B | 0.1 | Variable | j = 1 | [0, 0, 0, 0, 0] | π/2 | 5 | 0.01 |
C | 0.1 | i =1 | Variable | [0, 0, 0, 0, 0] | π/2 | 5 | 0.01 |
D | 0.1 | i =1 | j = 1 | Variable | π/2 | 5 | 0.01 |
E | 0.1 | i =1 | j = 1 | [0, 0, 0, 0, 0] | Variable | 5 | 0.01 |
F | 0.1 | i =1 | j = 1 | [0, 0, 0, 0, 0] | π/2 | Variable | 0.01 |
G | 0.1 | i =1 | j = 1 | [0, 0, 0, 0, 0] | π/2 | 5 | Variable |
将A组数据代入式(25)中,可以得到不同初始缺陷系数对复合材料板振动屈曲的影响,如图4所示。由Lcr-Ncr曲线可知:在应力波传播过程中,Ncr呈指数下降,分为敏感区和非敏感区。应力波在Lcr < 0.4 m区域传播时,Ncr的变化较陡峭,该区域为敏感区;应力波在Lcr > 0.4 m区域传播时,Ncr的变化趋于平缓,该区域为非敏感区,因此敏感分界点为0.4 m。当选取的初始缺陷系数增大时,临界荷载Ncr也随之增大。在敏感区,初始缺陷系数对临界荷载Ncr的影响较大,且随应力波传播呈减小趋势。此外,初始缺陷系数对非敏感区的影响较小。图4表明,初始几何缺陷系数越大,板越容易发生屈曲。
将B组数据代入式(25)中,可以得到不同x方向模态阶数对复合材料板振动屈曲的影响,如图5所示。由Lcr-Ncr曲线可知:当选取的x方向模态阶数增大时,临界荷载Ncr随之明显增大。在敏感区,x方向模态阶数对临界荷载Ncr的影响很大,并随应力波的传播不断减小,到达非敏感区之后影响较小并趋于稳定。图5表明,x方向模态阶数的增加会显著增大板的屈曲临界荷载。
将C组数据代入式(25)中,可以得到不同的y方向模态阶数对复合材料板振动屈曲的影响,如图6所示。由Lcr-Ncr曲线可知:当选取的y方向模态阶数增大时,临界荷载Ncr也随之变大。在应力波传播过程中,在敏感区y方向模态阶数对临界载荷基本没有影响,而在非敏感区有极小的影响。图6表明,y方向模态阶数的变化对板屈曲临界荷载基本没有影响。
将D组数据代入式(25)中,得到不同铺层角度对复合材料板振动屈曲的影响,如图7所示。由Lcr-Ncr曲线可知:在敏感区,不同的铺设角度对临界荷载Ncr的影响较大,且随应力波的传播不断减小,到达非敏感区后趋于稳定。图7表明,铺层角度小的单层板的层数越多,板的临界屈曲荷载越大,说明复合材料板的铺设角度直接影响板的屈曲临界荷载。
将E组数据代入式(25)中,得到不同初相位对复合材料板振动屈曲的影响,如图8所示。由Lcr-Ncr曲线可知:振型函数的初相位越大,对应的临界荷载越大。在敏感区,振型函数的初相位对临界荷载Ncr的影响较小,且随应力波的传播不断减小;到达非敏感区之后,影响趋于平缓。图8表明,振型函数的初相位越大,板的屈曲临界荷载越大。
将F组数据代入式(25)中,得到不同铺层层数对复合材料板振动屈曲的影响,如图9所示。由Lcr-Ncr曲线可知:当按照不同层数铺设时,敏感区的临界荷载Ncr的变化较大,且随应力波的传播不断减小;到达非敏感区之后变化较小并趋于平缓。图9表明,对于厚度固定且对称铺设的板,当铺设层数达到7时,其屈曲荷载随层数增加趋于稳定。
将G组数据代入式(25)中,得到不同铺设厚度对复合材料板振动屈曲的影响,如图10所示。由Lcr-Ncr曲线可知:板的铺设厚度越大,对应的临界荷载越大。在敏感区,不同的板厚对临界荷载Ncr的影响很大,且随应力波的传播不断减小,到达非敏感区后趋于稳定。图10表明,复合材料板的铺设厚度将直接决定板的屈曲临界荷载。
基于Kirchhoff薄板理论和Hamilton变分原理,建立了具有初始几何缺陷的四边简支复合材料板的振动控制方程。采用伽辽金法,选取符合边界条件的振型函数求解控制方程,得到屈曲临界载荷表达式。数值计算结果表明:应力波在未发生反射前的传播过程中,复合材料板的振动屈曲临界载荷随着临界长度的增大、铺设厚度的减小、初始几何缺陷系数的增大、振型函数初相位的减小而减小;复合材料板的各层铺层角度与荷载作用方向的夹角越小,屈曲临界载荷越大,当对称铺设层数达7层时,临界荷载趋于稳定。研究结果可为工程中复合材料板的结构设计与应用提供一定的参考。
[1] |
李红, 闫维嘉, 张禹, 等. 先进航空材料焊接过程热裂纹研究进展 [J]. 材料工程, 2022, 50(2): 50–61. doi: 10.11868/j.issn.1001-4381.2021.000676
LI H, YAN W J, ZHANG Y, et al. Research progress of hot crack in fusion welding of advanced aeronautical materials [J]. Journal of Materials Engineering, 2022, 50(2): 50–61. doi: 10.11868/j.issn.1001-4381.2021.000676
|
[2] |
ESPINOSA H D, RIM J E, BARTHELAT F, et al. Merger of structure and material in nacre and bone-Perspectives on de novo biomimetic materials [J]. Progress in Materials Science, 2009, 54(8): 1059–1100. doi: 10.1016/j.pmatsci.2009.05.001
|
[3] |
MIZUTA N, MATSUURA K, KIRIHARA S, et al. Titanium aluminide coating on titanium surface using three-dimensional microwelder [J]. Materials Science and Engineering: A, 2008, 492(1/2): 199–204. doi: 10.1016/j.msea.2008.03.028
|
[4] |
NOFAR M, HOSSEINI H R M, KOLAGAR-DAROONKOLAIE N. Fabrication of high wear resistant Al/Al3Ti metal matrix composite by in situ hot press method [J]. Materials & Design, 2009, 30(2): 280–286. doi: 10.1016/j.matdes.2008.04.071
|
[5] |
HARACH D J, VECCHIO K S. Microstructure evolution in metal-intermetallic laminate (MIL) composites synthesized by reactive foil sintering in air [J]. Metallurgical and Materials Transactions A, 2001, 32(6): 1493–1505. doi: 10.1007/s11661-001-0237-0
|
[6] |
PRICE R D, JIANG F C, KULIN R M, et al. Effects of ductile phase volume fraction on the mechanical properties of Ti-Al3Ti metal-intermetallic laminate (MIL) composites [J]. Materials Science and Engineering: A, 2011, 528(7/8): 3134–3146.
|
[7] |
LI T Z, JIANG F C, OLEVSKY E A, et al. Damage evolution in Ti6Al4V-Al3Ti metal-intermetallic laminate composites [J]. Materials Science and Engineering: A, 2007, 443(1/2): 1–15.
|
[8] |
ADHARAPURAPU R R, VECCHIO K S, JIANG F C, et al. Effects of ductile laminate thickness, volume fraction, and orientation on fatigue-crack propagation in Ti-Al3Ti metal-intermetallic laminate composites [J]. Metallurgical and Materials Transactions A, 2005, 36(6): 1595–1608. doi: 10.1007/s11661-005-0251-8
|
[9] |
张海广, 王瑜, 安连浩, 等. 冲击载荷下分支交错层状仿生复合材料动态断裂行为的实验研究和数值模拟 [J]. 高压物理学报, 2022, 36(1): 014101. doi: 10.11858/gywlxb.20210776
ZHANG H G, WANG Y, AN L H, et al. Experimental study and numerical simulation of dynamic fracture behavior of branch staggered laminated biomimetic composites under impact loading [J]. Chinese Journal of High Pressure Physics, 2022, 36(1): 014101. doi: 10.11858/gywlxb.20210776
|
[10] |
武晓东, 张海广, 王瑜, 等. 冲击载荷下仿贝壳珍珠层Voronoi结构的动态力学响应 [J]. 高压物理学报, 2020, 34(6): 064201. doi: 10.11858/gywlxb.20200559
WU X D, ZHANG H G, WANG Y, et al. Dynamic responses of nare-like Voronoi structure under impact loading [J]. Chinese Journal of High Pressure Physics, 2020, 34(6): 064201. doi: 10.11858/gywlxb.20200559
|
[11] |
WANG R Z, SUO Z, EVANS A G, et al. Deformation mechanisms in nacre [J]. Journal of Materials Research, 2001, 16(9): 2485–2493. doi: 10.1557/JMR.2001.0340
|
[12] |
YOURDKHANI M, PASINI D, BARTHELAT F. Multiscale mechanics and optimization of gastropod shells [J]. Journal of Bionic Engineering, 2011, 8(4): 357–368. doi: 10.1016/S1672-6529(11)60041-3
|
[13] |
ZAREMBA C M, BELCHER A M, FRITZ M, et al. Critical transitions in the biofabrication of abalone shells and flat pearls [J]. Chemistry of Materials, 1996, 8(3): 679–690. doi: 10.1021/cm9503285
|
[14] |
SARIKAYA M, GUNNISON K E, YASREBI M, et al. Mechanical property-microstructural relationships in abalone shell [J]. MRS Online Proceedings Library, 1989, 174(1): 109–116. doi: 10.1557/PROC-174-109
|
[15] |
JACKSON A P, VINCENT J F V, TURNER R M. The mechanical design of nacre [J]. Proceedings of the Royal Society B: Biological Sciences, 1988, 234(1277): 415–440. doi: 10.1098/rspb.1988.0056
|
[16] |
NUKALA P K V V, SIMUNOVIC S. A continuous damage random thresholds model for simulating the fracture behavior of nacre [J]. Biomaterials, 2005, 26(30): 6087–6098. doi: 10.1016/j.biomaterials.2005.03.013
|
[17] |
SHIN Y A, YIN S, LI X Y, et al. Nanotwin-governed toughening mechanism in hierarchically structured biological materials [J]. Nature Communications, 2016, 7(1): 10772. doi: 10.1038/ncomms10772
|
[18] |
GHAZLAN A, NGO T, TAN P, et al. Inspiration from nature’s body armours: a review of biological and bioinspired composites [J]. Composites Part B: Engineering, 2021, 205: 108513. doi: 10.1016/j.compositesb.2020.108513
|
[19] |
LI T Z, GRIGNON F, BENSON D J, et al. Modeling the elastic properties and damage evolution in Ti-Al3Ti metal-intermetallic laminate (MIL) composites [J]. Materials Science and Engineering: A, 2004, 374(1/2): 10–26. doi: 10.1016/j.msea.2003.09.074
|
[20] |
CAO Y, ZHU S F, GUO C H, et al. Numerical investigation of the ballistic performance of metal-intermetallic laminate composites [J]. Applied Composite Materials, 2015, 22(4): 437–456. doi: 10.1007/s10443-014-9416-1
|
[1] | HOU Zekai, LUAN Yunbo, LEI Keming, DONG Qianxi, NIU Tuyao, LI Yongcun. Hybrid Bionic Design Based on the Internal Structures of Nacre and Strombus Gigas Shell[J]. Chinese Journal of High Pressure Physics, 2024, 38(5): 054204. doi: 10.11858/gywlxb.20240724 |
[2] | HUANG Hongdong, YU Xiaohui, HE Duanwei. High Pressure Sintering of B-C-N-Ti Quaternary Superhard Composites[J]. Chinese Journal of High Pressure Physics, 2024, 38(4): 041102. doi: 10.11858/gywlxb.20240769 |
[3] | LI Teng, ZHANG Chenfan, DENG Qingtian, LI Xinbo, WEN Jinpeng. Optimized Design and Energy Absorption of TPU Material Based on Hierarchical Structure[J]. Chinese Journal of High Pressure Physics, 2022, 36(6): 064104. doi: 10.11858/gywlxb.20220542 |
[4] | WANG Haoyang, WEI Ying, QIAO Li. Simulation of Dynamic Crack Propagation in Superconducting Nb3Sn at Extreme Low Temperature[J]. Chinese Journal of High Pressure Physics, 2022, 36(3): 034201. doi: 10.11858/gywlxb.20210884 |
[5] | LI Mengshen, GE Tao, CHENG Wenpin, PEI Guoqing. A Comparative Study on the Finite Element Models for Projectiles Perforation into Reinforced Concrete Slabs[J]. Chinese Journal of High Pressure Physics, 2022, 36(2): 024206. doi: 10.11858/gywlxb.20210816 |
[6] | MENG Xiangsheng, WU Xiaodong, ZHANG Haiguang. Numerical Simulation on Interlaminar Fracture Toughness of 3D Printed Mortar Laminated Composites[J]. Chinese Journal of High Pressure Physics, 2020, 34(4): 044206. doi: 10.11858/gywlxb.20190827 |
[7] | WU Xiaodong, ZHANG Haiguang, WANG Yu, MENG Xiangsheng. Dynamic Responses of Nare-Like Voronoi Structure under Impact Loading[J]. Chinese Journal of High Pressure Physics, 2020, 34(6): 064201. doi: 10.11858/gywlxb.20200559 |
[8] | ZHENG Songlin. Advances in the Study of Dynamic Response of Crystalline Materials by Crystal Plasticity Finite Element Modeling[J]. Chinese Journal of High Pressure Physics, 2019, 33(3): 030108. doi: 10.11858/gywlxb.20190725 |
[9] | PENG Kefeng, PAN Hao, ZHAO Kai, ZHENG Zhijun, YU Jilin. Dynamic Compaction Behaviors of Copper Powders Using Multi-Particle Finite Element Method[J]. Chinese Journal of High Pressure Physics, 2019, 33(4): 044102. doi: 10.11858/gywlxb.20180665 |
[10] | ZHANG Pin-Liang, GONG Zi-Zheng, CAO Yan, LI Yu, XU Kun-Bo, MU Yong-Qiang, WU Qiang, SONG Guang-Ming. Macro- and Micro-Damage Behaviors of 47Zr45Ti5Al3V Alloy in Hypervelocity Impact[J]. Chinese Journal of High Pressure Physics, 2017, 31(3): 231-238. doi: 10.11858/gywlxb.2017.03.004 |
[11] | BAI Jing-Song, YU Ji-Dong, DAI Cheng-Da, WANG Yu, LIU Kun, LUO Guo-Qiang, SHEN Qiang, TAN Hua, WU Qiang, ZHANG Lian-Meng. Calculation Analysis of the Impact Melting and Resolidification Process for the Bismuth Using the Ti-Cu-W Pillow Flyer[J]. Chinese Journal of High Pressure Physics, 2015, 29(1): 15-22. doi: 10.11858/gywlxb.2015.01.003 |
[12] | LIU Jun, WANG Yan-Jin, FENG Qi-Jing, ZHANG Shu-Dao. Simulation Study on Effects of Material Properties of the Bulging Flyers Driven by Colliding Detonation Waves[J]. Chinese Journal of High Pressure Physics, 2015, 29(1): 63-68. doi: 10.11858/gywlxb.2015.01.011 |
[13] | ZHANG Wen-Kai, PENG Fang, GUO Zhen-Tang, GUAN Jun-Wei, LI Rong-Qi. 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 |
[14] | AN Pei, KOU Zi-Li, LI Zi-Yang, QIN Jia-Qian. Study on Thermal Stability of Ti2AlN-Al at High Pressure[J]. Chinese Journal of High Pressure Physics, 2010, 24(1): 71-75 . doi: 10.11858/gywlxb.2010.01.013 |
[15] | LI Qiang, YANG Chao, ZHAN Zai-Ji. Damage and Cracking of Zr41Ti14Cu12.5Ni10Be22.5 Bulk Amorphous Alloy under Planar Shock Compression[J]. Chinese Journal of High Pressure Physics, 2007, 21(1): 108-112 . doi: 10.11858/gywlxb.2007.01.018 |
[16] | WANG Wu-Xiao, JIE Wan-Qi, YUAN Sen. Chemical Reaction of CuO-Al System in Preparation of Al2O3/Cu Composite[J]. Chinese Journal of High Pressure Physics, 2004, 18(4): 353-358 . doi: 10.11858/gywlxb.2004.04.011 |
[17] | TIAN En-Ke, LONG Qi-Wei. Investigations of Complex Fractal Structures with Self-Similarity and Self-Affinity on Fractured Surfaces[J]. Chinese Journal of High Pressure Physics, 2001, 15(3): 178-185 . doi: 10.11858/gywlxb.2001.03.003 |
[18] | ZHAO Fang-Fang, LUO Jing-Run, TIAN Chang-Jin, JIN Zhou-Geng, HE Ying-Bo. The Crack Growth Process of Particulate Filled Polymer Monitored by Acoustic Emission[J]. Chinese Journal of High Pressure Physics, 2000, 14(3): 235-240 . doi: 10.11858/gywlxb.2000.03.014 |
[19] | HAN Jie, ZOU Guang-Tian, CUI Tian. Pressure-Induced Structural Phase Transition of Ba3Bi3Ti4NbO18[J]. Chinese Journal of High Pressure Physics, 1999, 13(1): 34-36 . doi: 10.11858/gywlxb.1999.01.006 |
[20] | DONG Yu-Bin, ZAHNG Wan-Jia, JING Fu-Qian, HAN Jun-Wan, CHEN Da-Nian, SU Lin-Xiang, FENG Jia-Bo. Numerical Analysis for Dynamic Damage Processes and LY-12 Aluminum Spallations[J]. Chinese Journal of High Pressure Physics, 1988, 2(4): 305-312 . doi: 10.11858/gywlxb.1988.04.003 |
Group | Initial defect coefficient | Order of mode | Laying angle/(°) | Initial phase | Number of layers laid | Thickness of the plate/m | |
x direction | y direction | ||||||
A | Variable | i = 1 | j = 1 | [0, 0, 0, 0, 0] | π/2 | 5 | 0.01 |
B | 0.1 | Variable | j = 1 | [0, 0, 0, 0, 0] | π/2 | 5 | 0.01 |
C | 0.1 | i =1 | Variable | [0, 0, 0, 0, 0] | π/2 | 5 | 0.01 |
D | 0.1 | i =1 | j = 1 | Variable | π/2 | 5 | 0.01 |
E | 0.1 | i =1 | j = 1 | [0, 0, 0, 0, 0] | Variable | 5 | 0.01 |
F | 0.1 | i =1 | j = 1 | [0, 0, 0, 0, 0] | π/2 | Variable | 0.01 |
G | 0.1 | i =1 | j = 1 | [0, 0, 0, 0, 0] | π/2 | 5 | Variable |
E1/GPa | E2/GPa | G12/GPa | μ12 | La/m | Lb/m |
140.0 | 8.6 | 5.0 | 0.35 | 0.60 | 0.50 |
Group | Initial defect coefficient | Order of mode | Laying angle/(°) | Initial phase | Number of layers laid | Thickness of the plate/m | |
x direction | y direction | ||||||
A | Variable | i = 1 | j = 1 | [0, 0, 0, 0, 0] | π/2 | 5 | 0.01 |
B | 0.1 | Variable | j = 1 | [0, 0, 0, 0, 0] | π/2 | 5 | 0.01 |
C | 0.1 | i =1 | Variable | [0, 0, 0, 0, 0] | π/2 | 5 | 0.01 |
D | 0.1 | i =1 | j = 1 | Variable | π/2 | 5 | 0.01 |
E | 0.1 | i =1 | j = 1 | [0, 0, 0, 0, 0] | Variable | 5 | 0.01 |
F | 0.1 | i =1 | j = 1 | [0, 0, 0, 0, 0] | π/2 | Variable | 0.01 |
G | 0.1 | i =1 | j = 1 | [0, 0, 0, 0, 0] | π/2 | 5 | Variable |