A simple four-node quadrilateral finite element for plates

A simple and accurate four-node quadrilateral finite element based on the Mindlin plate theory and Kirchhoff constraints is presented for general thin plate bending applications. The derivation of the element stiffness properties is straightforward, starting with a specified eight-node interpolation; usual discrete Kirchhoff (DK) constraints are employed to constrain out the four midside nodes of the element. The present resulting DK element passes patch tests with elements of arbitrary and even highly distorted mesh types. Numerical studies of the element convergence behaviours are undertaken for various plate bending problems so far investigated. It is indicated from comparative examples that fairly good convergence characteristics have been achieved.     

A discontinuous Galerkin method/HLLC solver for the Euler equations

This paper proposes a fully three‐dimensional non‐linear Euler methodology for solving aerodynamic and acoustic problems in the presence of strong shocks and rarefactions. It uses a discontinuous Galerkin method (DGM) within the element, and a Riemann solver (HLLC) at the boundaries to propagate rarefactions while preserving the entropy condition and capturing shocks with no spurious oscillations. This approach is thought to marry the best aspects of finite element and finite volume methods, achieving conservation while not requiring the solution of a large matrix. Examples in which shock and rarefaction waves are well captured are presented and the propagation of acoustic pulses is well demonstrated. Copyright © 2003 John Wiley & Sons, Ltd.     

Formation of β‐iPP in isotactic polypropylene/ethylene–propylene rubber blends: Effects of preparation method,composition, and thermal condition

The effects of preparation method, composition, and thermal condition on formation of β‐iPP in isotactic polypropylene/ethylene–propylene rubber (iPP/EPR) blends were studied using modulated differential scanning calorimeter (MDSC), wide angle X‐ray diffraction (WAXD), and phase contrast microscopy (PCM). It was found that the α‐iPP and β‐iPP can simultaneity form in the melt‐blended samples, whereas only α‐iPP exists in the solution‐blended samples. The results show that the formation of β‐iPP in the melt‐blended samples is related to the crystallization temperature and the β‐iPP generally diminishes and finally vanishes when the crystallization temperature moves far from 125 °C. The phenomena that the lower critical temperature of β‐iPP in iPP/EPR obviously increases to 114 °C and the upper critical temperature decreases to 134 °C indicate the narrowing of temperature interval, facilitating the formation of β‐iPP in iPP/EPR. Furthermore, it was found that the amount of β‐iPP in melt‐blended iPP/EPR samples is dependent on the composition and the maximum amount of β‐iPP formed when the composition of iPP/EPR blends is 85:15 in weight. The results through examining the effect of annealing for iPP/EPR samples at melt state indicate that this annealing may eliminate the susceptibility to β‐crystallization of iPP. However, only α‐iPP can be observed in solution‐blended samples subjected to annealing for different time. The PCM images demonstrate that an obvious phase‐separation happens in both melt‐blended and solution‐blended iPP/EPR samples, implying that compared with the disperse degree of EPR in iPP, the preparation method plays a dominant role in formation of β‐iPP. It is suggested that the origin of formation of β‐iPP results from the thermomechanical history of the EPR component in iPP/EPR. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 1704–1712, 2007     

Effect of nanotube functionalization on the properties of single‐walled carbon nanotube/polyurethane composites

A commercially available aliphatic thermoplastic polyurethane formulated with a methylene bis(cyclohexyl) diisocyanate hard segment and a poly(tetramethylene oxide) soft segment and chain‐extended with 1,4‐butanediol was dissolved in dimethylformamide and mixed with dispersed single‐walled carbon nanotubes. The properties of composites made with unfunctionalized nanotubes were compared with the properties of composites made with nanotubes functionalized to contain hydroxyl groups. Functionalization almost eliminated the conductivity of the tubes according to the conductivity of the composites above the percolation threshold. In most cases, functionalized and unfunctionalized tubes yielded composites with statistically identical mechanical properties. However, composites made with functionalized tubes did have a slightly higher modulus in the rubbery plateau region at higher nanotube fractions. Small‐angle X‐ray scattering patterns indicated that the dispersion reached a plateau in the unfunctionalized composites that was consistent with the plateau in the rubbery plateau region. The room‐temperature modulus and tensile strength increase was proportionally higher than almost all increases seen previously in thermoplastic polyurethanes; however, the increase was still an order of magnitude below what has been reported for the best nanotube–polymer systems. Nanotube addition increased the hard‐segment glass transition temperature slightly, whereas the soft‐segment glass transition was so diffuse that no conclusions could be drawn. Unfunctionalized tubes suppressed the crystallization of the hard segment; whereas functionalized tubes had no effect. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 490–501, 2007     

Eigenvalues of the antiplane-shear crack problem for a power-law material

