基于微分求积法求解非粘滞阻尼系统的时程响应

研究一种含有指数型非粘滞阻尼线性多自由度振动系统的时程分析问题。该非粘滞阻尼模型假设阻尼力与质点速度的时间历程相关,数学表述为质点速度与核函数的卷积。由于阻尼模型的改变,常用的数值积分方法（如Newmark-β法、Wilson-θ法）不能直接应用于这种非粘滞阻尼系统。基于一种无条件稳定的微分求积方法,给出了这种非粘滞阻...     

柔性梁斜碰撞问题的非线性动力学建模和实验研究

本文研究柔性梁点面斜碰撞问题。用Hertz接触模型处理法向撞击力,分别用Hertz切向接触模型和Coulomb摩擦力模型处理粘滞状态和滑动状态的摩擦力。从精确的应变与位移的关系出发,用绝对节点坐标法建立了柔性梁的动力学方程。为了准确地处理斜碰撞切向运动的复杂状态,提出滑动-粘滞切换的准则,在此基础上,设计了斜碰撞实验,数值对比了法向撞击力和法向速度的时间历程的仿真计算结果与实验结果,验证了Hertz理论在斜碰撞情况下的正确性。另一方面切向速度的实验与理论的结果对照表明滑动-粘滞切换准则的有效性。     

粘滞阻尼器减震结构非线性随机振动的时域显式降维迭代随机模拟法

粘滞阻尼器在大型复杂结构减震设计中应用广泛。由于粘滞阻尼器的非线性阻尼力特性,粘滞阻尼器减震结构非平稳随机地震反应分析是一个典型的局部非线性随机振动问题。利用减震结构动力响应时域显式表达式的降维列式优势,仅针对与粘滞阻尼器相关的局部自由度进行非线性迭代计算,提出了局部非线性随机振动问题的时域显式降维迭代随机模拟法,为设置粘滞阻尼器的大型复杂减震结构非线性地震反应分析提供一种高效的随机振动方法。以安装了四个纵桥向粘滞阻尼器的某主跨1200m悬索桥为工程实例,开展E2水准地震激励下的非线性随机振动分析。计算结果显示,设置阻尼器后,主梁的纵桥向位移得到明显控制,降幅达到80%,大桥的关键截面内力也有5%左右的降幅。     

用自制仪器精确测定变温液体的粘滞系数

介绍了用恒流法自制仪器精确测定液体的粘滞系数.该自制仪器,能有效地保证底部的密封性、液体水位恒定和可操作性,还增设一个玻璃管连通器,附加一根标尺和游标使液体水位的测量精确度提高.并在接毛细管输出液体的量筒处加上光电计使测量液体流速的精确度提高,采用该仪器测定了水和乙醇两种液体在变温条件下的粘滞系数,经实验比对,所测得结果与理论值能较好地吻合.同时设计了可调换不同管径的毛细管,能快速准确地测量多种变温液体的粘滞系数.     

用有限单元法研究金属粘塑性稳态流动问题

七十年代以来,用有限单元法研究金属塑性变形的方法有了很大发展,这种方法对那些目前只能用试验方法和工程近似方法解决的问题,提供了完全用计算机方法解决的途径。本文采用刚粘塑性有限单元法,求解薄板轧制时金属塑性流动问题。这一方法忽略金属的弹性性质。因此可以不考虑金属的弹性行为,并且不用考虑金属从弹性状态到塑性状态的过渡,这使求解过程大大简化。对于金属塑性加工问题,这样的近似是允许的。这一方法采用了粘塑性材料的本构关系,所以它不仅能够分析理想塑性材料和变形强化材料,而且可以分析包括变形强化性质的变形速度强化材料,这可以使塑性变形力学分析更接近真实情况。     

高速铁路连续梁桥的共振响应与振动控制研究

以液体粘滞阻尼器为振动控制外部装置,主桥采用欧拉伯努利梁,通过中国和谐号动车组CRH380AL、日本新干线Shinkansen700和欧洲高速列车HSLMA8三种不同类型的高速列车对比,模拟分析了高速列车作用下桥梁结构共振响应的影响因素,以及粘滞阻尼器的阻尼系数与安装位置对桥梁结构振动响应的减振效果。研究结果表明,(1)合理有效地布置列车荷载轴距,可使桥梁结构发生基频共振的列车时速在运营时速之外,避免了桥梁结构发生较大振动峰值响应,即桥梁结构的基频共振;(2)随着粘滞阻尼器阻尼系数的增大,桥梁的加速度峰值在列车不同时速下均在减少,对桥梁振动有着不同程度的减振效果;(3)通过合理安置液体粘滞阻尼器,可有效降低高速列车作用下桥梁结构的共振响应;(4)随着粘滞阻尼器与主梁的连接点位置逐渐远离支座,粘滞阻尼器的减振效果逐渐明显;(5)随着粘滞阻尼器与桥台的连接点位置逐渐靠近支座,粘滞阻尼器的减振效果略有提升,但不明显。     

