夯扩挤密碎石桩加固液化砂土地基的动力数值分析

砂土液化问题一直是土动力学与岩土地震工程研究领域的重要课题之一。基于南水北调中线某工程,通过现场和室内试验获取土体的物理力学参数,利用岩土数值分析软件FLAC3D对夯扩挤密碎石桩加固干渠液化砂土地基进行了动力数值分析。结果表明,由于夯扩挤密碎石桩的排水作用,干渠底部饱和砂土地基中的超静孔隙水压力和孔压比与加固前相比明显减小;干渠渠道底部饱和砂土中的监测曲线表明,随着地震荷载持续时间的增加,饱和砂土地基中超静孔隙水压力和孔压比峰值较加固前大幅值降低,且时程曲线达到峰值之后也由加固前的基本保持不变改为迅速消减降低;由于夯扩挤密碎石桩的排水和挤密作用,有效消除了干渠渠道底部以及渠堤坡面外侧平台至坡脚底部砂土层的液化现象,加固后干渠底部饱和砂土地基中没有液化现象产生。

一种应用于水动力数值分析的浮体湿表面网格生成技术

给出了一种浮体湿表面网格生成方法。利用累加弦长理论对型值点进行等弦长均匀化处理,构造通过均匀型值点的双三次B样条曲面,完成曲面描述,获得曲面网格。以立柱和横梁的网格为基础,利用直线与平面的相交关系,求出立柱与横梁交线上的型值点,依据无限插值理论生成交线周围过渡区的网格,并完成浮体整体网格的拼接。

基于正交设计和数值分析的夯扩挤密碎石桩加固液化砂土方案优化研究

砂土液化问题和处理一直是土动力学与岩土地震工程研究领域的重要课题之一。基于南水北调中线某工程液化砂土处理为例,应用正交设计方法,以桩径、桩长和桩间距3个参数作为正交试验的3个因素,每个因素取3个水平,设计了9组加固处理方案,进行了夯扩挤密碎石桩处理砂土液化动力数值分析,获得了不同加固方案下的饱和砂土中超静孔隙水压力和孔压比。采用极差分析法分析了桩径、桩长和桩间距对超静孔隙水压力和孔压比的影响规律,并做出了各因素对砂土超静孔隙水压力和孔压比影响的直观分析图,得到了夯扩挤密碎石桩的最优加固方案。

双拱初始缺陷海底管线水平向整体屈曲数值模拟分析

为了研究具有双拱反对称初始缺陷海底管线的整体屈曲特性,采用模态分析法将最可能出现的缺陷形态引入数值分析模型中。针对管线在高温高压作用下发生整体屈曲的动态变形特征,运用显式动力数值模拟方法研究了管线整体屈曲过程中水平向变形与轴向变形随温度和内压的变化规律,建立了在整体屈曲过程中屈曲管段与滑动管段轴力的变化过程与初始缺陷形态的关系。将数值模拟结果同经典解析解和室内模型实验结果进行对比,验证了本方法的可靠性。工程算例的分析结果表明,管线整体屈曲的发生是一个由低阶向高阶发展的过程,具有双拱缺陷的管线首先发生二阶模态的整体屈曲,而后过渡到四阶模态;管线整体屈曲的变形包括屈曲段的水平向变形和滑动管段的轴向缩进,其中水平变形释放了管壁内的轴力,轴力的释放量随初始缺陷尖锐程度的降低而增大;轴向缩进变形由于受到地基土的摩阻力使滑动管段内的轴力发生累积,轴力的累积量随初始缺陷的尖锐程度的降低而增加。以上研究成果对指导实际工程具有现实意义。

大跨连续刚构薄壁墩防船撞动力数值模拟

文章根据桥位航道等级及本项目航道通航条件影响评价报告确定船撞代表船型及水流速度,采用国内外五个常见的船撞力计算公式分别计算船撞力,并根据经验选取与实际情况较吻合的美国AASHTO公式数值作为桥墩的船撞力标准值。在中、高水位情况下,选取3°正向和25°侧向两个不同角度撞击桥墩的工况组合,采用大型有限元计算分析程序ANSYS,对桥墩分别进行有限元碰撞仿真计算。对计算结果进行综合分析,分别确定在未加防撞设施情况下和有柔性防撞消能设施情况下,桥墩自身抗力是否满足要求。

