A boundary-field equation method for a nonlinear exterior elasticity problem in the plane

This paper presents a new method for solving a nonlinear exterior boundary value problem arising in two-dimensional elasto-plasticity. The procedure is based on the introduction of a sufficiently large circle that divides the exterior domain into a bounded region and an unbounded one. This allows us to consider the Dirichlet-Neumann mapping on the circle, which provides an explicit formula for the stress in terms of the displacement by using an appropriate infinite Fourier series. In this way we can reduce the original problem to an equivalent nonlinear boundary value problem on the bounded domain with a natural boundary condition on the circle. Hence, the resulting weak formulation includes boundary and field terms, which yields the so called boundary-field equation method. Next, we employ the finite Fourier series to obtain a sequence of approximating nonlinear problems from which the actual Galerkin schemes are derived. Finally, we apply some tools from monotone operators to prove existence, uniqueness and approximation results, including Cea type error estimates for the corresponding discrete solutions.     

Numerical study of the parametric evolution of bifurcations in the three-dimensional ring problem of N bodies

We study the dynamics of a massless particle in an annular configuration of N bodies, N − 1 of which have equal masses m and are located in equal distances on a fictitious circle and one has mass βm and is located at the center of the circle. Our interest is focused on the bifurcation points from planar to three-dimensional families of symmetric periodic orbits in the above problem. We study numerically the evolution of these bifurcation points with respect to the variation of the mass parameter β. In particular we investigate the continuous evolution of bifurcation points for values of β from 2 up to 1000. The two distinct cases of the system’s behavior at β = 2 and 1000 are examined comparatively and various conclusions are drawn regarding the overall dynamical evolution of the three-dimensional system as the relative mass of the central body grows.     

高超声速飞行器横侧向失稳非线性分岔分析

针对滑翔式高超声速飞行器大攻角横侧向失稳问题,采用延拓算法和分岔理论,求解并分析了以俯仰舵偏为连续参数的稳态平衡分岔图和以副翼舵偏为连续参数的横侧向机动稳态平衡分岔图,对平衡分支的稳定性和突变点进行了分析,并给出了特征根拓扑结构变化.研究表明,高超声速飞行器存在极限分岔点、Hopf分岔点以及叉型分岔点,且从叉型分岔点延伸出多个平衡分支,引起横侧向的自滚转失稳;从Hopf分岔点延伸出极限环分支,该分支对应较为复杂的极限环运动,其中还包含倍周期分岔、花环分岔、极限环极限点分岔等复杂的分岔现象;在横侧向机动飞行情况下,模型存在横向操作偏离失稳问题,且存在多个不稳定的平衡点.研究结果为实现高超声速飞行器的稳定飞行和控制器的设计提供了极其重要的动力学信息.      

Absolute stabilization related to circle criterion: An LMI-based approach

A method is proposed for synthesizing output feedback controllers for nonlinear Lur' e systems . The problem of designing an output dynamic controller for uncertain-free systems and systems subject to multiplicative norm-bounded perturbations in the linear part were proposed respectively. The procedure is based on the use of the absolute stability, through the circle criterion, and a linear matrix inequalities (LAI) formulation. The controller existence conditions are given in terms of existence of suitable solutions to a set of parameter-dependent LMIs.     

Combining magnetic resonance measurements with numerical simulations – Extracting blood flow physiology information relevant to the investigation of intracranial aneurysms in the circle of Willis

Cerebral aneurysms in the region of the circle of Willis have an incidence of 3–6% in western populations and involve the risk of rupture with subsequent subarachnoidal bleeding. The patient specific blood flow patterns are of substantial importance for understanding the pathogenesis of the lesions and may eventually contribute to deciding on the most efficient treatment procedure for a specific patient.A non-invasive method for performing in vivo measurements on blood velocity is 4D phase-contrast magnetic resonance angiography (PC-MRA), on the basis of which a flow field with all its parameters can be simulated. We are using this approach to investigate the hemodynamic parameters in the circle of Willis and, by analyzing the values at common locations of aneurysms, trying to find potential parameters to predict the development of aneurysms. Methodologically, we are acquiring the artery geometry with 3D-time-of-flight magnetic resonance (TOF) measurements and the blood velocity in the feeding arteries with 4D PC-MRA measurements in a healthy volunteer. These measurements are combined with computational fluid dynamics (CFD) to describe detailed hemodynamic patterns within the circle of Willis.     

