基于Fourier级数的时变周期系数Riccati微分方程精细积分

结合Fourier级数展开方法,本文提出了基于精细积分的时变周期系数Riccati微分方程求解高效算法.首先,利用Fourier级数展开方法将周期系统表示成三角级数形式,在一个积分步内使用精细积分方法得到对应Hamilton系统状态转移矩阵的表达式.然后,通过Riccati变换的方法,得到含有状态转移矩阵的时变周期系数Riccati微分方程解的递推格式.本文方法充分利用了方程本身的周期性特点,文中的数值算例表明算法具有计算效率高、结果可靠等优势.     

输液管道临界流速的傅里叶级数法

傅里叶级数法被用于计算输液管道的临界流速，与有限元等数值法相比，更为简单可靠。     

Optimal control of material concentration using Fourier series

This paper presents an optimal control of the material concentration using Fourier series and finite element method. It is assumed that the optimal control value can be expanded into a Fourier series. The Fourier coefficient is identified to minimize the performance function and the optimal control value is determined. The Sakawa–Shindo algorithm is used for the minimization algorithm. The advection–diffusion equation and shallow water equation are used for the analysis of material concentration and water flow. The Crank–Nicolson scheme and finite element method using bubble function element with stabilized control parameter are employed as temporal and special discretization. Copyright © 2004 John Wiley & Sons, Ltd.     

Obtaining Fourier series graphically from large amplitude oscillatory shear loops

Large amplitude oscillatory shear (LAOS) flow has been used to characterize the nonlinear viscoelasticity of polymer melts and solutions. Results are frequently reported with shear stress versus strain loops, or with shear stress versus shear rate loops. A Fourier analysis of the stress response to LAOS is often desired for comparison with theory, or for quantitative comparison between resins. A method is presented which employs the discrete Fourier transform to obtain the Fourier series coefficients from LAOS loops.     

串联结构体系可靠度的二元泰勒级数展开

本文给出串联结构体系可靠度的二元泰勒级数展开表达式，并给出有关一维无穷积分的数值解法。当用于计算串联结构体系的可靠度时，与一元泰勒级数展开式相比，计算效率大大提高，精度足以满足工程应用要求。特别是如果不考虑级数展开式中的二次修正项，则使结构体系可靠度的计算变得非常简单。     

复合材料回转壳的有限元级数解

本文用有限元法和Fourier级数展开技术求解复合材料回转壳体在各种荷载作用下的弯曲问题,文中利用回转壳在几何上的轴对称性质,将各物理量在环向展开为Fourei级数,而在母线和壳厚方向分割单元,所采用的单元为6节点18自由度等参元,它考虑了剪切变形和挤压变形的影响,能计算厚度方向的挤压应力,数值算例表明,本文提出的单元性能优良,算法稳定收敛。     

Improving convergence of incremental harmonic balance method using homotopy analysis method

We have deduced incremental harmonic balance an iteration scheme in the （IHB） method using the harmonic balance plus the Newton-Raphson method. Since the convergence of the iteration is dependent upon the initial values in the iteration, the convergent region is greatly restricted for some cases. In this contribution, in order to enlarge the convergent region of the IHB method, we constructed the zeroth-order deformation equation using the homotopy analysis method, in which the IHB method is employed to solve the deformation equation with an embedding parameter as the active increment. Taking the Duffing and the van der Pol equations as examples, we obtained the highly accurate solutions. Importantly, the presented approach renders a convenient way to control and adjust the convergence.     

Nonlinear dynamic analysis of a quarter vehicle system with external periodic excitation

In this paper, a nonlinear dynamic model of a quarter vehicle with nonlinear spring and damping is established. The dynamic characteristics of the vehicle system with external periodic excitation are theoretically investigated by the incremental harmonic balance method and Newmark method, and the accuracy of the incremental harmonic balance method is verified by comparing with the result of Newmark method. The influences of the damping coefficient, excitation amplitude and excitation frequency on the dynamic responses are analyzed. The results show that the vibration behaviors of the vehicle system can be control by adjusting appropriately system parameters with the damping coefficient, excitation amplitude and excitation frequency. The multi-valued properties, spur-harmonic response and hardening type nonlinear behavior are revealed in the presented amplitude-frequency curves. With the changing parameters, the transformation of chaotic motion, quasi-periodic motion and periodic motion is also observed. The conclusions can provide some available evidences for the design and improvement of the vehicle system.     

