GENERAL SOLUTION OF THE OVERALL BENDING OF FLEXIBLE CIRCULAR RING SHELLS WITH MODERATELY SLENDER RATIO AND APPLICATIONS TO THE BELLOWS (Ⅳ)-CALCULATION FOR U-SHAPED BELLOWS

This is one of the applications of Part (I), in which the angular stiffness, and the corresponding stress distributions of U-shaped bellows were discussed. The bellows was divided into protruding sections, concave sections and ring plates for the calculation that the general solution (I) with its reduced form to ring plates were used respectively, but the continuity of the surface stresses and the meridian rotations at each joint of the sections were entirely satisfied. The resent results were compared with those of the slender ring shell solution proposed earlier by the authors, the standards of the Expansion Joint Manufacturers Association (EJMA), the experiment and the finite element method. It is shown that the governing equation and the general solution (I) are very effective. Contributed by HUANG Qian Biography: ZHU Wei-ping (1962-)     

断裂力学问题的杂交边界点方法

提出了一种求解断裂力学的新的边界类型无网格方法-杂交边界点法.以修正变分原理和移动最小二乘近似为基础,同时具有边界元法和无网格法的优良特性,求解时仅仅需要边界上离散点的信息.该文将杂交边界点方法应用到弹性断裂问题中,将移动最小二乘方法中的基函数扩充,能更好的模拟裂纹尖端应力场的奇异性,推导了求解断裂力学的杂交边界点法方程,与传统的元网格方法相比,文中方法具有后处理简单,计算精度高的优点.数值算例表明了该方法的稳定性和有效性.     

Zum Laval-Druckverhältnis bei verlustbehafteter Expansion

Übersicht Die Bestimmung des Laval-Punktes bei verlustbehafteter Expansion als Ort des Maximums der Massenstromdichte aus den hierfür ermittelten Versuchswerten ist bekanntlich wegen deren unvermeidlicher Streuung ziemlich ungenau. Hier wird versucht, die Laval-Punkte dadurch genauer zu ermitteln, daß auf die Kurven des Geschwindigkeitsbeiwertes abhängig vom Expansionsdruckverhältnis zurückgegriffen wird, da diese Werte durch Impulsverlustmessungen mit höherer Genauigkeit erhalten werden können.

---

Summary Determining the Laval point in expansions with energy-losses by finding the location of the maximum mass flow density from experimental data is rather inaccurate because of inherent dispersion of these empirically found data. This report attempts to locate Laval points more accurately by considering the velocity coefficient as a function of the pressure ratio, since values of this function may be determined to a high degree of accuracy by measuring the impulse loss.

     

Inflation in ω-field cosmoloty

In this paper, we shall apply the ω-field theory as first proposed by Yu¹³ to cosmology. Under the assumption that the spacetime geometry of the Universe is described by the Robertson-Walker metric and the matter tensor-consists only of theω-field, the Universe is found to follow a de Sitter Expansion. The horizon and flatness problems may thus be explained in a simple and natural way.     

Again discussing about singular perturbation of general boundary value problem for higher order elliptic equation containing two-parameter

In this paper using the method of The Two-Variable Expansion Procedure [11] we again discuss the construction of asymptotic expression of solution of general boundary value problem for higher order ellitptic equation containing two-parameter whose boundary condition is more general than [1]. We give asymptotic expression of solution as well as the estimation corresponding to the remainder term.     

AGAIN DISCUSSING ABOUT SINGULAR PERTURBATION OF GENERAL BOUNDARY VALUE PROBLEM FOR HIGHER ORDER ELLIPTIC EQUATION CONTAINING TWO-PARAMETER

In this paper using the method of "The Two-Variable Expansion Procedure" we again discuss the construction of asymptotic expression of solution of general boundary value problem for higher order ellitptic equation containing two-parameter whose boundary condition is more general than [1]. We give asymptotic expression of solution as well as the estimation corresponding to the remainder term.     

An exact element method for plane problem

In this paper, based on the step reduction method^[1] and exact analytic method^[2] anew method-exact element method for constructing finite element, is presented. Since the new method doesn ’t need the variational principle, it can be applied to solve non-positive and positive definite partial differential equations with arbitrary variable coefficient. By this method, a quadrilateral noncompatible element with 8 degrees of freedom is derived for the solution of plane problem. Since Jacobi ’s transformation is not applied, the present element may degenerate into a triangle element. It is convenient to use the element in engineering. In this paper, the convergence is proved. Numerical examples are given at the endof this paper, which indicate satisfactory results of stress and displacements can be obtained and have higher numerical precision in nodes.     

A comparison of three different methods for measuring both normal stress differences of viscoelastic liquids in torsional rheometers

