聚合物熔体三维挤出胀大的数值模拟

采用有限元方法分析K-BKZ本构方程描述的聚合物熔体的三维挤出胀大.对于本构方程中偏应力张量的计算,首先给出质点的运动轨迹,分段求出局部的变形梯度张量,再求出整体的变形梯度、Cauchy-Green应变张量和 Finger应变张量,沿轨迹采用分段高斯积分计算应力.把应力作为方程的右端项,给出迭代方法,求解非线性方程组.并根据自由面处的边界条件,迭代得出出口处自由面的最终位置.对轴对称流道和矩形流道进行分析计算,并把结果与二维分析和实验结果进行了比较,显示方法是可行的.     

空间狭缝流道粘弹流体流动的有限柱解法

空间狭缝流道在粘弹性聚合物成型加工中较为常见.针对流道特点,仅仅在流动平面内对速度采用形函数插值,在厚度方向采用傅里叶级数逼近流动分布函数,推导弱解形式的单元方程后,通过坐标变换得到整体坐标下的有限元方程系数矩阵,再集合成整体系数矩阵,从而建立了空间狭缝流动的有限柱解法.分别采用有限柱法和三维有限元对积分型K-BKZ本构模型粘弹流体在L型流道的流动进行求解,发现有限柱法与三维有限元的结果在整体上十分吻合.出口处流量分布的误差小于2%,流量的结果仅仅在流道收敛处略有差异,但差异仅局限于很小的区域.相比与三维有限元方法,有限柱法的单元数、计算时间和对内存需求大大减少.研究表明有限柱法是一种分析狭缝流动的简便有效的方法.     

基于流固耦合方法的水轮机蜗壳结构振动及损伤机理研究

蜗壳流道内的内水压力是引起外围混凝土发生损伤的主要原因.基于流固耦合理论,并引入混凝土弹塑性损伤模型,建立了流体与蜗壳结构耦合振动分析的理论框架,提出了一套水轮机流道内水体流动诱发蜗壳外围混凝土振动损伤的数值计算方法.首先基于有限体积法建立水轮机蜗壳流道流动的数值模型,同时采用有限元方法建立蜗壳结构固体区域的三维有限元模型;进而将流体区域边界上动水压力作为外荷载实时传递给固体区域边界进行三维有限元非线性损伤瞬态分析,实现了大型水轮机蜗壳结构中流体流场到固体应力、位移场的单向耦合三维数值分析.通过计算分析得到了水体流动诱发蜗壳外围混凝土振动的响应规律以及混凝土损伤的发展规律.     

三维非牛顿熔体前沿界面的Level Set/Ghost/SIMPLEC模拟

建立了三维粘性不可压非牛顿流体流动的控制方程,采用Level Set/Ghost/SIMPLEC方法模拟了注塑成型充模阶段的三维流动过程;追踪到了不同时刻的熔体前沿界面,预测并分析了流动过程中不同时刻的压力、速度等重要的流动特征参数,并与牛顿流体相应的流动特征参数做了对比.研究结果表明:Lovol Set/Ghost/SIMPLEC方法可以准确追踪非牛顿熔体前沿界面:幂律熔体在流动过程中的压差明显大于牛顿熔体的压差,沿横截面的速度分布也有明显的差别.     

基于黏弹性的微注射成型流动模拟

应用有限元方法研究了微注射成型中瞬态、可压缩、非牛顿熔体流动的黏弹性对流动前沿及流动平衡的影响。基于Phan-Thien-Tanner模型建立了熔体流动的本构方程,利用Hele-Shaw假设和简化建立了瞬态、可压缩、非牛顿熔体流动的连续性方程、动量方程、能量方程;为了有效地描述微注射成型的尺寸效应,采用了边界滑移和表面张力边界条件。通过分部积分和待定系数法导出了带有边界信息的变分方程和求解应力分量的半解析公式,构造了有限元离散求解及超松驰迭代算法。模拟结果表明:熔体的黏弹性对浇口附近的压力和后续的熔体流动前沿有重要影响;与黏性模型相比,黏弹性模型可以控制模拟压力的快速增长,减少不同型腔之间的充填差异,与短射实验结果也更吻合。     

基于ALE描述的压缩流动三维粘弹性数值模拟

塑料熔体压缩流动中在主流动平面及厚度方向均有速度变化,结合了剪切和拉伸两种流动特征。为了准确描述和模拟塑料熔体的压缩流动,本文基于粘弹性及ALE原理建立了熔体三维流动理论模型,构造了有限元求解的变分方程。为了避免整体求解计算量大、稳定性差的缺陷,提出了两重迭代解耦合算法分别求解耦合的连续方程、动量方程、本构方程、能量方程,开发了模拟程序。开展了等厚度板及变厚度板的注压成型实验及相应的数值模拟,结果表明:压缩过程中出现压力变化小于4.57%的平台现象;温度呈指数规律下降;塑料入口、流动末端第一法向应力差比平均值分别高1.73MPa、0.87MPa,变厚度区域第一法向应力差比平均值高1.16MPa。本文提出的理论模型和数值算法能够较好地表征压缩过程中熔体的压力、温度变化、应力演化。     

