The study on a kind of control system with nonlinear parabolic distributed parameters

The modelling of one kind of nonlinear parabolic distributed parameter control system with moving boundary, which had extensive applications was presented. Two methods were used to investigate the basic characteristics of the system: 1) transforming the system in the variable domain into that in the fixed domain; 2) transforming the distributed parameter system into the lumped parameter system. It is found that there are two critical values for the control variable: the larger one determines whether or not the boundary would move, while the smaller one determines whether or not the boundary would stop automatically. For one-dimensional system of planar, cylindrical and spherical cases the definite solution problem can be expressed as a unified form. By means of the computer simulation the open-loop control system and close-cycle feedback control system have been investigated. Numerical results agree well with theoretical results. The computer simulation shows that the system is well posed, stable, measurable and controllable. Foundation item: the National Natural Science Foundation of China (59936140) Biographies:ZHOU Jian-jun (1953-); XU Yan-hou (1933-)     

A realtime observatory for laboratory simulation of planetary flows

Motivated by the mid-latitude atmospheric circulation, we develop a system that uses observations from a differentially heated rotating annulus experiment to constrain a numerical simulation in real-time. The coupled physical-numerical system provides a tool to rapidly prototype new methods for state and parameter estimation, and facilitates the study of prediction, predictability, and transport of geophysical fluids where observations or numerical simulations would not independently suffice. A computer vision system is used to extract measurements from the physical simulation, which constrain the model-state of the MIT general circulation model in a hybrid data assimilation approach. Using a combination of parallelism, domain decomposition and an efficient scheme to select ensembles of model-states, we show that estimates that effectively track the fluid-state can be produced. To the best of our knowledge, this is the first realtime coupled system for this laboratory analog of planetary circulation.     

AN ARBITRARILY MIXED EXPLICIT-IMPLICIT ASYNCHRONOUS INTEGRATION ALGORITHM BASED ON UNIFORM DISCRETIZATION FORMAT

动力学问题的有限元分析需要在每一时步求解系统信息,相对于静力学问题,其计算量要大得多.因而,提高计算效率,节省计算工作量是动力学求解方法研究的主要内容.该文针对大型复杂动力学系统的高效求解问题,提出了一种基于Newmark离散格式的显式、隐式任意混合异步算法,根据整体系统不同局部的物理力学特性和求解精度要求,在空间域及时间域内对动力学系统方程进行多尺度求解.该方法根据显式、隐式算法固有的信息传递机制,采取动态的可变边界处理方法,避免了异步边界上的误差积累;并通过对整体系统能量平衡的校验,动态地确定和修正仿真计算时步,可以有效地预防不稳定性的产生和发展.数值算例表明：该算法能在保持较高的计算精度的同时,极大地降低计算资源消耗,因而具有一定的实用价值.     

统一格式的显式与隐式任意混合异步算法

动力学问题的有限元分析需要在每一时步求解系统信息,相对于静力学问题,其计算量要大得多.因而,提高计算效率,节省计算工作量是动力学求解方法研究的主要内容.该文针对大型复杂动力学系统的高效求解问题,提出了一种基于Newmark离散格式的显式、隐式任意混合异步算法,根据整体系统不同局部的物理力学特性和求解精度要求,在空间域及时间域内对动力学系统方程进行多尺度求解.该方法根据显式、隐式算法固有的信息传递机制,采取动态的可变边界处理方法,避免了异步边界上的误差积累;并通过对整体系统能量平衡的校验,动态地确定和修正仿真计算时步,可以有效地预防不稳定性的产生和发展.数值算例表明:该算法能在保持较高的计算精度的同时,极大地降低计算资源消耗,因而具有一定的实用价值.      

利用改进的ADI方法实现单机千万节点的油水两相渗流数值模拟

针对实际油藏的非均质分布特征及其复杂的边界条件,本文通过引入迭代参数的压缩因子和放大因子,对现有的交替方向迭代法(ADI)进行改进,提出一种适用于大规模油藏数值模拟的新算法．改进的ADI 方法计算精度可靠,且与现有的算法相比计算效率有所提高．更为关键的是,ADI 算法将求解三维压力方程的七对角矩阵分解为三个方向的三对角矩阵依次迭代求解,所需的运算存储量大幅减少,最大的计算规模有了大幅的提升．使用改进的ADI 方法,在单机上成功实现了千万节点的油水两相渗流数值模拟．计算实例表明,在同等单机硬件条件下,改进ADI 算法的最大计算规模是现有算法的1.7 倍以上．     