NUMERICAL PREDICTION OF AIRCRAFT HYDRAULIC SYSTEM BASED ON THERMAL DYNAMIC ANALYSIS

A mathematical model of principal elements of the aircraft hydraulic system is presented based on the heat transfer theory. The dynamic heat transfer process of the hydraulic oil and the pump shells within an aircraft hydraulic system are analyzed by the difference method. A kind of means for the prediction to variational trends of the aircraft hydraulic system temperature is provided during operation. The numerical prediction and simulation under the operational conditions are presented for ground trial running and the decelerated operation in flight. Computational results show that there is a good coincidence between the experimental data and the numerical predictions.

Magnesium composites with hybrid nano-reinforcements:3D simulation of dynamic tensile response at elevated temperatures

3D numerical simulations of dynamical tensile response of hybrid carbon nanotube(CNT)and SiC nanoparticle reinforced AZ91D magnesium(Mg)based composites considering interface cohesion over a temperature range from 25 to 300℃ were carried out using a 3D representative volume element(RVE)approach.The simulation predictions were compared with the experimental results.It is clearly shown that the overall dynamic tensile properties of the nanocomposites at different temperatures are improved when the total volume fraction and volume fraction ratio of hybrid CNTs to SiC nanoparticles increase.The overall maximum hybrid effect is achieved when the hybrid volume fraction ratio of CNTs to SiC nanoparticles is in the range from 7:3 to 8:2 under the condition of total volume fraction of 1.0%.The composites present positive strain rate hardening and temperature softening effects under dynamic loading at high temperatures.The simulation results are in good agreement with the experimental data.

Simulating Study of Dynamic Load Spectra Identification Method of Machinery in Cepstrum Domain

Based on the platform of Matlab and the theory of digital signal processing, we propose a method in the cepstrum domain for dynamic load spectra identification of machinery. We demonstrate that the dynamic load spectra can be identified from the response signal of the system, based on cepstra. An ARMA model is built based on the harmonic retrieval by high-order spectra. The coefficients of a Green function are determined and the window width can be estimated. Finally the effectiveness of the method is validated by simulation results.

Numerical modeling and dynamic analysis of the 2017 Xinmo landslide in Maoxian County, China

A catastrophic landslide occurred at Xinmo village in Maoxian County, Sichuan Province,China, on June 24, 2017. A 2.87×106 m3 rock mass collapsed and entrained the surface soil layer along the landslide path. Eighty-three people were killed or went missing and more than 103 houses were destroyed. In this paper, the geological conditions of the landslide are analyzed via field investigation and high-resolution imagery. The dynamic process and runout characteristics of the landslide are numerically analyzed using a depth-integrated continuum method and Mac Cormack-TVD finite difference algorithm.Computational results show that the evaluated area of the danger zone matchs well with the results of field investigation. It is worth noting that soil sprayed by the high-speed blast needs to be taken into account for such kind of large high-locality landslide. The maximum velocity is about 55 m/s, which is consistent with most cases. In addition, the potential danger zone of an unstable block is evaluated. The potential risk area evaluated by the efficient depthintegrated continuum method could play a significant role in disaster prevention and secondary hazard avoidance during rescue operations.

CFD Simulation of Fixed and Variable Pitch Vertical Axis Tidal Turbine

In this paper, hydrodynamic analysis of vertical axis tidal turbine （both fixed pitch ＆ variable pitch） is numerically analyzed. Two-dimensional numerical modeling ＆ simulation of the unsteady flow through the blades of the turbine is performed using ANSYS CFX, hereafter CFX, which is based on a Reynolds-Averaged Navier-Stokes （RANS） model. A transient simulation is done for fixed pitch and variable pitch vertical axis tidal turbine using a Shear Stress Transport turbulence （SST） scheme. Main hydrodynamic parameters like torque T, combined moment CM, coefficients of performance Cp and coefficient of torque Cr, etc. are investigated. The modeling and meshing of turbine rotor is performed in ICEM-CFD. Moreover, the difference in meshing schemes between fixed pitch and variable pitch is also mentioned. Mesh motion option is employed for variable pitch turbine. This article is one part of the ongoing research on tm＇bine design and developments. The numerical simulation results are validated with well reputed analytical results performed by Edinburgh Design Ltd. The article concludes with a parametric study of turbine performance, comparison between fixed and variable pitch operation for a four-bladed turbine. It is found that for variable pitch we get maximum Ce and peak power at smaller revolution per minute N and tip sped ratio 2.

