爆炸荷载作用下方形夹芯板动力学响应与优化设计数值分析

为研究方形蜂窝铝板在爆炸荷载作用下的动力学响应,基于LS-DYNA非线性有限元软件,建立了TNT炸药-前后面板-蜂窝夹芯-空气的三维有限元模型。采用ALE(任意的拉格朗日欧拉)多物质流固耦合算法分析了蜂窝铝板在冲击荷载作用下的变形机理、塑性变形、能量吸收以及结构的优化。数值模拟结果表明:随着面板厚度、核心高度的增加,蜂窝铝板在冲击荷载作用下的塑性变形明显减小,抵抗变形的能力增强;随着爆轰入射角度的增加,结构的破坏程度有所减小,入射角越大这种效果越发明显。对结构给定边长和受冲击面积以及面板厚度配合比、夹芯量纲为一的高度进行了局部的优化分析,为设计优质铝蜂窝板提供参考。     

染料敏化太阳能电池载流子传输的数值模拟

由于在染料敏化太阳能电池(dye-sensitized solar cell, DSSC)中存在染料弛豫、半导体薄膜中电子与氧化态染料分子发生反应和电子在电解质中与氧化态离子复合等不利反应,利用一个更完善的DSSC载流子传输模型对电池的光电性能进行模拟就显得非常重要。为此,本文基于由多重俘获理论建立的DSSC中的包括电子、染料阳离子、碘化物和三碘化物在内的载流子传输模型,数值模拟得到了不同TiO₂薄膜厚度、不同入射光强度与不同染料分子吸收系数下DSSC的J-V曲线。结果表明,随着TiO₂薄膜厚度的增加,太阳能电池的短路电流密度增大,开路电压减小,光电转换效率先增大后减小。当DSSC的TiO₂薄膜厚度为20 μm时,光电转换效率达到最大值7.41%,同时光电转换效率随入射光强度与染料分子吸收系数的增大均有一定程度提高,其中在吸收系数为4 500 cm^-1时,光电转换效率为6.73%。以上结果可以为改进DSSC的光电性能提供理论指导。     

水化硅酸钙纳米压痕分子动力学模拟方法优化

针对水化硅酸钙纳米压痕模型忽视了压头与基底之间相互作用的问题，由尺寸差异引起的金刚石压头难以计算的问题，以及Wittmann模型无法得到实际接触面积的问题，提出了新的模型与计算方法．结合分子动力学方法，采用金刚石压头-Wittmann模型基底的组合方式构建无定形态水化硅酸钙纳米压痕试验模型．在建模阶段，考虑到压头模型与基底模型粒子间尺寸差异，提出了等比例替换模型，通过公式推导并就不同尺寸模拟结果验证了等比例替换模型的可行性．在计算阶段，提出了局部前处理的弛豫方法进行模拟．确定最大荷载位置处的接触面积为546 nm2，进而求出水化硅酸钙模型硬度H为0.84 GPa、折合模量Er为30.52 GPa．并通过纳米压痕试验，验证了模拟结果的准确性，证明了模型的科学性，对今后水化硅酸钙(C-S-H)纳米层面的模拟具有重要借鉴意义．     

Hydrogeochemical modeling for groundwater management in arid and semiarid regions using MODFLOW and MT3DMS: A case study of the Jeffara of Medenine coastal aquifer,South-Eastern Tunisia

The study of water quality and the quantification of reserves and their variations according to natural and anthropogenic forcing is necessary to establish an adequate management plan for groundwater resources. For this purpose, a modeling approach is a useful tool that allows, after calibration phase and verification of simulation, and under different scenarios of forcing and operational changes, to estimate and control the groundwater quantity and quality. The main objective of this study is to collect all available data in a model that simulates the Jeffara of Medenine coastal aquifer system functioning. To achieve this goal, a conceptual model was constructed based on previous studies and hydrogeological investigations. The regional groundwater numerical flow model for the Jeffara aquifer was developed using MODFLOW working under steady-state and transient conditions. Groundwater elevations measured from the piezometric wells distributed throughout the study area in 1973 were selected as the target water levels for steady state (head) model calibration. A transient simulation was undertaken for the 42 years from 1973 to 2015. The historical transient model calibration was satisfactory, consistent with the continuous piezometric decline in response to the increase in groundwater abstraction. The developed numerical model was used to study the system's behavior over the next 35 years under various constraints. Two scenarios for potential groundwater extraction for the period 2015–2050 are presented. The predictive simulations show the effect of the increase of the exploitation on the piezometric levels. To study the phenomenon of salinization, which is one of the most severe and widespread groundwater contamination problems, especially in coastal regions, a solute transport model has been constructed by using MT3DMS software coupled with the groundwater flow model. The best calibration results are obtained when the connection with the overlying superficial aquifer is considered suggesting that groundwater contamination originates from this aquifer. Recommendations for water resource managers