黄土的流变特性模拟与研究

工程建筑活动与土体的应力、变形、强度具有紧密的相关性,使得对土体的变形、破坏规律的研究显得十分必要。在对某现场黄土土样剪切蠕变试验的基础上,本文以分数线性蠕变方程为理论模型,模拟了该土体的蠕变、强度和流变特性,得到了该土体的剪切模量G随时间t和随法向应力σn的变化规律。通过流变分析,得到了粘滞系数η对法向应力σn的变化率受剪应力τ影响的规律,即当剪应力值较小时,粘滞系数η随法向应力σn变化率较大,而当剪应力增大后,粘滞系数η随法向应力σn变化率变小。还得出了加载时间t对粘滞系数η随剪应力τ变化的影响规律,即每级加载一定时间后,粘滞系数η随剪应力τ变化规律,以及时间t对粘滞系数η变化率的影响。同时还从土的微观结构的角度进行了分析讨论。     

用边界积分方程法解具有自由边界的粘滞不可压缩流问题

文中提出用边界积分方程法解具有自由边界的粘滞不可压缩流问题。Navier—Stokes方程中的非线性项是通过线性化近似求解的。作为对这种方法的检验,计算了失重情况下容器内粘滞流体因表面张力变化(由温差引起)而产生的旋流,并将计算结果与用有限差分法所得结果作了比较。此外,还初步计算了在管内压力作用下液滴的形成过程。     

高速铁路连续梁桥的共振响应与振动控制研究

以液体粘滞阻尼器为振动控制外部装置，主桥采用欧拉伯努利梁，通过中国和谐号动车组CRH380AL、日本新干线Shinkansen700和欧洲高速列车HSLMA8三种不同类型的高速列车对比，模拟分析了高速列车作用下桥梁结构共振响应的影响因素，以及粘滞阻尼器的阻尼系数与安装位置对桥梁结构振动响应的减振效果。研究结果表明，（1）合理有效地布置列车荷载轴距，可使桥梁结构发生基频共振的列车时速在运营时速之外，避免了桥梁结构发生较大振动峰值响应，即桥梁结构的基频共振；（2）随着粘滞阻尼器阻尼系数的增大，桥梁的加速度峰值在列车不同时速下均在减少，对桥梁振动有着不同程度的减振效果；（3）通过合理安置液体粘滞阻尼器，可有效降低高速列车作用下桥梁结构的共振响应；（4）随着粘滞阻尼器与主梁的连接点位置逐渐远离支座，粘滞阻尼器的减振效果逐渐明显；（5）随着粘滞阻尼器与桥台的连接点位置逐渐靠近支座，粘滞阻尼器的减振效果略有提升，但不明显。     

基于Volterra积分方程Taylor展开的卷积型非粘滞阻尼系统响应的时程分析方法

考虑一个具有卷积型非粘滞阻尼特性的多自由度系统响应的时程分析问题。非粘滞阻尼模型假设阻尼力与质点速度的时间历程相关,数学表达式为阻尼力等于质点速度与某一核函数的卷积,该模型为常用的粘滞阻尼模型的一般化形式。以一种在特定区间内求解第二类Volterra方程的Taylor展开法为基础,对所分析时段中各时间点的响应函数逐步作Taylor展开,代入卷积核来消去运动方程的积分项,通过求解推导出的时变线性方程组完成对卷积型阻尼模型系统的时程响应分析。数值算例验证了该方法的有效性。该方法增大时间步长可以大幅减少计算量,计算精度有所下降。     

A domain perturbation study of steady flow in a cone-and-plate rheometer of non-ideal geometry

The method of domain perturbation developed by Joseph is used to calculate velocity and stress profiles in a slightly misaligned cone-and-plate rheometer where the cone is spinning and the plate is stationary. Results for a Newtonian fluid, a Criminale-Ericksen-Filbey fluid, an upper-convected Maxwell fluid, and a White-Metzner fluid are presented and compared with earlier results in which the cone is stationary and the plate is spinning (Dudgeon and Wedgewood, 1993). Streamlines calculated for the Newtonian fluid show a very small recirculation region near the stationary plate. Velocity and stress contours are symmetric around the plane of largest gap width. For the elastic fluids studied, streamlines are asymmetric. The fluid response lags where the fluid is dominated by memory effects. Much larger recirculation regions are calculated for fluids dominated by shear thinning. These recirculation regions contain a large fraction of the fluid in the apparatus and have the effect of changing the shape of the flow domain for the remaining fluid that rotates around the cone's axis. Elasticity also has a pronounced effect on the stress profile, indicating that the accuracy of the cone and plate may be compromised even for small mis-alignments.     

