Analysis of the dynamic response of prestressed concrete frame structures under blast load
-
摘要: 为了研究预应力混凝土(prestressed concrete, PC)框架结构的抗爆性能，利用有限元软件LS-DYNA对一栋3层2跨的大跨有/无黏结PC框架结构在不同比例距离的外部远爆荷载作用下的动力响应进行了分析。分析结果表明：混凝土预应力框架在地表远爆荷载作用下，最大层间位移角与前墙所受峰值反射超压近似成线性关系；有黏结混凝土预应力框架结构层间位移角相较于无黏结混凝土预应力框架更小，损伤分布更均匀，结构抗爆性能更好；基于分析结果，给出了不同比例距离对应的损伤状态，可用于对混凝土预应力框架结构进行爆炸损伤状态快速评估。Abstract: To study the anti-blast performance of prestressed concrete (PC) frame structure, the dynamic response of a 3-story 2-span long-span bonded/unbonded PC frame structure under the action of external remote blast loads of different scaled distances was analyzed using the finite element software LS-DYNA. In the different blast conditions, the type of explosion was a surface detonation, the explosion distance was 100 m, and the scaled distances were 2−8 m/kg1/3. Explosive loads on the building surface were calculated according to unified facilities criteria 3-340-02. The different types of prestress in the numerical model were realized by controlling the direction of coupling of concrete and prestressing tendons. To validate the numerical model, the results of the numerical simulations were compared with the experimental results. The blast resistance mechanism, story drift ratios, structural damage mode, and damage assessment of PC frames were analyzed under blast loads at different scaled distances. The analysis results show that the ground floor of the large-span PC frame is easy to damage under the action of the external remote blast load, and the ground floor columns can be strengthened to improve the overall structural blast resistance. In the PC frame, columns relative to the frame beams are weaker and prone to damage. The PC frame structure with a low degree of redundancy is easy to collapse. The maximum story drift ratio of the prestressed concrete frame under the surface blast load is approximately linearly related to the peak reflected overpressure applied to the front wall. Compared with the unbonded prestressed concrete frame, the bonded prestressed concrete frame structure has a smaller story drift ratio, more uniform damage distribution, and better structural blast resistance. Based on the analysis results, the damage states at different scaled distances was given, which can realize the rapid assessment of the explosion damage state of prestressed concrete frame structures.
-
Key words:
- blast load /
- dynamic response analysis /
- prestress /
- concrete frame
-
泡沫铝是由铝或者铝合金基体与内部孔洞组成的新型结构与功能一体化材料，内部存在大量孔洞。与传统的合金材料相比，泡沫铝具有更低的体积密度、更强的比强度和比刚度，并且具有良好的缓冲吸能、阻尼减振、电磁屏蔽和隔音降噪等性能。因此，在航空航天、武器装备、建筑结构和交通运输等领域，泡沫铝具有广阔的应用前景。在实际使用过程中，泡沫铝夹层结构（aluminum foam sandwich，AFS）是最常用的结构形式，理论模型[1-3]表明这类夹芯结构具有优良的抗冲击性能。泡沫铝的孔隙率一般为30%～90%，孔径为0.3～7.0 mm。孔隙率或孔径以及孔洞分布形式的改变，对AFS的学性能有显著影响，研究AFS是力学及材料学科的热点问题之一。郭亚周等[4]采用一级空气炮开展冲击实验，研究了各种形状弹体冲击作用下不同泡沫铝夹层结构的动力性能；张博一等[5]采用一级轻气炮驱动泡沫铝子弹，研究了泡沫铝子弹冲击下泡沫铝夹层结构的动态响应。空气炮是研究AFS抗冲击性能最常用的实验手段，但是实验成本较高，难以广泛应用，因此数值模拟成为目前的重要方法。邓旭辉等[6]采用数值方法，分析了双层泡沫铝夹芯板的抗冲击性能；张永康等[7]采用LS-DYNA进行数值计算，对不同构型泡沫铝夹层结构进行了冰雹冲击分析，研究了中面板位置对缓冲吸能性能的影响；夏志成等[8]设计了落锤实验对泡沫铝夹芯板进行研究，并应用数值模拟软件还原夹芯板冲击过程，计算了泡沫铝结构的吸能效果；Cui 等[9]采用工业CT扫描技术，重建了泡沫铝夹芯结构三维细观有限元模型，分析了不同构型弹道冲击下泡沫铝夹芯结构的失效模式，并在实验和数值观测结果的基础上，建立了预测弹丸残余速度的理论模型。在防爆结构领域，泡沫铝也引起相关学者的注意。苏兴亚等[10]采用弹道冲击摆系统，开展了爆炸载荷下分层梯度泡沫铝夹芯板的变形、失效模式和抗冲击性能的实验研究；王涛等[11]开展了不同爆炸载荷作用下泡沫铝夹芯板变形的实验研究，获得了泡沫铝夹芯结构背板中心点挠度与爆炸冲量之间的关系。
现有的有限元方法能有效计算泡沫铝夹层结构的力学性能，预测其准静态及动态响应，但是传统的有限元方法是从连续介质力学出发，用偏微分方程来描述材料的力学行为，在处理泡沫铝夹层结构界面开裂、芯材溃裂、背板材料裂纹扩展等非连续问题时存在瓶颈。近场动力学（peridynamic，PD）理论[12-14]是一种通过求解空间积分方程来描述材料力学行为的非局部数值方法，因它在解决裂纹扩展、材料破坏、复合材料渐进损伤以及多尺度问题上独特的优势，越来越多的学者将PD方法引入到各种领域解决传统方法难以处理的问题。杨娜娜等[15]采用近场动力学，对破片冲击下复合材料的损伤进行了模拟，分析了冲击速度、复合材料铺层方式等对结构损伤的影响；刘宁等[16]采用短程斥力模型，描述了碰撞过程，模拟了冲击压缩条件下单裂纹圆孔板动态破坏行为；熊伟鹏等[17]采用近场动力学方法，模拟了冰球碰撞过程，并与实验对比取得良好的效果。基于键的近场动力学方法不支持塑性材料本构模型，多用于脆性材料破坏问题。而常规态型近场动力学方法解耦了几何形状变形和体积变形，可以实现材料塑性变形[18]，因此适用于泡沫铝这类具有良好塑性的材料。目前，我国在近场动力学塑性模型方面的研究还较罕见。
本文中，为了明确孔隙率对泡沫铝夹层结构的抗冲击性能的影响，采用Monte-Carlo方法生成泡沫铝的随机模型，基于常规态型近场动力学理论，引入Mises屈服准则和线性各向同性强化模型来描述铝的塑性行为，对泡沫铝夹层结构进行冲击计算，分析孔隙率分别为0.7、0.6、0.5和0.4的泡沫铝夹层结构的缓冲效能以及破坏形态，以期为此类问题提供借鉴。
1. 常规态型近场动力学
1.1 运动方程
如图1所示，质点x与其半径为δ的邻域Hx内其他质点存在相互作用。当材料发生变形时，质点x与其邻域内的另一质点x'分别移动到y和y'，质点x'对质点x产生力密度t，力密度矢量t的方向与质点之间的连线重合，其大小与质点x'的变形以及质点之间的相对位移相关。同理，质点x'也受到其邻域Hx'内其他质点的影响，因此质点x也对质点x'产生力密度矢量t'，其方向与t平行，大小与质点x'的变形以及质点之间的相对位移相关。质点x的运动状态是在与邻域Hx内所有其他质点相互影响下产生，在邻域内对力密度进行积分，可以得到关于质点x的拉格朗日方程[12]：
\rho({\boldsymbol{x}}) \ddot{{\boldsymbol{u}}}({\boldsymbol{x}}, t)=\int_{H_{x}}\left({\boldsymbol{t}}\left({\boldsymbol{u}}^{\prime}-{\boldsymbol{u}}, {\boldsymbol{x}}^{\prime}-{\boldsymbol{x}}, t\right)-{\boldsymbol{t}}^{\prime}\left({\boldsymbol{u}}-{\boldsymbol{u}}^{\prime}, {\boldsymbol{x}}-{\boldsymbol{x}}^{\prime}, t\right)\right) {\rm d} H+{\boldsymbol{b}}({\boldsymbol{x}}, {{t}}) (1) 近场动力学的运动方程不涉及求解空间坐标的偏微分方程，取而代之的是通过求解材料内部质点影响范围内的空间积分来描述材料宏观变形。常用的数值计算方法是将材料离散成空间中均匀分布的一系列质点，质点坐标代表了其所在区域的体积，于是式(1)可写成如下离散形式[13]：
\rho\left({\boldsymbol{x}}_{k}\right) \ddot{{\boldsymbol{u}}}\left({\boldsymbol{x}}_{k}, t\right)=\sum_{j=1}^{N}\left[{\boldsymbol{t}}_{k, j}\left({\boldsymbol{u}}_{j}-{\boldsymbol{u}}_{k}, {\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}, t\right)-{\boldsymbol{t}}_{j,k}\left({\boldsymbol{u}}_{k}-{\boldsymbol{u}}_{j}, {\boldsymbol{x}}_{k}-{\boldsymbol{x}}_{j}, t\right)\right]V_{j}+{\boldsymbol{b}}{({\boldsymbol{x}}_k}, {{t}}) (2) 式中：xk和xj为质点坐标，并且xj在xk的影响范围内；N为质点xk影响范围内的质点总数；ρ(xk)为质点的密度；uk为质点xk的位移向量；
\ddot{{\boldsymbol{u}}}\left({\boldsymbol{x}}_{k}, t\right) 为质点xk的加速度矢量；b(xk, t)为作用在质点xk上的外体力密度矢量；Vj为质点xj的体积。1.2 弹脆性材料本构模型
质点间相互作用的力密度可以通过近场动力学应变能密度导出[12-14]：
t_{k, j}=2\left(a_{k}-a_{\mu}\right) \frac{\delta d}{\left|{\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}\right|} \varLambda_{k, j} \theta_{k}+2 \delta b s_{k, j} (3) 式中：θk和sk, j分别为质点xk的体积应变和质点xk、xj之间的伸长率。θk和sk, j分别定义为：
\theta_{k}=d \sum_{j=1}^{N} \delta s_{{k}, j} \varLambda_{{k}, j} V_{j} (4) s_{k, j}=\frac{\left|{\boldsymbol{y}}_{j}-{\boldsymbol{y}}_{k}\right|-\left|{\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}\right|}{\left|{\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}\right|} (5) 参数Λk, j定义为：
\varLambda_{k, j}=\frac{{\boldsymbol{y}}_{j}-{\boldsymbol{y}}_{k}}{\left|{\boldsymbol{y}}_{j}-{\boldsymbol{y}}_{k}\right|} \cdot \frac{{\boldsymbol{x}}_{j} -{\boldsymbol{x}}_{k}} {\left|{\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}\right|} (6) 通过考虑简单加载情况下的应变能密度，可以将式(3)中的近场动力学参数aκ、aμ、b和d与传统连续介质力学中的体积模量κ和剪切模量μ建立联系，推导出[14, 18]：
\left\{ \begin{array}{l} a_{\kappa}=\dfrac{1}{2} \kappa,\;\;\; a_{\mu}=\dfrac{5 \mu}{6}, \;\;\;b=\dfrac{15 \mu}{2 {\text{π}} \delta^{4}},\;\;\;\;\; d=\dfrac{9}{4 {\text{π}} \delta^{4}} \;\;\;\;\;\;\;\text { 3D } \\ a_{\kappa}=\kappa, \;\;\;a_{\mu}=2 \mu, \;\;\;b=\dfrac{6 \mu}{{\text{π}} h \delta^{5}}, \;\;\;\;\;d=\dfrac{2}{{\text{π}} h \delta^{3}}\;\;\;\;\;\;\;\;\;\;\;\;\mathrm{2D} \\ a_{\kappa}=0, \;\;\;a_{\mu}=0, \;\;\;b=\dfrac{E}{2 \delta^{3} A}, \;\;\;\;\; d=\dfrac{1}{2 \delta^{2} A} \;\;\;\;\;\;\;\;\;\;\;\;\text { 1D } \end{array}\right. (7) 近场动力学中，材料损伤失效是通过截断物质点之间的相互作用描述的，并且截断之后不可恢复，因此，在力密度和体积应变中引入一个与时间相关的标量函数
\varphi 进行修正，其定义为：\varphi\left({\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}, t\right)=\left\{ \begin{array}{lll} 1\;\;\;\;\;s_{k, j}{\text{＜}}s_{\text{c}} \; \;({完好}) \\ 0\;\;\;\;\;s_{k, j}{\text{＞}}s_{\text{c}} \;\; {(断裂 }) \end{array}\right. (8) 式中：sc为临界伸长率。当质点之间的伸长率超过临界伸长率，则质点间的相互作用永久性破坏，这意味着这2个质点之间的微势能因裂纹面产生而消失。于是，得到添加损伤修正的力密度和体积应变：
t_{k, j}=2\left(a_{\kappa}-a_{\mu}\right) \frac{\delta d}{\left|{\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}\right|} \varLambda_{k, j} \theta_{k}+2 \delta b \varphi\left({\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}, t\right) s_{k, j} (9) \theta_{k}=d \delta \sum_{j=1}^{N} \varphi\left({\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}, t\right) s_{k, j} \varLambda_{k, j} V_{j} (10) 对穿过裂纹面的微势能在区域内进行积分，与材料断裂的临界能量释放率建立联系，可以得到临界伸长率的表达式为[12-14]：
s_{\text{c}}=\left\{ \begin{array}{l} \sqrt{\dfrac{G_{\text{c}}}{\left[3\mu+\left(\dfrac{3}{4}\right)^4\left(\kappa-\dfrac{5\mu}{3}\right)\right]\delta}}\;\;\;\;\;\;\;\;\;\;{\rm{3D}}\\ \sqrt{\dfrac{G_{\text{c}}}{\left[\dfrac{6}{{\text{π}}}+\dfrac{16}{9{\text{π}}^2}(\kappa-2\mu) \right]\delta}}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{\rm{2D}} \end{array}\right. (11) 实际上，对于复杂的材料目前仍无法简单地计算其临界伸长率，此时便需要采用实验标定。通过统计质点影响范围内截断的相互作用占初始相互作用总数的比例，可表征局部的损伤程度：
\phi\left({\boldsymbol{x}}_{k}, t\right)=1-\frac{\displaystyle\sum_{j=1}^{N} \varphi\left({\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}, t\right) V_{j}}{\displaystyle\sum_{j=1}^{N} V_{j}} (12) 1.3 近场动力学塑性理论
塑性变形过程与加载路径相关，并且是不可逆的，通常需要采用增量的方式描述载荷与变形之间的关系。将材料变形分解成可恢复的弹性部分和不可恢复的塑性部分，于是可以将质点之间的伸长率增量Δsk, j分解为：
\Delta s_{k, j}=\Delta s_{k, j}^{{\rm{e}}}+\Delta s_{k, j}^{\rm p} (13) 式中：
\Delta s_{k, j}^{\rm e} 和\Delta s_{k, j}^{\rm p} 分别为质点之间的弹性伸长率增量和塑性伸长率增量。当采用步进的方式加载时，第n+1个加载步中质点xk和xj之间的总伸长率为：s_{k, j}^{(n+1)}=s_{k, j}^{(n)}+\Delta s_{k, j}^{(n+1)} =\left(s^{{\rm{e}}(n)}_{k,j}+\Delta s^{{\rm{e}}(n+1)}_{k,j}\right)+\left(\Delta s^{{\rm{p}}(n+1)}_{k,j}\right) (14) 式中：
s_{k, j}^{\rm e} 、s_{k_, j}^{\rm p} 和s_{k, j} 分别为质点之间的弹性伸长率、塑性伸长率和总伸长率。假设在第n+1个加载步中，材料仅存在弹性变形增量，塑性变形增量\Delta s_{k, j}^{{\rm{p}}(n+1)}=0 ，那么弹性伸长率的试算值为：s_{k, j({\rm{trial}})}^{{\rm{e}}(n+1)}=s_{k, j}^{(n+1)}-s_{k, j}^{{\rm p}(n)} (15) 此时，根据式(3)，可以得到力密度的试算值为：
t^{(n+1)}_{k,j({\rm{trial}})}=2(a_{\kappa} -a_{\mu})\frac{\delta d}{|{\boldsymbol{x}}_j-{\boldsymbol{x}}_k|}\varLambda_{k,j}\theta_k \left(s^{{\rm{e}}(n+1)}_{k,j({\rm{trial}})}\right)+2\delta bs^{{\rm{e}}(n+1)}_{k,j({\rm{trial}})} (16) 根据Mises屈服准则，当材料内一点的应力偏量第二不变量达到某一值时，该点便开始进入塑性状态，当考虑线性各向同性强化时，可以建立屈服函数[18]：
F_{k}^{(n+1)}=b \sum_{j=1}^{N} \frac{\delta}{\left|{\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}\right|}\left(s_{k, j}^{{\rm{e}}(n+1)}\left|{\boldsymbol{x}}_{j}-{\boldsymbol{x}}_{k}\right|\right)^{2} V_{j}-a_{\mu}\left(\theta_{k}^{{\rm{e}}(n+1)}\right)^{2}-G\left(\overline s_{k}^{{\rm{p}} (n+1)}\right) (17) 硬化函数定义为：
G\left(\bar{s}_{k}^{\rm p}\right)=\frac{\left(\sigma_{\rm y}+E_{\rm t} \bar{s}_{k}^{\rm p}\right)^{2}}{6 \mu} (18) 式中：
\bar{s}_{k}^{\rm p} 和Et分别为材料的等效塑性应变和切线模量。将弹性伸长率的试算值
s^{{\rm{e}}(n+1)}_{k,j({\rm{trial}})} 代入屈服函数，当F_{k}^{(n+1)} {\text{≤}} 0 时，该加载步位于屈服面内，塑性伸长率增量\Delta s_{k, j}^{{\rm p}(n+1)}=0 ，弹性伸长率的试算值s^{{\rm{e}}(n+1)}_{k,j({\rm{trial}})} 即为真实发生的弹性伸长率。否则，加载状态位于屈服面上方，\Delta s_{k, j}^{\mathrm{p} (n+1)} \neq 0 ，需要确定真实发生的弹性伸长率以使F_{k}^{(n+1)} 退回到屈服面内：\left\{\begin{array}{l}s^{{\rm e}(n+1)}_{k,j}=s^{{\rm e}(n+1)}_{k,j({\rm{trial}})}-\Delta s^{{\rm p}(n+1)}_{k,j}\\ \Delta s^{{\rm p}(n+1)}_{k,j}=C_k\left(\dfrac{1}{2\delta b}|{\boldsymbol{x}}_j-{\boldsymbol{x}}_k|t^{(n+1)}_{k,j({\rm{trial}})}+\dfrac{d}{b}(a_{\kappa}-a_{\mu})\varLambda_{k,j} \theta_k\left( s^{{\rm e}(n+1)}_{k,j({\rm{trial}})} \right) \right) \end{array} \right. (19) 式中：Ck为一个大于0的系数，使得将上式中的弹性伸长率代入屈服函数后
F_{k}^{(n+1)} \to 0 。Ck可以采用牛顿迭代法确定，最终得到此加载步材料各点的塑性伸长率和弹性伸长率分别为：\left\{\begin{array}{l}s^{{\rm p}(n+1)}_{k,j}=s^{{\rm p}(n+1)}_{k,j}-\Delta s^{{\rm p}(n+1)}_{k,j}\\ s^{{\rm e}(n+1)}_{k,j}=s^{(n+1)}_{k,j}- s^{{\rm p}(n+1)}_{k,j} \end{array} \right. (20) 2. 泡沫铝夹层结构模型
如图2所示，光学舱防护结构采用的是泡沫铝夹层复合材料，由上层铝蒙皮面板、中间泡沫铝夹层以及底层有机玻璃背板组成。泡沫铝的孔隙率P有0.7、0.6、0.5和0.4等4种，所选用泡沫铝主孔径为1～6 mm。泡沫铝夹层结构水平放置在刚性平台上，上方有一个直径为25 mm的刚性小球以20 m/s的速度向下冲击。为了便于建模并提高计算效率，将模型简化成二维结构，简化成二维的刚性小球质量为3.5×10−3 kg。采用图3所示的Monte-Carlo方法随机生成泡沫铝模型，分3级进行胞孔投放，第1级孔径为4～6 mm，第2级孔径为2～4 mm，第3级孔径为1～2 mm，先投放大孔径后投放小孔径，以达到提高投放效率并尽可能获得更高的孔隙率的目的。将模型离散成均匀分布的球形粒子，为确保数值计算精度，需保证最小胞孔孔径大于近场范围，确定粒子直径Δx=0.2 mm，近场半径δ=3Δx。最终建立的4种不同孔隙率泡沫铝夹层复合结构的近场动力学离散模型如图4所示。
图 2 破片冲击泡沫铝夹层结构模型Figure 2. Schematic diagram of an aluminum foam sandwich impacted by a fragment
图 3 Monte-Carlo方法生成泡沫铝算法流程Figure 3. Algorithm flow chart of the Monte-Carlo method to generate aluminum foam
图 4 泡沫铝夹层结构的近场动力学离散模型Figure 4. Discrete models of aluminum foam sandwiches for peridynamics泡沫铝夹层结构相邻层之间用环氧树脂进行粘接，计算中不考虑环氧树脂的质量，仅通过环氧树脂的力学性能定义相邻材料界面之间的近场动力学参数，并对材料界面及边界处的参数进行修正[12-14, 18]。实验中，以20 m/s速度冲击时，铝板蒙皮及泡沫铝夹层主要表现为塑性变形，并且未观察到明显破裂，因此在这种工况下可以采用近场动力学塑性模型描述铝材料。有机玻璃和环氧树脂表现为脆性破坏，因此用近场动力学弹脆性模型描述。计算中采用的各材料力学参数如表1所示。
表 1 材料参数Table 1. Material parameters材料 ρ/(kg∙m−3) E/GPa ν σy/MPa Et/MPa sc 铝 2700 69.5 0.33 127 586 有机玻璃 1190 3.6 0.4 0.047 环氧树脂 3 0.37 0.047 3. 数值计算结果及分析
3.1 泡沫铝夹层结构变形形态
图5显示了泡沫铝夹层结构的变形过程，云图颜色表示各点竖直方向位移uy的绝对值：在冲击作用下，先是正对冲击物的部位铝面板发生局部凹陷，凹陷深度不断增大，内部泡沫铝逐渐发生大范围变形，随着泡沫铝变形范围进一步增大，底部有机玻璃背板在上层泡沫铝带动下发生弯曲，结构呈现左右两端向上翘起的整体变形。分析可知，这是由于泡沫铝内部的大量孔洞能大幅削减应力从冲击部位向周围的扩散范围，因此在冲击初期局部变形集中在结构与冲击物的接触部位，随着泡沫铝的孔洞被压实，结构内部应力积累到一定程度才导致发生整体变形。这是泡沫铝能起到缓冲作用的重要原因之一。
此外，影响夹层结构抗冲击性能的另一个主要因素是泡沫铝芯材塑性变形耗能[8]。冲击过程中，随着上层的泡沫铝塑性变形范围增大，冲击物的动能逐渐被耗散，最终使得作用在背板上的冲击力极大缓解，从而达到保护底部背板的目的。由此可认为，冲击作用下泡沫铝芯材的塑性变形程度与范围越大，对背板的保护作用越强，则泡沫铝夹层结构抗冲击性能更好。图6为具有不同孔隙率芯材的夹层结构在冲击结束后的塑性变形形态，可以看出，结构上越靠近冲击物的部位塑性变形程度越大，塑性变形呈以冲击凹坑为中心的放射状分布，并且孔隙率越高，冲击形成的凹坑越深，塑性变形范围越大。
3.2 冲击物的速度与加速度分析
通过分析冲击物在冲击过程中的速度变化以及加速度变化过程，可以了解冲击过程中动能传递规律，间接了解泡沫铝芯材的缓冲作用。图7是冲击物冲击不同孔隙率泡沫铝夹层结构时的速度v变化曲线，孔隙率越高，速度变化曲线越平缓。此外，根据表2的分析结果，孔隙率越高，冲击物的反弹速度越低，冲击结束时冲击物具有更低的残余动能，这说明高孔隙率的泡沫铝通过塑性变形消耗了更多冲击能量。当孔隙率P从0.4提升到0.7时，泡沫铝对冲击物的动能吸收率从90%提高到99%。图8是冲击物冲击不同孔隙率泡沫铝夹层结构时的加速度a变化曲线，可见孔隙率越低，加速度变化曲线越陡峭、峰值越高，冲击碰撞过程越偏向刚性。通过以上分析，可以进一步明确，泡沫铝芯材孔隙率对夹层结构的抗冲击性能有显著影响，在一定范围内，孔隙率越高，缓冲效果越好，抗冲击性能越强。需要注意的是，上述规律在孔隙率进一步提高时将不再适用，虽然提高孔隙率能增强缓冲能力，但过大的孔隙率会造成泡沫铝芯材承载能力显著下降，无法充分发挥其防护作用。
表 2 冲击过程中冲击物的主要运动参数Table 2. Main motion parameters of impact object in the process of impact孔隙率 反弹速度/(m∙s−1) 残余动能/mJ 加速度峰值/(km∙s−2) 动能吸收率/% 0.7 0.68 8 47 99 0.6 2.91 15 53 97 0.5 4.72 39 66 94 0.4 6.08 65 85 90 采用图9所示的落锤实验系统获取冲击物加速度峰值，并与模拟计算结果对比，图10显示模拟计算值稍高于实验值，但是与孔隙率之间关系的变化规律是一致的。出现这一情况的主要原因是计算中将问题简化为二维模型导致的差异，此外，模拟计算中未考虑泡沫铝和铝蒙皮面板的微观裂纹扩展情况也是误差的重要来源。尽管存在这一偏差，模拟计算值与实验结果总体量级上吻合良好，并且表现出相同的规律，证明了近场动力学模拟分析结果与结论的有效性。
图 9 采用落锤实验获取冲击物加速度峰值Figure 9. The peak impact acceleration of the impactor obtained by a drop-weight experiment3.3 有机玻璃背板的裂纹范围
某光学舱采用泡沫铝作为防护结构是为了保护底层光学材料不被破坏，考察泡沫铝夹层结构抗冲击性能最直观的方式便是观察底部有机玻璃背板的破坏形态。近场动力学中采用临界伸长率预测脆性材料断裂问题的方法已被证明是有效的，因此，基于此方法计算出有机玻璃背板的裂纹起始与终止形态如图11所示。
如图11所示，冲击过程中，有机玻璃背板最先发生破坏的位置是中间正对冲击的部位，并逐渐向两侧扩展。随着泡沫铝孔隙率的降低，有机玻璃背板的裂纹扩展范围明显增大，说明高孔隙率的泡沫铝比低孔隙率的泡沫铝有更好的防护效果，进一步验证高孔隙率的泡沫铝夹层结构具有更好的抗冲击性能。此外，上下粘接界面也出现了一定损伤，结合3.1节可知，这是由于在冲击过程中上层铝面板首先发生塑性变形，随着铝面板下凹，铝面板和泡沫铝之间的粘接界面势必会因变形不协调而产生横向剪切力，这是导致上层粘接界面开裂的主要原因，随着泡沫铝芯材变形增加，引起结构发生整体弯曲变形，进而使泡沫铝和有机玻璃背板之间的粘接界面产生剪切和张拉应力，导致粘接界面出现损伤。从图11可以看出不同孔隙率的泡沫铝粘接界面损伤程度基本一致，说明孔隙率对粘接界面开裂程度的影响较小。
4. 结 论
在常规态型近场动力学理论中引入了Mises屈服准则和线性各向同性强化模型，建立了近场动力学塑性本构的数值计算框架。采用随机方法创建了泡沫铝结构的二维细观模型，并基于近场动力学计算程序模拟了冲击作用下不同孔隙率泡沫铝的塑性变形以及有机玻璃背板的裂纹扩展形态，验证了近场动力学解决泡沫铝夹层结构抗冲击问题的可行性，为相关领域问题提供了新的思路。探讨了在20 m/s的低速冲击作用下，孔隙率对泡沫铝夹层结构抗冲击性能的影响规律，得到以下结论。
(1)泡沫铝夹层结构良好的塑性变形能力是其发挥缓冲与防护作用的主要因素。一方面，冲击作用下泡沫铝结构先表现出局部变形，随着泡沫铝内部孔洞被压实，结构逐渐表现出整体弯曲变形，在这一变形过程中，泡沫铝内部孔洞能大幅削减应力从冲击部位向周围的扩散范围，起到缓冲作用；另一方面，泡沫铝芯材发生大范围塑性变形，冲击物的动能在这一过程中被耗散，使作用在背板上的冲击力被缓解。
(2)在一定范围内，泡沫铝芯材孔隙率越高，则夹层结构具有越好的抗冲击性能。孔隙率从0.4升高到0.7时，泡沫铝对冲击物的动能吸收率从90%提高到99%，冲击物加速度峰值从85 km/s2降低到47 km/s2。模拟结果与实验结果具有较好的一致性，验证了模拟结果的准确性和分析结论的有效性。
(3)近场动力学方法能有效预测有机玻璃背板的裂纹扩展以及粘接界面开裂过程，模拟结果表明，泡沫铝芯材孔隙率越高，有机玻璃背板的裂纹扩展范围越小，进一步说明本文所研究的某光学舱防护结构应尽量选择更高孔隙率的泡沫铝。
-
图 2 首层中柱尺寸和配筋
Figure 2. Size and reinforcement of the middle column of the first floor
图 5 大跨PC框架尺寸和配筋(单位：mm)
Figure 5. Dimensions and reinforcement details of large span PC frame (unit: mm)
图 10 工况1中结构的破坏模式
Figure 10. Collapse process of the structure under working condition 1
图 11 工况2中结构的倒塌过程
Figure 11. Collapse process of the structure under working condition 2
图 13 各工况层间位移角时程
Figure 13. Time histories of story drift ratio under each working condition
图 15 前墙峰值超压对层间位移角的影响
Figure 15. Effect of overpressure of the front wall on interlayer drift
表 1 钢筋材料参数
Table 1. Material properties of reinforcement
钢筋 面积/mm2 屈服强度/MPa 极限强度/MPa W0.5 3.22 441 513 D1 6.45 399 610 D5 32.20 449 513 
下载: 导出CSV
表 2 混凝土材料参数
Table 2. Material properties of concrete
龄期/d 弹性模量/GPa 密度/(kg·m−3) 抗压强度/MPa 28 28.7 2068 42.0 103～132 30.3 2068 44.7
下载: 导出CSV
表 3 各工况的爆炸参数
Table 3. Blast load parameters
工况 Z/(m·kg−1/3) W/kg pr,front/kPa ta/ms trf/ms pr,roof/kPa td//ms tof/ms pr,rear/kPa t'_{\rm{a}} /ms t'_{\rm{b}} /ms t'_{\rm{of}} /ms 1 2.0 125000 1056.0 84.8 119 91.80 115 206 90.30 149 188 269 2 3.0 37037 331.0 118.0 163 39.70 167 275 40.10 199 242 323 3 3.5 23323 223.0 132.0 180 28.50 184 301 30.20 218 261 341 4 4.0 15625 162.0 145.0 194 21.60 196 315 23.60 234 277 357 5 5.0 8000 101.0 165.0 215 13.70 219 337 15.80 259 303 375 6 6.0 4630 71.3 181.0 229 9.53 231 348 11.50 277 324 390 7 7.0 2916 54.7 193.0 238 7.01 242 358 8.72 290 338 398 8 8.0 1953 44.1 202.0 245 5.45 252 366 6.76 301 350 404
下载: 导出CSV
-
[1] 李砚召, 王肖钧, 张新乐, 等. 预应力混凝土结构抗爆性能试验研究 [J]. 实验力学, 2005, 20(2): 179–185. DOI: 10.3969/j.issn.1001-4888.2005.02.004.LI Y Z, WANG X J, ZHANG X L, et al. Test study on anti-detonation quality of prestressed concrete structure [J]. Journal of Experimental Mechanics, 2005, 20(2): 179–185. DOI: 10.3969/j.issn.1001-4888.2005.02.004. [2] 李砚召, 郭晓辉, 曹海. 预应力混凝土梁平面装药加载试验研究 [R]//中国国防科学技术报告. 2002. [3] 胡志坚, 张一峰, 俞文生, 等. 近场爆炸时预应力混凝土梁体抗爆分析 [J]. 中国公路学报, 2019, 32(3): 71–80. DOI: 10.19721/j.cnki.1001-7372.2019.03.008.HU Z J, ZHANG Y F, YU W S, et al. Anti-blast resistance analysis of prestressed concrete bridges under close-by blast [J]. China Journal of Highway and Transport, 2019, 32(3): 71–80. DOI: 10.19721/j.cnki.1001-7372.2019.03.008. [4] CHEN W S, HAO H, CHEN S Y. Numerical analysis of prestressed reinforced concrete beam subjected to blast loading [J]. Materials and Design, 2015, 65: 662–674. DOI: 10.1016/j.matdes.2014.09.033. [5] 宁显东. 有粘结预应力混凝土框架结构抗连续倒塌性能分析 [D]. 南宁: 广西大学, 2018: 33–51. DOI: 10.7666/d.Y3434506.NING X D. Analysis on progressive collapse performanceof bonded prestressed concrete frame construction [D]. Nanning: Guangxi University, 2018: 33–51. DOI: 10.7666/d.Y3434506. [6] 张鹏. 预应力混凝土框架结构的抗连续倒塌性能研究 [D]. 西安: 西安工业大学, 2016: 25–56. DOI: 10.7666/d.Y3078301.ZHANG P. Performance research of progressive collapse on prestressed concrete frame structure [D]. Xi’an: Xi’an Technological University, 2016: 25–56. DOI: 10.7666/d.Y3078301. [7] QIAN K, ZHANG X D, FU F, et al. Progressive collapse-resisting mechanisms of planar prestressed concrete frame [J]. ACI Structural Journal, 2019, 116(4): 77–90. DOI: 10.14359/51715567. [8] MALVAR L J, CRAWFORD J E, WESEVICH J W, et al. A plasticity concrete material model for DYNA3D [J]. International Journal of Impact Engineering, 1997, 19(9/10): 847–873. DOI: 10.1016/S0734-743X(97)00023-7. [9] LI J, HAO H. Influence of brittle shear damage on accuracy of the two-step method in prediction of structural response to blast loads [J]. International Journal of Impact Engineering, 2013, 54: 217–231. DOI: 10.1016/j.ijimpeng.2012.11.008. [10] WOODSON S C, BAYLOT J T. Quarter-scale building/column experiments [C]// Elgaaly M. Advanced Technology in Structural Engineering. Philadelphia, Pennsylvania, USA: American Society of Civil Engineers, 2000: 1–8. DOI: 10.1061/40492(2000)99. [11] BAYLOT J T, BEVINS T L. Effect of responding and failing structural components on the airblast pressures and loads on and inside of the structure [J]. Computers & Structures, 2007, 85(11): 891–910. DOI: 10.1016/j.compstruc.2007.01.001. [12] WOODSON S C, BAYLOT J T. Structural collapse: quarter-scale model experiments: SL-99-8[R]. US Army Corps of Engineers Engineer Research and Development Center, 1999. [13] US Army Corps of Engineers. Structures to resist the effects of accidental explosions: TM5-1300 [S]. Washington DC: Department of the Army, 1990. [14] SHI Y C, HAO H, LI Z X. Numerical derivation of pressure-impulse diagrams for prediction of RC column damage to blast loads [J]. International Journal of Impact Engineering, 2008, 35(11): 1213–1227. DOI: 10.1016/j.ijimpeng.2007.09.001. [15] 中华人民共和国住房和城乡建设部. 混凝土结构设计规范: GB 50010—2010 [S]. 北京: 中国建筑工业出版社, 2011.Ministry of Housing and Urban-Rural Development of the People’s Republic of China. Code for design of concrete structures: GB 50010—2010 [S]. Beijing: China Architecture & Building Press, 2011. [16] US Department of Defense. Structures to resist the effects of accidental explosions: UFC 3-340-02 [R]. Washington DC: The US Department of Army, 2008. [17] United States General Services Administration. Progressive collapse analysis and design guidelines for new federal office buildings and major modernization projects: GSA2003 [S]. Washington, DC: General Services Administration, 2003. [18] 中华人民共和国住房和城乡建设部, 中华人民共和国国家质量监督检验检疫总局. 建筑抗震设计规范: GB 50011—2010 [S]. 北京: 中国建筑工业出版社, 2010.Ministry of Housing and Urban-Rural Development of the People’s Republic of China, General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China. Code for seismic design of buildings: GB 50011—2010 [S]. Beijing: China Architecture & Building Press, 2010. -
下载:
计量
- 文章访问数: 260
- HTML全文浏览量: 145
- PDF下载量: 88
- 被引次数: 0