• The results of this study show that Groundwater resources of Jeffara of Medenine coastal aquifer in Tunisia are under immense pressure from multiple stresses.
• The water resources manager must consider the impact of economic and demographic development in groundwater management to avoid the intrusion of saline water.
• The results obtained presented some reference information that can serve as a basis for water resources planning.
• The model runs to provide information that managers can use to regulate and adequately control the Jeffara of Medenine water resources.

     

基于连续模型的人员疏散仿真分析

行人流连续模型直观地反映人群疏散过程中的疏散特征，本文基于行人流连续模型。研究行人在典型疏散场景下的疏散特征．在COMSOL中建立行人流连续模型及其方程，通过编写MATLAB代码，实现了连续模型及其循环求解框架．利用快速扫描法求解Eikonal方程得到背景场值，在每一步迭代循环中将背景场值作为模型的初始变量导入，调用COMSOL计算模块求解模型的瞬态控制方程．通过两个标准算例，重现了典型的行人流自组织现象，验证了连续模型的合理性．结果表明，本文的疏散仿真分析模型和计算程序是可靠的，疏散仿真分析可以为实际工程中的人员疏散方案的制定以及平面设计与安全布置等方面提供技术支撑．     

Strategies to cure numerical shock instability in the HLLEM Riemann solver

The HLLEM scheme is a popular contact and shear preserving approximate Riemann solver that is known to be plagued by various forms of numerical shock instability. In this paper, we clarify that the shock instability exhibited by this scheme is primarily triggered by the spurious activation of the antidiffusive terms present in the first and third Riemann flux components on the transverse interfaces adjoining the shock front due to numerical perturbations. These erroneously activated terms are shown to counteract the favorable damping mechanism provided by its inherent HLL-type diffusive terms, causing an unphysical variation of the conserved quantity ρu both along and across the numerical shock. To prevent this, two distinct strategies are proposed termed as S elective W ave M odification and A nti D iffusion C ontrol. The former focuses on enhancing the quantity of the favorable HLL-type dissipation available on these critical flux components by carefully increasing the magnitudes of certain nonlinear wave speed estimates, while the latter focuses on directly controlling the magnitude of these critical antidiffusive terms. A linear perturbation analysis is performed to gauge the effectiveness of these cures and to estimate a von Neumann–type stability bounds on the CFL number associated with their use. Results from a variety of classic shock instability test cases show that the proposed strategies are able to provide excellent shock stable solutions even on grids that are highly elongated across the shock front without compromising the accuracy on inviscid contact or shear dominated viscous flows.     

高速铁路连续箱型梁桥的动力响应与减振分析

根据列车具体的轴距和轴重,建立了和谐号动车组CRH380AL型列车简化模型;对高速铁路两跨连续梁桥采用多自由度欧拉伯努利梁单元进行主梁的模拟,并将液体黏滞阻尼器模拟为有限元阻尼单元;采用Newmark直接积分方法求解了高速列车作用下的连续梁桥运动方程,数值分析了列车车速以及液体黏滞阻尼器的阻尼系数对于高速铁路连续梁桥振动响应的影响。结果表明:黏滞阻尼器对于桥梁具有明显的减振效果,阻尼力不仅与阻尼系数有关还与列车时速有关;同一黏滞阻尼器条件下,桥梁的最大加速度并不随列车速度的增加而单调增加,而是在某些特定列车车速下桥梁的最大加速度出现了峰值,且随着黏滞阻尼器的阻尼系数增大,桥梁振动响应峰值处的最大加速度减幅不同;同一列车时速的条件下,桥梁的减振效果并不是随着阻尼系数的递增呈正比递增,而是随着阻尼系数的增大,阻尼器的减振效果增幅在减小。     

一种限制周期气固两相流直接数值模拟中整体质量流率的新方法

全周期边界条件的气固两相流系统是研究气固相互作用的一个常用计算模型。为平衡周期气固系统在流向上的受力,常对系统施加平衡系统总重力的固定压力梯度,但由于数值误差的影响,这种固定压力梯度的方法会使得气固系统整体不断加速进而导致模拟失真。为了降低系统质心加速度的影响,本文先后采用了前人提出的固定整体质量流率方法和本文提出的压力梯度后验修正方法,并对两种方法的模拟结果进行了对比。     

群集建筑中大跨裙摆屋盖风荷载特性的数值模拟研究

以苏州太平金融大厦为工程背景,针对其大跨裙摆屋盖的风荷载作用,首先采用RNG k-ε模型模拟分析了其平均风压分布规律,以及风向变化对屋盖表面风荷载体型系数的影响;其次,引入干扰因子IF,探讨了周边建筑对大跨裙摆屋盖风荷载的气动干扰作用。结果表明:0°风向下,走廊上方屋盖两侧区域出现“上吸下顶”的叠加作用;90°风向下屋盖北侧飘带末端区域受到狭道风效应出现正压集中现象;风向变化对大跨裙摆屋盖的风荷载体型系数分布影响较大;且周围建筑物对大跨裙摆屋盖的气动干扰效应明显,主要表现为风压“遮挡效应”,而局部区域表现为风压“放大效应”。