Multi-scale constitutive modeling of the small deformations of semi-crystalline polymers

A multi-scale constitutive model for the small deformations of semi-crystalline polymers such as high density Polyethylene is presented. Each macroscopic material point is supposed to be the center of a representative volume element which is an aggregate of randomly oriented composite inclusions. Each inclusion consists of a stack of parallel crystalline lamellae with their adjacent amorphous layers.Micro-mechanically based constitutive equations are developed for each phase. A viscoplastic model is used for the crystalline lamellae. A new nonlinear viscoelastic model for the amorphous phase behavior is proposed. The model takes into account the fact that the presence of crystallites confines the amorphous phase in extremely thin layers where the concentration of chain entanglements is very high. This gives rise to a stress contribution due to elastic distortion of the chains. It is shown that the introduction of chains’ elastic distortion can explain the viscoelastic behavior of crystalline polymers. The stress contribution from elastic stretching of the tie molecules linking the neighboring lamellae is also taken into account.Next, a constitutive model for a single inclusion considered as a laminated composite is proposed. The macroscopic stress-strain behavior for the whole RVE is found via a Sachs homogenization scheme (uniform stress throughout the material is assumed).Computational algorithms are developed based on fully implicit time-discretization schemes.     

Fluid-induced vibration of composite natural gas pipelines

Advancements in materials bonding techniques have led to the use of reinforced composite pipelines. The use of steel pipe with a fiber-reinforced composite over-wrap together has produced an exceptionally strong pipe with positive advantages in weight and corrosion resistively. Understanding the dynamic characteristics of this kind of sub-sea composite pipelines, which often accommodate axial flow of gas, and prediction of their response is of great interest.This paper presents a state-variable model developed for the analysis of fluid-induced vibration of composite pipeline systems. Simply supported, clamped and clamped-simply supported pipelines are investigated. The influence of fluid’s Poisson ratio, the ratio of pipe radius to pipe-wall thickness, laminate layup, the ratio of liquid mass density to pipe-wall mass density, the fluid velocity, initial tension and fluid pressure are all considered. The results of our proposed methodology are compared with those of finite element analysis, using ANSYS software.     

Composite power law holdup correlations in horizontal pipes

A wide range of experimental holdup data, from different sources, are analyzed based on a theoretical model proposed in this work to evaluate the holdup in horizontal pipes. 2276 gas–liquid flow experiments in horizontal pipelines with a wide range of operational conditions and fluid properties are included in the database. The experiments are classified by mixture Reynolds number ranges and composite analytical expressions for the relationship between the liquid holdup and no-slip liquid holdup vs. the gas–liquid volumetric flow rate are obtained by fitting the data with logistic dose curves. The Reynolds number appropriate to classify the experimental data for gas–liquid flows in horizontal pipes is based on the mixture velocity and the liquid kinematic viscosity. Composite power law holdup correlations for flows sorted by flow pattern are also obtained. Error estimates for the predicted vs. measured holdup correlations together with standard deviation for each correlation are presented. The accuracy of the correlations developed in this study is compared with the accuracy of 26 previous correlations and models in the literature. Our correlations predict the liquid holdup in horizontal pipes with much greater accuracy than those presented by previous authors.     

Modelling of Non-coalescing Liquid Drops in the Presence of Thermocapillary Convection

A fluid-dynamic model of two non-coalescing liquid drops of the same liquid, pressed against one another in the presence of thermocapillary convection, is proposed to correlate experimental results on the deformation of the drop surfaces, on the pressure distribution and on the thickness of the air film between the drops. The two-point boundary value problem for the Gauss–Laplace hydrostaticequation, subjected to the constant volume constraint, is solved by a fourth order Runge Kutta method to evaluatethe shapes of the drops for different values of the applied pressure jump across the surface exposed to a constant ambient pressure. The flow fields in the liquid drops and in the air layer are obtained by numerical solutions of the dynamic problem. Thenumerical results qualitatively agree with the experimental ones, explain why an air film could be created between the two drops and show that film thicknesses of some microns exist with excess pressures of the same order of magnitude of the pressure needed to deform the drops.     

On the motion of a fluid-saturated porous solid

Two-parameter structure model of a porous solid is proposed as an approximation of a real porous structure and the macroscopic mass and momentum balance equations are derived for such a medium filled with liquid. The approach presented leads to the equations of motion for a fluid-saturated porous medium with coupling terms via cross-mass couplings. The linear form of these equations is equivalent to the well-known Biot equations.     

EQUIVALENT DYNAMICS MODELING AND ANALYSIS OF LIQUID-FILLED SPACECRAFT WITH FUEL CONSUMPTION 1)

