压力诱导下CsGeBr<sub>3</sub>的结构相变

曲佳; 王弈铭; 王欣; 杨文革

doi:10.11858/gywlxb.20230769

压力诱导下CsGeBr₃的结构相变

doi: 10.11858/gywlxb.20230769

曲佳^{1, 2,},
王弈铭²,
王欣^1,*, ,,
杨文革^2,*, ,

1.
吉林大学物理学院，超硬材料国家重点实验室, 吉林长春　130012
2.
北京高压科学研究中心, 上海　201203

详细信息

作者简介:
曲　佳（1995－），女，博士研究生，主要从事卤化物钙钛矿的高压物性研究.E-mail：qujia19@mails.jlu.edu.cn

通讯作者:
王　欣（1973－），男，博士，教授，主要从事极端条件下材料的物性研究. E-mail：xin_wang@jlu.edu.cn

杨文革（1968－），男，博士，教授，主要从事高压材料的结构和物性研究. E-mail：yangwg@hpstar.ac.cn

中图分类号: O521.2
计量
- 文章访问数: 409
- HTML全文浏览量: 93
- PDF下载量: 47
出版历程
- 收稿日期: 2023-10-25
- 修回日期: 2023-12-15
- 录用日期: 2023-12-25
- 网络出版日期: 2024-07-19
- 刊出日期: 2024-09-29

Pressure-Induced Structural Phase Transition in Halide Perovskite CsGeBr₃

QU Jia^{1, 2
,},
WANG Yiming²,
WANG Xin^{1,*
, ,},
YANG Wenge^{2,*
, ,}

1.
State Key Laboratory of Superhard Materials, College of Physics, Jilin University, Changchun 130012, Jilin, China
2.
Center for High Pressure Science and Technology Advanced Research, Shanghai 201203, China

摘要

摘要: 近年来，压力下卤化物钙钛矿成为新的研究热点，呈现出许多优异的电学和光学等特性。高压下钙钛矿结构演变研究是所有物性研究的基石和重点。利用金刚石对顶砧压机，结合原位高压同步辐射X射线衍射、原位高压拉曼光谱、紫外-可见-近红外分光光度计测量技术和第一性原理计算，对全无机卤化物钙钛矿CsGeBr₃在高压下的结构演变进行了系统研究。结果表明：CsGeBr₃在常压下是菱方 $R3m$ 结构；在1 GPa时，CsGeBr₃发生菱方 $R3m$ 到立方 $Pm\overline{3}m$ 的结构相变；在更高的压力下保持立方结构；菱方 $R3m$ 到立方 $Pm\overline{3}m$ 的相变是可逆的。研究结果为进一步探索卤化物钙钛矿在压力下的性质、拓展其应用前景提供了重要的科学依据。
- 高压 /
- CsGeBr₃ /
- 结构相变 /
- 卤化物钙钛矿
Abstract: In recent years, pressure-induced physical properties of halide perovskites have attracted significant research interests due to their excellent optical and electronic properties. The study of the structural evolution of perovskite under compression is the foundation and key point of all physical property researches. In this paper, we systematically investigated the structural evolution of the all-inorganic halide perovskite CsGeBr₃ under compression using in situ high-pressure synchrotron X-ray diffraction, in situ high-pressure Raman spectroscopy, ultraviolet/visible/near-infrared spectrophotometry, and first-principles calculations. Our results show that CsGeBr₃ undergoes a reversible rhombohedral $R3m$ to cubic $Pm\overline{3}m$ structural phase transition at 1 GPa, and the cubic $Pm\overline{3}m$ phase maintains at higher pressures. This study provides important scientific basis for further exploration of the properties and applications of halide perovskites under compression.
- high pressure /
- CsGeBr₃ /
- structural phase transition /
- halide perovskite

HTML全文

高桩码头广泛应用于淤泥质海岸及河口地区，是我国码头的重要结构形式。水下爆炸是高桩码头战时遭受的主要威胁。开展水下爆炸作用下高桩码头毁伤效应研究，对于提高高桩码头的战时生存能力和保障能力、提高军用高桩码头设计水平，具有十分重要的意义。

国内混凝土结构水下抗爆研究的对象主要集中在大坝、桥梁、码头等重要的水上或者海上大型项目、平台和军事战略工程。基于三峡三期RCC围堰拆除工程，刘美山等^[1]、李裕春等^[2]通过试验研究了流场环境对混凝土屈服强度、结构破坏过程和毁坏程度的影响，为混凝土水下爆炸研究提供依据。赵根等^[3]通过混凝土墩浅水爆炸试验，分别研究了单、双装药对称及不对称设置同时起爆条件下混凝土墩的破坏情况，研究发现混凝土墩处于多轴应力状态。相较于实地试验，数值仿真方法效率更高。Georgin等^[4]研究了有限元混凝土模型在爆炸冲击下的应变率效应。李建阳等^[5]通过数值软件研究了混凝土立方体在不同工况水下爆炸荷载下的破坏过程，结果表明：爆炸近场区水深对混凝土损伤影响不大，爆炸远场区增加水深能削弱混凝土破坏。对于大型水工建筑，张社荣等^[6]、王高辉等^[7]利用数值模拟软件，考虑水深等因素，研究了混凝土坝的毁伤机理和破坏形态。闫秋实等^[8]采用数值仿真软件分析炸深和爆距对钢筋混凝土单桩的毁伤效应，得到了特定深度下单桩的抗爆区域。流场边界对水下爆炸冲击波和气泡脉动影响很大^[9]，考虑到高桩码头桩基处于复杂的流场边界条件，有必要分别对冲击波阶段、气泡脉冲阶段高桩码头桩基动态响应和毁伤效应进行研究。

考虑港池环境影响，采用LS-DYNA有限元软件开展水下爆炸下高桩码头毁伤效应研究，从冲击波传播和气泡脉动两个阶段分析高桩码头桩基受到的荷载特性及其动态响应，研究高桩码头的毁伤机理和破坏模式，为进一步开展高桩码头抗爆研究奠定基础。