Evaluation of the variable order of singularities by pseudocaustics

An experimental technique based on the optical method of pseudocaustics was developed for the evaluation of the order of singularity for singular elastic fields. According to this method, the mapping of an infinitesimal circle close to the singularity into an infinitesimal ellipse by means of the pseudocaustics may be used for evaluating directly the stress singularity. An experimental procedure and a numerical scheme were used for accurately defining the order of singularity. Results with variable values of stress singularities indicate the potentialities of the method.     

Dynamics of a quantized circle homeomorphism with quasiperiodic perturbation

We introduce the notion of a quantized circle homeomorphism that is a discontinuous mapping of an interval shift, which is widely used in modern digital radio electronics. For a two-dimensional dynamical system given by a triangular mapping, which is a quantized circle homeomorphism with quasiperiodic perturbation, we prove, under some assumptions, that there exist an invariant absorbing belt and a repulsive contour, study properties of these structures, and get estimates for their sizes. To make the exposition complete, we, first, study the corresponding problems for three less complicated systems, namely, a proper circle homeomorphism, a proper circle homeomorphism with quasiperiodic perturbation, and a quantized circle homeomorphism without perturbation.     

Natural oscillations of an ideal liquid in a thin toroidal tube rotating about the axis of symmetry

The natural oscillation spectrum of an ideal liquid in a rotating toroidal tube is considered. A recurrent procedure is presented which is applicable for any form of tube cross section and permits determination of spectrum points and eigenfunctions in the form of a series in the ratio of tube thickness to ring radius. The procedure is applied to the cases of rectangular and circular cross sections. In both cases the spectrum is computable and densely fills the real axis.     

Areal moments of inertia revisited: on the distinction between the principal directions

Three commonly used methods to determine the principal moments of inertia of a plane area and their directions are based on: (i) the stationarity condition for the axial moment of inertia, (ii) the eigenvalue analysis, and (iii) Mohr’s circle. In this paper we provide two new derivations, which are based on: (a) the matrix diagonalization and the invariant tensor properties, and (b) the conjugacy property of the moment of inertia vectors. A new general expression is derived which specifies the principal directions of inertia, as well as the directions of the maximum and minimum product of inertia. A comparative study of the five presented approaches is given, which is of interest from both conceptual and methodological points of view. The connection between the deviatoric part of the moment of inertia tensor and Land’s circle of inertia is also given. The presented analysis applies to any two-by-two symmetric second order tensor.     

随机响应面法最优概率配点数目分析

系统地研究随机响应面法采用线性无关原则选取概率配点的优越性,给出了基于线性无关原则选取概率配点的流程图,比较了基于回归方法和基于线性无关原则选取概率配点的优缺点。算例结果表明,基于回归方法选取概率配点时,配点数目应保证Hermite系数矩阵达到满秩,此时随机响应面法的计算精度才能得到保证,计算效率也远远高于传统的蒙特卡洛模拟方法。基于线性无关原则选取概率配点的随机响应面法在保证计算精度的同时,其计算效率远远高于基于回归方法选取概率配点的随机响应面法,它是结构可靠度分析一种有效的方法,尤其适用于极限状态方程不能用显式函数表达的复杂结构可靠度问题。研究成果为随机响应面法最优概率配点数目的确定奠定了一定的基础。     

2D and 3D separation of stresses using automated photoelasticity