基于傅里叶级数的汽轮机叶根轮槽结构设计

提出了一种基于傅里叶级数的枞树型叶根轮槽的优化设计方法。在叶根轮槽构型表征方面,采用傅里叶级数-直线法,与传统的圆弧-直线法相比,不但具有高保真度和完备性,还扩大了优化的搜索空间,减少了设计变量数目,有效降低计算代价。在优化设计方面,为避免以最小化峰值应力为目标函数引起的迭代震荡和不收敛现象,采用KS函数凝聚应力值,降低了优化问题的非线性程度,保证优化快速稳定收敛。以三齿枞树型叶根轮槽为实例,优化结果表明,傅里叶级数-直线法和KS凝聚法优化可有效降低汽轮机叶根轮槽的应力水平,提高汽轮机结构的可靠性。同时,提出的结构表征和优化设计方法具有较高的通用性和实际工程应用价值,可方便扩展至其他工程结构型线设计。     

频响函数计算的高精度级数展开法

把作者近期提出的用于频响函数计算的级数展开法进行了推广，使其计算精度进一步提高，并且能适用于自由－自由系统。数值示例表明，本文提出的方法是有效的。     

Connection between Volterra series and perturbation method in nonlinear systems analyses

This paper is concerned with the connection between the Volterra series and the regular perturbation method in nonlinear systems analyses. It is revealed for the first time that, for a forced polynomial nonlinear system, if its derived linear system is a damped dissipative system, the steady response obtained through the regular perturbation method is exactly identical to the response given by the Volterra series. On the other hand, if the derived linear system is an undamped conservative system, then the Volterra series is incapable of modeling the forced polynomial nonlinear system. Numerical examples are further presented to illustrate these points. The results provide a new criterion for quickly judging whether the Volterra series is applicable for modeling a given polynomial nonlinear system.     

Visco-elastic systems under both deterministic harmonic and random excitation

IntroductionForlinearviscoelasticsystemsunderbothadditiveandmultiplicativebroad_bandexcitationexcitations,Ariaratnam[1]studiedthestochasticstabilityofthesystembyusingthemethodofstochasticaveragingprocedure .Itwasshownthatthevisco_elasticforcecontributedtowarddamping ,hence ,stabilityofthesystem .However,thestiffnesseffectofthevisco_elasticcomponentwasnotfullyaccountedfor.FurthermoreAriaratnam[2 ]studiedthestochasticstabilityofthesystembutthemodelislinear.Inthetheoryofnonlinearrandomvibration…     

Improvement of Fourier Polarimetry for Applications in Tomographic Photoelasticity

The use of the Fourier Polarimetry method has been demonstrated to extract the three characteristic parameters in integrated photoelasticity. In contrast to the phase-stepping method, it has been shown that the Fourier method is more accurate. However, the Fourier method isn't very efficient as it requires that a minimum of nine intensity images be collected during a whole revolution of a polarizer while the phase-stepping method only needs six intensity images. In this paper the Fourier transformation is used to derive the expression for determination of the characteristic parameters. Four Fourier coefficients are clearly identified to calculate the three characteristic parameters. It is found that the angular rotation ratio could be set arbitrarily. The angular rotation ratio is optimized to satisfy the requirements of efficiency and proper data accuracy, which results in data collection about three times faster than the methods suggested by previous researchers. When comparing their performance in terms of efficiency and accuracy, the simulated and experimental results show that these angular rotation ratios have the same accuracy but the optimized angular rotation ratio is significantly faster. The sensitivity to noise is also investigated and further improvement of accuracy is suggested.     

A multi-scale modeling framework for instabilities of film/substrate systems

