流变学的前沿方向—液晶高分子流变学

简要叙述液晶高分子的结构,论述液晶分子流变学特征,液晶高分子流体的各类本构方程,其中包括作者发展的液晶高分子流体型本构方程。     

共转Oldroyd B流体剪切-拉伸流动研究

应用宏观理论的共转模型研究纺丝拉伸过程中测粘-拉伸流动的流变学，以探讨液晶高分子流体纺丝拉伸初步机制。由共转Oldroyd B流体模型的本构方程导出了测粘-拉伸流动的法向应力差和拉伸粘度的解析表达式，利用边界条件及运动方程列出了主拉伸速度与纺程的微分方程。应用Matlab软件绘制出拉伸粘度和主拉伸速度随其它参数变化的曲线，从而得到这些参数对纺丝拉伸流动流变学性能的影响，这对纺丝过程的工艺控制有一定的指导意义。     

非定常微分型黏弹性本构模型

在应力作用下, 材料的力学参数随着微观结构的变化而变化, 需要考虑参数的时间效应. 利用黏滞系数随时间变化的黏性元件, 构造出非定常Maxwell模型、非定常Kelvin模型和非定常Zener模型. 求解非定常模型的微分型本构方程得到它们的松弛模量、蠕变柔量和卸载方程. 结果表明, 可以把常见的经验松弛函数和经验蠕变函数视为非定常微分型本构模型.     

液晶高分子-各向异性流体挤出拉伸分岔研究

应用共转导数型本构方程研究了液晶高分子纺丝挤出过程的拉伸黏度,应用计算机符号运算软件 Maple得出解析表达式,拉伸黏度与拉伸率之间关系(随剪切速率变化)表明存在分岔现象,得出拉伸黏度显著高于相应的剪切黏度,解释了液晶高分子熔体挤出时不发生挤出胀大的物理机制．     

材料三维微结构表征及其晶体塑性有限元模拟

材料的力学性能,尤其是在有限变形下所呈现的宏观各向异性,是材料结构设计和服役寿命考虑的关键因素。由于宏观模型不能较好地反映材料微观结构(晶粒的形貌和取向等)对宏观塑性各向异性的影响,因此,本文建立了能实际反映晶粒形貌的三维Voronoi模型,并基于晶体塑性理论对铝合金在有限变形下的响应进行计算。首先,建立反映材料微结构的代表性体积单元RVE模型进行计算,并与实验结果进行对比验证。然后,以单向拉伸为例,分析了有限变形过程中试件的晶粒形貌和取向分布等微观因素对宏观各向异性演化的影响,并从材料和结构两个层面讨论了微观结构对宏观力学性能的影响。结果表明,本文模型能够反映微观结构对宏观力学性能的影响,为实际生产制造领域构件的力学性能提供可靠的预测。     

材料三维微结构表征及其晶体塑性有限元模拟

材料的力学性能，尤其是在有限变形下所呈现的宏观各向异性，是材料结构设计和服役寿命考虑的关键因素。由于宏观模型不能较好地反映材料微观结构（晶粒的形貌和取向等）对宏观塑性各向异性的影响，因此，本文建立了能实际反映晶粒形貌的三维Voronoi模型，并基于晶体塑性理论对铝合金在有限变形下的响应进行计算。首先，建立反映材料微结构的代表性体积单元RVE模型进行计算，并与实验结果进行对比验证。然后，以单向拉伸为例，分析了有限变形过程中试件的晶粒形貌和取向分布等微观因素对宏观各向异性演化的影响，并从材料和结构两个层面讨论了微观结构对宏观力学性能的影响。结果表明，本文模型能够反映微观结构对宏观力学性能的影响，为实际生产制造领域构件的力学性能提供可靠的预测。     

液晶高分子取向特征及其流变性质研究

研究液晶高分子取向特征及其流变性质，在向矢输运方程基础上，研究了管内Posieulle流动和纺丝拉伸流动中一维向矢部分，应用偏光显微镜研究了热致性和溶致性液晶高分子织构，显示出液晶态的各向异性特征，测定了溶致性液晶高分子羟基丙基纤维素（HPC）的流变特性。     

基于格子Boltzmann方法的粘弹性流体圆柱绕流分析

将光滑界面法引入到格子Boltzmann方法中分析粘弹性流体绕流问题，分别采用单松弛模型和对流扩散模型求解运动方程和Oldroyd-B本构方程，针对圆形和椭圆内部边界条件，给出连续界面插值函数，在此基础上，运用光滑界面法将内部边界转换为作用力项施加到演化方程中。首先分析圆柱绕流问题，给出不同材料参数情况下的流场分布和阻力系数计算结果，比较发现与宏观数值模拟结果相吻合。将模型拓展到绕椭圆流动中，分析椭圆形状和材料参数对粘弹性流体绕柱流的影响，发现随着椭圆长轴与短轴比值的增加和维森伯格数的增加，阻力系数逐渐下降，并且长短轴比对迭代收敛有较大影响。     

