拉延成形多目标序列响应面法优化设计方法

为了克服常规响应面法在整个设计空间进行逼近导致精度低和传统的单目标优化设计,只能针对其中的一个目标进行优化的缺陷,提出了一种多目标序列响应面法,优化设计拉延筋几何参数的方法. 该方法通过移动、缩放等方式在设计空间中不断更新兴趣域,在不同的兴趣域中将实验设计、响应面法和多目标粒子群优化算法相结合,获得了一组最小化起皱、拉裂缺陷的等效拉延阻力非劣解. 利用最小距离选解法从非劣解集中挑选出一组成形效果最好的解, 并以此解作为下一迭代步兴趣域的中心,直到收敛至一组最优的等效拉延约束阻力. 以最优等效拉延阻力为约束条件,利用等效拉延筋阻力模型结合遗传算法对拉延筋的真实几何参数进行设计. 优化设计的拉延筋几何参数提高了板料的成形性能. 算例表明,该方法具有较高的精度和较强的工程实用性.      

数值模拟界面流方法进展

用数值模拟方法研究界面流,对于工程设计和理论分析复杂流动都是非常重要的.单纯用现有的高分辨格式无法清晰定出流体界面,故此对界面必须进行特殊的处理,例如界面跟踪(front tracking)方法．本文对界面流的各种数值模拟方法进行了综述．这些方法包括锋面跟踪法、移动网格法、基於粒子法、边界积方法、连续对流格式、体积跟踪方法和水平集的方法．我们着重介绍了现阶段处理流体界面问题使用较多的方法:体积跟踪方法和水平集方法．同时我们也对各种方法的优缺点作了评述.      

两相流PIV粒子图像处理方法的研究

本文在单相PIV技术的基础上研究了两相流动PIV图像处理方法，采用摸板匹配法和灰度加权标定法对两相粒子进行了识别、区分和标定，采用灰度互相关法对区分后的单相粒子图像进行了处理，应用基于以上方法编制的Windows应用软件，首先对由美国Minnesota大学复杂流动实验室提供的两相流动粒子图片进行了处理，通过对比分析可见，应用本文所采用的方法能对两相粒子进行有效的识别和区分，然后以搅拌槽内液固两相流场为例对此方法进行了应用。     

两相流场粒子成像测速技术(PTV-PIV)初探

在单相ＰＴＶ－ＰＩＶ研究基础上，研究液固两相流场粒子成像图像数字测速技术．首次用ＰＴＶ技术得到直槽道中液固两相流场二维瞬时全场两相流速分布．初步探讨了两相ＰＴＶ－ＰＩＶ技术中与单相ＰＴＶ－ＰＩＶ测速技术的关键不同之处，提出两相ＰＴＶ－ＰＩＶ测速采样的两相相容性准则．为深入研究两相流场粒子成像测速技术奠定了基础．     

关于等效刚度法适用条件的讨论

关于等效刚度法适用条件的讨论黄义仿，贾汝民（兰州炼油厂职工大学，兰州７３００６０）本文利用势能驻值原理对文［１］中提出的方法前提条件的充分性和必要性进行了讨论，认为等效刚度法并不普遍适用．设文［１］中图１（ａ）和图１（ｂ）所示的两压杆处于临界状态．按...     

不可压黏性流问题的无网格法分析

不可压缩黏性流问题一般采用Navier-Stokes方程来描述，基于加权残值法，推导了问题的无网格伽辽金法（EFGM）离散Navier-Stokes方程，在时间域上采用分步方法计算，速度和压力由相互独立的方程以解耦的形式求解，并采用同阶移动最小二乘近似，在每一时间步中，对压力解和速度解采用了Newton-Raphson迭代法进行修正，最后将所得到的方法应用到剪切驱动空腔流问题中，验证了方法的有效性，且解的精度高、稳定性好。     

管壳式换热器内气液两相流冲刷垂直管速的阻力特性研究

