高温后高强高性能混凝土双轴压力学性能

利用大型静动真三轴试验机,进行了常温20${^\circ}$C以及200${^\circ}$C$\sim $ 600${^\circ}$C\,6个温度等级高温后高强高性能混凝土在7种应力比双轴压应力状态下的强度与变形试验.测得了双轴主压方向的静态强度、峰值应变与应力应变曲线,剖析了温度和应力比对单、双轴压强度与峰值应变发展趋势的影响规律性以及试件破坏形态. 试验结果表明:随温度的升高,高强高性能混凝土的单轴压减摩强度并不一定降低;双轴压强度相对于单轴压强度的提高倍数取决于应力比、不同温度等级后的高强高性能混凝土``脆硬性'. 提出了带有温度和应力比参数的Kupfer-Gerstle破坏准则公式.      

泡沫铝材料动态本构参数的实验确定

基于泡沫材料的动态刚性-线性硬化塑性-刚性卸载(D-R-LHP-R)模型,结合连续性方程,动量守恒方程及刚体的运动方程,得到了激波在泡沫材料中的量纲一消失位置X_s/L₀和动态屈服应力Y_i、激波波速c_p、冲击初始应变ε_i之间的如下关系式: $\frac{X_{\mathrm{s}}}{L_{0}}=\exp \left(-\frac{\rho_{0} c_{\mathrm{p}} v_{\mathrm{i}}}{Y}\right)=\exp \left(1-\frac{\sigma_{\mathrm{i}}}{Y}\right)=\exp \left(-\frac{\rho_{0} c_{\mathrm{p}}^{2} \varepsilon_{\mathrm{i}}}{Y}\right)$ 采用Taylor-Hopkinson装置进行实验,当直接测得泡沫铝试样密度ρ₀、边界初始应力σ_i、初始打击速度v_i、泡沫铝杆原长L₀及激波在泡沫铝杆中消失长度X_s后,利用方程式(a)可反演求得D-R-LHP-R模型下的泡沫铝动态应力应变曲线。最后通过与泡沫铝准静态实验数据对比,表明该泡沫铝是应变率敏感性材料。     

利用维氏和玻氏压头表征半导体材料断裂韧性

压痕法是测量材料断裂韧性 ($K_{\rm IC})$ 的常用方法之一, 如何根据不同的材料、不同的压头选择适合的公式, 是当前面临的一大问题. 因此,在不同载荷下对单晶硅 (111) 和碳化硅 (4H-SiC, 0001面) 这两种半导体材料进行了维氏微米硬度和玻氏纳米压痕实验, 对实验产生的裂纹长度$c$进行了统计分析, 并采用13个压痕公式计算材料的$K_{\rm IC}$, 开展了微米划痕实验, 验证压痕法评估半导体材料$K_{\rm IC}$的适用性. 研究结果表明: 为了消除维氏压痕实验产生的$c$的固有离散性, 需要多次测量取平均值; 裂纹长度与压痕尺寸的比值随压痕载荷的增大而增大; 材料的裂纹类型与载荷相关且低载荷下表现为巴氏裂纹, 高载荷下表现为中位裂纹; 与微米划痕实验得到的单晶硅和碳化硅材料的$K_{\rm IC}$平均值 (分别为0.96 MPa,$\cdot$,$\sqrt{\rm m}$和2.89 MPa,$\cdot$,$\sqrt{\rm m}$) 相比, 在同一压头下无法从13个公式中获得同时适用于单晶硅和碳化硅材料的压痕公式,但在同一材料下可以获得同时适用于维氏和玻氏压头的$K_{\rm IC}$计算公式; 基于中位裂纹系统发展而来的压痕公式更适合用于评估半导体材料的$K_{\rm IC}$, 且维氏压头下的$K_{\rm IC}$与玻氏压头下$K_{\rm IC}$的关系不是理论上的1.073倍, 应为1.13$\pm 压痕法是测量材料断裂韧性(K_(IC))的常用方法之一,如何根据不同的材料、不同的压头选择适合的公式,是当前面临的一大问题.因此,在不同载荷下对单晶硅(111)和碳化硅(4H-Si C, 0001面)这两种半导体材料进行了维氏微米硬度和玻氏纳米压痕实验,对实验产生的裂纹长度c进行了统计分析,并采用13个压痕公式计算材料的K_(IC),开展了微米划痕实验,验证压痕法评估半导体材料K_(IC)的适用性.研究结果表明:为了消除维氏压痕实验产生的c的固有离散性,需要多次测量取平均值;裂纹长度与压痕尺寸的比值随压痕载荷的增大而增大;材料的裂纹类型与载荷相关且低载荷下表现为巴氏裂纹,高载荷下表现为中位裂纹;与微米划痕实验得到的单晶硅和碳化硅材料的K_(IC)平均值(分别为0.96 MPa·m~(1/2)和2.89 MPa·m~(1/2))相比,在同一压头下无法从13个公式中获得同时适用于单晶硅和碳化硅材料的压痕公式,但在同一材料下可以获得同时适用于维氏和玻氏压头的K_(IC)计算公式;基于中位裂纹系统发展而来的压痕公式更适合用于评估半导体材料的K_(IC),且维氏压头下的K_(IC)与玻氏压头下K_(IC)的关系不是理论上的1.073倍,应为1.13±0.01.     

