专利名称：用于模拟塑性材料和气体的多相流的方法
技术领域：
本发明涉及用于模拟塑性材料和气体的多相流的计算机化方法，其有助于分析于腔室(Chamber)中混炼的例如未硫化橡胶、树脂等塑性材料的流动状态。
背景技术：
近年来，已经提出了各种计算机模拟法。例如，下述非专利文献I提出了用于计算在班伯里密炼机中混炼时在交联之前的塑性材料例如未硫化橡胶或树脂的流动状态的方法。[非专利文献I]“附聚物的分散混合的数值和实验研究”，V.Collin⑴'E.Peuvrel-Disdier(I)等人。非专利文献I的方法的前提是腔室完全被塑性材料充满(即填充率=100%)。然而，通过使用班伯里密炼机等难以使混炼在填充率为100%的条件下进行，并且已知随着填充率的改变，混炼的塑性材料的品质和状态也改变。尽管存在这些事实，但非专利文献I的方法不能通过将塑性材料的填充率设定为除100%以外的不同值来实施。的确，要使混炼在填充率100%的条件下进行几乎是不可能的，因此，非专利文献I的模拟结果不能通过实验来验证。另一方面，塑性材料的粘度很大程度上不同于空气的粘度，因此，如果通过流动计算(Flowcalculation)来模拟腔室中的塑性材料和空气两相,则计算倾向于变得不稳定。在通过利用空气的实际粘度进行流动计算的情况下，即使在不被打断的情况下完成计算，计算时间也变得非常长，并且需要高昂的计算成本。
发明内容
因此，研究出本发明来解决上述问题，并且本发明的目的在于提供一种用于模拟存在于腔室中的塑性材料和气体的多相流的计算机化方法，在所述方法中，可以在不大幅增加计算时间的条件下稳定地进行各种填充率下的多相流计算，并且由此可实现精确模拟。根据本发明，用于模拟在腔室中的塑性材料的流动状态的计算机化方法，包括:生成作为腔室的有限体积模型的腔室模型的步骤、定义塑性材料的材料模型的步骤、定义气体的气相模型的步骤、对在腔室模型中的材料模型和气相模型进行流动计算的步骤，其中，材料模型设定有塑性材料的粘度，并且气相模型设定有大于气体实际粘度的粘度。所述腔室可以是班伯里密炼机的混炼空间，在其中设置至少一个转子，并且塑性材料可以是待混炼的未硫化橡胶或树脂材料。在流动计算中，优选进行材料模型的剪切生热的计算，但是不进行气相模型的剪切生热的计算。在塑性材料模型和腔室模型之间的界面上，可对滑移边界条件进行定义。
图1是用于混炼塑性材料的班伯里密炼机一部分的示意性横截面图。图2是根据本发明的方法的一个例子的流程图。图3是腔室模型的透视图。图4是用于说明腔室模型的腔壁(Wall)和内部空间的啮合的横截面图。图5是单独显示腔室模型的功能部分的横截面图。图6是显示配置了材料模型和气相模型的腔室模型状态的横截面图。图7是流动计算程序的流程图。图8是显示通过根据本发明的模拟方法获得的塑性材料的混炼状态的图形。
具体实施例方式现在结合附图对本发明的实施方式加以详细说明。本发明的目的在于提供一种用于通过使用计算机(未图示)来评估设置于腔室中的塑性材料的流动状态的模拟方法。塑性材料，例如交联之前的未硫化橡胶、树脂或弹性体，具有粘度。然而，其并不限于这些材料。只要其具有稳定的流动状态，所述塑性材料可以是任意材料。在为交联之前未硫化橡胶的情况下，可将在80°C左右充分混炼的混合物的状态作为稳定的流动状态。所述腔室为基本上闭合的空间，在其中设置有塑性材料，并且所述塑性材在其中流动。所述腔室可以具有任意的构造(Configuration)。所述腔室的一个典型例子是班伯里密炼机的混炼空间。如图1所示，本实施方式中的腔室4形成于班伯里密炼机I的外壳2和设置于外壳2中的一对转子3之间，并且所述腔室4具有像数字8的截面形状。然而，所述腔室并不限于这样的构造。图2是根据本发明的模拟法的一个例子的流程图。*生成腔室模型的步骤Si。在该实施方式中，首先，计算机生成腔室模型5，其是腔壁和腔室4内部三维空间的有限体积模型。图3是腔室模型5的透视图。图4是其横截面图。腔室模型5具有由下述部分封闭而成的三维空间:由外壳2的内表面限定的外周面50、由两个旋转的转子3的外周面限定的内周面51、以及转子3沿轴向两侧的两个端面5s0该三维闭合空间被分(被离散)成三维单元(e)例如四面体单元、六面体单元和多面体单元。对于各单元(Element)，计算塑性材料(材料模型)的物理量例如压力、温度和/或速度。
所述外周面5ο和两个端面5s不被移动(Move)。然而，内周面5i根据转子3的旋转而移动，并且相应地，腔室模型5的构造被改变。在图5中显示的例子中，腔室模型5由四个功能部分组成:一对旋转部分5A及5B ；夹在它们之间的中间部分(In-between part) 5C ;以及包绕这些部分(5A、5B和5C)的外框架部分。各旋转部分5A/5B皆为筒状并且具有圆筒形外周面5Ao/5Bo以及对应于一个转子3外周面的内周面5i。将旋转部分5A和5B设置在外框架部分中，并且将它们限定为可环绕各中轴线Oa和Ob旋转，来显示由转子3的旋转引起的腔室体积构造的变化。与此相反，所述中间部分5C保持固定在旋转部分5A和5B之间，并且具有紧靠于各旋转部分5A和5B的两个凹面j。在所述凹面j和各圆筒形外周面5Ao和5Bo上，定义了作为滑移面的边界条件。这使得发生在旋转部分5A和5B之中的物理作用(力、热等)经过凹面j被传递至存在于中间部分5C之中的材料模型。外框架部分为筒状并且包绕旋转部分5A、5B和中间部分5C。两个轴末端均被两个端面5s封闭。在外框架部分和旋转部分5A与5B之间的界面上，以及在外框架部分和中间部分5C之间的界面上，定义了作为滑移面的边界条件。这使得产生于旋转部分5A和5B中的物理作用(力、热等)穿过它们之间的界面被传递至外框架部件外框架部分因转子的运行而受到相对较大的剪切力。因此，为了更详细地计算材料的速度等，将组成外框架部分的元件的尺寸制得比旋转部分5A和5B以及中间部分5C的尺寸更小。由此，可更详细地计算靠近腔室模型5内表面的材料模型的速度分布图等。定义材料模型的步骤s2*接着，通过计算机定义并且储存材料模型。材料模型是在腔室4中流动或移动的塑性材料的模型。在材料模型上，对塑性材料的物理特性例如剪切粘度、比热、热导率和比重进行定义并且储存于计算机中。如上所述，本实施方式中的塑性材料是处于塑化状态的橡胶混合物，因此，预先输入并且储存这样的橡胶混合物的物理特性。[剪切粘度]为了得到剪切粘度，首先，在多个温度条件下测定塑性材料的粘弹性(G^和G")。然后，根据Cox-Merz规则，将测得的粘弹性换算成剪切粘度。根据幂次法则，如上得到的剪切粘度可通过下述式(I)近似得到。式(I):η =m Y ' n_1其中，n:剪切粘度、 m:绝对温度的函数的系数、
Y /:剪切速度、以及η:系数。分析对象(塑性材料)的比热可以通过绝热连续加热法(以25 °C )得到，并且预先将测得的比热值输入并储存在电脑中。分析对象(塑性材料)的热导率可以通过热线法(以25°C)得到，并且预先将测得的值输入并储存在电脑中。定义气相模型的步骤S3*接着，通过计算机定义并且储存气相模型。所述气相模型是存在于腔室中的气体的模型。在本发明中，在腔室中塑性模型的填充率小于100%，因此，为了实现流动计算，未充满塑性模型的部分用气相模型充满。本实施例中的气相模型设定了粘度和比重。将气相模型的比重设定为气体比重的实际值。然而，通过计算机将气相模式的粘度设定为大于气体粘度实际值的值。当腔室中的塑性材料的填充率小于100%时，有必要对塑性材料和气体的多相流进行流动计算。 通常，气体和塑性材料之间的粘度差别非常大。因此，如果将因为速度引起的剪切生热纳入考虑来进行计算，则在气相模型和材料模型之间的界面上剪切生热会超常增加，并且结果是计算容易变得不稳定。根据本发明，因此，为了使计算稳定，在不会给计算结果带来负面影响的范围内，尽可能地提高气相模型的粘度。结果是，有可能模拟中增加时间间隔(时间增量)来加快计算速度。本发明的发明人仅改变气相模型的粘度，进行模拟来比较腔室模型5中的压力场。结果，发现如果气相模型的粘度超过实际粘度的10倍，则腔室模型5中的压力开始增加，并且给压力场带来负面影响。然而，如果气相模型的粘度低于实际粘度的5倍，则变得难以使计算稳定。因此，优选将气相模型的粘度设定在气体实际粘度的5 10倍的范围内。定义边界条件的步骤s4*接着，对流动计算所需的包括边界条件在内的各种条件进行定义。在腔室模型5的表面上，边界条件可能包括流速边界条件和温度边界条件。关于流速边界条件，可根据模拟的目的、所需的精度等对(a)壁表面无滑移条件，或(b )壁表面滑移条件进行定义。在壁表面非滑移条件的情况下，在腔室模型5的表面上的材料模型的流速始终是零。在壁表面滑移条件的情况下，在腔室模型5的表面上的塑性模型的流速可具有非零值、也可具有零值。在这种情况下，可将下述式(2) (Navier/ sLaw)应用于在材料模型和腔室模型5之间的界面上的滑移现象。
式(2):fs = -Fslip (vwall-vs) I Vs-Vwall IesllIri式(2)，:τ ω = Fslip [vslip-vwall]其中fs:在腔室模型5的表面上的材料模型的剪应力、Vslip:在腔室模型5的表面上的材料模型的滑移速度(切向速度)、Vwall:在腔室模型5的表面上的材料模型的切向速度、Fslip, eslip:材料特性(如果eslip=l，则为线性；如果0〈eslip〈l，则为幂次法则)。为了定义在腔室模型5的表面上的材料模型的滑移速度vslip，首先，得到初始阶段的速度场，然后由该速度场(或修正后的速度场)得到在腔室模型5的表面上的材料模型的剪切力，然后通过上述式(2)得到在腔室模型5的表面上的材料模型的滑移速度。在此情况下，如果需要，为了控制滑移速度的突然增加，优选通过将滑移速度乘以小于I的松弛系数，在常规的SIMPLE算法的循环之中使腔室模型5各表面的壁表面滑移速度逐渐收敛于真实值。 (见后述)关于温度边界条件，可定义为:(a)其中热量不会经过腔室模型5的表面从腔室模型5逃逸至外界的绝热条件，或(b)其中整个腔室模型5的表面皆具有恒温(例如50°C)的条件。进一步地，所述条件可能包括塑性模型的初始温度(在本实施方式中为20°C(293K))、由旋转部分5A和5B的转数代表的转子的转数、腔室模型5的表面的滑移比、塑性模型相对于腔室模型5的体积的填充率等。此外，所述条件可以包括流动计算的初始状态、用于计算的时间间隔、内部处理中的迭代数、计算的最大周期等。关于初始状态，例如参见图6中所示，可定义为气相模型的区域A，所述区域A在水平界面s的上侧，所述水平界面被定义为延伸穿过腔室模型5，以及液体模型的区域M，所述区域M在水平界面s的下侧。因此，通过改变界面s的水平面，可调节塑性材料(材料模型)的填充率。可根据模拟的目的任意地定义这些条件。进行流动计算的步骤s5*接着，在腔室模型5中设置材料模型和气相模型，并且基于上述条件进行流动计
笪
ο在流动计算中，计算限定材料模型运动状态的沿三个坐标轴方向(X、y、z)的分速度、以及限定材料模型内部状态的温度T和压力p。因此，待测的未知量为这5个物理量。本实施方式中的流动计算是通过使用用于不可压缩性流动的纳维尔斯托克斯方程来进行的，将各个材料模型和气相模型的密度设定为恒定值。在本发明中，在整个可能的温度范围内，将材料模型视作液体。因此，纳维尔斯托克斯方程，即质量守恒和能量方程的公式，必须作为流体的方程组来求解。
在本实施方式中，为了处理气体和塑性材料的多相流，采用用于计算自由界面的流动的VOF (流体体积)法。VOF法不直接计算两种流体之间的界面的运动。在VOF法中，自由界面通过对在各单元的体积范围内的各流体的填充率(体积分数)进行定义而表示。支配方程如下。[运动方程]在本实施方式中，将在腔室模型5中气相模型和材料模型汇流的两相流视为单相流。在此情况下，待求解的运动方程为沿三个坐标轴方向x、y和z的下述式(3)。通过VOF法对所述双相进行平均，并将其视为单相，结果是上述运动方程的求解变得可能。式(3):
权利要求
1.一种用于模拟腔室中塑性材料的流动状态的计算机化方法，包括:生成作为腔室的有限体积模型的腔室模型的步骤，定义塑性材料的材料模型的步骤、定义气体的气相模型的步骤、对处于腔室模型中的材料模型和气相模型进行流动计算的步骤，其中，所述材料模型设定有塑性材料的粘度，并且所述气相模型设定有比气体实际粘度更大的粘度。
2.如权利要求1的计算机化方法，其特征在于，腔室是设置有至少一个转子的班伯里密炼机的混炼空间，并且塑性材料是待混炼的未硫化橡胶或树脂材料。
3.如权利要求1或2的计算机化方法，其特征在于，在流动计算中，进行材料模型的剪切生热计算，但是不进行气相模型的剪切生热计算。
4.如权利要求1、2或3的计算机化方法，其特征在于，在塑性材料模型和腔室模型之间的界面上定义滑移边界条件。
全文摘要
一种用于模拟腔室中的塑性材料的流动状态的计算机化方法，其包括下述步骤生成作为腔室的有限体积模型的腔室模型；定义塑性材料的材料模型；定义气体的气相模型；以及对腔室模型中的材料模型和气相模型进行流动计算，其中，所述材料模型设定有塑性材料的粘度、并且所述气相模型设定有比气体实际粘度更大的粘度。
文档编号G06F17/50GK103218471SQ20131001850
公开日2013年7月24日 申请日期2013年1月18日 优先权日2012年1月23日
发明者角田昌也 申请人:住友橡胶工业株式会社