基于砼SHPB试验数值分析的HJC模型参数研究

为了研究混凝土HJC本构模型的参数取值及其参数对混凝土动态性能的影响,本文基于HJC模型原始参数,运用有限元软件ANSYS/LS-DYNA对混凝土SHPB试验进行了数值模拟,结果显示数值计算结果与试验结果存在一定的差异。据此根据模型中参数的物理意义确定出可能影响混凝土动态性能的关键参数,通过保持其他参数不变、改变关键参数的方法对混凝土SHPB试验进行数值分析,得到了各关键参数对混凝土动态性能的影响规律,并基于该分析结果对HJC模型的原始参数进行了修正。数值模拟结果表明,采用本文修正参数的计算结果与试验结果吻合良好,能较准确地反映混凝土的动态性能。     

高速流场气动光学RANS/DSMC混合数值模拟算法研究

提出了适用于高速流场气动光学数值模拟的RANS/DSMC混合算法.通过RANS对全局时均流场进行数值模拟,再对其中局部流场的脉动量采用DSMC进行数值模拟,以Maxwellian速度分布实现RANS宏观量信息向DSMC微观量信息的传递.采用超声速环境下尖劈模型对该混合算法进行校验,对比试验结果论证了算法的正确性.在计算...     

高导无氧铜圆柱-平板冲击实验及不同本构模型效果比较

利用一级气炮对高导无氧铜(OFHC)进行了圆柱以205 m/s速度冲击平板实验,并进行了数值模拟。用锰铜应力计测试了靶中应力随时间的变化,并进行了回收观测。采用Johnson-Cook(J-C)、Zerilli-Armstrong(Z-A)、Steinberg-Cochran-Guinan(S-C-G)3种本构模型对实验进行了数值模拟。实验结果与数值模拟结果比较表明:就峰值应力而言,采用J-C、Z-A及S-C-G本构模型的计算结果都比较接近实验;就圆柱变形而言,Z-A及S-C-G模型的计算较J-C模型结果更符合实验。然而,速度为500 m/s冲击实验的数值模拟结果表明:3种本构模型的计算结果差异明显。     

非结构混合网格化学非平衡粘性绕流数值模拟

在非结构混合网格上对化学非平衡粘性绕流进行了数值模拟。控制方程为考虑了化学非平衡效应的二维Navier-Stokes方程,化学动力学模型为7组元、7反应模型。控制方程中的对流项采用VanLeer逆风分裂格式处理,并应用MUSCL方法及Minmod限制器扩展到二阶精度,粘性项用中心差分格式处理。时间推进采用显式5步龙格-库塔方法。为了适应高超声速流场计算,对VanLeer通量分裂方法进行了改进,并引入了化学反应时间步长。对RAMC-II模型的飞行试验流场进行了数值模拟,计算结果与试验测量数据符合较好,并与参考文献中的数值模拟结果吻合。     

混凝土三维细观接触面模型数值模拟与CT试验验证

对混凝土三维随机骨料模型进行了改进,在骨料和砂浆之间加入了"面-面接触单元"模拟两者的接触特性,编制了相应的命令程序,称新模型为混凝土接触面模型.采用双折线损伤演化模型进行计算,对数值模拟结果进行了分析研究.从混凝土破坏过程图和荷载-位移曲线图两方面与CT试验结果进行了比较,表明试件破坏时裂纹的萌生、扩展过程与CT试验...     

混凝土K&C模型材料参数分析与模拟

混凝土K&C模型材料参数一般取国外文献中的原始数值,没有根据混凝土强度等级和单元尺寸的不同而作相应的调整．根据相关的试验研究成果,提出了一种确定K&C模型强度参数值的方法,并阐述了K&C模型损伤参数值的调整方法,使得数值计算结果更加合理．运用有限元显式动力分析软件ANSYS/LSDYNA,采用流固耦合方法模拟爆炸荷载作用下钢筋混凝土板的动态响应,混凝土K&C模型取本文确定的参数值,计算结果与试验结果吻合较好,从而验证了K&C模型材料参数取值的正确性．     

钨合金杆材料属性与侵彻性能的关系

通过对钨合金杆侵彻半无限厚铝合金靶问题的数值模拟,研究了高速侵彻中动能杆材料属性与侵彻性能的影响关系.计算采用的弹塑性流体本构模型中,用来描述材料属性的主要参数有:失效应变、屈服强度、剪切模量等.数值结果表明,失效应变是影响侵深的主要参数,而屈服强度和剪切模量对侵深的影响很小.另外,在此基础上,对着速在785 m/s～1924 m/s范围内的动能杆侵彻问题进行了数值模拟,计算结果与试验数据符合较好.     

筒状弹体垂直侵彻钢筋混凝土墙的数值模拟

采用ANSYS/LS-DYNA软件,对钢筋混凝土墙在受到筒状弹体侵彻撞击、钢筋约束力,及壁面摩擦力等综合作用下的力学效应进行了数值模拟,展现了弹体冲击作用下,钢筋混凝土靶体破碎、飞溅成坑、墙体整体滑移和背面层裂等现象,计算结果与试验宏观破坏现象吻合良好,说明材料模型和参数正确,模拟方法可行.     

