复杂加载条件下的砂土本构模型

试验表明,饱和砂土的应力应变关系具有显著的密度以及压力依存性,上述两点构成了描述砂土静力加载下变形特性无法忽视的因素. 此外,在循环加载等复杂加载作用下,砂土还会表现出明显的应力诱导各向异性以及相变转换特性. 基于在e--p空间中存在唯一的临界状态线这一基本假定,通过在e--p空间中引入当前状态点与临界状态线的距离R来作为反映密度与压力依存特性的状态参量, 将变相应力比以及峰值应力比表达为状态参量的指数函数,将上述应力比参量引入到统一硬化参量中可准确地反映初始状态下围压、密度 对于单调加载下应力应变关系的影响规律,能描述砂土剪缩、剪胀,应变软化、硬化等特性. 采用非相关联流动法则,p--q空间中采用水滴型屈服面,塑性势面为椭圆面,松砂在单调加载下的静态液化现象也可描述. 为反映循环加载下塑性体积应变的累积特性以及塑形偏应变的滞回特性,在循环加载下将状态参量R表达为应力比参量,并在硬化参数中引入描述应力诱导各向异性特性的旋转硬化部分,所提模型可有效地描述循环加载下剪切模量的衰减特性、刚度衰化性质、强度减小特性,在不排水约束作用下,则会产生往返活动性现象. 通过一系列的模型模拟与试验结果对比,验证了本构模型的有效性及适用性.      

应力诱发界面迁移下晶内孔洞的演化

随着微电子技术的迅猛发展, 集成电路中内连导线的失效问题引起广泛关注. 内连导线内部孔洞萌生、长大、漂移和失稳变形成狭长裂纹, 从而导致电路的开路失效. 这是内连导线失效的常见形式. 而界面迁移是导致微结构形态演化的主要机制之一. 本文基于界面迁移下微结构演化的经典理论和弱解描述, 建立了应力诱发界面迁移下微结构演化的有限单元法, 并验证了算法的可靠性. 对铜内连导线中晶内孔洞的演化进行了数值模拟, 详细分析了应力、线宽及形态比对晶内孔洞演化的影响. 研究结果表明, 椭圆形晶内孔洞存在生长和收缩两种演化分叉趋势. 通过大量数值分析得到了晶内孔洞演化的临界应力{\stackrel{\text{?}}{\sigma }}_c 、临界线宽{\stackrel{\text{?}}{h}}_c 和临界形态比{\beta }_c . 当\stackrel{\text{?}}{h}\mathrm{?}{\stackrel{\text{?}}{h}}_c, \stackrel{\text{?}}{\sigma }或\stackrel{\text{?}}{h} 时, 晶内孔洞会沿长轴长大; 反之, 晶内孔洞会收缩甚至愈合. 此外, 应力\beta 越大、线宽\stackrel{\text{?}}{\sigma } 越小或形态比\stackrel{\text{?}}{h} 越大, 晶内孔洞越易发生长大, 且孔洞面积增大速度越快; \beta 越小、\stackrel{\text{?}}{\sigma } 越大或\stackrel{\text{?}}{h} 越小, 晶内孔洞越易发生收缩, 且孔洞面积减小速度越快.      

功能梯度材料动力学问题的POD模型降阶分析

为了快速分析非均质材料结构在复杂载荷作用下的动态响应, 提出一种模型降阶方法, 只需计算结构在简单均质材料情况下的动力学问题, 进而用其计算结果对非均质材料结构进行分析. 首先, 采用结构内部任意一点处的材料参数值作为整个结构的材料参数, 利用有限元分析软件计算该均质材料结构在动态载荷作用下的位移场建立数据库, 该数据库包含计算模型各个节点(自由度为L)在某时间段内L个时刻的位移; 其次, 对数据库中的信息按照时间离散的特定方式组集成瞬像矩阵, 并利用特征正交分解方法对其进行分解, 得到该模型的H个特征正交基底, 选取其中能反应模型主要特征的个(其中\mathrm{~})作为一组最优基底, 通过这组基底建立模型的低阶离散控制方程; 最后, 求解低阶离散微分方程组, 得到功能梯度材料结构在复杂载荷作用下的位移场. 文中分别给出二维和三维算例, 比较了降阶模型和全阶模型计算结果, 验证了该方法的有效性, 并且计算效率能提高1{\mathrm{}}^{1)}2个数量级.      

一种考虑液相惯性的一维准静力固结模型