A procedure for the separation of full-field photoelastic images for use with an automated polariscope is described. Regions of background in the image are identified thus producing the boundary of the model. The shear difference method is used to calculate the components of stress along all raster lines in the image using photoelastic parameters at the boundary points to calculate the initial values of stress. Algorithms were also used to evaluate the stress components along raster lines which did not contain boundary points. A plastic template was used to evaluate the efficiency of the boundary routine. It was found that it was able to identify edges to within approximately one pixel on screen. The complete procedure for stress separation was evaluated using a stress frozen disc in compression and a turbine slot. The values of stress found using the automated polariscope with the stress-separation procedure were found to agree well with theory and with results determined using the method of Tardy compensation and manual analysis. The automated polariscope was also used to analyze three-dimensional stress components along arbitrary lines of a 3D model. A two-model slicing regime was used to analyze a strut subjected to a vertical load. This work was compared to results obtained by Frocht and Guernsey on an identical model machined from Fosterite and subjected to a higher load. Good agreement was found between the results for points away from the region of loading. Significant differences were found near to the load point, however. A finite element analysis of the same problem suggested that this was due to the effects of plasticity.     

Approximated Repair Methods for Outlier Strain Data from Hole-Drilling Residual Measurements

The current ASTM E 837 standard gives the standardized procedure for the evaluation of uniform and non uniform residual stresses, that is, stresses that do not vary and vary significantly with depth from the specimen surface, respectively. For non uniform stresses, the standard states that many small increments should be done in order to have a stable calculation of the residual stress profile. In addition, it states that irregularities as well as outlier strain points should be investigated and if necessary, the hole-drilling test should be repeated. In some applications outside the laboratory, the availability to repeat a test with outlier points is not possible. In these cases, the standard does not show the more appropriated way which should be followed to use only valid measured data (without outliers). For this reason, a stress profile corresponding to a shot peening test was simulated and one hole step was included in different hole depths as outlier point in order to evaluate the feasibility of some proposed ways of computation. These ways were: (a) following the ASTM procedure but replacing outlier strain points by linear interpolation of neighboring good ones, (b) following the ASTM procedure but replacing outlier strain points by a parabola interpolation of neighboring good ones and (c) using cumulative relaxation strain functions and only good measured points. Statistical criteria were also introduced and developed in order to identify outlier points. Results show a practical procedure to detect outlier points in experimental strain data.     

A simple mechanism for the polygonalization of railway wheels by wear

A mechanism is presented which explains the growing non-circularity of railway wheels by wear as a result of the interaction between the deviation from the circle and the rotational inertia of the wheels. The behaviour of the model is investigated by perturbation techniques, including the method of multiple time scales. Results conform with observations: The faster the train the quicker surge lower harmonics of the non-circularity. The procedure can be applied to study various other wear mechanisms even such where wear is influenced by locally variable material characteristics.     

Automated critical point identification for PIV data using multimodal particle swarm optimization

A hybrid sequential niche algorithm is used for the automated identification of critical points of velocity fields. This method combines an adaptive sequential niche technique with deterministic local optimization to detect critical points: focus, node and saddle points. A particle swarm algorithm performs a global search whereas vortex core identification functions compute the precise location as the extremum of the corresponding function. Once a critical point is found, a rectangular niche is constructed around the point. The particle swarm then proceeds to explore different regions of the velocity field. The process advances sequentially, avoiding areas near previously found critical points by blocking niches obtained from previous steps. The niche size is automatically adjusted each time a search enters inside an existing niche. Vortex core functions are used for critical point identification and calculating its precise location inside each niche. The procedure is validated on particle image velocimetry data obtained with two types of flows, an impinging jet flow and a flow downstream of a model building. The hybrid algorithm proved to be very efficient and robust for automated detection and identification of critical points. It can be used as a first step for studying the time‐dependent dynamic behavior of instantaneous velocity fields by tracking topological critical points. This is the first study that uses a multi‐modal particle swarm algorithm for critical point identification. Copyright © 2011 John Wiley & Sons, Ltd.     

Scattering by an anisotropic circle

The scattering by a circle is considered when the outside medium is isotropic and the inside medium is anisotropic (orthotropic). The problem is a scalar one and is phrased as a scattering problem for elastic waves with polarization out of the plane of the circle (SH wave), but the solution is with minor modifications valid also for scattering of electromagnetic waves. The equation inside the circle is first transformed to polar coordinates and it then explicitly contains the azimuthal angle through trigonometric functions. Making an expansion in a trigonometric series in the azimuthal coordinate then gives a coupled system of ordinary differential equations in the radial coordinate that is solved by power series expansions. With the solution inside the circle complete the scattering problem is solved essentially as in the classical case. Some numerical examples are given showing the influence of anisotropy, and it is noted that the effects of anisotropy are generally strong except at low frequencies where the dominating scattering only depends on the mean stiffness and not on the degree of anisotropy.     

