精密金属异型挤压塑性成形共形上限解及模腔优化

针对异型材精密挤压塑性成形及模腔建模理论的焦点课题，借助共形映射理论和金属塑性成形理论成果，将异型材截面域的金属挤压三维塑性流动转化为二维轴对称成形问题，建立了金属异型塑性流动的流函数、应变速度场和模具模腔等数学解析模型，应用能量上限原理，以含有对称轴叶形和六边形截面域的挤压型材模腔为例，求解了模腔优化参数，使优化模腔曲面的CAD／CAM一体化技术目标成为可能．     

基于注塑成型过程的液晶高分子Hele-Shaw流动模拟分析

针对注射模腔内液晶高分子Hele-Shaw流动进行模拟分析。模腔为规则的方形,温度为恒温300℃,以小分子液晶TIF理论为基础,建立内部流动分析方程组;采用这种方法预测整个注塑过程中任意时间点流动前沿的位置,随时掌握腔内压力、速度场的变化情况。通过模拟分析得知:在入射口处压力、速度场分布变化最大,向内部逐渐减小,且分布均匀,无畸变点;向矢场的变化则与模腔高度有关,靠近腔壁处分子向矢按照流动方向均匀排列,而中心层相对复杂,由X-Y平面内流动剪切速率决定;注射时间与初始压力、模腔高度的平方成反比。本文的结果为进一步分析注塑过程中非恒温下液晶高分子流动提供了参考。     

杯外径不相等的杯杯型轴对称复合挤压凸模单位压力

本文建立了杯外径不相等的杯杯型轴对称复合挤压的运动学容许速度场.在此基础上,用上限法分析了此类复合挤压的变形过程,确定了此类复合挤压变形不同阶段所需的上凸模平均单位压力.研究结果表明:作用于上凸模上的平均单位压力随着反挤变形程度的变化呈近似抛物线型;随着正挤变形程度的变化呈单调增函数.在挤压过程中,随着上凸模的逐渐压入,变形体的变形模式发生改变,即由早期阶段的稳态流动进入晚期阶段的非稳态流动.而正挤、反挤部分的变形程度组合是除剩余坯料相对厚度以外影响变形模式转换的主要因素.本文的理论分析结果可较为方便地运用于实际生产中的工艺设计.实验表明:理论分析结果与实际结果颇为一致.     

圆模拉拔挤压速度场的一种不同解法

文中对模面为圆的平面变形挤压拉拔设定了直角坐标系运动许可速度场,对此速度场采用曲线积分与变上限积分得到拉拔.挤压应力的上界解析解.     

全附体潜艇粘性流场的RANS模拟及场量和涡量的校验分析

采用相同拓扑结构和相近网格质量的4套网格和5种湍流模型,对全附体Suboff潜艇粘性流场进行RANS模拟,分析了网格密度、节点空间分布规律和湍流模型对计算精度的影响,详细校验了其力积分量、速度场量和涡量特征。结果表明:网格密度最大的G4网格(140万)计算精度最高,总阻力较实验值误差为0.723%,其采用SST湍流模型时最优。计算得到的压力系数和剪切应力系数分布均与实验值吻合很好;桨盘面速度等值线分布计算精度与文献相当,轴向相对速度0.9以上的计算半径稍大于实验值,其余半径与实验吻合较好;桨盘面上0.25倍半径处速度分量沿周向分布计算精度较文献高,轴向分量与实验值吻合较好,径向分量峰值稍小于实验值,但峰值所处周向位置与实验值一致。成功捕捉到了附体端面绕流诱导对旋涡、附体叶根截面下游处项链形涡对、尾翼端面尾缘上方附着涡蹄、附体马蹄涡系、尾翼截面通道流体挤压作用诱导涡以及桨盘面涡量汇集的潜艇涡量场特征,且围壳端面绕流诱导对旋涡沿流动方向持续稳定,不影响桨盘面涡量场,均与文献中由大涡模拟模拟得到的定性结论一致。研究表明,在网格密度较大、节点分布合理、网格质量较高、湍流模型选取适当和壁面函数使用有效的条件下,RANS模拟潜艇粘性流场的场量和涡量特征同样具有很高的计算精度,能够在工程应用中有力支撑新型艇型设计与性能分析。     

三维非牛顿熔体前沿界面的Level Set/Ghost/SIMPLEC模拟

建立了三维粘性不可压非牛顿流体流动的控制方程,采用Level Set/Ghost/SIMPLEC方法模拟了注塑成型充模阶段的三维流动过程;追踪到了不同时刻的熔体前沿界面,预测并分析了流动过程中不同时刻的压力、速度等重要的流动特征参数,并与牛顿流体相应的流动特征参数做了对比.研究结果表明:Lovol Set/Ghost/SIMPLEC方法可以准确追踪非牛顿熔体前沿界面:幂律熔体在流动过程中的压差明显大于牛顿熔体的压差,沿横截面的速度分布也有明显的差别.     

注塑模充模过程动态分析