变阻流区厚度的T型模具流道设计方法

热塑性塑料片材或流涎膜常用T型模具挤出成型,采用等歧管半径和等阻流区厚度的流道结构时,通常采用增加歧管半径或减小阻流区厚度的方法提高熔体出口流率沿流道宽度方向的均匀性,前者有利于降低挤出压力,但熔体停留时间显著增加,而后者有利于降低熔体停留时间,但挤出压力急剧增大。基于对流道中熔体流动的分析,采用沿流道宽度方向减小歧管半径的结构减小熔体停留时间,沿流道宽度方向增加阻流区厚度的结构降低挤出压力,在满足沿流道宽度方向熔体出口流率均匀的条件下,利用流变学理论推导了阻流区厚度沿流道宽度方向变化的微分方程,对该方程进行数值求解并拟合可得到阻流区厚度。与等歧管半径和等阻流区厚度的流道相比,变歧管半径和变阻流区厚度的流道可以显著降低挤出压力和熔体停留时间。同时,通过对歧管尺寸和形状的设计,在不显著增加挤出压力的情况下能够显著降低熔体停留时间,而通过阻流区长度和厚度的设计,在不增加熔体停留时间的情况下调整挤出压力以适应不同的成型要求。     

三维液体非线性晃动及其复杂现象

主要讨论三维液体非线性晃动问题 。将ALE（任意的拉格朗日-欧拉）运动学描述引人到Navier-Stokes方程的分步有限元计算格式中;在时间域上采用分步离散方法中的速度修正格式,利用Galerkin加权余量方法得到了系统的有限元离散方程;推导了考虑表面张力效应时有限元边界条件的弱积分形式;模拟了三维液体的非线性晃动问题,得到了一系列三维液体非线性晃动的复杂现象.进一步模拟了考虑表面张力效应以及在微重力环境下三维液体的非线性晃动,揭示了考虑表面张力效应以及在微重力环境下液体非线性晃动的重要特征.井将所得结论与现有的实验结果进行了比较.从而证实了该方法的有效性与正确性.     

可压缩多介质粘性流体的数值计算

将考虑热传导和粘性情况下的Navier Stokes方程描述的物理过程分解成3个子过程进行数值计算,即把整个流量计算分解成无粘性流量、粘性流量和热流量3部分,采用多介质流体高精度parabolic piecewise method（PPM）方法、二阶空间中心差方法和两步Rung-Kutta时间推进方法相结合进行数值计算。给出了激波管中Riemann问题和二维、三维Richtmyer-Meshkov界面不稳定性的Navier Stokes方程和Euler方程对比计算结果,显示了粘性对界面不稳定性的影响。     

叶轮机械流道中粘性可压缩流动的数值计算——(Ⅰ)静止叶列中轴对称流动的计算

本文提出一种计算叶轮机械流道中粘性可压缩定常层流问题的数值方法,通过在任意非正交曲线方向上的焓梯度方程、能量方程和熵方程之间的迭代计算,可以得到整个流道中粘性可压缩定常流动问题的数值解。本文首先在静止坐标系中进行分析和讨论,并描述了叶轮机械的静止叶列中轴对称流动的计算方法。计算实例表明,本方法的特点是简单明了,计算速度快,可以广泛地应用于工程设计之中。     

用PSM模型模拟聚合物熔体轴对称收缩流动

采用模拟黏弹流体挤出胀大的方法，计算了IUPAC-LDPE熔体经过4：1轴对称收缩流道的流动．计算的相对涡强度、入口校正和献中的结果基本一致，给出的流场也显示出计算结果是合理的．表明该方法能够适用于用积分型PSM模型表征的黏弹流体在收缩流道内的流动模拟．     

Modeling and simulation of three‐dimensional extrusion swelling of viscoelastic fluids with PTT,Giesekus and FENE‐P constitutive models

The investigation of the extrusion swelling mechanism of viscoelastic fluids has both scientific and industrial interest. However, it has been traditionally difficult to afford theoretical and experimental researches to this problem. The numerical methodology based on the penalty finite element method with a decoupled algorithm is presented in the study to simulate three‐dimensional extrusion swelling of viscoelastic fluids flowing through out of a circular die. The rheological responses of viscoelastic fluids are described by using three kinds of differential constitutive models including the Phan‐Thien Tanner model, the Giesekus model, and the finite extensible nonlinear elastic dumbbell with a Peterlin closure approximation model. A streamface‐streamline method is introduced to adjust the swelling free surface. The calculation stability is improved by using the discrete elastic‐viscous split stress algorithm with the inconsistent streamline‐upwind scheme. The essential flow characteristics of viscoelastic fluids are predicted by using the proposed numerical method, and the mechanism of swelling phenomenon is further discussed.Copyright © 2013 John Wiley & Sons, Ltd.     

NUMERICAL SIMULATION OF FLOW OF HIGHLY VISCOELASTIC FLOW IN THREE－DIMENSIONAL VARYING THICK SLIT CHANNEL

