Experimental and Numerical Analysis of Initial Plasticity in P91 Steel Small Punch Creep Samples

To date, the complex behaviour of small punch creep test (SPCT) specimens has not been completely understood, making the test hard to numerically model and the data difficult to interpret. This paper presents a novel numerical model able to generate results that match the experimental findings. For the first time, pre-strained uniaxial creep test data of a P91 steel at 600 ^°C have been implemented in a conveniently modified Liu and Murakami creep damage model in order to simulate the effects of the initial localised plasticity on the subsequent creep response of a small punch creep test specimen. Finite element (FE) results, in terms of creep displacement rate and time to failure, obtained by the modified Liu and Murakami model are in good agreement with experimental small punch creep test data. The rupture times obtained by the FE calculations which make use of the non-modified creep damage model are one order of magnitude shorter than those obtained by using the modified constitutive model. Although further investigation is needed, this novel approach has confirmed that the effects of initial localised plasticity, taking place in the early stages of small punch creep test, cannot be neglected. The new results, obtained by using the modified constitutive model, show a significant improvement with respect to those obtained by a ’state of the art’ creep damage constitutive model (the Liu and Murakami constitutive model) both in terms of minimum load-line displacement rate and time to rupture. The new modelling method will potentially lead to improved capability for SPCT data interpretation.     

Modifications for an explicit algebraic stress model

An extension of the explicit algebraic stress model, developed by Gatski and Speziale [Gatski TB, Speziale CG. On the explicit algebraic stress models for complex turbulent flows. Journal of Fluid Mechanics 1993; 254: 59–78] is proposed. The extension implicates some essential characteristics of second‐order closure models. The strain‐dependent coefficients are modified, resulting in an alleviation of the numerical instabilities involved in the model. A new near‐wall damping function f_μ in the eddy viscosity relation is introduced. To enhance dissipation in near‐wall regions, the model constant C_ϵ1 is modified and an extra positive source term is included in the dissipation equation. In addition, a realizable time scale is incorporated to remove the wall singularity. Computed results show that the modified Gatski–Speziale (MGS) model predictions are in better agreement with the direct numerical simulation (DNS) and experimental data than those of the original Gatski–Speziale (OGS) model. Copyright © 2001 John Wiley & Sons, Ltd.     

Visco-plastic constitutive modeling on Ohno-Wang kinematic hardening rule for uniaxial ratcheting behavior of Z2CND18.12N steel

Experimental results of monotonic uniaxial tensile tests at different strain rates and the reversed strain cycling test showed the characteristics of rate-dependence and cyclic hardening of Z2CND18.12N austenitic stainless steel at room temperature, respectively. Based on the Ohno-Wang kinematic hardening rule, a visco-plastic constitutive model incorporated with isotropic hardening was developed to describe the uniaxial ratcheting behavior of Z2CND18.12N steel under various stress-controlled loading conditions. Predicted results of the developed model agreed better with experimental results when the ratcheting strain level became higher, but the developed model overestimated the ratcheting deformation in other cases. A modified model was proposed to improve the prediction accuracy. In the modified model, the parameter m_i of the Ohno-Wang kinematic hardening rule was developed to evolve with the accumulated plastic strain. Simulation results of the modified model proved much better agreement with experiments.     

Free coating of viscoelastic and viscoplastic fluids onto a vertical surface

A modified four-constant Oldroyd model is used to predict the film thickness in free coating. The parameters are evaluated from previous steady shear data. The flow is assumed to be one of combined steady shearing and elongation. Despite the use of an integral momentum balance approach the agreement between theory and experiment is not good. However, the results seem to indicate that a more adequate characterisation of the coating fluid might produce better agreement.     

5A06铝合金的动态本构关系实验

运用材料试验机和分离式霍普金森压杆装置(SHPB)对3种不同加工及热处理状态的5A06铝合金在常温~500 C、应变率为10-4~103 s-1 下的力学行为进行了实验研究。基于Johnson-Cook (JC) 本构模型,通过实验数据拟合得到了每种状态下材料的本构模型参量。对Johnson-Cook本构模型中的应变率强化项作了修正,修正后的Johnson-Cook本构模型与实验数据基本吻合,从而确立了3种状态下5A06铝合金的动态本构关系。     

