注塑成型充模过程的温度场计算

对注塑成型过程中充填阶段的非等温效应的数值模拟进行了研究 ,模拟了非牛顿流体在任意形状型腔内的流动。这里 ,数学模型是基于非牛顿流体在非等温状态下的广义 Hele-Shaw流动 ,利用有限元 /有限差分混合数值方法求解流场中的压力和温度分布     

非等温黏弹性复杂流动的改进SPH方法模拟

非等温黏弹性流体广泛存在于自然界和工业生产中,准确预测黏弹性流体的非等温流动机理和复杂流变特性有着重要的应用价值.文章提出一种改进的光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法对非等温黏弹性复杂流动进行了数值模拟,其中流体的黏弹特性通过eXtended Pom-Pom本构模型来表征.为了提高模拟结果的精度,采用了一种核函数梯度的修正算法;为了灵活地施加边界条件,发展了边界粒子和虚拟粒子相联合的边界处理方法;为了消除流动过程中的拉伸不稳定性,施加了粒子迁移技术.运用改进SPH方法数值模拟了液滴撞击固壁和F型腔注塑成型问题,通过与Basilisk软件得到的结果进行比较验证了改进SPH方法求解非等温黏弹性流体的有效性.通过利用不同粒子初始间距进行计算,评价了改进SPH方法的数值收敛性.研究了非等温流动相较于等温流动的不同流动特征,深入分析了不同热流变参数对流动过程的影响.数值结果表明,文章提出的改进SPH方法可稳定、准确地描述非等温黏弹性复杂流动的传热机理、复杂流变特性和自由面变化特性.     

渐开线齿轮传动非牛顿润滑介质的线弹流数值分析研究

采用适合各种流变模型的广义Reynolds方程，通过数值联立求解非牛顿介质的线弹流润滑基本方程组，获得了渐开线齿轮啮合过程的油膜压力、膜厚、表面剪应力分布，并分析了啮合过程中非牛顿效应对齿轮传动最小油膜厚度的影响。在数值计算方向引入延拓方法，使表面煎应力迭代具有大范围收敛性。     

HIGH-ORDER DISCONTINUOUS GALERKIN SOLUTION OF N-S EQUATIONS ON HYBRID MESH

针对层流NS方程发展了混合网格上的高阶间断有限元方法,给出了物面边界高阶近似的具体步骤以及近物面弯曲单元的处理方法。对数值离散产生的非线性方程组采用牛顿迭代进行求解,每个牛顿循环采用预处理广义最小余量法求解产生的大型稀疏线性系统。使用该方法得到了典型算例的数值结果,并跟前人的计算结果进行了比较。计算结果表明,混合网格上应用高阶间断有限元方法求解黏性流动具有很好的应用前景。     

注塑成型三维流动的数值模拟

建立了非等温、粘性、不可压缩、非牛顿流体流动的控制方程。为了避免同时求解耦合的压力场、速度场,本文通过修改Galerkin方法的变分方程,导出了关于压力场的拟Poisson方程,用迭代法独立地求解连续性方程、动量方程,并进行速度一粘度迭代求出最终的压力场、速度场。由于直接使用Galerkin方法求解能量方程容易引起温度场的振荡,本文采用隐式格式及“上风”法离散能量方程,用超松驰迭代法求解温度场的代数方程组。比较了模拟结果与等温管道流动的解析解及法兰的实际注射结果,算例表明本文方法可以预测注射成型流动过程中的一些重要特征。与传统Galerkin方法相比,本文方法可以减少内存,提高数值方法的稳定性。     

非牛顿体通用模型线接触弹流润滑的数值分析

基于润滑剂在弹流润滑状态下表现为非牛顿体特性,根据弹流润滑理论,采用一种新的非牛顿体流变模型,建立了适用于非牛顿体的修正Ｒｅｙｎｏｌｄｓ方程,进行了等温弹流润滑的数值计算,并在等温解的基础地温度场分析。数值分析结果表明,由于滑滚比和模型参数对剪应力影响较大,因而在滑滚比和模型参数较大时应进行热弹流计算。通过温度场分析可证明：非牛顿体通过模型可用于等温弹流润滑和热弹流润滑计算。     

注塑模充模过程的数值分析

利用非等温、粘性广义Hele－Shaw流动理论,建立了三维薄壁型腔塑料熔体充填过程的数学模型．耦合利用有限元／有限差分法数值求解压力和能量方程,实现充模过程的动态模拟．     

螺旋槽液体润滑轴承膜压力的算子分裂法计算

采用算子分裂法（Ｏｐｅｒａｔｏｒ－ｓｐｌｉｔｔｉｎｇ ｍｅｔｈｏｄ）求解满足ＪＦＯ空泡压力边界条件和全油膜质量连续性的广义雷诺方程。润滑油膜流动由剪切流动和压力差流动两部分组成。先采用算子分裂法的迎风差分法求解剪切流动分量;再采用质量集中的有限元法求解压力差流动。结果表明：由算了分裂法得到的一维滑块轴承数值解与基于Ｅｌｏｒｄ算法的结果一致。同时还计算了人字型螺旋槽液体润滑轴承的压力分布、承载力和偏     

二维定常不可压缩粘性流动N-S方程的数值流形方法

将流形方法应用于定常不可压缩粘性流动N-S方程的直接数值求解,建立基于Galerkin加权余量法的N-S方程数值流形格式,有限覆盖系统采用混合覆盖形式,即速度分量取1阶和压力取0阶多项式覆盖函数,非线性流形方程组采用直接线性化交替迭代方法和Nowton-Raphson迭代方法进行求解.将混合覆盖的四节点矩形流形单元用于阶梯流和方腔驱动流动的数值算例,以较少单元获得的数值解与经典数值解十分吻合.数值实验证明,流形方法是求解定常不可压缩粘性流动N-S方程有效的高精度数值方法.     

基于黏弹性的微注射成型流动模拟

应用有限元方法研究了微注射成型中瞬态、可压缩、非牛顿熔体流动的黏弹性对流动前沿及流动平衡的影响。基于Phan-Thien-Tanner模型建立了熔体流动的本构方程,利用Hele-Shaw假设和简化建立了瞬态、可压缩、非牛顿熔体流动的连续性方程、动量方程、能量方程;为了有效地描述微注射成型的尺寸效应,采用了边界滑移和表面张力边界条件。通过分部积分和待定系数法导出了带有边界信息的变分方程和求解应力分量的半解析公式,构造了有限元离散求解及超松驰迭代算法。模拟结果表明:熔体的黏弹性对浇口附近的压力和后续的熔体流动前沿有重要影响;与黏性模型相比,黏弹性模型可以控制模拟压力的快速增长,减少不同型腔之间的充填差异,与短射实验结果也更吻合。