（杨伯源）（李勇）ＮＵＭＥＲＩＣＡＬＳＩＭＵＬＡＴＩＯＮＯＦＦＬＯＷＯＦＨＩＧＨＬＹＶＩＳＣＯＥＬＡＳＴＩＣＦＬＯＷＩＮＴＨＲＥＥ－ＤＩＭＥＮＳＩＯＮＡＬＶＡＲＹＩＮＧＴＨＩＣＫＳＬＩＴＣＨＡＮＮＥＬ￥ＹａｎｇＢｏｙｕａｎ；ＬｉＹｏｎｇ（Ｄｅｐａｒｔ...     

Adaptive finite element method for high-speed flow-structure interaction

An adaptive finite element method for high-speed flow-structure interaction is presented. The cell-centered finite element method is combined with an adaptive meshing technique to solve the Navier-Stokes equations for high-speed compressible flow behavior. The energy equation and the quasi-static structural equations for aerodynamically heated structures are solved by applying the Galerkin finite element method. The finite element formulation and computational procedure are described. Interactions between the high-speed flow, structural heat transfer, and deformation are studied by two applications of Mach 10 flow over an inclined plate, and Mach 4 flow in a channel. The project supported by the Thailand Research Fund (TRF)     

A simple method for simulating viscoelastic fluid flows in slit channel

A low-cost semi-analysis finite element technique, named the finite piece method (FPM) is presented in this article. It aims to solve three-dimensional (3D) viscoelastic slit flows. The viscoelastic stress of the fluid is modelled using an K-BKZ integral constitutive equation of the Wagner type. Picard iteration is used to solve non-linear equations. The FPM is tested on flow problems in both planar and contraction channels. The accuracy of the method is assessed by comparing flow distributions and pressure with results obtained by 3D finite element method (FEM). It shows that the solution accuracy is excellent and a substantial amount of computing time and memory requirement can be saved.     

Viscoelastic flow past confined objects using a micro–macro approach

This paper is concerned with the numerical prediction of viscoelastic flow past a cylinder in a channel and a sphere in a cylinder using molecular-based models. The basis of the numerical method employed is a micro–macro model in which the polymer dynamics is described by the evolution of an ensemble of Brownian configuration fields. The spectral element method is used to discretize the equations in space. Comparisons are made between the macroscopic simulations based on the Oldroyd B constitutive model and microscopic simulations based on Hookean dumbbells, and excellent agreement is found. The micro–macro approach can be used to simulate models, such as the finitely extensible nonlinear elastic (FENE) dumbbell model, which do not possess a closed-form constitutive equation. Numerical simulations are performed for the FENE model. The influence of the model parameters on the flow is described and, in particular, the dependence of the drag as a function of the Weissenberg number.     

Viscoelastic multi-mode simulations of wire-coating

A time-stepping finite element method is used to predict the viscoelastic stresses that arise in a tube-tooling wire-coating problem. The polymer melt HDPE is modelled by a multi-mode Phan-Thien/Tanner constitutive equation. Different flow geometries are considered to address optimisation of the process with respect to minimising the stress induced within the coating produced. The influence of the die itself and the various modes are considered. Relaxation times range for a three-mode model from 10⁻² to 10² s and for a seven-mode model from 10⁻³ to 10³ s. Typical Weissenberg numbers may range up to 10⁴. Three modes are sufficient to adequately describe the flow, and shorter/narrower draw-down regions are identified as being preferable. Once an adequate land length has been gathered, that has relaxed the flow stresses prior to draw-down, the actual details of die design are found to be inconsequential to the induced stresses in the delivered coatings.     

Numerical study of the rheology of rigid fillers suspended in long-chain branched polymer under planar extensional flow

We report a detailed numerical study of the rheology of two-dimensional rigid fillers suspended in branched polymer melt under planar extensional flow. The polymer melt is modelled using the pom-pom constitutive equation. The numerical method uses a finite element solution of the flow in a unit cell within the self-replicating lattice for planar extensional flow identified by Kraynik and Reinelt [Int. J. Multiphase flow 18 (1992) 1045]. This is implemented using a quotient space representation that maps all points in space onto points within the unit cell. We show that the Kraynik and Reinelt cell allows simulations of suspensions under planar extensional flow to be conducted to large strains in a truly periodic cell. The influence of both the pom-pom parameters and the particle area fraction on the rheology of the suspension are investigated. We find a reduction in the degree of extension-rate thickening with the increase of both particles concentration and Weissenberg numbers in agreement with published experimental and numerical findings on other viscoelastic models. This reduction is found to be due to flow disturbance created by the particles which disrupts the alignment of backbone tube segments with the extensional axis.     

形状记忆合金管接头空间轴对称有限元分析

本文采用形状记忆合金（ＳＭＡ）的三维本构方程和有限变形理论，考虑拉、压不同应力状态对相变点移动的规律，编制了ＳＭＡ轴对称大变形的有限元程序，与单向拉伸下解析所得的应力、应变曲线相比，证实程序的正确性．文末计算一ＳＭＡ管接头，并指出按空间轴对称计算的必要性．