在轨航天器贮腔内的液体可能表现出多种不同的运动模式, 主要包括液体相对于贮腔的整体性刚体运动、自由液面横向晃动、液体起旋后逐步发生明显的旋转晃动及液体自旋运动; 复合三自由度刚体摆晃动模型能够较为全面地描述这些液体运动模式, 同时为研究起旋阶段的液体晃动动力学问题提供了有效手段. 本文对非线性液体晃动刚体摆复合模型作进一步发展, 考虑模型等效参数随贮腔充液比的变化, 提出了变参数的刚体摆复合模型, 该模型适用于研究燃料消耗下非线性晃动类充液航天器大范围运动耦合动力学问题. 采用刚体摆复合模型对球形贮腔内的液体晃动进行等效后, 基于混合坐标意义下的拉格朗日方程推导了一类充液航天器轨道-姿态-晃动全耦合的动力学方程组, 并展开了充液航天器大角度三轴稳定姿态机动和零冲量轨道机动仿真以及航天器耦合动力学响应特性分析. 研究表明: 液体相对于贮腔的运动会造成航天器主刚体位置发生偏移, 当航天器在执行零冲量机动时, 燃料消耗会造成航天器的轨道平动速度无法收敛到零; 贮腔偏心布放时, 航天器在执行轨道机动过程中贮腔内液体易发生剧烈而且形式复杂的晃动行为, 进而可能造成航天器刚体运动的不稳定.     

燃料消耗下充液航天器等效动力学建模与分析

在轨航天器贮腔内的液体可能表现出多种不同的运动模式, 主要包括液体相对于贮腔的整体性刚体运动、自由液面横向晃动、液体起旋后逐步发生明显的旋转晃动及液体自旋运动; 复合三自由度刚体摆晃动模型能够较为全面地描述这些液体运动模式, 同时为研究起旋阶段的液体晃动动力学问题提供了有效手段. 本文对非线性液体晃动刚体摆复合模型作进一步发展, 考虑模型等效参数随贮腔充液比的变化, 提出了变参数的刚体摆复合模型, 该模型适用于研究燃料消耗下非线性晃动类充液航天器大范围运动耦合动力学问题. 采用刚体摆复合模型对球形贮腔内的液体晃动进行等效后, 基于混合坐标意义下的拉格朗日方程推导了一类充液航天器轨道-姿态-晃动全耦合的动力学方程组, 并展开了充液航天器大角度三轴稳定姿态机动和零冲量轨道机动仿真以及航天器耦合动力学响应特性分析. 研究表明: 液体相对于贮腔的运动会造成航天器主刚体位置发生偏移, 当航天器在执行零冲量机动时, 燃料消耗会造成航天器的轨道平动速度无法收敛到零; 贮腔偏心布放时, 航天器在执行轨道机动过程中贮腔内液体易发生剧烈而且形式复杂的晃动行为, 进而可能造成航天器刚体运动的不稳定.      

A thermomechanical model for the analysis of light alloy solidification in a composite mould

A coupled thermomechanical model to simulate light alloy solidification problems in permanent composite moulds is presented. This model is based on a general isotropic thermoelasto-plasticity theory and considers the different thermomechanical behaviours of each component of the mould as well as those of the solidifying material during its evolution from liquid to solid. To this end, plastic evolution equations, a phase-change variable and a specific free energy function are proposed in order to derive temperature-dependent material constitutive laws.The corresponding finite element formulation and the staggered scheme used to solve the coupled non-linear system of equations are also presented. Finally, the temperature and displacement predictions of the model are validated with laboratory measurements obtained during an experimental trial.     

非平整、多孔介质海底上波浪传播的复合方程

为了反映近岸区域实际存在的多孔介质海底效应，并且考虑到波浪在刚性海底上传播模型的 最新研究进展，运用Green第二恒等式建立了波浪在非平整、多孔介质海底上传播的复合方 程. 假设水深和多孔介质海底层厚度均由两种分量组成：慢变分量，其水平变化的长度尺度大于 表面波的波长；快变分量，其水平变化的长度尺度与表面波的波长等阶，但其振幅小于表面 波的振幅. 另外，多孔介质层下部边界的快变分量比水深的快变分量小1个量级. 针对水体层和多孔介质层，选择Green第二恒等式方法给出了波浪传播和渗透的复合方程， 它在交接面上满足压力和垂直渗透速度的连续性条件,可充分考虑波数变化的一般连续性， 并包含了某些著名的扩展型缓坡方程.     

求解非线性方程组的混沌粒子群算法及应用