Rheology of polymer blends: non-linear model for viscoelastic emulsions undergoing high deformation flows

Non-linear rheology of a mixture of two viscoelastic immiscible liquids undergoing high deformation flow was considered. Using Grmela's formalism (Grmela 1986, 1989, 1993a, b) and the coarse grained picture given by Onuki (1987) and Doi and Ohta (1991), we have derived a set of highly non-linear time-dependent transport equations that take into account a direct coupling between the rheology and morphology. Breakup, coalescence, and the high deformation of the interface were considered. Models of Doi and Ohta (1991), Lee and Park (1994), and Grmela et al. (1998) are recovered as special cases. The parameters of?the model were given a physical meaning in both shear and elongational flows and the predictions of the model were ¶compared to the predictions of the previous models on the basis of experimental results obtained on two model?blends PDM/PB polydimethylsiloxane/poly(1-butene) and PP/PS polypropylene/polystyrene blends.     

A New modelling strategy for phase-change heat transfer in turbulent interfacial two-phase flow

The paper presents a modelling strategy for phase-change heat transfer in turbulent interfacial two-phase flow. The computational framework is based on interface tracking ITM (level set approach), combined with large-scale prediction of turbulence, a new methodology known as Large-Eddy & Interface Simulation (LEIS), where super-grid scale turbulence and interfaces are directly solved, whereas the sub-scale parts are modelled. Because steady-state flow conditions are difficult to attain, recourse is made of the Very Large-Eddy Simulation (V-LES) instead of LES, where the flow-dependent cut-off filter is larger and independent from the grid. The computational approach is completed by a DNS-based interfacial phase-change heat transfer model built within the Surface Divergence (SD) theory. The original SD model is found to return better results when modified to account for scale separation, i.e. to segregate low-Re from high-Re number flow portions in the same flow. The model was first validated for an experiment involving a smooth to wavy turbulent, stratified steam-water flow in a 2D channel (Lim et al., 1984, Condensation measurement of horizontal concurrent steam-water flow, ASME J. Heat Transfer 106, 425–432.), revealing that the original SD model performs better for high interfacial shear rates. This screening phase also demonstrated that the most critical issue is the accurate prediction of the interfacial shear using ITM. The model was then applied successfully to predict condensing steam in the event of emergency core cooling in a Pressurized Water Reactor (PWR), where water is injected into the cold leg during a postulated loss-of-coolant-accident. The simulation results agree fairly well with the COSI (short for COndensation at Safety Injections) data (Janicot and Bestion, 1993, Condensation modelling for ECC injection, Nucl. Eng. Des. 145, 37–45).     

Prediction of thermal conductivity of composites with spherical fillers by successive embedding

In this study, a new model to predict the thermal conductivity of composites with spherical fillers is proposed. The original Eshelby model is extended to a finite filler volume fraction by successively embedding small filler volume fraction. The predicted results by the present model are compared with bounds such as parallel, series, and Hashin–Shtrikman models, results by modified Eshelby model, generalized self-consistent model, and effective medium theory, and the experimental results from the literature. It is found that the present model always lies between the bounds and shows better agreement with the experimental results than the other models for various filler volume fractions and thermal conductivity ratios.     

Influence of the particle size and volume fraction on micro-damage of the composites

The bulk and shear modulus of metal matrix composites with various volume fractions of particles are modified based on the Eshelby’s equivalent inclusion method combined with self-consistent scheme. By introducing the modified modulus, a new model, which can predict the particle size effects on the stress–strain relation under interfacial debonding damage between matrix and particles, is established. The results obtained from the present investigation show a better agreement with the experimental data.     

Stability of Periodic Solutions¶of Conservation Laws with Viscosity:¶Analysis of the Evans Function