经典Terzaghi一维固结理论不考虑孔隙流体惯性影响,且该理论在不同时期模型推导和表述结果差别较大,导致当前仍存在诸多困惑甚至认识混乱的现象. 在笔者前期研究大变形动力固结理论框架内,忽略固相惯性而重点考虑液相惯性影响,经过合理简化建立反映孔隙流体惯性的一维小变形固结波动模型. 该固结波模型具有频散和耗散特性. 采用分离变量法,可得到单面排水和瞬时加载条件下无量纲形式固结波解析解答. 算例分析结果表明:固结波发展规律受无量纲数D_{\mathrm{c}}变化影响而呈现不同性态;D_{\mathrm{c}}数值较大时固结波响应会出现阶跃和正负波动现象;当D_{\mathrm{c}}值较小时,可能出现Mandel-Cryer效应等特殊现象. 通过对早期和后期Terzaghi固结模型的分析和对比,初步探明Terzaghi固结理论模型内部的矛盾性,在普通土体坐标和固相体积坐标两种不同解读条件下,早期Terzaghi (1923,1925)模型可以分别诠释为具有小变形和大变形属性的不同固结模型. 在经典一维固结理论模型的不同诠释背景下,固结波模型也可以据此作出相应拓展和表述. 固结波理论揭示缩尺固结试验中土体物理力学参数与固结波响应两种因素之间存在一种不确定性矛盾,据此建议微观土力学研究重视尺度效应. 固结波模型的意义还在于,可为Terzaghi经典固结模型理论精度分析提供新的依据.      

高雷诺数下多柱绕流特性研究