针对非线性方程组的求解在工程上具有广泛的实际意义,经典的数值算法如牛顿法存在其收敛性依赖于初值而实际计算中初值难确定的问题,提出以混沌粒子群算法求解非线性方程。它通过将混沌搜索机制有机地引入粒子群算法,使每个粒子从混沌搜索机制与粒子群算法搜索机制中获得适当的搜索方向,以混沌变量的遍历性增强粒子的搜索性能与更全面地应用目标函数的信息,并反映到逐代更新的个体极值和群体极值中,可更有效地调整粒子的移向并最终获得最优解。测试结果表明这一尝试的有效性。最后将所提的方法用于建立复合材料结构的疲劳寿命与应力、温度、湿度的关系模型。     

General integral equations of thermoelasticity in micromechanics of composites with imperfectly bonded interfaces

One considers a linear thermoelastic composite medium, which consists of a homogeneous matrix containing a statistically inhomogeneous random set of heterogeneities with various interface effects and subjected to essentially inhomogeneous loading by the fields of the stresses, temperature, and body forces (e.g., for a centrifugal load). The general integral equations connecting the stress and strain fields in the point being considered and the surrounding points are obtained for the random and deterministic fields of inclusions. The method is based on a centering procedure of subtraction from both sides of a new initial integral equation their statistical averages obtained without any auxiliary assumptions such as the effective field hypothesis (EFH), which is implicitly exploited in the known centering methods. The new initial integral equation is presented in a general form of perturbations introduced by the heterogeneities and taking into account both the spring-layer model and coherent imperfect one. Some particular cases, asymptotic representations, and simplifications of proposed equations as well as a model example demonstrating the essence of two-step statistical average scheme are considered. General integral equations for the doubly and triply periodical structure composites are also obtained.     

Dynamics of liquid membranes. II: Adaptive finite difference methods

Two domain-adaptive finite difference methods are presented and applied to study the dynamic response of incompressible, inviscid, axisymmetric liquid membranes subject to imposed sinusoidal pressure oscillations. Both finite difference methods map the time-dependent physical domain whose downstream boundary is unknown onto a fixed computational domain. The location of the unknown time-dependent downstream boundary of the physical domain is determined from the continuity equation and results in an integrodifferential equation which is non-linearly coupled with the partial differential equations which govern the conservation of mass and linear momentum and the radius of the liquid membrane. One of the finite difference methods solves the non-conservative form of the governing equations by means of a block implicit iterative method. This method possesses the property that the Jacobian matrix of the convection fluxes has an eigenvalue of algebraic multiplicity equal to four and of geometric multiplicity equal to one. The second finite difference procedure also uses a block implicit iterative method, but the governing equations are written in conservation law form and contain an axial velocity which is the difference between the physical axial velocity and the grid speed. It is shown that these methods yield almost identical results and are more accurate than the non-adaptive techniques presented in Part I. It is also shown that the actual value of the pressure coefficient determined from linear analyses can be exceeded without affecting the stability and convergence of liquid membranes if the liquid membranes are subjected to sinusoidal pressure variations of sufficiently high frequencies.     

Experimental investigation on capillary force of composite wick structure by IR thermal imaging camera

A novel sintered–grooved composite wick structures has been developed for two-phase heat transfer devices. With ethanol as the working fluid, risen meniscus test is conducted to study the capillary force of wick structures. Infrared (IR) thermal imaging is used to identify and locate the liquid meniscus. The effects of sintered layer, V-grooves and powder size on capillary force are explored. The results show that the capillary force of composite wick structures is larger than that of grooved and sintered ones. Interaction wetting between groove and sintered powder happens during the liquid rise in composite wick, which provides an additional source of capillary force. It exhibits a variation of capillary force of composite wicks with different powder size due to the difference of open pore size and quantity in sintered porous matrix.     

Simulation of fiber reinforced composite materials mold filling process and mechanical properties analysis

A dynamic simulation of fiber reinforced composite materials mold filling process with double inlets is presented based on the gas–solid–liquid model proposed by Yang et al. [B.X. Yang, J. Ouyang, J. Tao, C.T. Liu, Modeling and simulation of fiber reinforced polymer mold filling process by level set method, CMES – Computer Modeling in Engineering and Sciences 63 (3) (2010) 191–222]. Numerical results show that the fibers far away from the melt interface are in skin-core-skin structure, while those near the interface are almost parallel to the arc of the interface. When the two streams of melts meet, the weld line will be formed, where the orientation of fibers is perpendicular to the flow direction. The orientation of fibers of the numerical result shows well agreement with the experimental results. Finally, the mechanical properties of fiber reinforced composite materials are analyzed. The composite materials with skin-core-skin structure are regarded as laminated orthogonal plywood and the elastic modulus, the shear modulus and Poisson’s ratio are predicted under different slenderness ratios and fiber volume fractions.     