A novel pressure sensor plate (normal stress sensor (NSS) from RheoSense, Inc.) was adapted to an Advanced Rheometrics Expansion System rheometer in order to measure the radial pressure profile for a standard viscoelastic fluid, a poly(isobutylene) solution, during cone–plate and parallel-plate shearing flows at room temperature. We observed in our previous experimental work that use of the NSS in cone-and-plate shearing flow is suitable for determining the first and second normal stress differences N ₁ and N ₂ of various complex fluids. This is true, in part, because the uniformity of the shear rate at small cone angles ensures the existence of a simple linear relationship between the pressure [i.e., the vertical diagonal component of the total stress tensor (Π₂₂)] and the logarithm of the radial position r (Christiansen and coworkers, Magda et al.). However, both normal stress differences can also be calculated from the radial pressure distribution measured in parallel-plate torsional flows. This approach has rarely been attempted, perhaps because of the additional complication that the shear rate value increases linearly with radial position. In this work, three different methods are used to investigate N ₁ and N ₂ as a function of shear rate in steady shear flow. These methods are: (1) pressure distribution cone–plate (PDCP) method, (2) pressure distribution parallel-plate (PDPP) method, and (3) total force cone–plate parallel-plate (TFCPPP) method. Good agreement was obtained between N ₁ and N ₂ values obtained from the PDCP and PDPP methods. However, the measured N ₁ values were 10–15% below the certified values for the standard poly(isobutylene) solution at higher shear rates. The TFCPPP method yielded N ₁ values that were in better agreement with the certified values but gave positive N ₂ values at most shear rates, in striking disagreement with published results for the standard poly(isobutylene) solution.

A Bubble Element for the Compressible Euler Equations

In this paper, we present a new Galerkin finite element method with bubble function for the compressible Euler equations. This method is derived from the scaled bubble element for the advection-diffusion problems developed by Simo and his colleagues, which is based on the equivalence between the Galerkin method employing piecewise linear interpolation with bubble functions and the Streamline-Upwind/Petrov Galerkin (SUPG) finite element method using P₁ approximation in the steady advection-diffusion problem. Simo and this author have applied this approach to transient advection-diffusion problems by using a special scaled bubble function called P-scaled bubble, which is designed to work in the transient advection-diffusion problems for any Peclet number from 0 to ∞. The method presented in this paper is an application of this p-scaled bubble element to a pure hyperbolic system.     

Nonstationary pandom vibpation analysis of linear elastic structures with finite element method

At present, the finite element method is an efficient method for analyzing structural dynamic problems. When the physical quantities such as displacements and stresses are resolved in the spectra and the dynamic matrices are obtained in spectral resolving form, the relative equations cannot be solved by the vibration mode resolving method as usual. For solving such problems, a general method is put forward in this paper. The excitations considered with respect to nonstationary processes are as follows: P(t)={P_i(t)},P_i(t)=a_i(t)P_i(t), a_i(t) is a time function already known. We make Fourier transformation for the discretized equations obtained by finite element method, and by utilizing the behaviour of orthogonal increment of spectral quantities in random process^[1], some formulas of relations about the spectra of excitation and response are derived. The cross power spectral denisty matrices of responses can be found by these formulas, then the structrual safety analysis can be made. When a_i(t)=l (i= 1,2,…n), the. method stated in this paper will be reduced to that which is used in the special case of stationary process.     

On the multiple scales perturbation method for difference equations

In the classical multiple scales perturbation method for ordinary difference equations (O Δ Es) as developed in 1977 by Hoppensteadt and Miranker, difference equations (describing the slow dynamics of the problem) are replaced at a certain moment in the perturbation procedure by ordinary differential equations (ODEs). Taking into account the possibly different behavior of the solutions of an O Δ E and of the solutions of a nearby ODE, one cannot always be sure that the constructed approximations by the Hoppensteadt–Miranker method indeed reflect the behavior of the exact solutions of the O Δ Es. For that reason, a version of the multiple scales perturbation method for O Δ Es will be presented and formulated in this paper completely in terms of difference equations. The goal of this paper is not only to present this method, but also to show how this method can be applied to regularly perturbed O Δ Es and to singularly perturbed, linear O Δ Es.     

A study ofJ-integral of the orthotropic composite material

The relation between J-integral near model I crack tip in the orthotropic plate and displacement derivative is derived in this paper. Meanwhile, the relation between stress intensity factor K _I and displacement is also given in this paper. With sticking film moire interferometry method, the three-point bending beam is tested, thus the values of J-integral and K _I can be obtained from the displacement field, and then the truth of relation formula between J-integral and K _I in the orthotropic composite materials is experimentally verified.     

APPLICATION OF THE DUAL RECIPROCITY BOUNDARY ELEMENT METHOD TO SOLUTION OF NONLINEAR DIFFERENTIAL EQUATION

In this paper the dual reciprocity boundary element method is employed to solve nonlinear differential equation ∇² u+u+ɛu ³ =b. Results obtained in an example have a good agreement with those by FEM and show the applicability and simplicity of dual reciprocity method(DRM)in solving nonlinear differential equations.     

Determination of the dynamic stress intensity factorsK I d andK II d,for a mixed-mode propagating crack

In this paper, the dynamic propagation problem of a mixed-mode crack was studied by means of the experimental method of caustics. The initial curve and caustic equations were derived under the mixed-mode dynamic condition. A multi-point measurement method for determining the dynamic stress intensity factors,K _I ^d , andK _II ^d , and the position of the crack tip was developed. Several other methods were adopted to check this method, and showed that it has a good precision. Finally, the dynamic propagating process of a mixed-mode crack in the three-point bending beam specimen was investigated with our method.     