Adaptive Boundary Conditions for Exterior Flow Problems

We consider the problem of solving numerically the stationary incompressible Navier–Stokes equations in an exterior domain in two dimensions. This corresponds to studying the stationary fluid flow past a body. The necessity to truncate for numerical purposes the infinite exterior domain to a finite domain leads to the problem of finding appropriate boundary conditions on the surface of the truncated domain. We solve this problem by providing a vector field describing the leading asymptotic behavior of the solution. This vector field is given in the form of an explicit expression depending on a real parameter. We show that this parameter can be determined from the total drag exerted on the body. Using this fact we set up a self-consistent numerical scheme that determines the parameter, and hence the boundary conditions and the drag, as part of the solution process. We compare the values of the drag obtained with our adaptive scheme with the results from using traditional constant boundary conditions. Computational times are typically reduced by several orders of magnitude.     

基于可拓决策法的车辆自适应避撞控制方法研究

车辆的避撞控制可以有效避免或缓解车辆的碰撞事故,是自动驾驶汽车的关键控制技术之一.各种交通条件、不确定的道路附着系数以及复杂的液压制动执行系统都会降低避撞控制的有效性.因此,本文提出了一种基于可拓决策法的自适应避撞控制,该控制方法对路面附着系数具有自适应性,并能够精确控制制动系统的制动液压力.首先,设计了滑模观测器来估计轮胎纵向力,并基于观测得到的轮胎纵向力,进一步提出带遗忘因子的递推最小二乘法估计道路附着系数.其次,基于递推最小二乘法的估计值,提出了基于路面附着系数自适应调节的自适应避撞控制方法,该方法基于可拓决策法的先决判定决定当前时刻应采用何种避撞控制策略,即采用可拓决策方法判断进行点刹预警制动、全制动或者不制动.再次,通过对执行系统-电控液压制动系统进行精确的液压控制,实现主动减速控制效果.最后,采用软件联合仿真手段对上述方法进行了验证,结果表明所提出的算法在避撞行驶工况中具有良好的避撞效果.     

Locally similar solutions for hydromagnetic and thermal slip flow boundary layers over a flat plate with variable fluid properties and convective surface boundary condition

This paper presents heat transfer process in a two-dimensional steady hydromagnetic convective flow of an electrically conducting fluid over a flat plate with partial slip at the surface of the boundary subjected to the convective surface heat flux at the boundary. The analysis accounts for both temperature-dependent viscosity and temperature dependent thermal conductivity. The local similarity equations are derived and solved numerically using the Nachtsheim-Swigert iteration procedure. Results for the dimensionless velocity, temperature and ambient Prandtl number within the boundary layer are displayed graphically delineating the effect of various parameters characterizing the flow. The results show that momentum boundary layer thickness significantly depends on the surface convection parameter, Hartmann number and on the sign of the variable viscosity parameter. The results also show that plate surface temperature is higher when there is no slip at the plate compared to its presence. For both slip and no-slip cases surface temperature of the plate can be controlled by controlling the strength of the applied magnetic field. In modelling the thermal boundary layer flow with variable viscosity and variable thermal conductivity, the Prandtl number must be treated as a variable irrespective of flow conditions whether there is slip or no-slip at the boundary to obtain realistic results.     

二维含裂纹结构与声耦合问题研究

应用分形有限元方法结合边界元方法研究了二维含裂纹结构和声耦合问题.采用二级分形有限元方法对含裂纹的弹性结构体进行离散处理,这样可以使得自由度数大大地减少;无限大外域声场的计算使用边界元方法,可以自动满足无穷远辐射条件.数值仿真算例结果表明:结构声耦合系统的共振频率随着裂纹深度的增加而下降;裂纹附近的声场所受的影响较为明显.     

圆锥壳自由振动传递函数解