1. 数值模型

1.1 高桩码头及港池环境

以某高桩梁板式码头为例，港池水深17 m，水底为饱和黏土。整个码头由标准结构段连接。桩、梁、板均采用钢筋混凝土结构，上部面板宽26.5 m，包含7个桩基，其中3个直桩和2对斜叉桩，直桩与斜桩间距4.5 m，斜桩间距1.7 m，直桩和斜桩均为截面600 mm × 600 mm的方形桩，斜桩与横梁截面的夹角为8°，与纵梁截面的夹角为10°。单个排架跨度为7.5 m，一般12个排架之间设置变形缝，忽略码头系泊、防护和前沿设备影响，沿海岸线方向截取3个排架，如图1所示，开展水下爆炸作用下高桩码头毁伤效应研究。

图 1 码头模型三视图（单位：cm）

Figure 1. Three views of wharf model （Unit：cm）

下载: 全尺寸图片幻灯片

1.2 有限元模型

高桩码头计算区域如图2所示。炸药采用Mark80系列中使用最广泛的Mark82通用炸弹，内部装有192磅（87 kg）Tritonal高爆炸药，等效TNT当量117.45 kg，炸药位于两个直桩之间。考虑到有限元模型的对称性和计算效率，选取炸药附近桩基和上部结构对称建立模型，横向选取3跨，纵向选取2跨。有限元模型如图3所示。该模型包括混凝土、钢筋、空气、水、土、炸药6种材料，其中：混凝土和钢筋采用拉格朗日单元，空气、水、炸药、黏土采用欧拉单元。通过ALE_COUPLING_NODAL_CONSTRAINT定义混凝土与钢筋之间的黏结，采用ALE算法实现结构大变形和流固耦合作用，通过CONSTRAINED_LAGRANGE_IN_SOLID定义结构与流场之间的耦合作用，流域四周边界采用环境单元（Ambient），通过LOAD_BODY_Z和INITIAL_STRESS_DEPTH初始化重力场和静水压力。通过关键字INITIAL_VOLUME_FRACTION_GEOMETRY和INITIAL_DETONATION实现球形装药和中心起爆。以炸药中垂面为对称面，建立1/2有限元模型，对称面上施加单向约束，码头上部结构边界处施加单向约束。近场冲击波属于高频荷载，为防止冲击波衰减过快和耦合泄露，在流固耦合区域采用加密网格，流场网格尺寸为9 cm，桩基网格尺寸为10 cm，确保准确捕捉冲击波与桩基的耦合作用；同时为了提高计算效率，在远离炸药和桩基区域采用渐变过渡网格^[10]。

图 2 高桩码头计算区域

Figure 2. Calculation region of high-piled wharf

下载: 全尺寸图片幻灯片

图 3 有限元模型

Figure 3. Finite element model

下载: 全尺寸图片幻灯片

采用正确的材料参数对数值模拟结果十分重要。混凝土采用Concrete_Damage_Rel3模型模拟，此时只需要提供混凝土的无侧限抗压强度、泊松比以及密度即可自动生成混凝土参数^[11]。该模型引入初始屈服面、极限强度面和残余强度面3种强度面，并考虑偏应力不变量对强度破坏面的影响。在Mat_Concrete_Damage_Rel3混凝土材料模型中：当应力达到初始屈服面但未达到极限强度面时，通过初始屈服面和极限强度面的线性插值表示；当应力达到极限强度面但未达到残余强度面时，通过极限强度面和残余强度面的线性插值表示

$\Delta {\sigma _{\rm{f}}} = \eta \Delta {\sigma _{\rm{m}}}{\rm{ + (1 - }}\eta {\rm{)}}\Delta {\sigma _{\rm{y}}} \;\;\;\;\; \lambda \leqslant {\lambda _{\rm{m}}}({\text{后继屈服面}})$

(1)

$\Delta {\sigma _{\rm{h}}} = \eta \Delta {\sigma _{\rm{m}}}{\rm{ + (1 - }}\eta {\rm{)}}\Delta {\sigma _{\rm{r}}}\;\;\;\;\;\lambda > {\lambda _{\rm{m}}}\;({\text{后继软化面}})$

(2)

式中：λ为损伤变量，是等效塑性应变的函数；λ_m为损伤转折点，是强化段和软化段的边界；η为损伤变量λ的函数，λ < λ_m为强化段，由零增至1，λ > λ_m为软化段，由1减至零。Mat_Concrete_Damage_Rel3混凝土模型能够输出混凝土的比例损伤变量 ${\delta}$ （Scaled damage factor）

$\delta {\rm{ = }}2\lambda /(\lambda + {\lambda \rm_m})$

(3)

当δ = 0时，混凝土处于弹性阶段；当0 < δ < 1时，混凝土屈服，进入塑性阶段；当1 < δ < 2时，混凝土进入软化阶段。

码头桩基混凝土强度等级为C50，上部梁板混凝土强度等级为C40，密度为2 550 kg/m³，泊松比为0.2。通过动载增大系数DIF曲线设定混凝土应变率效应，混凝土抗压动载增大系数CDIF和抗拉动载增大系数TDIF分别为^[12-13]