Experimental Study of Kinetic Properties of Pyrolysis for Conveyor Belt in Coal Mine

The combustion of conveyor belt is a leading factor of mine fire. In this paper, the pyrolysis properties of ordinary conveyor belt and fire-resistant belt were studied experimentally with thermo-gravimetric analysis and derivative thermo-gravimetric analysis, and the curves of pyrolysis properties were achieved. On this basis, the activation energy and reaction order of pyrolysis were obtained in combination with theoretical analysis, aiming to provide data for further numerical simulation and simulating experiment of mine fire.

Dynamics and Thermodynamics of Porous HMX-like Material Under Shock

Strong shock may induce complex processes in porous materials. We use the newly developed materialpoint-method to simulate such processes in an HMX-like material. To pick out relevant information, morphological characterization is used to treat with the temperature map. Via the Minkowski funetional analysis the dynamics and thermodynamics of the shock wave reaction on porous HMX-like material are studied. The geometrical and topological properties of the ＂hot-spots＂ are revealed. Numerical results indicate that, shocks in porous materials are not simple jump states as classically viewed, but rather are a complex sequence of compressions and rarefactions. They cover a broad spectrum of states. We can use coarse-grained description to the wave series. A threshold value of temperature presents a Turing pattern dynamical procedure. A higher porosity is generally preferred when the energetic material needs a higher temperature for initiation. The technique of data analysis can be used to other physical quantities, for example, density, particle velocity, some specific stress, etc. From a series of studies along the line, one may get a large quantity of information for desiring the fabrication of material and choosing shock strength according to what needed is scattered or connected ＂hot-spots＂. PACS numbers： 05.70.Ln, 05 Key words： porous material 70.-a, 05.40.-a, 62.50.Ef shock wave, Minkowski functionals

Resistance and Seakeeping Numerical Performance Analyses of a Semi-Small Waterplane Area Twin Hull at Medium to High Speeds

The hydrodynamic analysis of a new semi-small waterplane area twin hull （SWATH） suitable for various applications such as small and medium size passenger ferries is presented. This may be an attractive crossover configuration resulting from the merging of two classical shapes： a conventional SWATH and a fast catamaran. The final hull design exhibits a wedge-like waterline shape with the maximum beam at the stem; the hull ends with a very narrow entrance angle, has a prominent bulbous bow typical of SWATH vessels, and features full stern to arrange waterjet propellers. Our analysis aims to perform a preliminary assessment of the hydrodynamic performance of a hull with such a complex shape both in terms of resistance of the hull in calm water and seakeeping capability in regular head waves and compare the performance with that of a conventional SWATH. The analysis is performed using a boundary element method that was preliminarily validated on a conventional SWATH vessel.

Using High-Performance Computing for the Modeling of High-Speed Dynamics

Since computers and software have spread into all fields of industry, extensive efforts are currently made in order to improve the safety by applying certain numerical solutions. For many engineering problems involving shock and impact, there is no single ideal numerical method that can reproduce the various regimes of a problem. An approach wherein different techniques may be applied within a single numerical analysis can provide the ＂best＂ solution in terms of accuracy and efficiency. This paper presents a set of numerical simulations of ballistic tests which analyze the effects of soda lime glass laminates. The goal is to find an appropriate solver technique for simulating brittle materials and thereby improve bullet-proof glass to meet current challenges. To have the correct material model available is not enough. In this work, the main solver technologies are compared to create a perfect simulation model for soda lime glass laminates. The calculation should match ballistic trials and be used as the basis for further studies. These numerical simulations are performed with the nonlinear dynamic analysis computer code ANSYS AUTODYN.

Numerical simulation study of gob air leakage field and gas migration regularity in downlink ventilation

Aiming at the issue that mass of gas emission from mining gob and the gas exceeded in working face, gob air leakage field and gas migration regularity in downlink ventilation was studied. In consideration of the influence of natural wind pressure to analyze the stope face differential pressure, gob air leakage field distribution and gas migration regularity theoretically. Established a two-dimensional physical model with one source and one doab, and applied computational fluid dynamics analysis software Fluent to do numerical simulation, analyzed and contrasted to the areas of gob air leakage on size and gas emission from gob to working face on strength when using the downlink ventilation and uplink ventilation. When applied downward ventilation in stope face, the air leakage field of gob nearly working face, and the air leakage intensity were smaller than uplink, this can effectively reduce the gas emission from gob to working face; when used downlink ventilation, the air leakage airflow carry the lower amount of gas to doab than uplink ventilation, and more easily to mix the gas, reduced the possibility of gas accumulation in upper comer and the stratified flows, it can provide protection to mine with safe and effective production.