本文在线性弹性理论基础上，给出了一种求解圆锥薄壳自由振动的渐进传递函数解法，壳体的三个位移分量，外力和边界条件首先沿环向展开的Ｆｏｕｒｉｅｒ级数，然后关于时间变量进行Ｌａｐｌａｃｅ变换，这样就将壳体的控制方程化为一系列含复参数ｓ的变系数常微分方程组，通过定义状态变量。得到了壳体动力学问题的状态空间控制微分方程，引入一小参数，并利用摄动技术就可以得到微分方程的渐进传递函数解，将各于锥段的解进行综合，     

物理参数直接识别法对于分布系统的应用

本文将物理参数直接识别法用于均匀分布系统等效离散模型的物理参数识别，数值模拟结果表明，识别得到的离散模型的动力特性较好地接近分布系统的动力特性，离散系统物理参数识别的有关结论对分布系统仍然成立，数值模拟与实验结果还表明，均匀分布系统时域响应测量点应等距离布置。     

Effect of uncertainty on the bifurcation behavior of pitching airfoil stall flutter

In this paper the effect of system parametric uncertainty on the stall flutter bifurcation behavior of a pitching airfoil is studied. The aerodynamic moment on the two-dimensional rigid airfoil with nonlinear torsional stiffness is computed using the ONERA dynamic stall model. The pitch natural frequency, a cubic structural nonlinearity parameter, and the structural equilibrium angle are assumed to be uncertain. The effect on the amplitude of the response, the bifurcation of the probability distribution, and the flutter boundary is considered. It is demonstrated that the system parametric uncertainty results already in 5% probability of pitching stall flutter at a 12.5% earlier position than the point where a deterministic analysis would predict unstable behavior. Probabilistic collocation is found to be more efficient than the Galerkin polynomial chaos method and Monte Carlo simulation for modeling uncertainty in the post-bifurcation domain.     

弹性力学问题的局部边界积分方程方法

提出了弹性力学平面问题的局部边界积分方程方法。这种方法是一种无网格方法,它采用移动最小二乘近似试函数,且只包含中心在所考虑节点的局部边界上的边界积分。它易于施加本质边界条件。所得系统矩阵是一个带状稀疏矩阵。它组合了伽辽金有限元法、整体边界元法和无单元伽辽金法的优点。该方法可以容易推广到求解非线性问题以及非均匀介质的力学问题。计算了两个弹性力学平面问题的例子,给出了位移和能量的索波列夫模,所得计算结果证明：该方法是一种具有收敛快、精度高、简便有效的通用方法。     

Shear deformable bars of doubly symmetrical cross section under nonlinear nonuniform torsional vibrations—application to torsional postbuckling configurations and primary resonance excitations

In this paper a boundary element method is developed for the nonuniform torsional vibration problem of bars of arbitrary doubly symmetric constant cross section, taking into account the effects of geometrical nonlinearity (finite displacement—small strain theory) and secondary twisting moment deformation. The bar is subjected to arbitrarily distributed or concentrated conservative dynamic twisting and warping moments along its length, while its edges are subjected to the most general axial and torsional (twisting and warping) boundary conditions. The resulting coupling effect between twisting and axial displacement components is also considered and a constant along the bar compressive axial load is induced so as to investigate the dynamic response at the (torsional) postbuckled state. The bar is assumed to be adequately laterally supported so that it does not exhibit any flexural or flexural–torsional behavior. A coupled nonlinear initial boundary value problem with respect to the variable along the bar angle of twist and to an independent warping parameter is formulated. The resulting equations are further combined to yield a single partial differential equation with respect to the angle of twist. The problem is numerically solved employing the Analog Equation Method (AEM), a BEM based method, leading to a system of nonlinear Differential–Algebraic Equations (DAE). The main purpose of the present contribution is twofold: (i) comparison of both the governing differential equations and the numerical results of linear or nonlinear free or forced vibrations of bars ignoring or taking into account the secondary twisting moment deformation effect (STMDE) and (ii) numerical investigation of linear or nonlinear free vibrations of bars at torsional postbuckling configurations. Numerical results are worked out to illustrate the method, demonstrate its efficiency and wherever possible its accuracy.

     

时变参数时滞减振控制研究