This paper discusses the problem of finding the eigenvalue spectrum in determining the stress and strain fields at the tip of an antiplane-shear crack in a power-law material. It is shown that the perturbation method provides an analytical dependence of the eigenvalue on the material nonlinearity parameter and the eigenvalue of the linear problem. Thus, it is possible to find the entire spectrum of eigenvalues and not only the eigenvalue of the Hutchinson-Rice-Rosengren problem. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 1, pp. 173–180, January–February, 2008.     

Identification of nonstationary axisymmetric load distributed along a cylindrical shell

The paper outlines a procedure to identify the space-and time-dependent external nonstationary load acting on a closed circular cylindrical shell of medium thickness. Time-dependent deflections at several points of the shell are used as input data to solve the inverse problem. Examples of numerical identification of various nonstationary loads, including moving ones are presented. The relationship between the external load and the stress-strain state of the shell is described by the Volterra equation of the first kind. The identification problem is solved using Tikhonov's regularization method and Apartsin's h-regularization method __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 7, pp. 91–100, July 2008.     

Stability analysis of abnormal multiplication of plankton using parameter identification technique

This paper presents the mathematical approach for the abnormal multiplication of plankton. An abnormal multiplication can be expressed as an unstable problem and the stability of the system is investigated by introducing eigenvalues of a mathematical equation. The stability of the system can be judged by an eigenvalue based on the Lyapunov's stability theory. In this paper, the Arnoldi‐QR method is used to obtain eigenvalues and eigenvectors of the system. The mode superposition method is employed to create spatial distribution needed to analyse the stability. To obtain the objective eigenvalue, the parameter identification technique is employed. The finite element method is used for the discretization in space. Lake Kasumigaura, which is located in Ibaraki Prefecture in Japan, is selected and actual data in 1975, 1976, 1991 and 2000 are used in order to investigate the stability of the specified lake in Japan. Copyright © 2004 John Wiley & Sons, Ltd.     

Indirect boundary element method for unsteady linearized flow over prolate and oblate spheroids and hemispheroidal protuberances

The indirect boundary element method was used to study the hydrodynamics of oscillatory viscous flow over prolate and oblate spheroids, and over hemispheroidal bodies hinged to a plate. Analytic techniques, such as spheroidal coordinates, method of images, and series representations, were used to make the numerical methods more efficient. A novel method for computing the hydrodynamic torque was used, since for oscillatory flow the torque cannot be computed directly from the weightings. Instead, a Green's function for torque was derived to compute the torque indirectly from the weightings. For full spheroids, the method was checked by comparing the results to exact solutions at low and high frequencies, and to results computed using the singularity method. For hemispheroids hinged to a plate, the method for low frequencies was checked by comparing the results to previous results, and to exact solutions at high frequencies. Copyright © 2004 John Wiley & Sons, Ltd.     

A projection method for solving incompressible viscous flows on domains with moving boundaries

In this paper, a projection method is presented for solving the flow problems in domains with moving boundaries. In order to track the movement of the domain boundaries, arbitrary‐Lagrangian–Eulerian (ALE) co‐ordinates are used. The unsteady incompressible Navier–Stokes equations on the ALE co‐ordinates are solved by using a projection method developed in this paper. This projection method is based on the Bell's Godunov‐projection method. However, substantial changes are made so that this algorithm is capable of solving the ALE form of incompressible Navier–Stokes equations. Multi‐block structured grids are used to discretize the flow domains. The grid velocity is not explicitly computed; instead the volume change is used to account for the effect of grid movement. A new method is also proposed to compute the freestream capturing metrics so that the geometric conservation law (GCL) can be satisfied exactly in this algorithm. This projection method is also parallelized so that the state of the art high performance computers can be used to match the computation cost associated with the moving grid calculations. Several test cases are solved to verify the performance of this moving‐grid projection method. Copyright © 2004 John Wiley Sons, Ltd.     

Operator‐splitting method for the analysis of cavitation in liquid‐lubricated herringbone grooved journal bearings

This paper presents an operator‐splitting method (OSM) for the solution of the universal Reynolds equation. Jakobsson–Floberg–Olsson (JFO) pressure conditions are used to study cavitation in liquid‐lubricated journal bearings. The shear flow component of the oil film is first solved by a modified upwind finite difference method. The solution of the pressure gradient flow component is computed by the Galerkin finite element method. Present OSM solutions for slider bearings are in good agreement with available analytical and experimental results. OSM is then applied to herringbone grooved journal bearings. The film pressure, cavitation areas, load capacity and attitude angle are obtained with JFO pressure conditions. The calculated load capacities are in agreement with available experimental data. However, a detailed comparison of the present results with those predicted using Reynolds pressure conditions shows some differences. The numerical results showed that the load capacity and the critical mass of the journal (linear stability indicator) are higher and the attitude angle is lower than those predicted by Reynolds pressure conditions for cases of high eccentricities. Copyright © 2004 John Wiley & Sons, Ltd.