纤维悬浮流中的方向分布函数和张量研究.pdfVIP

纤维悬浮流中的方向分布函数及张量研究‘ 林建忠 张凌新 吴法理 (浙江大学力学系，杭州，310027)(浙江菲达机电集团公司，浙江诸暨) 摘要本文叙述了纤维悬浮流的应用背景及．论述了纤维方向分布函数和方向张量在描述悬 浮流中的重要性，给出了分布函数和方向张量的表述形式和求解方法．并将其应用于楔形流 场的计算． 关键词纤维悬浮流，分布函数，方向张量，数值模拟 1．纤维悬浮流的应用背景 纤维包含在流体中的纤维悬浮流，是多相流和非牛顿流理论研究中的难点问题，同时， 它还具有广泛的应用背景。在复合材料中，短纤维复合材料是最常见的，它的成型和加工就 与纤维悬浮流的动力学特性密切相关，加工过程中由流动诱导的纤维取向、纤维对于流动的 影响直接决定了成品的质量。目前一种新型的人工木质材料．是通过将原木碎化成木纤维， 然后对木纤维在高压下进行重组．形成质地坚硬的材料，这种材料可以取代现有的一些金属， 而纤维的重组过程与气纤的相互作用有关。纺织工业中大部分纺织品的成型与加工都是短纤 维分组成束、成片的过程，该过程是典型的气纤流动，纤维在气流中的指向、均匀度、绕卷 度、扭曲度等直接决定了成品的质地、强度和手感。此外，纤维悬浮流还与造纸过程、陶瓷 加工等有着密切的关系。 2．纤维方向分布函数及张量的作用 纤维的方向分布函数是描述纤维在流场中主轴取向的参数．而纤维在流场中取向直接 决定了材料的特性。取向的差异导致材料的力学性质不是各向同性．通常顺着纤维主轴的方 向具有最大的强度和硬度，而其它方向则相对弱些。因此，为了更好地制造和使用材料，就 必须了解纤维的取向，进而通过人工控制使得材料符合使用的要求。 描述纤维取向的最简单的方法是使用一个标度，通常是纤维与一个坐标轴之间的夹角： 对于三维情况，是在球坐标下使用两个角度量西和口来刻画取向状态。但是，这种简明的描 述方法却不适应于确定流变属性。为了使纤维的取向描述与流变属性之间的关系更直接一 些，必须引进方向分布函数。 方向分布函数定义为在给定时间、位置下，纤维处在某一角度的概率。这种定义既能 够使它与方向描述紧密相连，同时又能在研究流变特性的应力构造方程中起到重要作用，如 应力的表述为：f=一，7，户+cf9，啊?PPPPdP，其中妒是分布函数。 方向分布函数也存在缺陷．它用来描述简单的流动如平面流比较方便，但对于三维流 动以及复杂边界的流场，计算方向分布函数显得过于复杂，并且刻画流变特性不够简洁。 为了克服分布函数的缺陷．可咀引入方向张量。方向张量最初由分布函数间接积分而得 +国家杰出青年科学基金资助项目(№ 到，仅作为方向状态对应力张量影响的中间量。随着对方向张量研究的深入，现在它独立地 用来描述纤维的取向分布特性。虽然它看上去比较复杂，但在描述流体的流变特性时却有效 而清晰，因此它比方向分布函数更为优越。 求解分布函数，用的是分布函数的控制方程，即F0kker—Planck方程：而求解方向张 量，用的是由Fokker—P1anck方程经过适当变形后得到的张量的控制方程。这个方程并不封 闭，即二阶的控制方程含有四阶的张量．四阶的控制方程含有六阶的张量……。 3．分布函数的求解 分布函数伊晓，卢，，)定义为时刻r、位置‘处，在方向户上存在粒子的可能性．即概率。 因此在d芦范围内，存在粒子的数目可以简单地写为妒伍，户．f舻。可见，分布函数既能方 便地描述方向状态，又能直接同粒子的数目联系在一起。 分布函数满足如下的归一化和周期性条件： e‘p船l矽=1，妒船)=妒移+疗) 其中，庐是纤维轴线与流场对称轴线的夹角· 分布函数的控制方程由Fokker-Planck方程描述…： 鲁=一品∽)新 ad”“ cz，～ 以上是一维波动方程，乒是纤维角度≯和位置的函数，对单根纤维，有以下表达式… 净警嚣+警+挚删～扣嘞胁z妒’驴酗2 2’4……’’ (3) 其中户，是应变率张量分量，国F是涡状态张量，它们表示为： 2塑堂+蚓 o [％】= 3