On unsteady unidirectional flows of a second grade fluid

Some properties of unsteady unidirectional flows of a fluid of second grade are considered for flows produced by the sudden application of a constant pressure gradient or by the impulsive motion of one or two boundaries. Exact analytical solutions for these flows are obtained and the results are compared with those of a Newtonian fluid. It is found that the stress at the initial time on the stationary boundary for flows generated by the impulsive motion of a boundary is infinite for a Newtonian fluid and is finite for a second grade fluid. Furthermore, it is shown that initially the stress on the stationary boundary, for flows started from rest by sudden application of a constant pressure gradient is zero for a Newtonian fluid and is not zero for a fluid of second grade. The required time to attain the asymptotic value of a second grade fluid is longer than that for a Newtonian fluid. It should be mentioned that the expressions for the flow properties, such as velocity, obtained by the Laplace transform method are exactly the same as the ones obtained for the Couette and Poiseuille flows and those which are constructed by the Fourier method. The solution of the governing equation for flows such as the flow over a plane wall and the Couette flow is in a series form which is slowly convergent for small values of time. To overcome the difficulty in the calculation of the value of the velocity for small values of time, a practical method is given. The other property of unsteady flows of a second grade fluid is that the no-slip boundary condition is sufficient for unsteady flows, but it is not sufficient for steady flows so that an additional condition is needed. In order to discuss the properties of unsteady unidirectional flows of a second grade fluid, some illustrative examples are given.     

具有弹性支承的圆管在内外部流激发下动力特性及实验研究

本文建立了具有弹性支承的圆管在内外部流激发下的力学模型.推导了内部流与静止外部流作用下圆管的耦联振动方程.提出了确定弹性系数的方法.采用振型叠加法分析圆管动力特性问题.对内部流与静止外部流情况下圆管固有频率进行了计算和测量,计算值与实验值吻合较好.另外,对内外流同时激发下圆管的固有频率进行了测量,得到若干对工程实际有用的结论.     

Strongly oblique interactions between internal solitary waves with the same mode

I.IntroductionZabuskyandKruskal(l965)foundthattwoKdVsolitarywavesofthesamemodekeeptheiroriginalshapesandspeedsafterstronginteractions,andcalledthesewavessolitons.However,solitaryx"avessometimestravelintwodimensionalspace,otherthaninonedimensionalspace.Mil…     

Numerical simulation of time‐dependent hydrodynamic removal of a contaminated fluid from a cavity

The time‐dependent hydrodynamic removal of a contaminated fluid from a rectangular cavity on the floor of a duct is analysed numerically. Laminar duct flows are considered for Reynolds numbers of 50 and 1600 where the characteristic length is the duct height. Two cases are considered where: (1) the fluid density in the cavity is the same as that for the duct fluid and (2) the cavity fluid has a higher density than the duct fluid but the two fluids are miscible. The flow is solved by a numerical solution of the time‐dependent Navier–Stokes equations. Attention is focused on the convective transport of contaminated fluid out from the cavity and the effect of duct flow velocity profile on the cleaning process. Passive markers are used in the numerical simulation for the purpose of identifying the contaminated cavity fluid. The results show that the flow patterns in the cavity are influenced by the type of duct flow. From a cleaning perspective, the results suggest that it is easier for the duct flow to penetrate a cavity and to remove contaminated cavity fluid when the duct flow is of the Poiseuille type and the aspect ratio is large. Copyright © 2003 John Wiley & Sons, Ltd.     

A Novel Relative Permeability Model Based on Mixture Theory Approach Accounting for Solid–Fluid and Fluid–Fluid Interactions

