Boundary Element Method (BEM) Analysis for Galvanic Corrosion of Hot Dip Galvanized Steel Immersed in Seawater

A numerical analysis of galvanic corrosion of hot-dip galvanized steel immersed in seawater was presented. The analysis was based on the boundary element methods (BEMs) coupled with Newton-Raphson iterative technique to treat the nonlinear boundary conditions, which were determined by the experimental polarization curves. Results showed that galvanic current density concentrates on the boundary of steel substrate and zinc coating, and the sacrificial protection of zinc coating to steel substrate results in overprotection of steel cathode. Not only oxygen reduction but also hydrogen reduction could occur as cathode reactions, which probably led up to the adsorption and absorption of hydrogen atoms. Flat galvanized steel tensile sample shows a brittle behavior similar to hydrogen embrittlement according to the SSRT (show strain rate test) in seawater.

A meso elastoplastic constitutive model for polycrystalline metals based on equivalent slip systems with latent hardening

A meso material model for polycrystalline metals is proposed, in which the tiny slip systems distributing randomly between crystal slices in micro-grains or on grain boundaries are replaced by macro equivalent slip systems determined by the work-conjugate principle. The elastoplastic constitutive equation of this model is formulated for the active hardening, latent hardening and Bauschinger effect to predict macro elastoplastic stress-strain responses of polycrystalline metals under complex loading conditions. The influence of the material property parameters on size and shape of the subsequent yield surfaces is numerically investigated to demonstrate the fundamental features of the proposed material model. The derived constitutive equation is proved accurate and efficient in numerical analysis. Compared with the self-consistent theories with crystal grains as their basic components, the present theory is much simpler in mathematical treatment.

Application of enthalpy method for moving boundary problem in laser surface melting

A numerical analysis was carried out to study the moving boundary problem in the physical process of pulsed Nd-YAG laser surface melting prior to vaporization. The enthalpy method was applied to solve this two-phase axisymmetrical melting problem Computational results of temperature fields were obtained, which provide useful information to practical laser treatment processing. The validity of enthalpy method in solving such problems is presented.

激波在带扩容室的变截面管中传播和消波效果的数值分析

在利用实验结果验证程序的基础上，对激波在具有复杂结构管道内的传播特性进行了数值模拟，比较了不同入射激波条件和特征尺寸对扩容室消波效应的影响，并考虑了二维(轴对称、平面)模型和三维模型的影响。计算结果表明，同一尺寸的扩容室对不同马赫数人射激波的消波效果是相近的，扩容室特征尺寸不同则消波效果差别很大。

高超声速喷管非平衡尺度效应的数值分析

采用轴对称热化学非平衡全N-S方程,数值分析了高超声速喷管流场中非平衡引起的尺度效应.并给出计算条件下喷管轴线上的物理量和化学组元分布.计算结果表明,在喉部下游不远处流动趋近于化学冻结和热力冻结状态.流场(特别是核心无粘流区)的特性不仅取决于流动尺度的相对值,也取决于其绝对值.

不同热化学非平衡模型对高超声速喷管流场影响的数值分析

采用轴对称NS方程，数值研究不同热化学非平衡模型对高超声速喷管流场的影响，包括：（1）不同组元数的高温空气模型（5组元、7组元、11组元）的比较；（2）热力非平衡（双温度）的化学动力过程与热力平衡（单温度）的化学动力过程的比较。计算结果表明，高焓风洞实验条件下喷管流场处于热力和化学都是非平衡的状态。在计算条件下，数值模拟以采用7组元或11组元的热化学非平衡模型为宜。

Wave pattern characteristics of a two-phase nozzle flow by shock propagation

In this paper, the wave pattern characteristics of shock-induced two-phase nozzle Hows with the quiescent or moving dusty gas ahead of the incident-shock front has been investigated by using high-resolution numerical method. As compared with the corresponding results in single-phase nozzle flows of the pure gas, obvious differences between these two kinds of flows can be obtained.

Dynamic damage and failure in viscoplastic materials

In this paper, a dynamic damage model in ductile solids under the application of a dynamic mean tensile stress is developed. The proposed model considers void nucleation and growth as parts of the damage process under intense dynamic loading (strain rates epsilon greater than or equal to 10(3) s(-1)). The evolution equation of the ductile void has the closed form, in which work-hardening behavior, rate-dependent contribution and inertial effects are taken into account. Meanwhile, a plate impact test is performed for simulating the dynamic fracture process in LY12 aluminum alloy. The damage model is incorporated in a hydrodynamic computer code, to simulate the first few stress reverberations in the target as it spalls and postimpact porosity in the specimen. Fair agreement between computed and experimental results is obtained. Numerical analysis shows that the influence of inertial resistance on the initial void growth in the case of high loading rate can not be neglected. It is also indicated that the dynamic growth of voids is highly sensitive to the strain rates.

Spall damage in aluminum-alloy

A void growth relations for ductile porous materials under intense dynamic general loading condition is presented. The mathematical model includes the influence of inertial effects, material rate sensitivity, as well as the contribution of void surface energy and material work-hardening. Numerical analysis shows that inertia appears to resist the growth of voids. The inertial effects increase quickly with the loading rates. The theoretical analysis suggests that the inertial effects cannot be neglected at high loading rates. Plate-impact tests of aluminum alloy are performed with light gas gun. The processes of dynamic damage in aluminum alloy are successfully simulated with a finite-difference dynamic code in which the theoretical model presented in this paper is incorporated.