本文对垂直布置的管壳式换热器中，气液两相液冲刷管束时的流动阻力进行了实验研究。实验段采用管壳式换热器模型，其圆形外壳内的管束由４９根管子组成，并沿长度方向用三块折流板将管束分成四个冲刷流程，当气液两相流自下而上冲刷管束时，分别测量摩擦阻力和局部阻力，并用分相流动模型进行分析，得到良好的计算关联式。     

对估算金属切口强度的能量法与Neuber法的比较

估算金属切口的强度时,人们一般采用以Neuber准则为理论基础的Neuber法,或以Moski及Glinka的等效能量密度法为理论基础的能量法．本文以弹塑性幂硬化材料为研究对象,将两种方法得到的结果与实验结果进行了比较,发现Neuber法总是低估切口强度,而大多数情况下,能量法却高估切口强度,并且总的来说,能量法的精度高于Neuber法的精度;能量法的可靠性高于Neuber法的可靠性．     

多自由度Duffing系统受多相位平稳随机激励的响应

本文用虚拟激励法结合等效线性法来求解受多相位平稳随机激励的多自由度Ｄｕｆｆｉｎｇ系统，过程简洁，计算高效，给出了在该激励下多自由度Ｄｕｆｆｉｎｇ系统的数值模拟算法，二者进行了比较，得到了比较满意的结果。     

一类电化学反应的差分解法

对一类电化学反应中的粒子扩散问题提出了一种差分解法。在一个实际例子的数值计算中，得到了电极表面的特质浓度降为零所需要的时间。     

Creeping motion of a slip spherical particle in a circular cylindrical pore

A combined analytical–numerical study for the creeping flow caused by a spherical fluid or solid particle with a slip-flow surface translating in a viscous fluid along the centerline of a circular cylindrical pore is presented. To solve the axisymmetric Stokes equations for the fluid velocity field, a general solution is constructed from the superposition of the fundamental solutions in both cylindrical and spherical coordinate systems. The boundary conditions are enforced first at the pore wall by the Fourier transforms and then on the particle surface by a collocation technique. Numerical results for the hydrodynamic drag force acting on the particle are obtained with good convergence for various values of the relative viscosity or slip coefficient of the particle, the slip parameter of the pore wall, and the ratio of radii of the particle and pore. For the motion of a fluid sphere along the axis of a cylindrical pore, our drag results are in good agreement with the available solutions in the literature. As expected, the boundary-corrected drag force for all cases is a monotonic increasing function of the ratio of particle-to-pore radii, and approaches infinity in the limit. Except for the case that the cylindrical pore is hardly slip and the value of the ratio of particle-to-pore radii is close to unity, the drag force exerted on the particle increases monotonically with an increase in its relative viscosity or with a decrease in its slip coefficient for a constant ratio of radii. In a comparison for the pore shape effect on the axial translation of a slip sphere, it is found that the particle in a circular cylindrical pore in general acquires a lower hydrodynamic drag than in a spherical cavity, but this trend can be reversed for the case of highly slippery particles and pore walls.     

Slow motions of a circular cylinder experiencing slip near a plane wall

A theoretical study is presented for the two-dimensional creeping flow caused by a long circular cylindrical particle translating and rotating in a viscous fluid near a large plane wall parallel to its axis. The fluid is allowed to slip at the surface of the particle. The Stokes equations for the fluid velocity field are solved in the quasi-steady limit using cylindrical bipolar coordinates. Semi-analytical solutions for the drag force and torque acting on the particle by the fluid are obtained for various values of the slip coefficient associated with the particle surface and of the relative separation distance between the particle and the wall. The results indicate that the translation and rotation of the confined cylinder are not coupled with each other. For the motion of a no-slip cylinder near a plane wall, our hydrodynamic drag force and torque results reduce to the closed-form solutions available in the literature. The boundary-corrected drag force and torque acting on the particle decrease with an increase in the slip coefficient for an otherwise specified condition. The plane wall exerts the greatest drag on the particle when its migration occurs normal to it, and the least in the case of motion parallel to it. The enhancement in the hydrodynamic drag force and torque on a translating and rotating particle caused by a nearby plane wall is much more significant for a cylinder than for a sphere.     