Sorting-free Hill-based stability analysis of periodic solutions through Koopman analysis

In this paper, we aim to study nonlinear time-periodic systems using the Koopman operator, which provides a way to approximate the dynamics of a nonlinear system by a linear time-invariant system of higher order. We propose for the considered system class a specific choice of Koopman basis functions combining the Taylor and Fourier bases. This basis allows to recover all equations necessary to perform the harmonic balance method as well as the Hill analysis directly from the linear lifted dynamics. The key idea of this paper is using this lifted dynamics to formulate a new method to obtain stability information from the Hill matrix. The error-prone and computationally intense task known by sorting, which means identifying the best subset of approximate Floquet exponents from all available candidates, is circumvented in the proposed method. The Mathieu equation and an n-DOF generalization are used to exemplify these findings.

     

A closed form solution for stress intensity factors of shear modes for 3-D finite bodies with cracks by energy release rate method

In this paper, a new analytical-engineering method of closed form solution for stress intensity factors of shear modes for 3-D finite bodies with cracks is derived by the time saving energy release rate method. Hence a complete series of useful results of stress intensity factorsK _II andK _III can be obtained. And the results provided by this method are in good agreement with some obtained by other methods.     

New numerical method for partial differential equations. 1: Application to the diffusion equation

In this paper a new, highly accurate method called PH is presented for the numerical integration of partial differential equations. The method is applied for the solution of the one-dimensional diffusion equation. Upon integrating the equation within a subdomain of space and time using the prismoidal approximation, a three-point implicit scheme is obtained with a truncation error of order O(k⁴, h⁶), where k and h represent the time and space steps respectively. The method is stable under the condition s = αk/h² ? S(δ), where the function S(δ) increases as the parameter δ decreases from 1/12 to negative values. In practice the method behaves as unconditionally stable upon choosing an appropriate value for δ. A new formula is also adopted for the implementation of a Neumann boundary condition, introducing a truncation error of order O(h⁴). Numerical solutions are obtained incorporating Dirichlet and Neumann boundary conditions. The results prove that our method is far more accurate than any other-implicit or explicit method.     

Local field correction PIV: on the increase of accuracy of digital PIV systems

Cross-correlation Particle Image Velocimetry (PIV) has become a well known and widely used experimental technique. It has been already documented that difficulties arise resolving velocity structures smaller than the interrogation window. This is caused by signal averaging over this window. A new cross-correlation PIV method that eliminates this restriction is presented. The new method brings some other enhancements, such as the ability to deal with large velocity gradients, seeding density inhomogeneities, and high dispersion in the brightness of the particles. The final result is a method with a remarkable capability for accurately resolving small scale structures in the flow, down to a few times the mean distance between particles. When compared to particle tracking velocimetry, the new method is capable of obtaining measurements at high seeding density concentrations. Therefore, better overall performance is obtained, especially with the limited resolutions of video CCDs. In this paper, the new method is described and its performance is evaluated and compared to traditional PIV systems using synthetic images. Application to real PIV data are included and the results discussed. Received: 9 March 1998 / Accepted: 25 August 1998     

The Laplace transform method for Burgers' equation

The Laplace transform method (LTM) is introduced to solve Burgers' equation. Because of the nonlinear term in Burgers' equation, one cannot directly apply the LTM. Increment linearization technique is introduced to deal with the situation. This is a key idea in this paper. The increment linearization technique is the following: In time level t, we divide the solution u(x, t) into two parts: u(x, t^k) and w(x, t), t^k⩽t⩽t^k+1, and obtain a time‐dependent linear partial differential equation (PDE) for w(x, t). For this PDE, the LTM is applied to eliminate time dependency. The subsequent boundary value problem is solved by rational collocation method on transformed Chebyshev points. To face the well‐known computational challenge represented by the numerical inversion of the Laplace transform, Talbot's method is applied, consisting of numerically integrating the Bromwich integral on a special contour by means of trapezoidal or midpoint rules. Numerical experiments illustrate that the present method is effective and competitive. Copyright © 2009 John Wiley & Sons, Ltd.     

THE METHOD OF THE RECIPROCAL THEOREM OF FORCED VIBRATION FOR THE ELASTIC THIN RECTANGULAR PLATES(Ⅰ)-RECTANGULAR PLATES WITH FOUR CLAMPED EDGES AND WITH THREE CLAMPED EDGES

In this paper the method of the reciprocal theorem (MRT) is extended to solve the steady state responses of rectangular plates under harmonic disturbing forces. A series of the closed solutions of rectangular plates with various boundary conditions are given and the tables and figures which have practical value are provided. MRT is a simple, convenient and general method for solving the steady state responses of rectangular plates under various harmonic disturbing forces. The paper contains three parts: (I) rectangular plates with four clamped edges and with three clamped edges; (II) rectangular plates with two adjacent clamped edges; (III) cantilever plates. We are going to publish them one after another.