采用改进的延迟分离涡方法数值模拟了高雷诺数下的柱体绕流,包括单圆柱绕流、单方柱绕流、串列双圆柱绕流和串列双方柱绕流,研究了不同雷诺数下圆柱绕流与方柱绕流的水动力特性. 计算结果与实验数据及其他文献的数值计算结果吻合良好,研究表明,单方柱绕流在范围内未出现类似于单圆柱绕流的阻力危机现象,其平均阻力系数\overline{C_d}、升力系数均方根{C\mathrm{\text{'}}}_l及斯特劳哈尔数\mathit{St}维持在一定范围内波动. 串列双圆柱绕流与串列双方柱绕流中,均选取L/D=2.0, 2.5, 3.0, 3.5和4.0这五中间距比进行计算. 串列双圆柱绕流中,当\mathit{Re}=2.2\times 10^4时,在内存在一临界间距比(L_c/D)使得L_c/D前后上下游圆柱的升阻力系数发生跳跃性变化,且当时,下游圆柱的阻力系数为负数. 而当\mathit{Re}=3.0\times 10^6时,则不存在临界间距比,且下游圆柱的阻力系数始终为正数. 串列双方柱绕流在\mathit{Re}=1.6\times 10^4和\mathit{Re}=1.0\times 10^6两种工况下的临界间距比分别处于和区间内,且当时,两个雷诺数下的下游方柱阻力系数均为负数.      

扩展腔对方波型微混合器混合性能的影响研究

微混合器凭借节约试剂、混合强度高和易于集成等优点,在材料合成、医药制备和生化检测等领域中具有广泛的应用. 为了进一步提高混合性能,保证混合过程的安全性及生化反应结果的准确性,设计了一种带扩展腔的新型方波型微混合器. 在综合考虑混合强度和压降的前提下,通过实验研究和数值模拟分析了窄缝宽度、窄缝长度和扩展腔高度对微混合器混合性能的影响并得到了不同雷诺数\mathit{Re}条件下的最优结构参数. 与方波型微混合器的混合性能进行比较,发现\mathit{Re}=20时,带扩展腔的方波型微混合器的混合强度更高,其中\mathit{Re}时两者混合强度相差最多,可达12%. 在相同\mathit{Re}下,带扩展腔的方波型微混合器的压降要低于方波型微混合器. 对带扩展腔的方波型微混合器进行内部流场分析,发现扩展腔结构能在流体层流状态的基础上引入涡流,使通道中流体的流动状态发生改变、对流增强,进而混合性能提高.      

错列角度对双圆柱涡激振动影响的数值模拟研究

为研究错列角度\alpha对双圆柱涡激振动问题的影响,采用自主研发的基于CIP (constrained interpolation profile)方法的数值模型,对雷诺数\mathit{Re}=100、错列角度\alpha =0^{°}\mathrm{~}90^{°} (间隔15 #x00B0;)的等直径双圆柱涡激振动问题进行数值模拟. 模型在笛卡尔网格系统下建立,采用具有三阶精度的 CIP 方法求解 N-S (Navier--Stokes)方程,采用浸入边界法处理流--固耦合问题,避免了任意拉欧方法下的网格畸变和重叠动网格技术中的大量信息交换问题,保证了模型的计算效率. 重点分析不同错列角度\alpha上下游圆柱的升阻力系数、位移响应、涡脱频率和尾涡模态等. 结果表明:折合速度U_r=2.0\mathrm{~}3.0时,上下游圆柱升阻力随错列角度的增大基本呈单调增大的趋势;U_{\mathrm{r}}=5.0\mathrm{~}8.0时,随错列角度的增大,上下游圆柱阻力变化较小,升力呈“上凸”趋势,在\alpha =15^{°}\mathrm{~}30^{°}取得最大值;U_r=10.0\mathrm{~}13.0时,随错列角度的增大,上下游圆柱阻力变化较小,升力呈“下凹”趋势,在\alpha =30^{°}\mathrm{~}45^{°}取得最小值,且柱体横流向振幅和升力没有明显的对应关系. 最后,结合尾涡模态对以上规律的成因进行分析. 研究结果可为相关海洋工程设计提供参考.      

PMMA膨胀环动态拉伸碎裂实验研究

在强动载作用下, 脆性材料的碎裂问题是一个重要的研究课题, 而脆性材料在冲击拉伸载荷下的力学行为的实验研究相对较匮乏. 提出了一种动态拉伸断(碎)裂的液压膨胀环实验技术, 可用于准脆性/脆性材料的动态拉伸. 利用该技术对有机玻璃(PMMA)圆环试件进行了不同膨胀速度下的动态碎裂实验研究. 从回收碎片的断口形貌和碎片内部残余裂纹观察可知试件的破碎由环向拉伸应力造成, 碎片断口处发出的稀疏波会将周围的拉伸应力卸载, 从而抑制其他裂纹的进一步发展. 利用超高速相机记录了试件的膨胀碎裂过程, 利用DISAR激光速度干涉仪获得了试件外表面粒子的径向膨胀速度历史, 通过试件上的应变片获得了试件的应变历史和断裂应变. 实验结果表明: 在拉伸应变率150\mathrm{~}500s^{-1}范围, 材料的动态断裂应变低于准静态加载下的断裂应变, 体现出“动脆”现象; 随着加载应变率的提高, PMMA 材料的碎片尺寸减小; 无量纲化的PMMA圆环的平均碎片尺寸介于韧性碎裂模型和脆性碎裂模型的预测数值之间, 反映出材料的准脆性特性.      

Singularly perturbed nonlinear second order elliptic equation

In this paper, we study singular perturbation problems of some semi-linear second order elliptic equations with nonlinear boundary value conditions: where ε is a small positive parameter and u/ l is a directional derivative, which lies on an oblique vector (x,ε). We have given a construction of the asymptotic solutions and proof of their asymptotic correctness, which is based on the principle of contraction mapping.     

碱类型对地聚合物混凝土应变率效应影响的对比研究

采用经波形整形技术改进后的?100 mm SHPB实验装置开展了强度等级均为C30的NS激发矿渣粉煤灰基地聚合物混凝土(NaOH and sodium silicate-activated slag and fly ash based geopolymer concrete,NSSFGC)和NN激发矿渣粉煤灰基地聚合物混凝土(NaOH and Na₂CO₃-activated slag and fly ash based geopolymer concrete,NNSFGC)的动态压缩实验,并对比分析了在冲击荷载作用下的应变率效应。结果表明:NSSFGC和NNSFGC的峰值应力均随应变率的增加而越大。这表明,GC为应变率敏感材料;在动态压缩状态下,NSSFGC和NNSFGC的应变率敏感阈值分别为51.82、28.89 s-1;GC的应变率敏感性强于普通混凝土;NN型激发剂更有利于发挥GC的整体强度特性。由此可知,NNSFGC的应变率敏感性明显强于NSSFGC。     

大应变率范围内AM80镁合金的变形行为及组织演变

采用INSTRON准静态压缩试验机和分离式霍普金森压杆装置，研究固溶态AM80镁合金在室温准静态和冲击载荷下的变形行为及组织演变。准静态载荷下，流变应力随应变率（3×10-5~4×10-1 s-1）的升高逐渐降低，表现为负应变率敏感性；冲击载荷下，流变应力随应变率（7.00×102~5.20×103 s-1）的升高而升高，呈现出明显的正应变率敏感性。冲击载荷下AM80镁合金的变形机制以基面滑移和孪生为主，大量细小致密的形变孪生以及适量非基面滑移的启动是AM80镁合金在冲击载荷下流变应力明显高于准静态载荷的重要原因。此外，随应变率的升高，AM80镁合金变形的均匀性明显增强，当应变速率升至3.65×103 s-1时，冲击变形所引起的局部绝热温升软化大于应变硬化与应变速率硬化的总和，部分晶粒产生了明显的动态回复，使得孪晶密度和变形均匀性反而降低。     

A DISCUSSION ON “REPRESENTING GENERAL SOLUTION OF EQUATIONS IN THEORY OF ELASTICITY BY HARMONIC FUNCTIONS”

But,generally,it is impossible to determine P,Q,R,φ,ψ,A,B from(4)unless the domain ofelasticity is convex in x-direction,y-direction,and z-direction.This problem was resolved byEubanks and Sternberg in 1956.And,they gave examples which illustrat that solutions of Eq.(4)may not exist in nonconvex regions.So the solution stated in[1]is only effective in three-directconvex region,for instance,sphere,cube and so on.And it may not hold true if it is structures whichare commonly used in engineering such as the sheel and the I-shaped beam.It is well known that in 1932 and 1934,Papkovich and Neuber proposed a general solution forEq.(1).It has the form     

Experimental and numerical study of turbulent flow and heat transfer inside hexagonal duct

Flow and heat transfer characteristics in transition and turbulent regions are studied experimentally and numerically in a horizontal smooth regular hexagonal duct under constant wall temperature boundary condition covering a range of Reynolds number from 2.3 × 10³ to 52 × 10³. Two types of k-omega (standard and shear stress transport (SST)) and three types of k-ε (standard, renormalization (RNG), and realizable) turbulence model are employed for transition and turbulent regions, respectively. Both average and fully developed Darcy friction factor and Nusselt number are presented as a function of Reynolds number. It is seen that k-omega SST and k-ε realizable turbulence models gave the best agreement with the experimental data in transition and turbulent regions, respectively. All the experimental results are correlated within an accuracy of ±13 % and ±7 % for Nusselt number and Darcy friction factor, respectively. Results obtained in this study are compared with circular duct results using hydraulic diameter.     

Influence of disodium hydrogen phosphate dodecahydrate on hydrothermal formation of hemihydrate calcium sulfate whiskers

The influence of Na₂HPO₄·12H₂O on the hydrothermal formation of hemihydrate calcium sulfate (CaSO₄·0.5H₂O) whiskers from dihydrate calcium sulfate (CaSO₄·2H₂O) at 135 °C was investigated. Experimental results indicate that the addition of phosphorus accelerates the hydrothermal conversion of CaSO₄·2H₂O to CaSO₄·0.5H₂O via the formation of Ca₃(PO₄)₂ and produces CaSO₄·0.5H₂O whiskers with thinner diameters and shorter lengths. Compared with the blank experiment without Na₂HPO₄·12H₂O, the existence of minor amounts (8.65 × 10⁻⁴–4.36 × 10⁻³ mol/L) of Na₂HPO₄·12H₂O led to a decrease in the diameter of CaSO₄·0.5H₂O whiskers from 1.0–10.0 to 0.5–2.0 μm and lengths from 70–300 to 50–200 μm.     

Heat transfer characteristics of staggered wing-shaped tubes bundle at different angles of attack

An experimental and numerical study has been conducted to clarify heat transfer characteristics and effectiveness of a cross-flow heat exchanger employing staggered wing-shaped tubes at different angels of attack. The water-side Re_w and the air-side Re_a were at 5 × 10² and at from 1.8 × 10³ to 9.7 × 10³, respectively. The tubes arrangements were employed with various angles of attack θ_1,2,3 from 0° to 330° at the considered Re_a range. Correlation of Nu, St, as well as the heat transfer per unit pumping power (ε) in terms of Re_a and design parameters for the studied bundle were presented. The temperature fields around the staggered wing-shaped tubes bundle were predicted by using commercial CFD FLUENT 6.3.26 software package. Results indicated that the heat transfer increased with the angle of attack in the range from 0° to 45°, while the opposite was true for angles of attack from 135° to 180°. The best thermal performance and hence the efficiency η of studied bundle occurred at the lowest Re_a and/or zero angle of attack. Comparisons between the experimental and numerical results of the present study and those, previously, obtained for similar available studies showed good agreements.     

FLOW PERIODICITY AND ACOUSTIC RESONANCE IN PARALLEL TRIANGLE TUBE BUNDLES

Parallel triangle tube arrays with pitch ratios in the range of 1·2–4·2 have been tested at Reynolds numbers up to 9×10⁴to investigate the vorticity shedding and acoustic resonance mechanisms. Three different components of flow periodicities have been observed. The flow periodicity with the highest Strouhal number (S₃) is the weakest component. It is caused by a shear-layer instability in small pitch ratio arrays and at low Reynolds numbers. The S₂component is associated with small-scale vortex shedding at the first row. The third component has the lowest Strouhal number (S₁) and is the strongest. It is generated by large-scale, alternating vortex shedding at deeper rows, and it becomes dominant at all rows at high Reynolds numbers. For tube arrays with pitch ratios less than 3·4, the onset of acoustic resonances could not be related to the natural flow periodicities mentioned above. This behaviour is in contrast with that of normal triangle arrays, but similar to the acoustic behaviour of in-line arrays. Strouhal number charts for the natural flow periodicities and for the onset of acoustic resonances have been developed.     

MoS2基复合薄膜制备及其结构与摩擦学性能研究

利用多靶磁控溅射法制备MoS_2基系列复合薄膜,通过扫描电子显微镜、X-ray衍射、拉曼光谱对薄膜微观形貌及晶体结构进行表征,利用纳米压入表征薄膜硬度及弹性模量,通过微动摩擦磨损试验对比分析了该系列薄膜在大气环境下的润滑性能.研究结果表明:MoS_2薄膜中复合C、Ti元素可有效抑制薄膜生长过程中柱状结构的形成,薄膜更为致密;复合薄膜的摩擦系数及磨损率显著降低.其中,MoS_2+C+Ti三元复合薄膜硬度为7.43 GPa,其弹性模量及弹性恢复量分别为98.1 GPa和61.7%.在大气环境(RH 35%~50%)下,法向载荷10 N时MoS_2+C+Ti薄膜具有最低的磨损率4.67×10–17 m~3/(N·m),该薄膜在不同载荷的微动摩擦试验中均具有最好的承载力.     

Poly(N-isopropylmethacrylamide-co-4-acrylamidobenzo-18-crown-6) microgels with expanded networks for excellent adsorption of lead(II) ions

Lead(II) (Pb²⁺) ions are toxic heavy metal ions that can accumulate in the human body through water and cause severe health problems, including neurotoxicity, nephrotoxicity, hematological toxicity and even genotoxicity effects. To remove Pb²⁺ selectively and effectively from aqueous solutions, we develop a novel type of Pb²⁺-recognizable microgels with excellent adsorption capacity to Pb²⁺, which are fabricated from 4-acrylamidobenzo-18-crown-6 (B18C6Am) and N-isopropylmethacrylamide (NIPMAM) monomers by precipitation copolymerization method. The prepared poly(N-isopropylmethacrylamide-co-4-acrylamidobenzo-18-crown-6) (PNMB) microgels exhibit expanded structures, because the electron-donating methyl groups at α-carbon could descend the polarity of CO in the NIPMAM monomers and thus weaken the polymer segment···segment interactions. The expanded structures of PNMB microgels are beneficial for adsorption of Pb²⁺ due to the low steric hindrance in the polymeric networks. The Pb²⁺ adsorption isotherms of PNMB microgels are consistent with the Langmuir model for monolayer adsorption. The results indicate that the prepared Pb²⁺-recognizable PNMB microgels are highly promising for the selective removal of lead(II) ions from aqueous solutions.     

Dynamical deformation of a flat liquid–liquid interface

The passage of solid spheres through a liquid–liquid interface was experimentally investigated using a high-speed video and PIV (particle image velocimetry) system. Experiments were conducted in a square Plexiglas column of 0.1 m. The Newtonian Emkarox (HV45 50 and 65% wt) aqueous solutions were employed for the dense phase, while different silicone oils of different viscosity ranging from 10 to 100 mPa s were used as light phase. Experimental results quantitatively reveal the effect of the sphere’s size, interfacial tension and viscosity of both phases on the retaining time and the height of the liquid entrained behind the sphere. These data were combined with our previous results concerning the passage of a rising bubble through a liquid–liquid interface in order to propose a general relationship for the interface breakthrough for the wide range of Mo ₁/Mo ₂ ∈ [2 × 10⁻⁵–5 × 10⁴] and Re ₁/Re ₂ ∈ [2 × 10⁻³–5 × 10²].