时滞动力吸振器对谐波激励有着良好的减振控制效果,但对随机激励的减振控制效果却并不明显,具体表现为时滞动力吸振器对随机激励的减振控制效果与被动吸振器几乎相同。针对上述问题,本文提出了一种新的时变参数时滞减振控制方法。在原有时滞减振控制方法的基础上,首先将时滞增益系数由定值形式变为时间函数形式,然后通过时变优化得到多组时滞控制参数并使其以一定时间周期循环作用于减振控制过程,通过这种方法进一步改善了时滞动力吸振器减振性能。本文最后以二自由度时滞动力吸振器减振模型为例,以主系统的振动响应为仿真对象,运用精细积分法求解了具有时变时滞参数的时滞动力学方程,以此得到了在谐波激励和随机激励作用下主系统振动的时域仿真结果。研究结果表明,在时变参数时滞动力吸振器的控制下,主系统无论是受谐波激励作用还是受随机激励作用,其振动位移、振动速度和振动加速度均比在定值参数时滞动力吸振器控制下时有大幅的减少,时滞动力吸振器的减振性能有了明显的改善。      

Control of a liquid–liquid heat exchanger

 The analytical solutions of equations describing the dynamics of distributed parameter systems are usually complicated in form and derivations and inconvenient to use for simulation and control system design. The liquid–liquid counter flow heat exchanger is an example of these disturbed parameter systems. An analytical solution of the dynamics of a symmetrically operated counter flow heat exchanger in the form of transfer function matrix is investigated in open-loop and close-loop conditions. The resulting non-linear model was linearized using perturbation approach. A feed-forward path controller to counteract at any disturbances in the boundary temperature and a non-interactive controller to decouple the outputs were implemented. Received on 22 February 2001 / Published online: 29 November 2001     

Numerical simulation of trajectory and deformation of bubble in tip vortex

According to the behaviors of a bubble in the ship wake flow,the numerical simulation is divided into two stages,quasi-spherical motion and non-spherical motion,based on whether the bubble is captured by the vortex or not.The one-way coupled particle tracking method(PTM) and the boundary element method(BEM) are adopted to simulate these two stages,respectively.Meanwhile,the initial condition of the second stage is taken as the output of the first one,and the entire simulation is connected and completed.Based on the numerical results and the published experimental data,the cavitation inception is studied,and the wake bubble is tracked.Besides,the split of the bubble captured by the vortex and the following sub-bubbles are simulated,including motion,deformation,and collapse.The results provide some insights into the control on wake bubbles and optimization of the wake flow.     

Finite deflection dynamics of elastic beams

Solutions are obtained for the problem of an infinite elastic beam subjected to essentially constant velocity boundary conditions at one point of the beam. The effects of finite deflections, normal force, rotatory inertia and shear deformation are included. The equations of the problem are converted into non-dimensional form and a perturbation approach is used to obtain a consistent approximation. Numerical solutions are obtained for the bending moment, shear force and the normal force for different velocities of impact. It is shown that the solution to the problem depends on a combined geometrical and material parameter which does not vary significantly for compact sections and a loading parameter which determines the amplitude of the response. Finally the linear Timoshenko beam theory is shown to predict the bending moment and shear force extremely well even when the deflections are large enough to cause appreciable stretching of the centroidal axis.     

Unsteady conjugate problem of a dissipative fluid in a horizontal channel with a periodic variation temperature

The unsteady two-dimensional transient heat transfer problem referring to a fully laminar flow developing in a parallel-plane channel exposed to a periodic variation surface temperature with distance is numerically studied. The effects of channel thickness, Péclet number, wall-to-fluid conductivity ratio, thermal diffusivity ratio, angular frequency and the viscous dissipation parameter are determined in the solutions. The non-linear equations are discretized by means an implicit finite difference scheme and the electric analogy to the resulting system is applied to convert these equations into a network-electrical model that was solved using a computer code (electric circuits simulator). In this scheme, only spatial discretization is necessary, while time remains as a real continuous variable, and its programming does not require manipulation of the sophisticated mathematical software that is inherent in other numerical methods. The network simulation method, which satisfies the conservation law for the heat flux variable and the uniqueness law for temperature, also permits the direct visualization of the local and/or integrated transport variables at any point or section of the medium.