基于D稳定域和ITAE准则的主动队列管理算法

主动队列管理（active queue management，简称AQM）是网络拥塞控制的研究热点之一，其中的关键问题是如何设计反馈控制策略．提出一种新的基于D稳定域和时间乘以误差绝对值乘积积分（integral of time-weighted absolute error,简称ITAE）性能准则的比例-积分-微分（proportional-integral-differential，简称PID）优化设计方法（简称DITAE-PID），并用于AQM控制器的设计，控制闭环系统的理想动态性能．首先在复平面上设定一组理想的D稳定域，然后以ITAE为目标函数，通过数值优化算法求出控制器的参数，使得闭环系统的所有特征根都在D稳定域内，以降低排队延时，提高有效吞吐量．对比仿真实验结果表明孩算法能够预先探测和控制拥塞，有较好的鲁棒性，链路利用率更高，丢包率更小，平均队列长度更趋于期望值，同时，趋于期望队列长度的时间更短，其综合性能明显优于典型的随机早期探测（random early detection，简称RED）和比例-积分（proportional-integral，简称PI）算法．

Infrared heater arrays for warming ecosystem field plots

There is a need for methodology to warm open-field plots in order to study the likely effects of global warming on ecosystems in the future. Herein, we describe the development of arrays of more powerful and efficient infrared heaters with ceramic heating elements. By tilting the heaters at 45 degrees from horizontal and combining six of them in a hexagonal array, good uniformity of warming was achieved across 3-m-diameter plots. Moreover, there do not appear to be obstacles (other than financial) to scaling to larger plots. The efficiency [eta(h) (%); thermal radiation out per electrical energy in] of these heaters was higher than that of the heaters used in most previous infrared heater experiments and can be described by: eta(h) = 10 + 25exp(-0.17 u), where u is wind speed at 2 m height (m s(-1)). Graphs are presented to estimate operating costs from degrees of warming, two types of plant canopy, and site windiness. Four such arrays were deployed over plots of grass at Haibei, Qinghai, China and another at Cheyenne, Wyoming, USA, along with corresponding reference plots with dummy heaters. Proportional integral derivative systems with infrared thermometers to sense canopy temperatures of the heated and reference plots were used to control the heater outputs. Over month-long periods at both sites, about 75% of canopy temperature observations were within 0.5 degrees C of the set-point temperature differences between heated and reference plots. Electrical power consumption per 3-m-diameter plot averaged 58 and 80 kW h day(-1) for Haibei and Cheyenne, respectively. However, the desired temperature differences were set lower at Haibei (1.2 degrees C daytime, 1.7 degrees C night) than Cheyenne (1.5 degrees C daytime, 3.0 degrees C night), and Cheyenne is a windier site. Thus, we conclude that these hexagonal arrays of ceramic infrared heaters can be a successful temperature free-air-controlled enhancement (T-FACE) system for warming ecosystem field plots.

基于自调整神经元的无人直升机航向控制方法

直升机航向动力学包含输入非线性、时变参数和主-尾旋翼之间的强耦合,传统的比例积分微分(Proportional integral differential,PID)方法很难达到良好的控制性能。基于以上原因,通过把自调整神经元与滑模控制相结合,提出一种能够解决带有输入非线性的航向自适应控制方法。与常规自适应控制相比,用滑模条件代替误差函数作为目标函数,使控制器在保证闭环稳定性的同时,能够进一步使跟踪误差满足期望精度。证明了该方法的稳定性,针对实际模型直升机试验平台航向动力学模型的仿真结果,以及与传统PID方法的比较都表明了该方法的有效性。

Complex potentials and singular integral equation for curve crack problem in antiplane elasticity

Four types of the fundamental complex potential in antiplane elasticity are introduced: (a) a point dislocation, (b) a concentrated force, (c) a dislocation doublet and (d) a concentrated force doublet. It is proven that if the axis of the concentrated force doublet is perpendicular to the direction of the dislocation doublet, the relevant complex potentials are equivalent. Using the obtained complex potentials, a singular integral equation for the curve crack problem is introduced. Some particular features of the obtained singular integral equation are discussed, and numerical solutions and examples are given.