Nonclassical conservation laws with viscosity arising in multiphase fluid and solid mechanics exhibit a rich variety of traveling-wave phenomena, including homoclinic (pulse-type) and periodic solutions along with the standard heteroclinic (shock, or front-type) solutions. Here, we investigate stability of periodic traveling waves within the abstract Evans-function framework established by R. A. Gardner. Our main result is to derive a useful stability index analogous to that developed by Gardner and Zumbrun in the traveling-front or -pulse context, giving necessary conditions for stability with respect to initial perturbations that are periodic on the same period T as the traveling wave; moreover, we show that the periodic-stability index has an interpretation analogous to that of the traveling-front or -pulse index in terms of well-posedness of an associated Riemann problem for an inviscid medium, now to be interpreted as allowing a wider class of measure-valued solutionsor, alternatively, in terms of existence and nonsingularity of a local “mass map” from perturbation mass to potential time-asymptotic T-periodic states. A closely related calculation yields also a complementary long-wave stability criterion necessary for stability with respect to periodic perturbations of arbitrarily large period NT, N → ∞. We augment these analytical results with numerical investigations analogous to those carried out by Brin in the traveling-front or -pulse case, approximating the spectrum of the linearized operator about the wave.The stability index and long-wave stability criterion are explicitly evaluable in the same planar, Hamiltonian cases as is the index of Gardner and Zumbrun, and together yield rigorous results of instability similar to those obtained previously for pulse-type solutions; this is established through a novel dichotomy asserting that the two criteria are in certain cases logically exclusive. In particular, we obtain results bearing on the nature and mechanism for formation of highly oscillatory Turing-like patterns observed numerically by Frid and Liu and ?ani? and Peters in models of multiphase flow. Specifically, for the van der Waals model considered by Frid and Liu, we show instability of all periodic waves such that the period increases with amplitude in the one-parameter family of nearby periodic orbits, and in particular of large- and small-amplitude waves; for the standard, double-well potential, this yields instability of all periodic waves.Likewise, for a quadratic-flux model like that considered by ?ani? and Peters, we show instability of large-amplitude waves of the type lying near observed patterns, and of all small-amplitude waves; our numerical results give evidence that intermediate-amplitude waves are unstable as well. These results give support for an alternative mechanism for pattern formation conjectured by Azevedo, Marchesin, Plohr, and Zumbrun, not involving periodic waves.     

Modeling study of three-phase low liquid loading flow in horizontal pipes

A theoretical study is conducted to model the flow characteristics of three-phase stratified wavy flow in horizontal pipelines with a focus on the low liquid loading condition, which is commonly observed in wet gas pipelines. The model predictions are compared to the experimental data of Karami et al. (2016a, b). These experiments were conducted with water or 51 wt% of MEG in the aqueous phase, and inlet aqueous phase fraction values from 0 to100%.Modeling of three-phase flow can be described as a combination of two-phase gas-liquid flow modeling, and a liquid phase oil-water mixing modeling. A mechanistic model is proposed to predict flow characteristics of three-phase stratified wavy flow in pipeline. For the gas-liquid interactions, Watson's (1989) combined momentum balance equation derivation was applied. However, the calculation procedure was reversed, and the wave celerity was assumed as an input, while interfacial friction factor was one of the model's outputs. The liquid-liquid interactions were modeled using a simple energy balance equation and shift in liquid phase center of gravity calculations. The liquid phases can be separated, partially mixed, or fully mixed. The bottom aqueous film velocity was calculated using the law of the wall formulation, and was used to calculate the flowing aqueous phase fraction.The model predictions of different flow characteristics for two and/or three-phase flows were compared with available experimental data. The pressure gradient, wave amplitude, and aqueous phase fraction predictions were in good agreement with the experimental data. However, the liquid holdup predictions were slightly under-predicted by the model. Overall, an acceptable agreement was observed for all cases.Most of the common multiphase stratified flow models are developed with the assumption of steady-state conditions and with constant interfacial friction factor value. This study proposes a novel method to model stratified flow. The predictions are in acceptable agreement with experimental data conducted under stratified wavy flow pattern conditions.     

Nanocrystalline aluminum and iron: Mechanical behavior at quasi-static and high strain rates,and constitutive modeling