Time-harmonic BEM for 2-D piezoelectricity applied to eigenvalue problems

We will derive the fundamental generalized displacement solution, using the Radon transform, and present the direct formulation of the time-harmonic boundary element method (BEM) for the two-dimensional general piezoelectric solids. The fundamental solution consists of the static singular and the dynamics regular parts; the former, evaluated analytically, is the fundamental solution for the static problem and the latter is given by a line integral along the unit circle. The static BEM is a component of the time-harmonic BEM, which is formulated following the physical interpretation of Somigliana’s identity in terms of the fundamental generalized line force and dislocation solutions obtained through the Stroh–Lekhnitskii (SL) formalism. The time-harmonic BEM is obtained by adding the boundary integrals for the dynamic regular part which, from the original double integral representation over the boundary element and the unit circle, are reduced to simple line integrals along the unit circle.The BEM will be applied to the determination of the eigen frequencies of piezoelectric resonators. The eigenvalue problem deals with full non-symmetric complex-valued matrices whose components depend non-linearly on the frequency. A comparative study will be made of non-linear eigenvalue solvers: QZ algorithm and the implicitly restarted Arnoldi method (IRAM). The FEM results whose accuracy is well established serve as the basis of the comparison. It is found that the IRAM is faster and has more control over the solution procedure than the QZ algorithm. The use of the time-harmonic fundamental solution provides a clean boundary only formulation of the BEM and, when applied to the eigenvalue problems with IRAM, provides eigen frequencies accurate enough to be used for industrial applications. It supersedes the dual reciprocity BEM and challenges to replace the FEM designed for the eigenvalue problems for piezoelectricity.     

鞍马全旋运动的动力学

本文从经典力学观点讨论人体在鞍马上的全旋运动,理想化的全旋运动要求人体以不变的章动角作匀速进动,同时绕纵轴朝相反方向自旋,以保证在运动过程中始终面向前方。文中规定一种带特殊完整约束条件的刚体规则进动作为这种运动的理想模式,并利用人体的多刚体模型分析鞍马上的直体全旋及托马斯全旋的运动规律,导出肩关节支撑力及肌肉控制力矩的计算公式。     

Dynamic instability of vibrating carbon nanotubes near small layers of graphite sheets based on nonlocal continuum elasticity

This article presents a new asymptotic method to predict dynamic pull-in instability of nonlocal clamped–clamped carbon nanotubes (CNTs) near graphite sheets. Nonlinear governing equations of carbon nanotubes actuated by an electric field are derived. With due allowance for the van der Waals effects, the pull-in instability and the natural frequency–amplitude relationship are investigated by a powerful analytical method, namely, the parameter expansion method. It is demonstrated that retaining two terms in series expansions is sufficient to produce an acceptable solution. The obtained results from numerical methods verify the strength of the analytical procedure. The qualitative analysis of system dynamics shows that the equilibrium points of the autonomous system include center points and unstable saddle points. The phase portraits of the carbon nanotube actuator exhibit periodic and homoclinic orbits.     

Singular perturbation theory for homoclinic orbits in a class of near-integrable Hamiltonian systems

This paper describes a new type of orbits homoclinic to resonance bands in a class of near-integrable Hamiltonian systems. It presents a constructive method for establishing whether small conservative perturbations of a family of heteroclinic orbits that connect pairs of points on a circle of equilibria will yield transverse homoclinic connections between periodic orbits in the resonance band resulting from the perturbation. In any given example, this method may be used to prove the existence of such transverse homoclinic orbits, as well as to determine their precise shape, their asymptotic behavior, and their possible bifurcations. The method is a combination of the Melnikov method and geometric singular perturbation theory for ordinary differential equations.