混凝土双K断裂参数计算的半解析有限元法

混凝土裂缝扩展的双$K$断裂准则，用于描述混凝土结构裂缝的起裂、稳定扩 展和失稳断裂. 其相应的双$K$断裂参数(起裂断裂韧度$K_{\rm IC}^{\rm ini} $和失 稳断裂韧度$K_{\rm IC}^{\rm un}$)一般通过简便的试验和基于虚拟裂缝扩展粘 聚力的解析方法确定. 利用平面扇形域哈 密顿体系的方程，通过分离变量法及共轭辛本征函数向量展开法，以解析的方法推导出基于混 凝土虚拟裂缝扩展线性粘聚力模型的平面裂缝解析元列式. 将该解析元与有限元相结合，构成 半解析的有限元法，可求解任意结构几何形状的混凝土平面裂缝双$K$断裂参数的计算问题. 数值计算结果表明半解析有限元法对该类问题的求解是十分有效的.     

等速上仰翼型动态失速现象研究

翼型大迎角绕流的静态失速将造成升力突降和气动性能急剧恶化，但利用非定常运动所产生 的动态失速效应，可以大大地延缓气流分离和失速现象的发生. 采用Rogers发 展的双时间步Roe格式，求解拟压缩性修正不可压N-S方程. 数值模拟了低雷诺数 ($Re=4.8 \times 10^{4}$)条件下NACA0015翼型作等速上仰($\alpha =0^{\circ} \sim 60^{\circ}$)的动态失速过程，同Walker的试验结果比 较，验证了计算结果的正确性. 研究了该过程中主涡、二次涡和三次涡的发展，升 力系数随攻角变化，以及不同上仰速度对动态失速效应所造成的影响.     

内埋武器高速风洞弹射投放模型试验关键技术研究

针对新一代战斗机超声速内埋武器弹射投放分离安全性问题,采用高速风洞投放实验技术研究内埋武器从开式武器舱弹射投放分离动态运动过程,风洞投放模型试验过程中采用除垂直加速度不足外,其余全部运动严格相似的轻模型相似设计方法,并针对轻模型法垂直加速度不足所导致的投放垂直位移偏离实物位移问题,采用一种简单易行的公式修正法进行补偿,试验给出了不同初始弹射投放分离条件下,内埋武器从载机投放分离后运动轨迹与姿态角随分离时间的变化规律,试验马赫数$Ma = 1.5$.研究结果表明:初始投放分离角速度对内埋武器投放分离后的运动轨迹及姿态角有较大的影响,当初始投放分离角速度$\omega _{z0}^s = 0^\circ/{\rm s}$时,内埋导弹出舱后先向下运动远离载机的流场干扰区,之后逐渐向载机方向抬升靠近并最终碰撞载机,高速风洞投放试验结果是不安全的,但经过公式修正后投放试验结果比较乐观,垂直方向运动仍然一直下降远离载机,这说明采用高速风洞投放试验得出的导弹不安全投放分离对真实载机来说不一定会出现,高速风洞投放试验结果比较保守. 当初始投放分离角速度$\omega _{z0}^s= 15^\circ/{\rm s}$和$\omega _{z0}^s = 30^\circ/{\rm s}$时,内埋导弹投放分离后运动趋势几乎一致, 均没出现向载机靠近的现象,内埋导弹具有一定的初始投放分离角速度有利于内埋武器的安全分离.      

基于微米力学实验的页岩I型断裂韧度表征