The responses of nanocrystalline aluminum powder of different grain sizes, prepared by ball milling and consolidated into bulk specimens by hot pressing, were determined under quasi-static and dynamic compression. The experiments demonstrated that the reduction in grain size resulted in several-fold increase in hardness and strength; the responses of nanocrystalline aluminum was found to be strain-rate-dependent. Using these measurements, Khan, Huang and Liang (KHL) viscoplastic model was modified by including a bi-linear Hall–Petch type relation to correlate with the response of nanocrystalline aluminum, including the variation of work hardening with grain size. The modified constitutive equation gives results very close to the experimental observations on nanocrystalline aluminum. In addition, the response of nanocrystalline iron, previously published, was also correlated with proposed model. Both correlations of strain-rate-dependent responses for different grain sizes were in good agreement with the experimental results over a wide range of grain sizes (micrometer to nanometer) and strain rates.     

弹体贯穿钢筋混凝土数值模拟

详细描述了依据损伤原理建立的连续损伤模型,并对该模型进行了改进。在LS DYNA程序用户自定义材料模型中加入改进的连续损伤模型,并对弹体侵彻钢筋混凝土的穿孔过程进行了数值模拟,其结果与实验结果相吻合,模型可以用于钢筋混凝土的动态破坏预报。     

Constitutive modeling of shape memory polymer based self-healing syntactic foam

In a previous study, it was found that the shape memory functionality of a shape memory polymer based syntactic foam can be utilized to self-seal impact damage repeatedly, efficiently, and almost autonomously [Li G., John M., 2008. A self-healing smart syntactic foam under multiple impacts. Comp. Sci. Technol. 68(15–16), 3337–3343]. The purpose of this study is to develop a thermodynamics based constitutive model to predict the thermomechanical behavior of the smart foam. First, based on DMA tests and FTIR tests, the foam is perceived as a three-phase composite with interfacial transition zone (interphase) coated microballoons dispersed in the shape memory polymer (SMP) matrix; for simplicity, it is assumed to be an equivalent two-phase composite by dispersing elastic microballoons into an equivalent SMP matrix. Second, the equivalent SMP matrix is phenomenologically assumed to consist of an active (rubbery) phase and a frozen (glassy) phase following Liu et al. [Liu, Y., Gall, K., Dunn, M.L., Greenberg, A.R., Diani J., 2006. Thermomechanics of shape memory polymers: uniaxial experiments and constitutive modeling. Int. J. Plasticity 22, 279–313]. The phase transition between these two phases is through the change of the volume fraction of each phase and it captures the thermomechanical behavior of the foam. The time rate effect is also considered by using rheological models. With some parameters determined by additional experimental testing, the prediction by this model is in good agreement with the 1D test result found in the literature. Parametric studies are also conducted using the constitutive model, which provide guidance for future design of this novel self-healing syntactic foam and a class of light-weight composite sandwich structures.     

The application of 2D base force element method (BFEM) to geometrically non-linear analysis

Based on the concept of the base forces by Gao, a new finite element method – the base force element method (BFEM) on complementary energy principle for two-dimensional geometrically non-linear problems is presented. A 4-mid-node plane element model of the BFEM for geometrically non-linear problem is derived by assuming that the stress is uniformly distributed on each sides of a plane element. The explicit formulations of the control equations for the BFEM are derived using the modified complementary energy principle. The BFEM is naturally universal for small displacement and large displacement problems. A number of example problems are solved using the BFEM and the results are compared with corresponding analytical solutions and those obtained from the standard displacement finite element method. A good agreement of the results, and better performance of the BFEM, compared to the displacement model, in the large displacement and large rotation calculations, is observed.     

理想八面体点阵结构力学模型优化及力学性能测试分析

利用光固化立体成型技术(SLA)制备了不同相对密度的光敏树脂基体理想八面体结构,并对其进行紫外 固化处理,通过理论计算、实验测试以及数值模拟的方式对八面体点阵结构的力学性能进行了分析．实验结果 表明八面体点阵结构的压缩模量以及抗压强度随着相对密度的增加而上升,相对密度从 11.3 %提升至 27.9 %, 弹性模量上升了 3.5倍左右,抗压强度上升 2.6倍左右．利用 ABAQUS有限元软件对实验过程进行了模拟,数 值模拟结果与实验结果具有良好的一致性．在理论计算中,对 Deshpande 等提出的八面体点阵结构计算模型 （DFA模型） 进行了修正,计算结果显示修正后的计算结果更接近实验值与模拟值．     