Spatial pattern formation in stiff thin films on soft substrates is investigated from a multi-scale point of view based on a technique of slowly varying Fourier coefficients. A general macroscopic modeling framework is developed and then a simplified macroscopic model is derived. The model incorporates Asymptotic Numerical Method (ANM) as a robust path-following technique to trace the post-buckling evolution path and to predict secondary bifurcations. The proposed multi-scale finite element framework allows sinusoidal and square checkerboard patterns as well as their bifurcation portraits to be described from a quantitative standpoint. Moreover, it provides an efficient way to compute large-scale instability problems with a significant reduction of computational cost compared to full models.     

On a Discrete Fourier Series Solution to Static Problem for Conical Shells of Circumferentially Varying Thickness

An approach is developed to solve the two-dimensional boundary-value problems of the stress-strain state of conical shells with circumferentially varying thickness. The approach employs discrete Fourier series to separate variables and make the problem one-dimensional. The one-dimensional boundary-value problem is solved by the stable discrete-orthogonalization method. The results obtained are presented as plots and tables __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 9, pp. 26–37, September 2005.     

Boundary element elastic stress analysis of 3D generally anisotropic solids using fundamental solutions based on Fourier series

The authors have very recently proposed an efficient, accurate alternative scheme to numerically evaluate etc. Green’s function, U(x), and its derivatives for three-dimensional, general anisotropic elasticity. These quantities are necessary items in the formulation of the boundary element method (BEM). The scheme is based on the double Fourier series representation of the explicit, exact, algebraic solution derived by Ting and Lee (1997) [Ting, T.C.T., Lee, V.G., 1997. The three-dimensional elastostic Green’s function for general anisotropic linear elastic solid. Q. J. Mech. Appl. Math. 50, 407–426] expressed in terms of Stroh’s eigenvalues. By taking advantage of some its characteristics, the formulations in this double Fourier series approach are revised and several of the analytical expressions are re-arranged in the present study. This results in significantly fewer terms to be summed in the series thereby improving further the efficiency for evaluating the Green’s function and its derivatives. These revised Fourier series representations of U(x) and its derivatives are employed in a BEM formulation for three-dimensional linear elastostatics. Some numerical examples are presented to demonstrate the veracity of the implementation and its applicability to the elastic stress analysis of generally anisotropic solids. The results are compared with known solutions in the literature where possible, and with those obtained using the commercial finite element code ANSYS. Excellent agreement is obtained in all cases.     

扇形截面杆扭转的应力分析

柱体扭转的基本方程为非齐次偏微分方程,在极坐标系下,利用分离变量法及傅立叶级数展开法,求出了扭转应力函数,进一步即可计算出扇形截面杆在外力偶作用下,扭转角和横截面上剪应力的精确解答。这种方法为精确解法,在各种机械及其他工程设备中,对受扭转作用的扇形截面杆设计,有一定实用价值。     

A novel finite point method for flow simulation

A novel finite point method is developed to simulate flow problems. The mashes in the traditional numerical methods are supplanted by the distribution of points in the calculation domain. A local interpolation based on the properties of Taylor series expansion is used to construct an approximation for unknown functions and their derivatives. An upwind‐dominated scheme is proposed to efficiently handle the non‐linear convection. Comparison with the finite difference solutions for the two‐dimensional driven cavity flow and the experimental results for flow around a cylinder shows that the present method is capable of satisfactorily predicting the flow separation characteristic. The present algorithm is simple and flexible for complex geometric boundary. The influence of the point distribution on computation time and accuracy of results is included. Copyright © 2002 John Wiley & Sons, Ltd.     

A new auxiliary equation method for finding travelling wave solutions to KdV equation

In this paper, a new auxiliary equation method is used to find exact travelling wave solutions to the （1＋1）-dimensional KdV equation. Some exact travelling wave solu- tions with parameters have been obtained, which cover the existing solutions. Compared to other methods, the presented method is more direct, more concise, more effective, and easier for calculations. In addition, it can be used to solve other nonlinear evolution equations in mathematical physics.