张量函数的表示理论──本构方程统一不变性研究（续上期）

张量函数的表示理论──本构方程统一不变性研究（续上期）郑泉水（Ｑ－ＳＺｈｅｎｇ）清华大学工程力学系，北京１０００８４４矢量和二阶张量的二维所有种类各向异性张量函数如前所述，理论和应用力学的研究者对平面问题一直怀有浓厚的兴趣．三维空间各向异性物理行为通...     

具有体膨胀和剪切效应的结构陶瓷相变塑性细观本构模型:Ⅰ.非比例加载历史

本文在对结构陶瓷的四方至单斜(t→m)马氏体相变进行细观力学、热力学和微观机制分析的基础上,导出了在非比例加载条件下考虑材料的体膨胀和剪切效应的相变塑性细观本构模型。作者首次采用 Mori-Tanaka 方法以自洽的方式导出了材料构元的 Helmho-ltz 自由能及余能函数的解析表达式,它是外加宏观应力(或应变)、温度、相变夹杂体积分数以及夹杂内平均相变应变的函数,其中夹杂体积分数和平均相变应变为描述材料构元微结构变化的内变量。最后按 Hill-Rice 本构理论框架导出相变塑性屈服面方程及增量本构关系。     

Constitutive equation of co-rotational derivative type for anisotropic-viscoelastic fluid

A constitutive equation theory of Oldroyd fluid B type, i.e. the co-rotational derivative type, is developed for the anisotropic-viscoelastic fluid of liquid crystalline (LC) polymer. Analyzing the influence of the orientational motion on the material behavior and neglecting the influence, the constitutive equation is applied to a simple case for the hydrodynamic motion when the orientational contribution is neglected in it and the anisotropic relaxation, retardation times and anisotropic viscosities are introduced to describe the macroscopic behavior of the anisotropic LC polymer fluid. Using the equation for the shear flow of LC polymer fluid, the analytical expressions of the apparent viscosity and the normal stress differences are given which are in a good agreement with the experimental results of Baek et al. For the fiber spinning flow of the fluid, the analytical expression of the extensional viscosity is given. The project supported by the National Natural Science Foundation of China (19832050 and 10372100)     

An unsymmetric constitutive equation for anisotropic viscoelastic fluid

A continuum constitutive theory of corotational derivative type is developed for the anisotropic viscoelastic fluid–liquid crystalline (LC) polymers. A concept of anisotropic viscoelastic simple fluid is introduced. The stress tensor instead of the velocity gradient tensor D in the classic Leslie–Ericksen theory is described by the first Rivlin–Ericksen tensor A and a spin tensor W measured with respect to a co-rotational coordinate system. A model LCP-H on this theory is proposed and the characteristic unsymmetric behaviour of the shear stress is predicted for LC polymer liquids. Two shear stresses thereby in shear flow of LC polymer liquids lead to internal vortex flow and rotational flow. The conclusion could be of theoretical meaning for the modern liquid crystalline display technology. By using the equation, extrusion–extensional flows of the fluid are studied for fiber spinning of LC polymer melts, the elongational viscosity vs. extension rate with variation of shear rate is given in figures. A considerable increase of elongational viscosity and bifurcation behaviour are observed when the orientational motion of the director vector is considered. The contraction of extrudate of LC polymer melts is caused by the high elongational viscosity. For anisotropic viscoelastic fluids, an important advance has been made in the investigation on the constitutive equation on the basis of which a series of new anisotropic non-Newtonian fluid problems can be addressed. The project supported by the National Natural Science Foundation of China (10372100, 19832050) (Key project). The English text was polished by Yunming Chen.     

Spectral element method for viscoelastic flows in a planar contraction channel

A new algorithm, which combines the spectral element method with elastic viscous splitting stress (EVSS) method, has been developed for viscoelastic fluid flows in a planar contraction channel. The system of spectral element approximations to the velocity, pressure, extra stress and the rate of deformation variables is solved by a preconditioned conjugate gradient method based on the Uzawa iteration procedure. The numerical approach is implemented on a planar four‐to‐one contraction channel for a fluid governed by an Oldroyd‐B constitutive equation. The behaviour of the Oldroyd‐B fluids in the contraction channel is investigated with various Weissenberg numbers. It is shown that numerical solutions obtained here agree well with experimental measurements and other numerical predictions. Copyright © 2003 John Wiley & Sons, Ltd.     

A VARIATIONAL ANALYTICAL APPROACH TO TIME DEPENDENT FLOW OF OLDROYD B FLUID IN CIRCULAR TUBE