混凝土骨料对高速侵彻弹体质量侵蚀的影响分析

高速侵彻时，弹靶之间发生强烈的局部作用，引起弹体头部发生质量侵蚀，从而影响弹体的侵彻性能。在弹体侵彻过程中，混凝土中的骨料对弹体的质量侵蚀有显著影响。本文通过对高速侵彻混凝土弹体的质量侵蚀实验数据进行分析，进一步分析讨论了混凝土骨料对弹体质量侵蚀的影响。将混凝土靶体视为骨料和砂浆基质混合的二相复合材料，引入混凝土骨料的体积分数χ和骨料的剪切强度τ₁代替骨料的莫氏硬度H，给出无量纲骨料修正因子β，建立了修正的弹体质量损失工程模型。模型预测结果与现有的实验数据符合得很好，更准确地表征了混凝土骨料对弹体质量损失的影响。     

Removal of shallow shell restriction from Touratier kinematic model

A mathematical model to quantify the variation of the displacement field through the thickness of a laminated shell has been proposed previously by Beakou and Touratier (1993). Transverse shear deformations were taken into account and thickness correction factors were not required but the analysis was restricted to shallow shells (i.e. the principal radii of shell curvature were both assumed to be large relative to the shell thickness). In this technical note the later restriction is removed by replacing the Cartesian coordinate form of elasticity tensor with the more general curvilinear form. The modified laminate level model coefficients are derived from the expressions for the curvilinear transverse shear stress components by: applying zero transverse shear stress boundary conditions at the top and bottom of the shell and enforcing interlayer layer continuity at each internal lamina interface. The purpose of this model enhancement is to facilitate the development of a degenerate finite element type that can be used to compute the deformation of a non-shallow laminated shell.     

跨音速有旋流动正问题的赝势函数变分有限元法

本文运用赝势函数变分有限元方法数值模拟了绕翼型的跨音速有旋流动。在含有激波的跨音速有旋流动中，势函数已不存在，但为了保留势函数模型在求解方面的优越性，上海大学的刘高联引入了一个通用函数一赝势函数，可以看出该赝势函数保持了势函数的所有好处，又突破了流动有势的限制，是势函数对有旋流动的一个自然的、物理上相容、数学上求解简便的推广，进一步地，刘高联还得到了赝势函数的变分原理族，为变分有限元法求解有旋流动打下了基础。另一方面，为了提高数值求解的收敛性和有效地捕捉流场中的激波，本文还采用了“人工密度”办法。绕翼型的跨音速有旋流动的计算实践证明了赝势函数的有效性。     

The Chaboche hardening rule: A re-evaluation of calibration procedures and a modified rule with an evolving material parameter

This study first re-examines key existing calibration procedures for the Chaboche decomposed nonlinear kinematic hardening rule. Special attention has been paid to the extent to which the predictions of calibrated models represent the physics of phenomena observed in strain- and stress-controlled cyclic tests. It has been shown and discussed that newer methods suffer from some weaknesses and may yield nonphysical results. Second, with the aim of improving the accuracy of its ratcheting strain predictions in unsymmetric stress-controlled tests, a modified calibration procedure and a new Chaboche hardening rule with an evolving material parameter is introduced. The former improves the accuracy of predicted ratcheting strains at the initial cycles, while the latter is responsible for the better agreement of the results at the higher cycles. The new model assumes there exists a linear relationship between one of the coefficient of the model and its corresponding effective/equivalent backstress. Finally, the predictions of the new model and those of the original Chaboche rule are compared with the available experimental results from the literature. It has been shown that, regardless of the level of mean stress of unsymmetric stress-controlled test, the ratcheting rates of the higher cycles are satisfactorily predicted by the proposed modified model.