液浮转子微陀螺筒间流动阻力矩参数影响

建立简化计算模型,将一种具有空间指数收敛特性的谱元-Fourier法引入复杂的液浮转子微陀螺筒间流动的直接数值模拟。对于柱坐标系下不可压缩Navier-Stokes方程,在子午面的径向和轴向采用Galerkin谱元法,进行标准化的Gauss-Lobatto-Legendre节点展开,在周向采用Fourier级数展开法。设定基准算例,计算分析了高径比、旋转速度和黏性系数等参数对内筒侧面和内筒端面阻力矩的影响。结果表明,两种阻力矩都随这些参数的增大而增大,且对总阻力矩的贡献量级相当。     

Evaluation of the drag force by integrating the energy dissipation rate in stokes flow for 2D domains using the FEM

The total drag force on the surface of a body, which is the sum of the form drag and the skin friction drag in a 2D domain, is numerically evaluated by integrating the energy dissipation rate in the whole domain for an incompressible Stokes fluid. The finite element method is used to calculate both the energy dissipation rate in the whole domain as well as the drag on the boundary of the body. The evaluation of the drag and the energy dissipation rate are post-processing operations which are carried out after the velocity field and the pressure field for the flow over a particular profile have been obtained. The results obtained for the flow over three different but constant area profiles—a circle, an ellipse and a cross-section of a prolate spheroid—with uniform inlet velocity are presented and it is shown that the total drag force times the velocity is equal to the total energy dissipation rate in the entire finite flow domain. Hence, by calculating the energy dissipation rate in the domain with unit velocity specified at the far-field boundary enclosing the domain, the drag force on the boundary of the body can be obtained.     

Slow Motion of a Porous Cylindrical Shell in a concentric cylindrical cavity

This paper concerns the Slow Motion of a Porous Cylindrical Shell in a concentric cylindrical cavity using particle-in-cell method. The Brinkman’s equation in the porous region and the Stokes equation for clear fluid in their stream function formulations are used. The hydrodynamic drag force acting on each porous cylindrical particle in a cell and permeability of membrane built up by cylindrical particles with a porous shell are evaluated. Four known boundary conditions on the hypothetical surface are considered and compared: Happel’s, Kuwabara’s, Kvashnin’s and Cunningham’s (Mehta-Morse’s condition). Some previous results for hydrodynamic drag force and dimensionless hydrodynamic permeability have been verified. Variation of the drag coefficient and dimensionless hydrodynamic permeability with permeability parameter σ, particle volume fraction γ has been studied and some new results are reported. The flow patterns through the regions have been analyzed by stream lines. Effect of particle volume fraction γ and permeability parameter σ on flow pattern is also discussed. In our opinion, these results will have significant contributions in studying, Stokes flow through cylindrical swarms.     

Effect of heated underwater boundary layer of a prolate spheroid on viscous drag

 A prolate spheroid submerged in water can be heated for decreasing the viscous drag because of a decrease in viscosity with increasing temperature. The heated boundary layer experiment, based on this principle, was carried out in a circulating water channel of NCKU, and the viscous drag of a spheroid, with a five-to-one ratio of length to mid-diameter, was also measured by means of wake surveys. The difference of total-head between wake and undisturbed region was measured by two total-head tubes, and the water speed in the wake area was measured by an electromagnetic flowmeter. The results indicate that the viscous drag of the model decreases with surface heating, and the decrease in viscous drag of the model is 13% when the surface temperature is 17 °C above the ambient water temperature. The velocity gradient in the wake region and the total-heating differential readings between undisturbed and wake region also decrease with surface heating due to the delayed laminar-turbulent transition in the boundary layer. The values of the coefficient of viscous drag obtained without surface overheating are found to be in agreement, for the range of Froude numbers investigated, with the results obtained from previous experiments. Received: 18 June 1996/Accepted: 29 April 1997     

Slow Motion of a Porous Sphere Translating Along the Axis of a Circular Cylindrical Pore Subject to a Stress Jump Condition

