对本构方程书写形式的探讨.pdf

对本构方程书写形式的探讨 张 平1 董贤勇2 张绍东2 李兆敏2 (1．石油大学(华东)储建学院，东营257061； 7061) 2．石油大学(华东)石油工程学院，东营25 摘要 笔者在阅读非牛顿流体圆管和环空流动有关文献和书籍时，发现有些文章的推导 结论错误，其错误产生原因在于对运动方程和非牛顿流体本构方程中出现的切应力出现了概 念上的混淆。为此，笔者在分析流体运动方程中各项物理意义的基础上．提出了一种统一的书 写本构方程的方法，使之用于研究非牛顿流体圆管和环空管流层流流动规律。结果表明，该方 法物理意义明确，可得出有关流动的正确解，从而克服了以往方法中的不足。 关键词 非牛顿流体，本构方程，圆管流动，环空流动，统一方法 l 前 言 表述物质受力与运动响应之间的关系的方程称为本构方程。建立本构方程有多种途径，大 体可以分为经验方法、半经验半理论方法和理论方法。就目前的情况来看，完全采用理论的方 法构造出的本构方程一般都比较复杂，很难用来解决工程问题。所以通常采用经验的方法，结 合一定的理论分析，构造出适合工程应用的一些本构方程Ⅲ。 流体的本构方程是研究流动问题的前提条件，其对流动问题的解具有实质性的影响。 2 问题的提出 在研究非牛顿流体圆管和环空管流层流流动问题时，最关键的就是确定其速度分布。其中 最常用的方法就是运用基本方程(动量方程和相应流体的本构方程结合)求解。在这两种方程 中，把其连接在一起的是切应力。 但是文献[1．2，3]中介绍非牛顿流体本构方程时，都没有明确指出其中的切应力的方向， 仅仅表示出了其大小；而动量方程中切应力r，。是一个具有明确物理意义的矢量，具有方向性， 表现存在正负。在将两者统一起来的过程中容易产生问题。 3问题的解决 正因为本构方程中的切应力与动量方程中的切应力不统一，所以就有必要使其统一。这里 笔者提出了一种方法，该方法简单明了，可以使其一致，并能准确地求解其在圆管和环空中的 速度分布。 非牛顿流体本构方程虽然形式多种多样，但是可以统一的表示成牛顿流体的模式，其中仉 ·168· 称为表观黏度，即通过黏度计能够测出的黏度，肯定是一个正数。 r=仇)， 这样剪切应力的符号就应该与剪切速率的符号保持一致，也可以表达为，两者的方向是相同 的。 3．1 以宾汉流体为例 其本构方程‘13按照式(1)的形式可以写为 y r=仉y=(等机j 对于给定的塑性流体来说，，r。和叩，分别是正的常数。为了保证其表观黏度仉为一个正数， 则对仉中的y取绝对值。 即一(号机户j。=(贵饥P (2) 然后，按照圆管和环空中的剪切速率的实际情况去掉绝对值符号，得到r，：的准确形式。 在圆管中夕一誓o，r，：=一r。+巩夕； 在环空中，内俱。夕一害o，l—r。+7，7； 外倾。尹=警orr=一一r。+咿 3．2以卡森流体为例 其本构方程：33按照式(1)的形式可以写为 …卜协+6彬l 为了保证其表观黏度仉为～个正数，而且根号下的数要保持为正数，则对叩。中的7取绝 对值。 (4) 即r=【在+6卜净k=[悟+6卜 然后，按照圆管和环空中的剪切速率的y实际情况去掉绝对值符号，得到r。的准确形式。 +‘12夕一～(厄+√二丽一 在…，一鲁o，r一『压 ·169· 在环空中，内侧穸=誊o，r。= 6卜一(厄+、丽)2；