基于力热耦合材料模型的Steven试验数值模拟方法

在力热耦合材料模型中,增加炸药自热放能模型,建立了Steven试验的力-热-化耦合的数值模拟方法。数值计算模型中,应力应变关系采用双线性硬化弹塑性模型,炸药受力后的热作用采用各向同性热材料模型,炸药的化学反应采用Arrhenius反应率函数,同时还考虑了升温和熔化对材料力学、热学性能的影响。针对标准Steven试验,通过数值分析得到了靶板的变形情况和炸药点火的速度阈值,将计算结果与实验数据进行了比较,两者符合较好。表明该方法可以较好地模拟Steven试验,而且与以往的分析模型和方法相比,本文的方法不需要增加经验性的点火准则和判据,具有更广泛的适用性,可以为研究低速撞击条件下炸药的力热响应和局域化点火问题提供参考。     

一种直降式滑橇减震器性能研究

为验证无人机伞降回收系统的直降式滑橇着陆装置的动力学性能,对其着陆缓冲过程进行数学模拟构造了该装置的动力学模型,通过数值求解得出了该装置缓冲过程的动态特性。对其着陆压缩过程的动力学模型进行了简化,计算结果可用于减震器参数的确定。落震试验结果表明:动力学模型能准确模拟着陆装置的缓冲过程;以5m/s着陆时过载仅为6.8g;本文所采用动力学模型的数值求解结果与实验结果一致,误差小于5%。     

用代数应力模型模化型腔充填过程中的紊流流动

针对型腔充填过程中的紊流流动，用代数应力模型研究紊流流动现象．成功地将代数应力模型引入ＰＨＯＥＮＩＣＳ软件中，完成了恒温流体充填过程三维时均速度和自由表面的计算模拟．与实验结果相比说明，本文基于ＰＨＯＥＮＩＣＳ软件开发的充填模拟具备了对复杂型腔充填过程的数值计算．     

水模拟连铸结晶器中流速值的氢气泡测试方法

为了获得连铸结晶器中钢液的流速值，以验证数学模拟结果，本文用相似准则原理设计的水模拟板坯连铸结晶，将氢气泡测试技术首先用于其中流速的定量测定，利用氢气泡团的时间－脉线组合图谱所提供的时间－位移信息，测定了不同拉速下水模拟结晶器中若干点的速度及分布，对发泡阴极丝的定位方法和测量位置选择进行了描述，总结了试验中的经验。     

Modeling of bubble behaviors and size distribution in a slab continuous casting mold

Population balance equations combined with Eulerian–Eulerian two-phase model are employed to predict the polydispersed bubbly flow inside the slab continuous-casting mold. The class method, realized by the MUltiple-SIze- Group (MUSIG) model, alongside with suitable bubble breakage and coalescence kernels is adopted. A two-way momentum transfer mechanism model combines the bubble induced turbulence model and various interfacial forces including drag, lift, virtual mass, wall lubrication, and turbulent dispersion are incorporated in the model. A 1/4th scaled water model of the slab continuous-casting mold was built to measure and investigate the bubble behavior and size distribution. A high speed video system was used to visualize the bubble behavior, and a digital image processing technique was used to measure the mean bubble diameter along the width of the mold. Predictions by previous mono-size model and MUSIG model are compared and validated against experimental data obtained from the water model. Effects of the water flow rate and gas flow rate on the mean bubble size were also investigated. Close agreements by MUSIG model were achieved for the gas volume fraction, liquid flow pattern, bubble breakage and coalescence, and local bubble Sauter mean diameter against observations and measurements of water model experiments.     

Convected level set method for the numerical simulation of fluid buckling

‘Fluid buckling’ is a phenomenon characterized mainly by the existence of fluid toroidal oscillations during flow. It appears when a high viscosity fluid flows vertically against a flat surface and may occur in industrial applications, as in injection molding of a propergol in complex‐shaped cavities. These coiling or folding oscillations appear during the mold filling stage, leading to air entrapment. To understand and to model this free surface flow problem, a convected level set method is proposed. First, a sinus filter is applied to the distance function to get a smooth truncation far from the interface. Second, the reinitialization is embedded in the transport equation model, avoiding it as a separate step during calculation. In order to validate the method, numerical results are presented on classical interface capturing benchmarks. Finally, results are shown on two‐dimensional and three‐dimensional viscous jet buckling problems. Copyright © 2010 John Wiley & Sons, Ltd.     

Multiphase flow with a simplified model for oil entrapment