页岩断裂韧度($K_{IC})$是页岩气储层水力压裂设计的基础参数之一,由于组成的非均质性,常规宏观力学测量方法存在制样困难、力学解释参数不连续、精度偏低等问题. 如何及时获取页岩的断裂特性,确保安全高效的工程施工,是当前面临的一大问题. 因此,提出了基于微米力学实验的页岩Ⅰ型断裂韧度分析方法,可用于页岩微裂纹起裂、发育直至形成宏观裂纹的机理研究,进行页岩宏观Ⅰ型断裂韧度预测. 基于页岩多尺度组成分析,开展了维氏压头和玻氏压头的页岩微米力学实验,分析了页岩残余压痕与压头间的相似关系、有效测试载荷以及压头参数的优化与选择. 分析了不同压入载荷下的页岩细观断裂韧度分布特征,开展了宏观巴西圆盘实验,验证页岩微米力学测试方法的适用性. 研究结果表明,在有效载荷范围内的页岩细观Ⅰ型断裂韧度波动性较小,当压入载荷过大时,由于岩样压痕区域出现局部剥落导致断裂韧度测量值偏小. 与宏观实验的比对分析显示,微米力学实验的$K_{IC}$平均值为0.86 MPa$\cdot \sqrt{m}$,直槽切缝巴西圆盘实验得到的$K_{IC}$平均值为0.92 MPa$\cdot \sqrt{m}$,两类方法的统计平均值较为接近,页岩局部组成的非均质性使得微米力学测量结果较宏观测试更为分散. 研究结果可用于页岩宏观Ⅰ型断裂韧度预测,为有效解决页岩气储层水力压裂参数评价提供新的思路和方法.      

不同壁面取向下超疏水平面直轨道上的气泡滑移

利用特定几何分布的超疏水表面实现气泡定向输运在矿物浮选和生物孵化等领域具有广阔的应用前景, 对平面直线超疏水轨道而言, 其壁面取向是相关工程结构的关键参数, 但超疏水壁面取向对倾斜壁面气泡滑移的影响尚不明确. 本文采用高速阴影成像系统研究了不同壁面取向($-90^\circ\leqslant \beta \leqslant 90^\circ$)及轨道倾角($45^\circ\leqslant \alpha \leqslant 75^\circ$)下, 气泡($D_{eq}=2.4$ mm, $Re=500$ $\sim$ 700, $We=7$ $\sim$ 13)在轨道宽度为2 mm的超疏水直线轨道上的运动特性. 气泡在轨道上的滑移近似为匀速, 形状为具有多脊的半子弹型. 根据气液界面波动程度的不同, 滑移气泡可分为波动型和稳定型, 稳定型气泡只在较小倾角且较大方位角时出现($45^\circ\leqslant \alpha < 70^\circ$, $| \beta | \geqslant 45^\circ$). 根据倾角不同, 滑移速度关于$\beta $有2种变化规律: 当$\alpha \leqslant 65^\circ$, 气泡滑移速度近似为关于$\beta =0^\circ$ 的单峰分布($\beta =0^\circ$时, 气泡滑移速度最大); 当$\alpha \geqslant 70^\circ$, 气泡滑移速度在不同的方位角下基本保持稳定. 气泡的最大滑移速度可达0.66 m/s ($\beta =0^\circ$, $\alpha =70^\circ$), 远大于相同尺度的自由上升气泡($\approx0.25$ m/s), 这主要是壁面浸润性分布和惯性力的耦合效应所致. 轨道取向(方位角$\beta )$及轨道倾角($\alpha )$通过改变气泡沿轨道方向的驱动力和气泡迎风面积影响气泡的滑移速度和气液界面稳定性.      

开裂孔隙材料渗透率的细观力学模型研究

采用细观力学方法对含随机裂纹网络的孔隙材料渗透性进行研究.开裂孔隙材料渗透性的影响因素包括裂纹网络的密度、连通度、裂纹的开度以及孔隙材料基体渗透性.对于不连通的裂纹网络,该文采用已有的相互作用直推法(interaction direct derivative,IDD)的理论框架,引入裂纹的密度$\rho$和裂纹开度比$b$,提出了裂纹夹杂$\!$-$\!$-$\!$基体两相复合材料渗透率的IDD理论解.对于部分连通裂纹网络,考虑局部裂纹团内部各个裂纹对有效渗透率的相互放大作用,引入裂纹网络的连通度$f$,定义与连通度相关的水平裂纹密度$\rho^{h}$,按照增量法将表征连通特征的水平裂纹嵌入有效基体中,以此方式来考虑裂纹夹杂间的相互搭接,提出了考虑裂纹连通特征的扩展IDD理论解,分别考虑了基体材料渗透率$K_{m}$、裂纹密度$\rho $、裂纹开度比$b$以及与连通度$f$相关的$\rho ^{\rm h}$.最后通过对有限区域内含随机裂纹网络孔隙材料渗透过程的有限元模拟分别验证了不连通和部分连通裂纹网络扩展IDD模型的适用性:(1)当裂纹不连通时,由于基体对流体渗透的阻隔作用,裂纹的开度对有效渗透率影响不大;(2)当裂纹部分连通时,裂纹密度分别小于1.1(无关联裂纹网络,分形维数为2.0)、1.2(关联裂纹网络,分形维数为1.75)时,扩展IDD模型能够很好地估计开裂孔隙材料的有效渗透率,但是随着裂纹进一步扩展,最大裂纹团主导作用凸显,扩展IDD模型不再适用.      