${\rm CDIF} = \frac{{{f\rm_{dc}}}}{{{f\rm_c}}} = \left\{ {\begin{array}{*{20}{l}} 1&{\dot \varepsilon < {{\dot \varepsilon }\rm_{stat}}}\\ {{{\left( {\dfrac{{\dot \varepsilon }}{{{{\dot \varepsilon }\rm_{stat}}}}} \right)}^{1.026\alpha }}}&{{{\dot \varepsilon }\rm_{stat}} \leqslant \dot \varepsilon \leqslant 30 \; {\rm s^{ - 1}}}\\ {\gamma {{\left( {\dfrac{{\dot \varepsilon }}{{{{\dot \varepsilon }_{\rm{stat}}}}}} \right)}^{1/3}}}&{30\;{\rm s^{ - 1}} < \dot \varepsilon < 300\;{\rm s^{ - 1}}} \end{array}} \right.$

(4)

${\rm{TDIF}} = \frac{{{f\rm_{dt}}}}{{{f\rm_t}}} = \left\{ {\begin{array}{*{20}{l}} {{{\left( {\dfrac{{\dot \varepsilon }}{{{{\dot \varepsilon }\rm_{stat}}}}} \right)}^\beta }} & {{{\dot \varepsilon }\rm_{stat}} < \dot \varepsilon \leqslant 1\;{\rm s^{ - 1}}}\\ {\xi {{\left( {\dfrac{{\dot \varepsilon }}{\dot \varepsilon}_{\rm{stat}}} \right)}^{1/3}}} & {1\;{\rm s^{ - 1}} < \dot \varepsilon < 160\;{\rm s^{ - 1}}} \end{array}} \right.$

(5)

式中：f_dc、f_dt分别为动态压缩和拉伸强度，f_c、f_t分别为准静态压缩和拉伸强度， ${\dot \varepsilon }$ 为应变率， ${{\dot \varepsilon}_{\rm{stat}}}$ 为准静态应变率，α、β、γ、ξ为参数。

钢筋采用Mat_Plastic_Kinematic弹塑性模型描述，采用HRB335强度，泊松比为0.3，应变率参数C和P分别取40.4和5。假设空气、水、炸药均为均匀连续，空气、水、炸药、黏土分别采用线性多项式状态方程、Grüneisen状态方程、标准JWL状态方程和线弹性模型描述。表1^[14]详细列出了材料参数，其中：C₀～C₆为线性多项式状态方程参数，E为初始单位质量内能，c、S₁～S₃为Grüneisen状态方程参数，Γ为Grüneisen常数，A、B、ω、R₁、R₂为标准JWL状态方程参数，E为弹性模量，G为剪切模量。

表 1 材料参数

Table 1. Material parameters

Material	ρ/(kg·m⁻³)	C₀, C₁	C₂, C₃	C₄	C₅	C₆	E/(J·kg⁻¹)
Air	1.29	0	0	0.4	0.4	0	2.5×10⁵

Material	ρ/(kg·m⁻³)	c	S₁	S₂	S₃	Γ
Water	1 000	1 480	2.56	−1.986	0.226 8	0.5

Material	ρ/(kg·m⁻³)	A	B	ω	R₁	R₂	E₀/(GJ·m⁻³)
Explosive	1 630	3.74 × 10¹¹	7.33 × 10⁹	0.3	4.15	0.95	7

Material	ρ/(kg·m⁻³)	E/MPa	G/MPa
Soil	1 800	16	8

下载: 导出CSV

| 显示表格

2. 计算结果及分析

2.1 有效性验证

Cole^[15]根据大量试验数据得到自由场水下爆炸的冲击波峰值、气泡最大半径和周期的经验公式

${p_{\max }} = 52.16{\left(\frac{{{W^{1/3}}}}{S}\right)^{1.13}}$

(6)

${R\rm_m} = 3.30{\left(\frac{W}{{H + 10.33}}\right)^{1/3}}$

(7)

$T = 2.08\frac{{{W^{1/3}}}}{{{{(H + 10.33)}^{5/6}}}}$

(8)

式中：W为炸药装量，kg；S为测点距离，m；p_max为冲击波峰压，MPa；H为炸药水深，m；T为第1次气泡脉动周期，s；R_m为气泡最大半径，m。复杂边界条件下，气泡很难保持球形运动，为衡量气泡体积变化，定义气泡最大等效半径： ${R\rm_m} = {\left(\dfrac{{3{V\rm_m}}}{{4\text{π} }}\right)^{1/3}}$ ，其中V_m为气泡最大体积，m³。

两个直桩之间的流场受冲击波反射、绕射的影响较小，在测距S为 2、4、6、8 m处选取测点，对比自由场冲击波峰值压力的经验公式和数值结果，如表2所示，其中p_max-free为自由场冲击波峰压，p_max-straight为直桩附近冲击波峰压，p_max-oblique为斜桩附近冲击波峰压。对比发现，数值仿真结果与经验公式计算结果符合较好。选取直桩、斜桩附近流场测点，测点水深与炸药深度一致，发现桩基迎爆面由于冲击波反射叠加，峰值压力较自由场略微增强，桩基背面由于冲击波绕射，峰值压力衰减。气泡受桩基阻碍，气泡周期和气泡最大半径均较自由场经验结果减小。

表 2 数值模拟与理论结果比较

Table 2. Comparison of the numerical and theoretical results

Method	p_max-free/MPa				p_max-straight/MPa		p_max-oblique/MPa		Bubble
Method	2 m	4 m	6 m	8 m	Blast face	Back blast face	Blast face	Back blast face	T/s	R_m/m
Theoretical	144.15	65.86	41.65	30.09	68.09	57.90	33.83	31.11	0.88	6.07
Numerical	153.21	68.21	42.27	27.34	72.35	21.76	33.84	16.06	0.62	5.29
Error/%	6.29	3.57	1.49	9.14	6.25	−62.42	0.03	−48.38	−29.55	−12.85

下载: 导出CSV

| 显示表格

2.2 冲击波阶段

炸药爆炸后首先作用桩基的是水下爆炸冲击波。爆炸冲击波具有峰值大、持续时间短的特点。考虑码头模型的对称性，在炸药附近直桩和斜桩上沿高度方向设置7个典型测点，见图4，直桩和斜桩上测点由水底向上依次为Z_h₁～Z_h₇和X_h₁～X_h₇。

图 4 炸药及测点位置示意图

Figure 4. Schematic of explosive and the measure points

下载: 全尺寸图片幻灯片

图5显示了直桩上测点的速度（v）和加速度（a）响应曲线。可见，直桩中部测点Z_h4的速度、加速度响应最大，并由桩身中部向桩顶、桩底两端逐渐减小。桩身迎爆面经历冲击波叠加、局部空化和空化闭合3个阶段：桩基在冲击作用下迅速加速达到第1个速度峰值；冲击波叠加导致迎爆面出现局部空化，空化的发生截断了压力载荷，桩基失去了加载载荷并开始减速；随着空化闭合，出现第2个速度峰值，然后逐渐减速，因此桩中测点Z_h4的速度、加速度曲线呈现明显的二次加载现象。桩中测点Z_h4的加速度在t = 2.02 ms到达正向峰值，在t = 2.34 ms达到反向峰值，桩基迎爆面空化和冲击波绕射共同导致桩身产生反向加速度，冲击波阶段后桩基受滞后流和反作用力共同影响，桩身加速度呈现周期性振荡。斜桩的速度、加速度响应如图6所示。可见，其响应规律与直桩基本一致，只是当冲击波作用在斜桩斜侧面时，冲击波反射较弱，空化现象不明显。通过对比桩身深水区、浅水区测点的加速度、速度时程曲线可以发现：桩基浅水区测点Z_h5、Z_h6的加速度响应峰值较相同爆炸距离的深水区测点Z_h3、Z_h2分别增加了53.51%和91.25%，因此水下爆炸作用下桩基浅水区的冲击响应较深水区更加剧烈。

图 5 直桩的速度和加速度响应

Figure 5. Velocity and acceleration of straight pile

下载: 全尺寸图片幻灯片

图 6 斜桩的速度和加速度响应

Figure 6. Velocity and acceleration of oblique pile

下载: 全尺寸图片幻灯片

2.3 气泡脉动阶段

图7为Mark82战斗部装药量水下爆炸时，高桩码头的破坏过程以及气泡脉动过程。气泡第1次脉动周期为0.62 s，最大气泡半径为5.29 m，气泡在t = 0.07 s时膨胀至直桩内侧并挤压直桩，在冲击波和滞后流的作用下，直桩和斜桩均产生明显外向变形，桩顶、桩中以及桩基与横梁连接处出现明显损伤。气泡在t = 0.29 s时达到最大体积并绕过直桩，此时4个桩基、桩顶受拉侧混凝土失效，桩中迎爆面和背面均产生较大损伤。气泡在达到最大体积后开始收缩，桩基受气泡Bjerknes力“吸引”向内侧产生变形，桩中和桩顶的毁伤范围和破坏程度进一步增加。气泡在t = 0.62 s时产生垂直水面向上的射流，射流击穿气泡上表面形成环状气泡，桩基附近气泡因受Bjerknes力吸引而收缩缓慢，形成尾状长条气泡。气泡第2次脉动以环状气泡开始膨胀，并在t = 1.00 s时由于气泡多位置不同速度坍塌，导致气泡溃散。研究表明，气泡第1次脉动后，剩余能量的破坏作用只有初始能量的7%左右^[16]，故气泡第2次脉动阶段并未对码头造成明显的毁伤效应，高桩码头的毁伤现象较第1次气泡脉动结束时刻基本维持不变，因此可以判断水下爆炸作用下高桩码头的毁伤现象在气泡第1次脉动结束阶段已经基本完成。

图 7 高桩码头毁伤和气泡脉动过程（Mark82）

Figure 7. Bubble pulse process and the damage process of the wharf (Mark82)

下载: 全尺寸图片幻灯片

图8、图9为水下爆炸作用下直桩和斜桩的速度及位移（x）响应曲线。桩基在气泡脉动过程中随气泡膨胀、收缩发生周期性往复运动，在冲击波初始阶段和气泡射流时刻桩身受到初始冲击波和气泡脉冲载荷，桩基速度出现两次突变，直桩和斜桩均在气泡膨胀末期和气泡第1次脉动结束阶段达到最大位移，分别为0.58 m（−0.30 m）和0.12 m（−0.09 m），直桩距离炸药较近，其速度、加速度响应峰值更大。对比桩基深水区、浅水区测点的速度响应发现：气泡膨胀初期，桩基中部的速度、位移响应最大；在气泡射流时，桩基浅水区测点Z_h5的峰值速度超过了桩基中部测点Z_h4。这是由于气泡在第1次脉动过程上浮了约4 m，此时气泡深度与测点Z_h5基本一致，因此桩基浅水区的速度、位移响应随着气泡不断上浮而逐渐增强，其毁伤效应大于桩基深水区。

图 8 直桩的速度和位移响应

Figure 8. Velocity and displacement of straight pile

下载: 全尺寸图片幻灯片

图 9 斜桩的速度和位移响应

Figure 9. Velocity and displacement of oblique pile

下载: 全尺寸图片幻灯片

2.4 高桩码头破坏过程和毁伤机理

图10显示了高桩码头在冲击波和气泡脉动作用下的毁伤形貌。在冲击波阶段，高桩码头损伤较轻微，桩基的中部、顶部受拉区混凝土进入塑性阶段，码头上部结构仍处于弹性阶段。在气泡脉动阶段，高桩码头桩基的损伤范围和损伤程度明显增大，横、纵梁连接处也出现不同程度的损伤。桩基的顶部和中部毁伤最为严重，桩基顶部混凝土由于固接横梁造成内部应力波反射叠加，连接处混凝土处于高应力状态，达到极限应变而失效；由于桩基在气泡作用下进行往复运动，桩基的迎爆面和背面处于周期性压、拉应力状态，因此近场直桩中部两侧混凝土均在达到极限应变后失效。高桩码头上部结构的损伤程度小，横梁与纵梁连接处出现轻微损伤，面板没有产生损伤。高桩码头的破坏形态在第1次脉动结束时已经基本完成。在冲击波阶段，桩基受到的冲击载荷峰值压力较大，但是由于荷载持续时间短，冲击波绕射导致桩基处于围压状态，因此桩基的毁伤程度较弱；在气泡脉动阶段，由于气泡膨胀挤压和气泡收缩吸引桩基产生往复变形，桩中、桩顶混凝土达到极限应变而失效，钢筋外露。

图 10 冲击波和气泡脉动作用下高桩码头的毁伤

Figure 10. Damage of high-piled wharf subject to shock wave and bubble impulse

下载: 全尺寸图片幻灯片

3. 炸药当量对毁伤效应的影响

3.1 高桩码头毁伤现象

3.1.1 Mark81

Mark81内部装有100磅（45.4 kg）Tritonal高爆炸药，等效61.29 kg TNT。图11为Mark81炸药引爆后高桩码头毁伤以及气泡脉动过程。冲击波阶段，直桩受拉侧混凝土进入塑性阶段，斜桩距离炸药稍远，仍处于弹性阶段。气泡初始呈现球形膨胀，在t = 0.25 s时达到最大气泡半径，此时桩基损伤程度和损伤范围增大，桩顶与横梁连接处产生一定损伤。气泡受桩基Bjerknes效应影响，在桩基附近形成尾状长条气泡，同时桩也受到气泡“吸引”出现向内侧变形，气泡在t = 0.5 s时完成第1次气泡脉动。高桩码头的整体毁伤程度较弱，直桩中部和顶部受拉侧部分混凝土单元失效。

图 11 高桩码头毁伤以及气泡脉动过程（Mark81）

Figure 11. Bubble pulse process and the damage process of the wharf (Mark81)

下载: 全尺寸图片幻灯片

3.1.2 Mark83

Mark83内部装有445磅（202 kg）Tritonal高爆炸药，等效272.70 kg TNT。图12为Mark83炸药引爆后高桩码头毁伤以及气泡脉动过程。冲击波阶段，直桩中部混凝土进入塑性阶段，斜桩受拉侧进入塑性阶段。初始时气泡呈现球形膨胀，在t = 0.30 s时达到最大气泡半径，由于桩身阻碍气泡膨胀，气泡水平方向稍扁并在桩身两侧绕过，由于桩身近水面气穴效应，水面向下凹陷引起气泡上表面凹陷。此时直桩的桩顶与横梁连接处混凝土完全失效，桩中背爆面混凝土失效，炸药两侧直桩已经失去承载力。斜桩的顶部与横梁连接处混凝土损伤严重，码头面板的损伤区域由中心向两侧发展，横梁、纵梁与面板连接处也出现局部损伤。气泡在收缩过程中受到桩基的Bjerknes效应，导致气泡在桩身附近的收缩速度较慢，气泡上下表面的收缩速度快，气泡多部位塌陷。在t = 0.71 s时气泡体积收缩至最小，气泡射流方向紊乱，气泡破裂成多个气泡并重新开始膨胀。高桩码头在气泡收缩的过程中，桩基受气泡“吸引”向内侧发生变形，桩基和码头面板的毁伤范围稍微增大。

图 12 高桩码头毁伤过程以及气泡脉动过程（Mark83）

Figure 12. Bubble pulse process and the damage process of the wharf (Mark83)

下载: 全尺寸图片幻灯片

Mark80系列战斗部装药量水下爆炸时，高桩码头在冲击波阶段的毁伤程度很弱，毁伤作用在第1次气泡脉动结束时基本完成，毁伤区域主要集中在桩基中部、顶部、底部以及横、纵梁连接处，其中桩顶与横梁连接处为最薄弱部位，码头面板在Mark81、Mark82工况下基本没有损伤，在Mark83工况下码头面板和横、纵梁均出现大面积损伤。高桩码头在气泡第1次脉动膨胀阶段的损伤积累较快，在气泡第1次脉动收缩阶段损伤缓慢增加并趋于稳定，在随后的气泡脉动过程中高桩码头的毁伤现象和毁伤范围基本维持不变。

3.2 高桩码头毁伤效应评估

对高桩码头的毁伤效应进行评估。考虑到高桩码头的桩基破坏最严重，桩中最大位移可以在一定程度上反映高桩码头的毁伤效应，然而当桩基呈剪切破坏时，桩中最大位移无法完全反映其承载性能的损伤效应，进而低估码头的破坏程度，因此根据高桩码头的竖向剩余承载力定义高桩码头损伤评估系数D^[16]

$D = 1 - \frac{{{p\rm_r}}}{{{p_0}}}$

(9)

式中：p_r为高桩码头在水下爆炸作用下的剩余承载力，p₀为高桩码头未施加荷载时的初始承载力。D 为0～0.2时，表示轻度损伤；D 为 0.2～0.5时，表示中度损伤；D为 0.5～0.8时，表示重度损伤；D为 0.8～1.0时，表示完全损伤。

高桩码头有限元模型计算结束后，分两步进行重启动计算：第1步读取高桩码头在气泡第1次脉动结束时刻的毁伤状态，删除流场载荷，待码头达到准静态；第2步在高桩码头上部1 mm放置刚性板（尺寸为1 500 cm× 1 100 cm× 50 cm，略大于码头面板尺寸），对其进行对称建模，采用Rigid刚体模型，以恒定速度（v₀= 0.01 m/s）的垂直竖向速度加载在码头模型上部，如图13所示，可以分别求得高桩码头在完好状态下的承载力和水下爆炸作用后高桩码头的剩余承载力，从而得到各工况下高桩码头的损伤评估系数。

图 13 重启动模型

Figure 13. Restart model

下载: 全尺寸图片幻灯片

通过爆轰因子反映水下爆炸冲击作用^[17]

${F} = \dfrac{{{W^{1/3}}}}{d}$

(10)

式中：W为炸药装药量，kg；d为炸药距离桩基的距离，m。表3为不同炸药当量作用下高桩码头的毁伤效应及评估，其中x_max为桩基最大变形，t_xmax为桩基变形最大时刻。

表 3 高桩码头毁伤评估

Table 3. Damage assessment of high-piled wharf

W/kg	F	x_max/cm	t_xmax/s	Damage phenomenon	D	Damage level
61.29	1.14	29.4	0.12	The middle and top part of the pile are slightly damaged. The bottom of the pile is slightly damaged. The connection of beams and piles is slightly damaged. The longitudinal beam and the wharf panel are not damaged.	0.06	Slight
117.45	1.42	57.7	0.15	The middle part and the blasting face in top part of pile are damaged, and the reinforcement is exposed. The bottom of the pile is slightly damaged. The connection of beams and piles is slightly damaged.	0.49	Moderate
272.70	1.88	113.9	0.25	The concrete at the top and middle of the pile are completely damaged, and the reinforcement is exposed. the bottom of the pile, longitudinal beam, transverse beamand panel are all damaged. The connection of beams and piles is damaged.	0.76	Severe

下载: 导出CSV

| 显示表格

4. 结　论

基于数值仿真方法建立水下爆炸高桩码头耦合模型，对水下爆炸作用下高桩码头的破坏过程和损伤机理进行了分析，探析了不同炸药当量高桩码头的破坏模式，评估了水下爆炸作用下高桩码头的毁伤效应，主要得出以下结论：

（1）冲击波阶段，由于冲击波持续时间短并发生绕射，桩基混凝土处于围压状态，冲击波对桩基的毁伤作用较弱；高桩码头在气泡第1次膨胀阶段迅速发生毁伤，并在气泡收缩阶段逐渐完成，在后续气泡脉动过程中高桩码头的毁伤现象基本维持不变。

（2）气泡脉动阶段，桩基随气泡膨胀收缩产生周期性往复变形，桩基中部两侧处于反复拉、压状态，桩顶与横梁连接处混凝土局部断裂，桩基整体呈现弯剪破坏。

（3）水下爆炸作用下，高桩码头的毁伤部位主要集中在桩基中部、顶部以及桩、梁连接处，其中桩基与横梁连接处为抗爆最薄弱部位，桩基浅水区的毁伤程度大于桩基深水区。

（4）随着炸药当量的增加，桩基中部和顶部的毁伤程度和毁伤范围增加，码头横、纵梁连接处以及码头面板相继产生不同程度损伤。

（5）在Mark81、Mark82、Mark83战斗部炸药当量的作用下，高桩码头的损伤评估分别为：轻度损伤、中度损伤和重度损伤，在Mark83战斗部炸药当量作用下炸药近场桩基顶部和中部混凝土完全失效，导致局部桩基丧失承载力。

图 1 (a) CsGeBr₃的SEM和EDS图像，(b) 常压下CsGeBr₃的XRD谱与PDF标准卡的对比结果

Figure 1. (a) SEM and EDS images of CsGeBr₃; (b) XRD pattern of CsGeBr₃ under ambient pressure compared with PDF card

下载: 全尺寸图片幻灯片

图 2 (a) CsGeBr₃的原位XRD谱随压力的变化，(b) 不同压力下CsGeBr₃的XRD环，(c) 不同压力下CsGeBr₃的XRD谱的Rietveld精修结果

Figure 2. (a) In situ XRD patterns of CsGeBr₃ as a function of pressure; (b) XRD rings of CsGeBr₃ at different pressures; (c) results of Rietveld refinement of XRD patterns of CsGeBr₃ at different pressures

下载: 全尺寸图片幻灯片

图 3 (a) CsGeBr₃的晶胞体积随压力的变化（低压相和高压相的p-V曲线均用Birch-Murnaghan状态方程拟），(b) CsGeBr₃的α、β和d_Ge-Br/(a/2)随压力的变化，(c) CsGeBr₃在常压和高压下的晶体结构

Figure 3. (a) Formula unit cell volume of CsGeBr₃ as a function of pressure (The p-V curves for low-pressure phase and high-pressure phase are both fitted with the Birch-Murnaghan equation of state); (b) α, β and d_Ge-Br/(a/2) of CsGeBr₃ at various pressures; (c) crystal structure of CsGeBr₃ under ambient pressure and high-pressure

下载: 全尺寸图片幻灯片

图 4 (a) 不同压力下CsGeBr₃的原位拉曼光谱，(b) 30～250 cm⁻¹范围内拉曼峰位随压力的变化，(c) 单晶CsGeBr₃的原位紫外-可见-近红外吸收谱，(d) CsGeBr₃在不同压力下的光学带隙（插图为常压下的Tauc图），(e) 加压过程中CsGeBr₃的光学影像

Figure 4. (a) In situ Raman spectra of CsGeBr₃ at diffrent pressures; (b) pressure dependence of the Raman peak positions in the wavenumber range of 30–250 cm⁻¹; (c) in situ ultraviolet-visible-near infrared (UV-Vis-NIR) absorption spectra of single crystal CsGeBr₃ under compression; (d) pressure dependence of the bandgaps of CsGeBr₃ (The illustration shows the Tauc plot for ambient pressure); (e) optical photos of CsGeBr₃ during compression

下载: 全尺寸图片幻灯片

图 5 (a) 三阶Birch-Murnaghan状态方程拟合的E-V曲线，(b) $R3m$ 和 $Pm\overline{3}m$ 相的焓差ΔH与压力的变化关系

Figure 5. (a) E-V curve fitted by third-order Birch-Murnaghan EOS; (b) relative enthalpies ΔH of $R3m$ and $Pm\overline{3}m$ phases vs. pressure

下载: 全尺寸图片幻灯片

表 1 不同压力下 CsGeBr₃的XRD谱的Rietveld精修结果

Table 1. Rietveld refinement results of XRD patterns for CsGeBr₃ at different pressures

Pressure/GPa	Crystal system	Space group	a/Å	α/(°)		β/(°)
0.80	Rhombohedral	R3m	5.544 43(12)	89.137 6(33)		94.018 1(30)
1.43	Cubic	${ Pm\overline{3}m }$	5.462 94(19)	90		90
3.01	Cubic	${Pm\overline{3}m }$	5.358 09(14)	90		90

Pressure/GPa	d_Ge-Br/Å	$d'_{{\mathrm{Ge\text{-}Br}}}$ /Å	V/Å³	Z	R_p/%	R_wp/%
0.80	2.585 19(6)	2.968 27(6)	170.382(11)	1	1.16	1.89
1.43	2.731 47(10)		163.034(17)	1	1.53	2.45
3.01	2.679 05(7)		153.826(12)	1	1.12	1.82

下载: 导出CSV

参考文献(45)

[1]	SAPAROV B, MITZI D B. Organic-inorganic perovskites: structural versatility for functional materials design [J]. Chemical Reviews, 2016, 116(7): 4558–4596. doi: 10.1021/acs.chemrev.5b00715
[2]	MITZI D B, CHONDROUDIS K, KAGAN C R. Organic-inorganic electronics [J]. IBM Journal of Research and Development, 2001, 45(1): 29–45. doi: 10.1147/rd.451.0029
[3]	LI W, WANG Z M, DESCHLER F, et al. Chemically diverse and multifunctional hybrid organic-inorganic perovskites [J]. Nature Reviews Materials, 2017, 2(3): 16099. doi: 10.1038/natrevmats.2016.99
[4]	MANSER J S, CHRISTIANS J A, KAMAT P V. Intriguing optoelectronic properties of metal halide perovskites [J]. Chemical Reviews, 2016, 116(21): 12956–13008. doi: 10.1021/acs.chemrev.6b00136
[5]	GOLDSCHMIDT V M. Die gesetze der krystallochemie [J]. Naturwissenschaften, 1926, 14(21): 477–485. doi: 10.1007/BF01507527
[6]	SHAMSI J, URBAN A S, IMRAN M, et al. Metal halide perovskite nanocrystals: synthesis, post-synthesis modifications, and their optical properties [J]. Chemical Reviews, 2019, 119(5): 3296–3348. doi: 10.1021/acs.chemrev.8b00644
[7]	FU Y P, ZHU H M, CHEN J, et al. Metal halide perovskite nanostructures for optoelectronic applications and the study of physical properties [J]. Nature Reviews Materials, 2019, 4(3): 169–188. doi: 10.1038/s41578-019-0080-9
[8]	CORREA-BAENA J P, SALIBA M, BUONASSISI T, et al. Promises and challenges of perovskite solar cells [J]. Science, 2017, 358(6364): 739–744. doi: 10.1126/science.aam6323
[9]	PETRUS M L, SCHLIPF J, LI C, et al. Capturing the sun: a review of the challenges and perspectives of perovskite solar cells [J]. Advanced Energy Materials, 2017, 7(16): 1700264. doi: 10.1002/aenm.201700264
[10]	LIN K B, XING J, QUAN L N, et al. Perovskite light-emitting diodes with external quantum efficiency exceeding 20 percent [J]. Nature, 2018, 562(7726): 245–248. doi: 10.1038/s41586-018-0575-3
[11]	WANG N N, CHENG L, GE R, et al. Perovskite light-emitting diodes based on solution-processed self-organized multiple quantum wells [J]. Nature Photonics, 2016, 10(11): 699–704. doi: 10.1038/nphoton.2016.185
[12]	SAIDAMINOV M I, ABDELHADY A L, MURALI B, et al. High-quality bulk hybrid perovskite single crystals within minutes by inverse temperature crystallization [J]. Nature Communications, 2015, 6: 7586. doi: 10.1038/ncomms8586
[13]	XING J, YAN F, ZHAO Y W, et al. High-efficiency light-emitting diodes of organometal halide perovskite amorphous nanoparticles [J]. ACS Nano, 2016, 10(7): 6623–6630. doi: 10.1021/acsnano.6b01540
[14]	AHMADI M, WU T, HU B. A review on organic-inorganic halide perovskite photodetectors: device engineering and fundamental physics [J]. Advanced Materials, 2017, 29(41): 1605242. doi: 10.1002/adma.201605242
[15]	DOU L T, YANG Y M, YOU J B, et al. Solution-processed hybrid perovskite photodetectors with high detectivity [J]. Nature Communications, 2014, 5: 5404. doi: 10.1038/ncomms6404
[16]	ZHOU Y X, HUANG Y Y, XU X L, et al. Nonlinear optical properties of halide perovskites and their applications [J]. Applied Physics Reviews, 2020, 7(4): 041313. doi: 10.1063/5.0025400
[17]	PARK N G. Perovskite solar cells: an emerging photovoltaic technology [J]. Materials Today, 2015, 18(2): 65–72. doi: 10.1016/j.mattod.2014.07.007
[18]	SONG J Z, LI J H, LI X M, et al. Quantum dot light-emitting diodes based on inorganic perovskite cesium lead halides (CsPbX₃) [J]. Advanced Materials, 2015, 27(44): 7162–7167. doi: 10.1002/adma.201502567
[19]	WANG Y M, LI M T, PEI C Y, et al. Critical current density and vortex phase diagram in the superconductor Sn_0.55In_0.45Te [J]. Physical Review B, 2022, 106(5): 054506. doi: 10.1103/PhysRevB.106.054506
[20]	LIU X Q, LI M T, ZHANG Q, et al. Pressure engineering promising transparent oxides with large conductivity enhancement and strong thermal stability [J]. Advanced Science, 2022, 9(31): 2202973. doi: 10.1002/advs.202202973
[21]	LI M T, ZHANG D J, HAN J, et al. Pressure-tuning structural and electronic transitions in semimetal CoSb [J]. Physical Review B, 2021, 104(5): 054511. doi: 10.1103/PhysRevB.104.054511
[22]	LI N N, FAN F R, SUN F, et al. Pressure-enhanced interplay between lattice, spin, and charge in the mixed perovskite La₂FeMnO₆ [J]. Physical Review B, 2019, 99(19): 195115. doi: 10.1103/PhysRevB.99.195115
[23]	LIU X Q, JIANG P, WANG Y M, et al. T_c up to 23.6 K and robust superconductivity in the transition metal δ-Ti phase at megabar pressure [J]. Physical Review B, 2022, 105(22): 224511. doi: 10.1103/PhysRevB.105.224511
[24]	YAN L M, DING C, LI M T, et al. Modulating charge-density wave order and superconductivity from two alternative stacked monolayers in a bulk 4Hb-TaSe₂ heterostructure via pressure [J]. Nano Letters, 2023, 23(6): 2121–2128. doi: 10.1021/acs.nanolett.2c04385
[25]	QU J, YAN L M, LIU H, et al. Pressure-induced structural phase transition in corundum-related class Cu₃TeO₆ [J]. High Pressure Research, 2021, 41(3): 318–327. doi: 10.1080/08957959.2021.1975699
[26]	LI N N, ZHANG Q, WANG Y G, et al. Perspective on the pressure-driven evolution of the lattice and electronic structure in perovskite and double perovskite [J]. Applied Physics Letters, 2020, 117(8): 080502. doi: 10.1063/5.0014947
[27]	LÜ X J, WANG Y G, STOUMPOS C C, et al. Enhanced structural stability and photo responsiveness of CH₃NH₃SnI₃ perovskite via pressure-induced amorphization and recrystallization [J]. Advanced Materials, 2016, 28(39): 8663–8668. doi: 10.1002/adma.201600771
[28]	LIN J, CHEN H, GAO Y, et al. Pressure-induced semiconductor-to-metal phase transition of a charge-ordered indium halide perovskite [J]. Proceedings of the National Academy of Sciences of the United States of America, 2019, 116(47): 23404–23409.
[29]	ASFIA M B, JAMAN S, RASHID M A. Pressure induced band gap shifting from ultra-violet to visible region of RbSrCl₃ perovskite [J]. Materials Research Express, 2022, 9(9): 095902. doi: 10.1088/2053-1591/ac8f88
[30]	LÜ X J, STOUMPOS C, HU Q Y, et al. Regulating off-centering distortion maximizes photoluminescence in halide perovskites [J]. National Science Review, 2021, 8(9): 288. doi: 10.1093/nsr/nwaa288
[31]	JAFFE A, LIN Y, MAO W L, et al. Pressure-induced metallization of the halide perovskite (CH₃NH₃)PbI₃ [J]. Journal of the American Chemical Society, 2017, 139(12): 4330–4333. doi: 10.1021/jacs.7b01162
[32]	TANG L C, HUANG J Y, CHANG C S, et al. New infrared nonlinear optical crystal CsGeBr₃: synthesis, structure and powder second-harmonic generation properties [J]. Journal of Physics: Condensed Matter, 2005, 17(46): 7275–7286. doi: 10.1088/0953-8984/17/46/011
[33]	HUANG L Y, LAMBRECHT W R L. Electronic band structure trends of perovskite halides: beyond Pb and Sn to Ge and Si [J]. Physical Review B, 2016, 93(19): 195211. doi: 10.1103/PhysRevB.93.195211
[34]	ZHANG Q Q, MUSHAHALI H, DUAN H M, et al. The linear and nonlinear optical response of CsGeX₃ (X = Cl, Br, and I): the finite field and first-principles investigation [J]. Optik, 2019, 179: 89–98. doi: 10.1016/j.ijleo.2018.10.159
[35]	LIN Z G, TANG L C, CHOU C P. Study on mid-IR NLO crystals CsGe(Br_xCl_{1– x})₃ [J]. Optical Materials, 2008, 31(1): 28–34. doi: 10.1016/j.optmat.2008.01.004
[36]	LIN Z G, TANG L C, CHOU C P. Characterization and properties of novel infrared nonlinear optical crystal CsGe(Br_xCl_{1– x})₃ [J]. Inorganic Chemistry, 2008, 47(7): 2362–2367. doi: 10.1021/ic7011777
[37]	HUANG L Y, LAMBRECHT W R L. Vibrational spectra and nonlinear optical coefficients of rhombohedral CsGeX₃ halide compounds with X = I, Br, Cl [J]. Physical Review B, 2016, 94(11): 115202. doi: 10.1103/PhysRevB.94.115202
[38]	SEO D K, GUPTA N, WHANGBO M H, et al. Pressure-induced changes in the structure and band gap of CsGeX₃ (X = Cl, Br) studied by electronic band structure calculations [J]. Inorganic Chemistry, 1998, 37(3): 407–410. doi: 10.1021/ic970659e
[39]	SCHWARZ U, WAGNER F, SYASSEN K, et al. Effect of pressure on the optical-absorption edges of CsGeBr₃ and CsGeCl₃ [J]. Physical Review B, 1996, 53(19): 12545–12548. doi: 10.1103/PhysRevB.53.12545
[40]	MAO H K, XU J, BELL P M. Calibration of the ruby pressure gauge to 800 kbar under quasi-hydrostatic conditions [J]. Journal of Geophysical Research, 1986, 91(B5): 4673–4676. doi: 10.1029/JB091iB05p04673
[41]	PRESCHER C, PRAKAPENKA V B. DIOPTAS: a program for reduction of two-dimensional X-ray diffraction data and data exploration [J]. High Pressure Research, 2015, 35(3): 223–230. doi: 10.1080/08957959.2015.1059835
[42]	LARSON A C, VON DREELE R B. General structure analysis system (GSAS): No. LAUR 86–748 [R]. Los Alamos: Los Alamos National Laboratory, 2004.
[43]	KOHN W, SHAM L J. Self-consistent equations including exchange and correlation effects [J]. Physical Review, 1965, 140(4A): A1133–A1138. doi: 10.1103/PhysRev.140.A1133
[44]	BLÖCHL P E. Projector augmented-wave method [J]. Physical Review B, 1994, 50(24): 17953–17979. doi: 10.1103/PhysRevB.50.17953
[45]	KRESSE G, FURTHMÜLLER J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set [J]. Physical Review B, 1996, 54(16): 11169–11186. doi: 10.1103/PhysRevB.54.11169