A computationally simple procedure is described to model effects of oil entrapment on three-phase permeability-saturation-capillary pressure relations. The model requires knowledge of airwater saturation-capillary pressure relations, which are assumed to be nonhysteretic and are characterized by Van Genuchten's parametric model; scaling factors equal to the ratio of water surface tension to oil surface tension and to oil-water interfacial tension; and the maximum oil (also referred to as nonwetting liquid in a three-phase medium) saturation which would occur following water flooding of oil saturated soil. Trapped nonwetting liquid saturation is predicted as a function of present oil-water and air-oil capillary pressures and minimum historical water saturation since the occurrence of oil at a given location using an empirically-based algorithm. Oil relative permeability is predicted as a simple function of apparent water saturation (sum of actual water saturation and trapped oil saturation) and free oil saturation (difference between total oil and trapped oil saturation), and water relative permeability is treated as a unique function of actual water saturation. The proposed method was implemented in a two-dimensional finite-element simulator for three-phase flow and component transport, MOFAT. The fluid entrapment model requires minimal additional computational effort and computer storage and is numerically robust. The applicability of the model is illustrated by a number of hypothetical one- and two-dimensional simulations involving infiltration and redistribution with changes in water-table elevations. Results of the simulations indicate that the fraction of a hydrocarbon spill that becomes trapped under given boundary conditions increases as a nonlinear function of the maximum trapped nonwetting liquid saturation. Dense organic liquid plumes may exhibit more pronounced effects of entrapment due to the more dynamic nature of flow, even under static water table conditions. Disregarding nonwetting fluid entrapment may lead to significant errors in predictions of immiscible plume migration.     

连铸结晶器卷渣过程的实验设计及数值模拟

连铸过程结晶器内卷渣过程是冶金流体力学课程教学中的重点和难点,在课堂讲授过程中学生难以掌握复杂的流体理论知识。结合四川省钒钛资源重点实验室现有的冶金过程结晶器模拟装置,通过在冶金流体力学课教学中引入结晶器内卷渣实验和数值模拟,让学生观察到结晶器内钢液流动,激发学生学习兴趣,改善教学效果。通过实验过程、理论知识及数值模拟分析,满足冶金流体力学课程的教学需求。     

Unsteady RANS simulation of oscillating mould flows

Mould flow oscillations are of major importance for the performance of the continuous casting process. They are suspected to promote entrainment of slag and other unwanted secondary phases into the melt pool. These oscillating turbulent flows are investigated by means of numerical simulations. The numerical model is based on the equation of continuity and the unsteady Reynolds averaged Navier–Stokes equations. The system of flow equations is closed by a Reynolds stress turbulence model in combination with non‐equilibrium wall functions. The unsteady simulation resolves low‐frequency oscillations of the flow field. These frequencies and numerically resolved mean values are in agreement with results of corresponding model experiments. The proposed model should be advantageous in order to investigate the mechanisms of the oscillations and the process of slag entrainment in more detail. Copyright © 2006 John Wiley & Sons, Ltd.     

Analysis of Composite Molding Using the Homogenization Method

In this paper, we present an application of the homogenization method to the analysis of Resin Transfer Molding (RTM) and Structural Reaction Injection Molding (SRIM). RTM and SRIM are relatively new molding processes for manufacturing continuous fiber reinforced polymer composites. First, the mold flow is analyzed. In the molding process, the resin experiences significant temperature changes as it fills the mold and forms a free boundary with air as it pushes out the air. In addition, the flow domain is the mold cavity packed with fiber perform, which is a porous medium. Here, the homogenization method is used to model the non-isothermal flow through porous media with free boundaries. A computer program is developed which is capable of simulating a three-dimensional mold flow using the finite element approximation. An example is provided for a three-dimensional part. Then, an analysis of the residual stress developed in the curing stage is given. The curing stage starts when the mold is completely filled and it involves chemical reaction and large temperature variation. In curing, the resin part undergoes larger volume shrinkage than the fiber part, and the residual stresses are developed due to this volume mismatch. In some cases, these stresses are large enough to cause micro-cracking and to exhaust the strength of the material. Here, a brief discussion of the application of the homogenization method to a residual stress analysis is given and one example is provided.     

Flowing simulation of injection molded parts with micro-channel简

In the micro-molding of component with a micro-sized channel, the ability for polymer melt to flowing into the micro-channel in a macro-sized part is a big challenge. The multidimensional flow behaviors are included in the injection molding the macro-component with a micro-channel. In this case, a simplified model is used to analyze the flow behaviors of the macro-sized part within a micro-channel. The flow behaviors in the macro-cavity are estimated by using the finite element and finite difference methods. The influence of the injection rate, micro-channel size, heat transfer coefficient, and mold temperature on the flowing distance is investigated based on the non-isothermal analytic method. The results show that an increase in the radius of the micro-channel and mold temperature can improve effectively the flowing distance in the micro-channel.     

The development and influence of gas bubbles in phreatic aquifers under natural flow conditions

In a phreatic aquifer, bubbles may result from the entrapment of air during groundwater recharge and/or bacterial metabolism. The calculated critical depth of about 1 m at which bubbles are most likely to be found in a granular aquifer, coincides with the depth of 0.60 m of an almost stagnant water layer (specific discharge 1 × 10^-6 cm sec^-1) found at the water table region under natural flow conditions. Bubbles clog pores and therefore reduce the hydraulic conductivity without significantly reducing the volumetric water content. Stagnation at the water table region results since prevailing pressures (in the order of 10^-1 atmospheres) are not sufficiently large to move bubbles through porous media in a water environment.