注塑成型是利用型腔模制造理想制品主要的成型加工方式 ,塑料熔体的流动行为将直接影响着最终塑件的质量 ,塑料熔体在三维薄壁型腔内的流动属于带有运动边界的粘性不可压缩流体的流动 ,本文针对塑料注塑成型特点 ,经过量纲分析和引入合理而必要的假设 ,得到了适合于充模分析的数学模型。控制方程的求解主要包括三个阶段 :压力场、温度场和流动前沿位置的确定。数值求解采用有限元法求解压力场、有限差分法求解温度场、并利用控制体积法跟踪熔体前沿 ,实现了充模过程的动态模拟     

变形-损伤耦合的起塑性恒压轴对称充模胀形的有限元模拟

采用大变形刚粘塑性有限元法模拟超塑性恒压轴对称充模胀形过程、分析了模具几何参数及材料参数对胀形过程中材料的流变行为、胀形制许厚度分布和成形时间的影响规律。给出了质点的流动轨迹、不同时刻制件的剖面形状及应力、应变分布；基于修正的Gurson粘塑性势推导了内部空洞体积分数累积增大模型并据此进行了变形-损伤耦合计算．     

利用格子Boltzmann方法数值模拟Rayleigh-Benard对流

采用格子Boltzmann方法对较大Rayleigh数范围下的二维Rayleigh-Benard对流进行了模拟研究.引入能量分布函数,利用该能量分布函数与粒子速度分布函数耦合来求解一个热流场,能量分布函数与粒子速度分布函数和Boltzmann方程构成了一个新的双分布格子Boltzmann模型.在考虑密度随温度变化的情况下,进行数值模拟,得到了Rayleigh-Benard对流速度、温度随时间的变化规律、系统的流线和等温线分布及平均Nusselt数与Rayleigh数的之间的关系,与相关文献数据进行了对比,模拟结果非常吻合,证明了改进的双分布格子Boltzmann模型的有效性.     

冲击载荷作用下柱壳链中的弹性波传播简化模型及其解析解

柱壳链能引起波形的弥散,具备操控波形的潜力。建立了柱壳链结构的等效连续介质模型和细观有限元模型,研究了质量块冲击作用下柱壳链中的弹性应力波传播过程及其几何弥散特性。基于考虑横向惯性修正的Rayleigh-Love波动方程,建立了柱壳链在质量块冲击下的控制方程,采用Laplace变换及其逆变换获得了位移场、速度场和应变场的解析解,所得结果与细观有限元模拟结果较好吻合。结果表明,在冲击过程中应变和速度峰值均逐渐减小,应变峰值、振荡幅度和波形前沿宽度与泊松比和惯性半径相关,泊松比和惯性半径越大,应变峰值越小,应变分布振荡越剧烈,波形前沿宽度越宽。     

大学物理《力学》中的自由度分析

对《力学》中的物体自由度进行多方面分析,以深化教学、提高学生正 确分析物理问题的能力.使用实际教学分析的研究方法,在《力学》范围内讨论自由度与坐标、 自由与约束的关系并得以下结论: (1) 同一物体的自由度随其所在的``空间'不同而不同, 不因坐标系的选取不同而 异, 在同类参考系中不因参考系的动静而有别;(2)自由度遵循叠加原理. 讨论了质点系的总自由度及相关计算问题,并指出研究《力学》中自由度的意义.     

Development and design of new methods of measurement of state of strain of structures

The present paper deals with development and design of new methods utilizing Wiedemann's effect for determination of state of strain in building structures. Wiedemann's effect and some features of torsional strain of magnetic field are the basis of new experimental method. Especially the point electromagnetic strain gages using the effect of pure torsion of electromagnetic field to enable universal examination. For strain-gage measurements, almost all physical quantities are used which can be related to the variation in length of the structures. From the electric strain measurements, the most commonly used methods are the measurements by resonance-wire strain gages or by electric-resistance strain gages. In this paper, electromagnetic strain gages are discussed using the Wiedemann effect, and the author describes some new measuring equipment and his own suggestions and methods based on an application of this effect.     

Calculation of aerodynamic characteristics of arrays of profiles of arbitrary shape

It is well known that the problem on nonseparating potential flow of an incompressible fluid about an array of profiles reduces to an integral equation for a certain real function, determined on the contours of the profiles of the array. As such a function one can take, as was done, for instance, in [1–5], the relative velocity of the fluid on the profiles of the array. For arrays of profiles of arbitrary shape it is necessary to solve the corresponding integral equation numerically. In the particular examples of the calculation of aerodynamic arrays that are available [1–3] the numerical methods used were based on the approximate evaluation of contour integrals by rectangle formulas. As investigations showed, sizeable errors arose thereby in the approximate solution obtained, these being especially significant in the case of curved profiles of relatively small bulk. In the present paper a method for the numerical solution of the integral equation obtained in [5] is proposed. The method is based on the replacement of a profile of the array with an inscribed N polygon, the length of whose sides is of the order N^–1 and whose internal angles are close to . Convergence with increasing N of the numerical solution to an exact solution of the integral equations at the reference points is demonstrated. Examples of the calculation are given.Novosibirsk. Translated from Izvestiya Akademii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 2, pp. 105–112, March–April, 1972.