A novel model is presented for estimating steady-state co- and counter-current relative permeabilities analytically derived from macroscopic momentum equations originating from mixture theory accounting for fluid–fluid (momentum transfer) and solid–fluid interactions (friction). The full model is developed in two stages: first as a general model based on a two-fluid Stokes formulation and second with further specification of solid–fluid and fluid–fluid interaction terms referred to as \(R_{{i}}\) (i =  water, oil) and R, respectively, for developing analytical expressions for generalized relative permeability functions. The analytical expressions give a direct link between experimental observable quantities (end point and shape of the relative permeability curves) versus water saturation and model input variables (fluid viscosities, solid–fluid/fluid–fluid interactions strength and water and oil saturation exponents). The general two-phase model is obeying Onsager’s reciprocal law stating that the cross-mobility terms \(\lambda _\mathrm{wo}\) and \(\lambda _\mathrm{ow}\) are equal (requires the fluid–fluid interaction term R to be symmetrical with respect to momentum transfer). The fully developed model is further tested by comparing its predictions with experimental data for co- and counter-current relative permeabilities. Experimental data indicate that counter-current relative permeabilities are significantly lower than corresponding co-current curves which is captured well by the proposed model. Fluid–fluid interaction will impact the shape of the relative permeabilities. In particular, the model shows that an inflection point can occur on the relative permeability curve when the fluid–fluid interaction coefficient \(I>0\) which is not captured by standard Corey formulation. Further, the model predicts that fluid–fluid interaction can affect the relative permeability end points. The model is also accounting for the observed experimental behavior that the water-to-oil relative permeability ratio \(\hat{{k}}_{\mathrm{rw}} /\hat{{\mathrm{k}}}_{\mathrm{ro}} \) is decreasing for increasing oil-to-water viscosity ratio. Hence, the fully developed model looks like a promising tool for analyzing, understanding and interpretation of relative permeability data in terms of the physical processes involved through the solid–fluid interaction terms \(R_{{i}}\) and the fluid–fluid interaction term R.     

Force acting on a particle in unsteady flow of a pseudoplastic fluid

The accelerated flow of a pseudoplastic fluid around a quiescent sphere at Reynolds numbers Re = 0–200 and dimensionless acceleration Ga = 10–10⁴ is studied by numerical simulation. It is shown that the analytical expression of the added mass force for an ideal fluid is appropriate for a pseudoplastic fluid. An expression for calculating the hereditary Basset force for a pseudoplastic fluid is proposed.     

Viscous fluid damping in a laterally oscillating finger of a comb-drive micro-resonator based on micro-polar fluid theory

Viscous damping is a dominant source of energy dissipation in laterally oscillating micro-structures. In micro-resonators in which the characteristic dimensions are compa-rable to the dimensions of the fluid molecules, the assumption of the continuum fluid theory is no longer justified and the use of micro-polar fluid theory is indispensable. In this paper a mathematical model was presented in order to predict the viscous fluid damping in a laterally oscillating finger of a micro-resonator considering micro-polar fluid theory. The coupled governing partial differential equations of motion for the vibration of the finger and the micro-polar fluid field have been derived. Considering spin and no-spin boundary conditions, the related shape functions for the fluid field were presented. The obtained governing differential equations with time varying boundary conditions have been trans-formed to an enhanced form with homogenous boundary conditions and have been discretized using a Galerkin-based reduced order model. The effects of physical properties of the micro-polar fluid and geometrical parameters of the oscillat-ing structure on the damping ratio of the system have been investigated.     

Criteria of negative wake generation behind a cylinder

It was demonstrated by simulation in our previous study that both the normal stress and its gradient are responsible for the negative wake generation (overshoot in the axial velocity) and streamline shifting. Extensional properties of the fluids dominate the generation of the negative wake, while other factors strengthen or weaken the formation of velocity overshoot. In this study, the criteria for the negative wake generation are discussed in detail for various fluid models, including the PTT, the FENE-CR, the FENE-P, and the Giesekus models. With the FENE-CR fluid, it is easier to generate negative wake than with the FENE-P fluid. This confirms that the constant shear viscosity FENE-CR fluid enhances the velocity overshoot, and that the shear-thinning viscosity FENE-P fluid delays the negative wake generation. The Giesekus fluid has a similar behaviour to the PTT fluid with regarding to the critical conditions of negative wake generation when appropriate fluid parameters are selected. The mechanism of wall proximity in enhancing the negative wake generation is also demonstrated with the analysis for the first time.     

Direct numerical simulation of particle–fluid systems by combining time-driven hard-sphere model and lattice Boltzmann method

A coupled numerical method for the direct numerical simulation of particle–fluid systems is formulated and implemented, resolving an order of magnitude smaller than particle size. The particle motion is described by the time-driven hard-sphere model, while the hydrodynamic equations governing fluid flow are solved by the lattice Boltzmann method (LBM). Particle–fluid coupling is realized by an immersed boundary method (IBM), which considers the effect of boundary on surrounding fluid as a restoring force added to the governing equations of the fluid. The proposed scheme is validated in the classical flow-around-cylinder simulations, and preliminary application of this scheme to fluidization is reported, demonstrating it to be a promising computational strategy for better understanding complex behavior in particle–fluid systems.