Influence of inertial and thermal effects on the dynamic growth of voids in porous ductile materials

The influence of inertial, thermal and rate - sensitive effects on the void growth at high strain rate in a thermal - viscoplastic solid is investigated by means of a theoretical model presented in the present paper. Numerical analysis of the model suggests that inertial, thermal and rate - sensitive effects are three major factors which greatly influence the behavior of void growth in the high strain rate case. Comparison of the mathematical model proposed in the present work and Johnson's model shows that if the temperature - dependence is considered, material viscosity eta can take the experimentally measured values.

A new formulated method of a quasi-compatible finite-element and its application in fracture-mechanics

Based on the local properties of a singular field, the displacement pattern of an isoparametric element is improved and a new formulated method of a quasi-compatible finite element is proposed in this paper. This method can be used to solve various engineering problems containing singular distribution, especially, the singular field existing at the tip of cracks. The singular quasi-compatible element (SQCE) is constructed. The characteristics of the elements are analysed from various angles and many examples of calculations are performed. The results show that this method has many significant advantages, by which, the numerical analysis of brittle fracture problems can be solved.

纤维束张紧力缠绕复合材料飞轮初应力的三维数值分析

利用平面应力型全弹性模型的思想(即将纤维束张紧力缠绕看成多层复合材料薄环连续过盈装配的过程),建立了三维纤维束张紧力缠绕复合材料飞轮初应力分析模型,并给出了基于面-面接触算法求解张紧力缠绕复合材料飞轮初应力的三维数值方法.算例分析表明,三维数值分析得到的飞轮的环向初应力及径向初始压应力(数值)均略低于平面应力模型的结果,且这种差距随着飞轮轴向长度的增加而缓慢增大;三维分析证实了平面应力模型关于张紧力缠绕复合材料飞轮的初应力分析有足够的精度.最后给出了三维模型轴向效应的表征方法.

线性组合边界中煤层气井的不稳定压力数值分析

文章建立并求解了不同线性组合边界条件下的煤层气井不定常渗流模型，得到了井底压力和压力导数的理论曲线，以及压力场分布图。在不稳定流动的早期，压力和压力导数曲线都是斜率为1的直线段。当流体流动进入稳定流动状态时，在压力导数曲线上出现一个值为0．5的水平直线段。在理论曲线的晚期，理论曲线与边界性质存在密切关系。对于全封闭煤层气藏，所有压力和压力导数曲线是斜率为1的直线段。对于存在定压边界煤层气藏，压力双对数曲线的晚期是水平直线段，压力导数双对数曲线的晚期都是急剧下掉的曲线。但是混合边界的压力导数存在上翘现象，上翘程度和早晚与封闭边界的数目多少以及距离大小有关。压力场中，在靠近封闭边界的等压线垂直于边界，在靠近定压边界的等压线平行于边界。在邻近定压和封闭边界交合处，定压边界起主导作用。所得到的理论曲线可以直接应用于煤层气田的试井测试资料分析，压力场分布图对了解煤层气藏的开发状况有积极的意义。

Nonlinear Dynamic Response of Tension Leg Platform

Tension Leg Platform (TLP) is a typical compliant offshore structure for oil exploitation in deep water. Most of the existing mathematical models for analyzing the dynamic response of TLP are based on explicit or implicit assumptions that displacements (translations and rotations) are small magnitude. Herein a theoretical method for analyzing the nonlinear dynamic behavior of TLP with finite displacement is developed, in which multifold nonlinearities are taken into account, i.e. finite displacement, coupling of the six degrees of freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force. Using this theoretical model, we perform the numerical analysis of dynamic response of a representative TLP. The comparison between the degenerative linear solution of the proposed nonlinear model and the published one shows good agreements. Furthermore, numerical results are presented which illustrate that nonlinearities exert a distinct influence on the dynamic responses of the TLP.

Vortex dislocations in wake-type flow induced by spanwise disturbances

Vortex dislocations in wake-type flow induced by three types of spanwise disturbances superimposed on an upstream velocity profile are investigated by direct numerical simulations. Three distinct modes of vortex dislocations and flow transitions have been found. A local spanwise exponential decay disturbance leads to the appearance of a twisted chainlike mode of vortex dislocation. A stepped spanwise disturbance causes a streamwise periodic spotlike mode of vortex dislocation. A spanwise sinusoidal wavy disturbance with a moderate waviness causes a strong unsteadiness of wake behavior. This unsteadiness starts with a systematic periodic mode of vortex dislocation in the spanwise direction followed by the spanwise vortex shedding suppressed completely with increased time and the near wake becoming a steady shear flow. Characteristics of these modes of vortex dislocation and complex vortex linkages over the dislocation, as well as the corresponding dynamic processes related to the appearance of dislocations, are described by examining the variations of vortex lines and vorticity distribution. The nature of the vortex dislocation is demonstrated by the substantial vorticity modification of the spanwise vortex from the original spanwise direction to streamwise and vertical directions, accompanied by the appearance of noticeable vortex branching and complex vortex linking, all of which are produced at the locations with the biggest phase difference or with a frequency discontinuity between shedding cells. The effect of vortex dislocation on flow transition, either to an unsteady irregular vortex flow or suppression of the Kaacutermaacuten vortex shedding making the wake flow steady state, is analyzed. Distinct similarities are found in the mechanism and main flow phenomena between the present numerical results obtained in wake-type flows and the experimental-numerical results of cylinder wakes reported in previous studies.