Crack Tip Field and J-Integral with Strain Gradient Effect

The mode I plane strain crack tip field with strain gradient effects is presented in this paper based on a simplified strain gradient theory within the framework proposed by Acharya and Bassani. The theory retains the essential structure of the incremental version of the conventional J_2 deformation theory No higher-order stress is introduced and no extra boundary value conditions beyond the conventional ones are required. The strain gradient effects are considered in the constitutive relation only through the instantaneous tangent modulus. The strain gradient measures are included into the tangent modulus as internal parameters. Therefore the boundary value problem is the same as that in the conventional theory Two typical crack Problems are studied: (a) the crack tip field under the small scale yielding condition induced by a linear elastic mode-I K-field and (b) the complete field for a compact tension specimen. The calculated results clearly show that the stress level near the crack tip with strain gradient effects is considerable higher than that in the classical theory The singularity of the strain field near the crack tip is nearly equal to the square-root singularity and the singularity of the stress field is slightly greater than it. Consequently, the J-integral is no longer path independent and increases monotonically as the radius of the calculated circular contour decreases.

Quasi-Dammann Grating With Proportional Intensity Array Spots

A quasi-Dammann grating is proposed to generate array spots with proportional-intensity orders in the far field. To describe the performance of the grating, the uniformities of the array spots are redefined. A two-dimensional even-sampling encode scheme is adopted to design the quasi-Dammann grating. Numerical solutions of the binary-phase quasi-Dammann grating with proportional-intensity orders are given. The experimental results with a third-order quasi-Dammann grating, which has an intensity proportion of 3:2:1 from zero order to second order, are presented. (C) 2008 Optical Society of America

Stress intensity factors calculation in anti-plane fracture problem by orthogonal integral extraction method based on femol

For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions (or between the derivatives of eigenfunctions) of anti-plane problem is exploited. We developed for the first time two sets of radius-independent orthogonal integrals for extraction of stress intensity factors (SIFs), so any order SIF can be extracted based on a certain known solution of displacement (an analytic result or a numerical result). Many numerical examples based on the finite element method of lines (FEMOL) show that the present method is very powerful and efficient.

The mixed-type integral equations for transient elastodynamic plane crack problems

In the present paper, by use of the boundary integral equation method and the techniques of Green fundamental solution and singularity analysis, the dynamic infinite plane crack problem is investigated. For the first time, the problem is reduced to solving a system of mixed-typed integral equations in Laplace transform domain. The equations consist of ordinary boundary integral equations along the outer boundary and Cauchy singular integral equations along the crack line. The equations obtained are strictly proved to be equivalent with the dual integral equations obtained by Sih in the special case of dynamic Griffith crack problem. The mixed-type integral equations can be solved by combining the numerical method of singular integral equation with the ordinary boundary element method. Further use the numerical method for Laplace transform, several typical examples are calculated and their dynamic stress intensity factors are obtained. The results show that the method proposed is successful and can be used to solve more complicated problems.

Study of three elastic-plastic constitutive models by non-proportional finite deformations of OFHC copper

Experimental stress-strain data of OFHC copper first under torsion to 13% and then under torsion-tension to about 10% are used to study the characteristics of three elastic-plastic constitutive models: Chaboche's super-positional nonlinear model, Dafalias and Popov's two surface model and Watanabe and Atluri's version of the endochronic model. The three models, originally oriented for infinitesimal deformation, have been extended for finite deformation. The results show (a) the Mises-type yield surface used in the three models brings about significant departure of the predictions from the experimental data; (b) Chaboche's and Dafalias' models are easier than Watanabe and Atluri's model in determining the material parameters in them, and (c) Chaboche's and Watanabe & Atluri's models produce almost the same prediction to the data, while Dafalias' model cannot accurately predict the plastic deformations when a loading path changes in its direction. Copyright (C) 1996 Elsevier Science Ltd