ＡＶＡＲＩＡＴＩＯＮＡＬＡＮＡＬＹＴＩＣＡＬＡＰＰＲＯＡＣＨＴＯＴＩＭＥＤＥＰＥＮＤＥＮＴＦＬＯＷＯＦＯＬＤＲＯＹＤＢＦＬＵＩＤＩＮＣＩＲＣＵＬＡＲＴＵＢＥＨａｎｓｈｉ－ｆａｎｓ（韩式方）（Ｒｅｓ．Ｌａｂ．Ｎｏｎ－ＮｅｗｔｏｎｉａｎＦｌｕｉｄＭｅｃｈ...     

A variational analytical approach to time dependent flow of oldroyd B fluid in circular tube

In the present investigation the time dependent flow of an Oldroyd fluid B in a horizontal cylindrical pipe is stuided by the variational analytical approach developed by author. The time dependent problem is mathematically reduced to a partial differential equation of third order. Using the improved variational approach due to Kantorovich the partial differential equation can be reduced to a system of ordinary differential equations for different approximations. The ordinary differential equations are solved by the method of the Laplace transform which is led to an analytical form of the solutions. Project supported by TWAS and Chinese Academy of Sciences and the National Science Foundation of China     

Coupling of finite element method and discontinuous Galerkin method to simulate viscoelastic flows

In this paper, a numerical method, which is about the coupling of continuous and discontinuous Galerkin method based on the splitting scheme, is presented for the calculation of viscoelastic flows of the Oldroyd‐B fluid. The momentum equation is discretized in time by using the Adams‐Bashforth second‐order algorithm, and then decoupled via the splitting approach. Considering the Oldroyd‐B constitutive equation, the second‐order Runge‐Kutta approach is selected to complete the temporal discretization. As for the spatial discretizations, the fundamental purpose is to make the best of finite element method (FEM) and discontinuous Galerkin (DG) method to handle different types of equations. Specifically speaking, for the subequations, FEM is chosen to treat the Poisson and Helmholtz equations, and DG is employed to deal with the nonlinear convective term. In addition, because of the hyperbolic nature, DG is also utilized to discretize the Oldroyd‐B constitutive equation spatially. This coupled method avoids resorting to extra stabilization technique occurred in standard FEM framework even for moderately high values of Weissenberg number and also reduces the complexity compared with unified DG scheme. The Oldroyd‐B model is applied to investigate several typical and challenging benchmarks, such as the 4:1 planar contraction flow and the lid‐driven cavity flow, with a wide range of Weissenberg number to illustrate the feasibility, robustness, and validity of our coupled method.     

Squeeze film flow of ideal elastic liquids

An exact solution is presented for the squeeze film flow of an Oldroyd B. fluid. The solution demonstrates that the flow kinematics is similar to the Newtonian (or Maxwellian) one. Theoretical predictions for constant velocity squeezing are compared to experimental observation for well characterized non-shear thinning elastic fluids. It is shown both theoretically and experimentally that the effect of elasticity in a constant velocity squeeze film flow is to always reduce the load relative to the inelastic (Newtonian) prediction and that this load reduction falls between the upper and lower asymptote prediction by the exact solution for the Oldroyd B fluid. The upper load asymptote is given by the Stefan solution for the viscosity of the polymer solution and the lower asymptote is given by the Stefan solution for the viscosity of the solvent. Experimental observations agree with the theoretical prediction for the Oldroyd B fluid at low shear rates where it is shown that the steady and dynamic flow properties of the test fluids used in the experimental program are well represented by the Oldroyd B constitutive equation. With the exception of the work of Lee et al. [6] for constant load squeezing of a Maxwell fluid, this work represents one of the few cases where experimental observation of large effects due to elasticity are indeed predicted with a constitutive equation which actually describes the steady and dynamic shear properties of the fluids used in the experimental program.     

On non-linear magnetohydrodynamic problems of an Oldroyd 6-constant fluid

In this work we develop a mathematical model to predict the velocity profile for an unidirectional, incompressible and steady flow of an Oldroyd 6-constant fluid. The fluid is electrically conducting by a transverse magnetic field. The developed governing equation is non-linear. This equation is solved analytically to obtain the general solution. The governing non-linear equation is also solved numerically subject to appropriate boundary conditions (three cases of typical plane shearing flows) by an iterative technique with the finite-difference discretizations. A parametric study of the physical parameters involved in the problems such as the applied magnetic field and the material constants is conducted. The obtained results are illustrated graphically to show salient features of the solutions. Numerical results show that the applied magnetic field tends to reduce the flow velocity. Depending on the choice of the material parameters, the fluid exhibits shear-thickening or shear-thinning behaviours.