微型热驱动回路的脉动及位差对热性能的影响

以甲醇为工质，采用高速数据采集系统测定了微型热驱动回路在不同运行参数下的压力 及温度脉动，其脉动周期及脉动幅度随蒸发段热流密度的增加而减小. 实验发现，在蒸发段 热流密度较低的情况下，蒸气管中是泡状流或弹状流交替存在，而在蒸发段热流密度较高时， 蒸气管中为环状流. 就位差对热性能的影响进行了详细的实验研究，并在冷凝器空气自 然对流和强迫对流情况下，以加热块温度90${^\circ}$C为上限，得出微通道蒸发器和冷凝 器在不同位差下的最大蒸发段热流密度. 通过对实验现象的观察及分析，以期开发出适用于 未来电子产品高功率需求的微型化电子冷却器.     

完整岩石拉压应力状态下的张裂破坏准则

一点的应力状态可由3个主应力$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$来表示, 当规定主应力以压为正时, 沿最大主应力$\sigma_{1}$方向将产生收缩变形, 若中间主应力$\sigma_{2}$和最小主应力$\sigma_{3}$都远小于$\sigma_{1}$, 则沿$\sigma_{2}$和$\sigma_{3}$方向会产生横向扩张变形, 当横向扩张变形达到一定极限时, 将会在平行于$\sigma _{1}$的方向产生张裂破坏. 如何建立这种张裂破坏的强度准则目前尚缺乏研究, 最大拉应变理论(第二强度理论)有时被用来解释张裂破坏, 但最大拉应变理论难以应用于三向受力状态. 本文分别用$\varepsilon_{1}$, $\varepsilon_{2}$表示最大张应变和次大张应变, 则最大拉应变理论认为当$\varepsilon_{1}$达到单向拉伸屈服应变时, 材料将产生破坏. 而本文将根据$\varepsilon_{1}+\varepsilon_{2}$之和达到极限值$\varepsilon_u$来建立张裂破坏准则. 可以证明$\varepsilon_{1} +\varepsilon_{2}$所表示的是$\sigma_{1}$主平面的面积增长率. 当$\sigma_{3}<\sigma_{2} \ll \sigma_{1}$时, 大部分岩石都具有脆性破坏的特点, 所以可将破坏前的岩石视为满足广义胡克定律的线弹性材料, 这样用$\varepsilon_{1}$, $\varepsilon_{2}$表示的强度准则可通过$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$来表示. 在这个过程中还可考虑岩石在拉伸和压缩时具有不同弹性参数和强度的特点, 并可通过单向拉伸和单向压缩的破坏状态来确定$\varepsilon_u$. 不管$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$是压应力, 还是拉应力, 或者$\sigma_{1}$, $\sigma_{2}$, $\sigma_{3}$中有拉有压的情形, 基于$\varepsilon_{1} +\varepsilon_{2} =\varepsilon_u$都可建立相应的强度准则. 所建立的准则可以反映中间应力$\sigma_{2}$对强度的影响规律, 通过建立的强度准则还可以证明: 静水拉力能引起屈服, 而静水压力不能产生屈服; 压缩破坏能使塑性体积增大, 其结果比Mohr-Coulomb准则更能反映实际情形. 并通过拉压应力状态下的试验数据验证了所建立的强度准则, 所得理论计算结果和已有的试验数据吻合得很好. 通过提出的强度准则和圆盘劈裂的试验结果, 可获得更为可靠的岩石单轴抗拉强度.      

ANALYSIS OF MECHANICAL PROPERTIES OF CEMENT CLINKER UNDER DIFFERENT TEMPERATURES BASED ON HERTZIAN INDENTATION METHOD1)