Computational analysis of the flow of pseudoplastic power-law fluids in a microchannel plate

The flow of pseudoplastic power-law fluids with different flow indexes at a microchannel plate was studied using computational fluid dynamic simulation.The velocity distribution along the microchannel plate and especially in the microchannel slits,flow pattern along the outlet arc and the pressure drop through the whole of microchannel plate were investigated at different power-law flow indexes.The results showed that the velocity profile in the microchannel slits for low flow index fluids was similar to the plug flow and had uniform pattern.Also the power-law fluids with lower flow indexes had lower stagnation zones near the outlet of the microchannel plate.The pressure drop through the microchannel plate showed huge differences between the fluids.The most interesting result was that the pressure drops for power-law fluids were very smaller than that of Newtonian fluids.In addition,the heat transfer of the fluids through the microchannel with different channel numbers in a wide range of Reynolds number was investigated.For power-law fluid with flow index(n=0.4),the Nusselt number increases continuously as the number of channels increases.The results highlight the potential use of using pseudoplastic fluids in the microheat exchangers which can lower the pressure drop and increase the heat transfer efficiency.

边坡动力响应规律研究

通过大量动力数值分析对边坡动力响应规律进行研究,发现了边坡动力响应的位移、速度、加速度三量放大系数等值线在边坡剖面上分布的节律性特点.对于一定的岩土体材料,当边坡高度在一定范围内时,边坡动力响应的位移、速度、加速度三量在铅直方向会出现放大作用,而当边坡高度超过这个范围时,边坡动力响应的三量分布会出现铅直方向的节律性特点.边坡边缘部位对振动的反应幅值较之内部存在放大现象,随着水平深度的进一步增加,位移、速度、加速度三量在水平方向也出现节律性的变化.边坡坡度的变化对位移、速度、加速度放大系数等值线图的分布形式有明显影响,坡度决定了三量分布的等值线方向和极值区的走向.边坡的动力响应图像可以用两个指标刻画,一个是等值线的走向,一个是平行坡面方向的极值个数.

Roadway dynamic response analysis under mining rockburst condition

In order to determine how a roadway withstands a momentum wave and determinethe extent of damage to rock surrounding the roadway under different force wavepeak impacts,the roadway dynamic response state was analysed using numerical simulationmethod.The roadway's critical peak force wave and fracture region under dynamicwave action were put forward.It is concluded that the method has practical value to roadwaysupport and rockburst prevention.

Hydrodynamics numerical investigation of hoistable masts for underwater vehicles

Using the unsteady incompressible Navier-Stokes equation as the governing equation, the large eddy simulation (LES) model is implemented to investigate the shedding of vortices, the flow pattern of turbulence, the unsteady pressure fluctuation and the time history of the lift coefficient and drag coefficient of hoistable masts with various mast shapes and various arrangements in this paper. Combining the FFT, combined time-frequency transform and wavelet power spectrum analysis, the characteristics of unsteady pressure can be obtained in both time and frequency domain. It shows that the main frequency of pressure fluctuation is near the frequency of vortex shedding in time domain using the FFT method. It can be inferred from the combined time-frequency transform that the unsteady pressure fluctuation has obviously the peak value and the second peak value in time domain. It could indicate that the fluctuation power varies from the fluctuation frequency through the power spectrum analysis. By the data analysis, it shows that the vortex shedding is the dominant cause of the periodically pressure fluctuation. And the interaction pattern of wake and interplay between wake and the walls of masts under different arrangements are also discussed in this paper.

Numerical Simulation and Kinetic Analysis of Turbine Sail

This paper firstly introduces the structure and working principle of turbine sail. Numerical model of a turbine sail is established with Gambit software. The aerodynamic characteristics of the turbine sail are described with RNG k-e turbulence model and the numerical simulation is carded out with Fluent software. The influence of sail＇s structure is analyzed including plate, separation type and height/width ratio. The lift coefficients and drag coefficients of the simulated turbine sail are calculated under different rotation angles, suction intensity and separation plate position. The calculated results are compared with the wind tunnel experimental results, which verifies the feasibility of the numerical results and establishes a foundation for the optimal design of turbine sails.