简单剪切流场中的长纤维形态变化的计算机模拟研究.pdfVIP

简单剪切流场中的长纤维形态变化的计算机模拟 王 刚r 许元泽1 周持兴2 摘 要：悬浮在简单剪切流场中的纤维会旋转、平移，当相对流场强度增大到某种程 度时，纤维表现出柔性，在流场中周期性的弯曲变形和重新伸直。如果流场强度足够大， 纤维会出现周期性的断裂，并最终导致一定时间内纤维长度存在双峰分布：原始长度和某 令l痞界值。 关键词：长纤维悬浮体系粒子模型甘算机模拟断裂 纤维悬浮体系存在于广泛的二[业应用中，按照纤维之间的相互碰撞机率可分为极稀 溶液、半浓溶液和浓溶液。熔体的流变性能、“纤维一熔体”及“纤维一纤维”阕的楣互作用 影镌体系的流变性麓，后者影响纾维的形态。一定范围内，纤维增强塑料的强度隧纾维长 度的增加而增大。经过加工，纤维长度会出现不同程度的减小。研究表明，挤出产品中纤 维的长度与原料中的纤维以及加入方式(直接或预处理等)等关系并不明显，而是与熔体 的黏度，剪切流场的强度和加工时间等有关【l’2]。“纤维～纤维”，“纤维一机械”和“纤维 一聚合物”之闻的耀互俸用是影噙制品中纤维长发酶关键蠢索。纤维悬浮体系是一释复 杂的多相体系，具有高度的菲线性力学性质，对其流变性质的实验测定和理论描述均十分 困难，尤其是长纤维除了旋转和取向，还存在着复杂的变形和可能的断裂。粒子模型是一 种有效的研究悬浮体系的手段，可以很方便地研究纤维的形态变化、纤维之间以及纤维与 流体之闻的攘互作孀等，由予每个粒子的运动均能褥裂求解，因此整令体系的微残结构和 宏观性能等也麓褥到很好的溪测。本文采蔫改进的粒子模型，研究了无限大简单剪切流 场中一个包含多根长纤维的参考单元，重点研究了长纤维在流动过程中的形态变化，包括 取向、变形和断裂等。 *联系夫麓余：王躐，募t薅±后，童要获事复杂溅落瀛褒学疆窕。E．mail：跫篓!毪基鱼堑!!：竺婆 ．--——172·--—— 1 模型描述 1．1 粒子模型 将纤维看作是一系列短刚性单元(以下简称短棒)通过铰接方式构成，每个短棒单元 的长度由一个经验公式确定[3]，其受力分析可以通过将短棒离散成更小的球单元来完成， 如图1所示。分析每个短棒受到的流体作用力和力矩，以及短棒内部的拉力、扭矩(弯曲 和扭转)，并遵从牛顿定律和纤维完整性原则，假定纤维末端自由运动作为边界条件，可以 得到一组方程。给定一个初始条件，然后采用差分的方法，在逐渐增加的时间步长计算每 个短棒的速度、位置、方向等参数，从而可以得到纤维的运动规律。 尼 一r 图1 短棒链纤维和球链短 rod--chainfiberandabead--chainrod Fig．1 1．2参考单元 考察无限大简单剪切流场中的一个立方体周期性单元，边长为L，体积V—L3，其中 L，就可以改变体系的浓度乒一Nna2l／L3。 1．3纤维间的相互作用 假设每根纤维内部不存在相互作用，考察流场中分别属于两根不同纤维的两个单元 面的相互距离h。足够小时，由流体润滑理论得到单元歹施加给单元i的作用力： ，。一一K(h“)^d咒d (1) 其中竹。是两个短棒之间的法向矢量，定义为当歹单元向正方向运动时，h“减小；阻尼常数 和狭缝的变化速度分别表示为： Kcho，一编丢帆一参“ ㈤ 1．4纤维断裂临界条件 研究表明，纤维弯曲断裂是主要因素，通常所需应力不到拉伸断裂的四分之一。判断 一】73— 涟 弯曲断裂的经验公式为：“=E1％，其中t-b是临界弯曲曲率半径，巩是纤雏的最大强度，E 是弯瞳模量。在我们的模型中，经过与短棒单元的几何与力学推导，．-f以变形为： 细嫩。sP(巫乒)2) ㈤ 为简化起见，断裂时发生的能量耗散忽略不计，断裂后生成的新端点可自由运动。 一。