Cauchy singular integral equation method for transient antiplane dynamic problems

In this paper, by use of the boundary integral equation method and the techniques of Green basic solution and singularity analysis, the dynamic problem of antiplane is investigated. The problem is reduced to solving a Cauchy singular integral equation in Laplace transform space. This equation is strictly proved to be equivalent to the dual integral equations obtained by Sih [Mechanics of Fracture, Vol. 4. Noordhoff, Leyden (1977)]. On this basis, the dynamic influence between two parallel cracks is also investigated. By use of the high precision numerical method for the singular integral equation and Laplace numerical inversion, the dynamic stress intensity factors of several typical problems are calculated in this paper. The related numerical results are compared to be consistent with those of Sih. It shows that the method of this paper is successful and can be used to solve more complicated problems. Copyright (C) 1996 Elsevier Science Ltd

A new type of boundary integral-equation for plane problems of elasticity including rotational forces

Based on the idea proposed by Hu [Scientia Sinica Series A XXX, 385-390 (1987)], a new type of boundary integral equation for plane problems of elasticity including rotational forces is derived and its boundary element formulation is presented. Numerical results for a rotating hollow disk are given to demonstrate the accuracy of the new type of boundary integral equation.

Numerical-Value Perturbation Method with Application of Solving Convective-Diffusion Integral Equation

The convective--diffusion equation is of primary importance in such fields as fluid dynamics and heat transfer hi the numerical methods solving the convective-diffusion equation, the finite volume method can use conveniently diversified grids (structured and unstructured grids) and is suitable for very complex geometry The disadvantage of FV methods compared to the finite difference method is that FV-methods of order higher than second are more difficult to develop in three-dimensional cases. The second-order central scheme (2cs) offers a good compromise among accuracy, simplicity and efficiency, however, it will produce oscillatory solutions when the grid Reynolds numbers are large and then very fine grids are required to obtain accurate solution. The simplest first-order upwind (IUW) scheme satisfies the convective boundedness criteria, however. Its numerical diffusion is large. The power-law scheme, QMCK and second-order upwind (2UW) schemes are also often used in some commercial codes. Their numerical accurate are roughly consistent with that of ZCS. Therefore, it is meaningful to offer higher-accurate three point FV scheme. In this paper, the numerical-value perturbational method suggested by Zhi Gao is used to develop an upwind and mixed FV scheme using any higher-order interpolation and second-order integration approximations, which is called perturbational finite volume (PFV) scheme. The PFV scheme uses the least nodes similar to the standard three-point schemes, namely, the number of the nodes needed equals to unity plus the face-number of the control volume. For instanc6, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D problems, 2~Dand 3-D flow model equations. Comparing with other standard three-point schemes, The PFV scheme has much smaller numerical diffusion than the first-order upwind (IUW) scheme, its numerical accuracy are also higher than the second-order central scheme (2CS), the power-law scheme (PLS), the QUICK scheme and the second-order upwind(ZUW) scheme.

Diffraction integral formulas of the pulsed wave field in the temporal domain

To resolve the diffraction problems of the pulsed wave field directly in the temporal domain, we extend the Rayleigh diffraction integrals to the temporal domain and then discuss the approximation condition of this diffraction formula. (C) 1997 Optical Society of America.

The accuracy of Kirchhoff's approximation in describing the far field speckles produced by random self-affine fractal surfaces

Based on the rigorous formulation of integral equations for the propagations of light waves at the medium interface, we carry out the numerical solutions of the random light field scattered from self-affine fractal surface samples. The light intensities produced by the same surface samples are also calculated in Kirchhoff's approximation, and their comparisons with the corresponding rigorous results show directly the degree of the accuracy of the approximation. It is indicated that Kirchhoff's approximation is of good accuracy for random surfaces with small roughness value w and large roughness exponent alpha. For random surfaces with larger w and smaller alpha, the approximation results in considerable errors, and detailed calculations show that the inaccuracy comes from the simplification that the transmitted light field is proportional to the incident field and from the neglect of light field derivative at the interface.