The coupled flow problem of an incompressible axisymmetrical quasisteady motion of a porous sphere translating in a viscous fluid along the axis of a circular cylindrical pore is discussed using a combined analytical–numerical technique. At the fluid–porous interface, the stress jump boundary condition for the tangential stress along with continuity of normal stress and velocity components are employed. The flow through the porous particle is governed by the Brinkman model and the flow in the outside porous region is governed by Stokes equations. A general solution for the field equations in the clear region is constructed from the superposition of the fundamental solutions in both cylindrical and spherical coordinate systems. The boundary conditions are satisfied first at the cylindrical pore wall by the Fourier transforms and then on the surface of the porous particle by a collocation method. The collocation solutions for the normalized hydrodynamic drag force exerted by the clear fluid on the porous particle is calculated with good convergence for various values of the ratio of radii of the porous sphere and pore, the stress jump coefficient, and a coefficient that is proportional to the permeability. The shape effect of the cylindrical pore on the axial translation of the porous sphere is compared with that of the particle in a spherical cavity; it found that the porous particle in a circular cylindrical pore in general attains a lower hydrodynamic drag than in a spherical envelope.     

超疏水沟槽表面通气减阻实验研究

减阻是解决航行体提速和增程的主要技术途径之一, 对缓解日益严峻的能源危机极为重要. 在重力式管道实验系统中, 测试给出了湍流状态下不同通气速率时减阻率随雷诺数及沟槽无量纲间距的变化规律和气膜铺展状态, 对比分析了单纯超疏水表面与超疏水沟槽表面上通气时减阻效果的差异.实验板材质为无色亚克力, 沟槽结构采用机械方法加工, 并在表面喷涂超疏水涂层. 结果表明, 持续通气能解决超疏水沟槽表面气膜层流失问题, 实现气膜层长时间稳定维持; 恒定雷诺数下, 随通气速率增大, 超疏水沟槽表面气膜铺展更趋均匀, 减阻率上升; 由于通气速率影响气膜横向扩展能力, 致使恒定通气速率下, 减阻率随雷诺数的变化呈现两种模式; 在固定雷诺数及通气速率时, 减阻率随沟槽尺寸的扩大先增后减, $S^{+}\approx 76$时减阻率最大. 分析其原因在于, 沟槽结构增大沾湿面积的同时, 显著提升了通气状态下超疏水表面气膜层的稳定性, 因而展示出与超疏水表面和沟槽表面均不相同的减阻规律, 且效果更佳.      

Minimization of the drag of blunt edges of projecting elements of a hypersonic object

An approximate method is given for calculating the pressure on the critical line of a cylindrical edge of an element projecting above the surface of a hypersonic object. Correlations of the pressure distribution in the interference regions are proposed on the basis of the calculations and analysis of the experimental data. Variational problems of finding the form of a leading edge having minimum drag are set up and solved on the basis of the pressure calculation model obtained. The solutions are sought with the use of Pontryagin's maximum principle. The effect on the minimum drag profile of Mach numbers and blunting of the body is analyzed. The advantage of the optimal form of the edge in comparison with a corresponding rectilinear edge is shown.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 3–9, May–June, 1976.     

基于微纳米压入法提取无铅焊料合金弹塑性力学参数的研究

本文通过微纳米压入法结合数值模拟研究了无铅焊料合金SnAg3.5 的弹塑性力学性能，分别采用圆柱形压头及两种不同锥角压头对无铅焊料合金进行压入测试：基于圆柱形压头测试过程中接触面积恒定的特点得到了无铅焊料的弹性模量，进一步采用塑性应变梯度理论对两种锥角压头的测试结果予以修正并通过数值模拟反分析得到相应的特征应力值，同时基于压入特征塑性应变与压头锥角的关系式得到两种不同锥角压头下的特征应变值，在此基础上经求解方程组得到焊料合金的初始屈服应力与应变强化因子，进而得到了焊料合金的幂强化弹塑性本构关系．该方法剔除了压入尺度效应的影响并保证了所得本构关系的唯一性，给出了一种通过原位压入测试表征金属材料弹塑性力学性能的有效方法．