Three-dimensional nonlinear vibrations of composite beams — I. Equations of motion

Newton's second law is used to develop the nonlinear equations describing the extensional-flexural-flexural-torsional vibrations of slewing or rotating metallic and composite beams. Three consecutive Euler angles are used to relate the deformed and undeformed states. Because the twisting-related Euler angle is not an independent Lagrangian coordinate, twisting curvature is used to define the twist angle, and the resulting equations of motion are symmetric and independent of the rotation sequence of the Euler angles. The equations of motion are valid for extensional, inextensional, uniform and nonuniform, metallic and composite beams. The equations contain structural coupling terms and quadratic and cubic nonlinearities due to curvature and inertia. Some comparisons with other derivations are made, and the characteristics of the modeling are addressed. The second part of the paper will present a nonlinear analysis of a symmetric angle-ply graphite-epoxy beam exhibiting bending-twisting coupling and a two-to-one internal resonance.     

On correct account of finite rotations in finite plasticity theory

The non-uniqueness of the trantition from nonobjective constitutive relations to objective ones with the use of the principle of material frame-indifference (PMFI) is shown. To eliminate it, the concept of finite strain without rotations (FSWR) for a given material type and each strain component (elastic, plastic) is introduced. In FSWR the rotation is excluded with respect to the natural preferred configuration for a given material. Considered are a simple solid, a liquid, a monocrystal, a polycrystal and a composite. The proecedure is proposed for consistent generalization of known infinitesimal relations for finite strains and rotations. The structure of constitutive relations is derived for anisotropic elasto-plastic mono- and polycrystalline materials.     

Mechanical effects of circular liquid inclusions inside soft matrix: role of internal pressure change and surface tension

The mechanical effects of dilute liquid inclusions on the solid-liquid composite are explored,based on an analytical circular inclusion model incorporating the internal pressure change of the liquid and the surface tension of the interface.Several simple explicit dependences of the stress field and effective stiffness on the bulk modulus and the size of the liquid,the surface tension,and Poisson’s ratio of the matrix are derived.The results show that the stresses in the matrix are reduced,and the stiffness of the solid-liquid composite is enhanced with the consideration of either the surface tension or the internal pressure change.Particularly,the effective Young’s modulus predicted by the present model for either soft or stiff matrices agrees well with the known experimental data.In addition,according to the theoretical results,it is possible to stiffen a soft solid by pressured gas with the presence of the surface tension of the solid-gas interface.     

The Riemann problem for a hyperbolic model of two-phase flow in conservative form

This article is to continue the present author's work (International Journal of Computational Fluid Dynamics (2009) 23 (9), 623–641) on studying the structure of solutions of the Riemann problem for a system of three conservation laws governing two-phase flows. While existing solutions are limited and found quite recently for the Baer and Nunziato equations, this article presents the first instance of an exact solution of the Riemann problem for two-phase flow in gas–liquid mixture. To demonstrate the structure of the solution, we use a hyperbolic conservative model with mechanical equilibrium and without velocity equilibrium. The Riemann problem solution for the model equations comprises a set of elementary waves, rarefaction and discontinuous waves of various types. In particular, such a solution treats both the wave structure and the intermediate states of the two-phase gas–liquid mixture. The resulting exact Riemann solver is fully non-linear, direct and complete. On this basis then, we use locally the exact Riemann solver for the two-phase flow in gas–liquid mixture within the framework of finite volume upwind Godunov methods. In order to demonstrate the effectiveness and accuracy of the proposed solver, we consider a series of test problems selected from the open literature and compare the exact and numerical results by using upwind Godunov methods, showing excellent oscillation-free results in two-phase fluid flow problems.     

Dynamics and Stability of a Liquid Triaxial Ellipsoid

The problem of the dynamics of a liquid triaxial ellipsoid confined inside a uniformly rotating rigid shell is formulated and solved within the framework of the Poincare-Joukovski problem. The system of three non-linear first-order differential equations has kinetic energy and circulation integrals representing two ellipsoids with displaced centers in parameter space. All the steady-state and time-dependent solutions are studied; time-dependent solutions exist in four zones, where they are represented by elliptic Jacobi functions, and on three zone boundary interval (represented in terms of elementary functions). On the three other boundary intervals there exist liquid steady-state ellipsoids for which the vortex vector either coincides with one of the principal axes of the ellipsoid or lies in one of its principal planes of symmetry.