本文通过高温球压法得到各种脆性和准脆性材料的表面压痕应力与应变之间的关系曲线。通过压痕应力-应变曲线的分析既可比较方便地确定出材料的压痕弹性模量、剪切模量和布氏硬度,又可比较不同温度下水泥熟料的变形性能。在不同温度(25 ${^\circ}\!$C$\sim $1400 ${^\circ}\!$C)处理下,球压应力松弛试验载荷松弛,在载荷峰值为100 N时,随着温度的升高,水泥熟料载荷松弛更明显。随着温度从500 ${^\circ}\!$C升高到1400 ${^\circ}\!$C,载荷松弛非常明显,尤其温度高于1200 ${^\circ}\!$C,水泥熟料样品内部的硅酸三钙(Ca$_3$SiO$_5$, 简称C3S)分解以及有部分液相的出现引起的应力松弛现象最为明显,在1275 ${^\circ}\!$C时熟料基本上已经软化,载荷急速松弛,所以认为1275 ${^\circ}\!$C为熟料的脆延转化温度。通过水泥熟料高温球压松弛试验可以确定水泥熟料在二次加热过程中的脆-延转化温度,测定熟料弹性模量和抗压强度急剧变化的温度范围。研究水泥熟料在不同温度下的力学行为和力学特性,探索提高粉磨效率的新途径,实现高温下的低能耗粉碎。     

The Dielectric Permittivity of Crystals in the Reduced Hartree–Fock Approximation

In a recent article (Cancès et al. in Commun Math Phys 281:129–177, 2008), we have rigorously derived, by means of bulk limit arguments, a new variational model to describe the electronic ground state of insulating or semiconducting crystals in the presence of local defects. In this so-called reduced Hartree–Fock model, the ground state electronic density matrix is decomposed as ${\gamma = \gamma^0_{\rm per} + Q_{\nu,\varepsilon_{\rm F}}}$ , where ${\gamma^0_{\rm per}}$ is the ground state density matrix of the host crystal and ${Q_{\nu,\varepsilon_{\rm F}}}$ the modification of the electronic density matrix generated by a modification ν of the nuclear charge of the host crystal, the Fermi level ε _F being kept fixed. The purpose of the present article is twofold. First, we study in more detail the mathematical properties of the density matrix ${Q_{\nu,\varepsilon_{\rm F}}}$ (which is known to be a self-adjoint Hilbert–Schmidt operator on ${L^2(\mathbb{R}^3)}$ ). We show in particular that if ${\int_{\mathbb{R}^3}\,\nu \neq 0, Q_{\nu,\varepsilon_{\rm F}}}$ is not trace-class. Moreover, the associated density of charge is not in ${L^1(\mathbb{R}^3)}$ if the crystal exhibits anisotropic dielectric properties. These results are obtained by analyzing, for a small defect ν, the linear and nonlinear terms of the resolvent expansion of ${Q_{\nu,\varepsilon_{\rm F}}}$ . Second, we show that, after an appropriate rescaling, the potential generated by the microscopic total charge (nuclear plus electronic contributions) of the crystal in the presence of the defect converges to a homogenized electrostatic potential solution to a Poisson equation involving the macroscopic dielectric permittivity of the crystal. This provides an alternative (and rigorous) derivation of the Adler–Wiser formula.     

Lower Semicontinuity for Integral Functionals in the Space of Functions of Bounded Deformation Via Rigidity and Young Measures

We establish a general weak* lower semicontinuity result in the space BD(Ω) of functions of bounded deformation for functionals of the form

A Nonlinear Elliptic PDE with Two Sobolev–Hardy Critical Exponents

In this paper, we consider the following PDE involving two Sobolev–Hardy critical exponents,

Incompressible 3D Navier–Stokes Equations as a Limit of a Nonlinear Parabolic System

In this paper, we consider the Cauchy problem for a nonlinear parabolic system ${u^\epsilon_t - \Delta u^\epsilon + u^\epsilon \cdot \nabla u^\epsilon + \frac{1}{2}u^\epsilon\, {\rm div}\, u^\epsilon - \frac{1}{\epsilon}\nabla\, {\rm div}\, u^\epsilon = 0}$ in ${\mathbb {R}^3 \times (0,\infty)}$ with initial data in Lebesgue spaces ${L^2(\mathbb {R}^3)}$ or ${L^3(\mathbb {R}^3)}$ . We analyze the convergence of its solutions to a solution of the incompressible Navier?CStokes system as ${\epsilon \to 0}$ .