New eight node serendipity quadrilateral plate bending element for thin and moderately thick plates using Integrated Force Method

A new 8-node serendipity quadrilateral plate bending element (MQP8) based on the Mindlin-Reissner theory for the analysis of thin and moderately thick plate bending problems using Integrated Force Method is presented in this paper. The performance of this new element (MQP8) is studied for accuracy and convergence by analyzing many standard benchmark plate bending problems. This new element MQP8 performs excellent in both thin and moderately thick plate bending situations. And also this element is free from spurious/zero energy modes and free from shear locking problem.

New twelve node serendipity quadrilateral plate bending element based on Mindlin-Reissner theory using Integrated Force Method

The Integrated Force Method (IFM) is a novel matrix formulation developed for analyzing the civil, mechanical and aerospace engineering structures. In this method all independent/internal forces are treated as unknown variables which are calculated by simultaneously imposing equations of equilibrium and compatibility conditions. This paper presents a new 12-node serendipity quadrilateral plate bending element MQP12 for the analysis of thin and thick plate problems using IFM. The Mindlin-Reissner plate theory has been employed in the formulation which accounts the effect of shear deformation. The performance of this new element with respect to accuracy and convergence is studied by analyzing many standard benchmark plate bending problems. The results of the new element MQP12 are compared with those of displacement-based 12-node plate bending elements available in the literature. The results are also compared with exact solutions. The new element MQP12 is free from shear locking and performs excellent for both thin and moderately thick plate bending situations.

Closed form solutions for element equilibrium and flexibility matrices of eight node rectangular plate bending element using integrated force method

Closed form solutions for equilibrium and flexibility matrices of the Mindlin-Reissner theory based eight-node rectangular plate bending element (MRP8) using integrated Force Method (IFM) are presented in this paper. Though these closed form solutions of equilibrium and flexibility matrices are applicable to plate bending problems with square/rectangular boundaries, they reduce the computational time significantly and give more exact solutions. Presented closed form solutions are validated by solving large number of standard square/rectangular plate bending benchmark problems for deflections and moments and the results are compared with those of similar displacement-based eight-node quadrilateral plate bending elements available in the literature. The results are also compared with the exact solutions.

Flow over a cylinder with a hinged-splitter plate

Previous work on rigid splitter plates in the wake of a bluff body has shown that the primary vortex shedding can be suppressed for sufficiently long splitter plates. In the present work, we study the problem of a hinged-splitter plate in the wake of a circular cylinder. The splitter plate can rotate about the hinge at the base of the cylinder due to the unsteady fluid forces acting on it, and hence the communication between the two sides of the wake is not totally disrupted as in the rigid splitter plate case. In our study, we investigate this problem in the limit where the stiffness and internal damping associated with the hinge are negligible, and the mass ratio of the splitter plate is small. The experiments show that the splitter plate oscillations increase with Reynolds numbers at low values of Re, and are found to reach a saturation amplitude level at higher Re, Re>4000. This type of saturation amplitude level that appears to continue indefinitely with Re, appears to be related to the fact that there is no structural restoring force, and has been observed previously for transversely oscillating cylinders with no restorin force. In the present case, the saturation tip amplitude level can be tip to 0.45D, where D is the cylinder diameter. For this hinged-rigid splitter plate case, it is found that the splitter plate length to cylinder diameter ratio (L/D) is crucial in determining the character and magnitude of the oscillations. For small splitter plate length (L/D <= 3.0), the oscillations appear to be nearly periodic with tip amplitudes of about 0.45D nearly independent of L/D. The nondiinensional oscillation frequencies (fD/U) on the other hand are found to continuously vary with L/D from fD/U approximate to 0.2 at L/D = 1 to fD/U approximate to 0.1 at L/D = 3. As the splitter plate length is further increased beyond L/D >= 4.0, the character of the splitter plate oscillations suddenly changes. The oscillations become aperiodic with much smaller amplitudes. In this long splitter plate regime, the spectra of the oscillations become broadband, and are reminiscent of the change in character of the wake oscillations seen in the earlier fixed-rigid splitter plate case for L/D >= 5.0. In the present case of the hinged-splitter plate, the sudden transition seen as the splitter plate length (L/D) is increased from 3 to 4 may be attributed to the fact that the wake vortices are no longer able to synchronize with the plate motions for larger splitter plate lengths. Hence, as observed in other vortex-induced vibration problems, the oscillations becomeaperiodic and the amplitude reduces dramatically.

Analysis of continuum using boundary compatibility conditions of integrated force method

Static and vibration problems of an indeterminate continuum are traditionally analyzed by the stiffness method. The force method is more or less non-existent for such problems. This situation is primarily due to the incomplete state of development of the compatibility conditions which are essential for the analysis of indeterminate structures by the flexibility method. The understanding of the Compatibility Conditions (CC) has been substantially augmented. Based on the understanding of CC, a novel formulation termed the Integrated Force Method (IFM) has been established. In this paper IFM has been extended for the static and vibration analyses of a continuum. The IFM analysis is illustrated taking three examples: 1. (1) rectangular plate in flexure 2. (2) analysis of a cantilevered dam 3. (3) free vibration analysis of a beam. From the examples solved it is observed that the force response of an indeterminate continuum with mixed boundary conditions can be generated by IFM without any reference to displacements in the field or on the boundary. Displacements if required can be calculated by back substitution.

MHD flow over a moving plate in a rotating fluid with magnetic field, Hall currents and free stream velocity

The non-similar boundary layer flow of a viscous incompressible electrically conducting fluid over a moving surface in a rotating fluid, in the presence of a magnetic field, Hall currents and the free stream velocity has been studied. The parabolic partial differential equations governing the flow are solved numerically using an implicit finite-difference scheme. The Coriolis force induces overshoot in the velocity profile of the primary flow and the magnetic field reduces/removes the velocity overshoot. The local skin friction coefficient for the primary flow increases with the magnetic field, but the skin friction coefficient for the secondary flow reduces it. Also the local skin friction coefficients for the primary and secondary flows are reduced due to the Hall currents. The effects of the magnetic field, Hall currents and the wall velocity, on the skin friction coefficients for the primary and secondary flows increase with the Coriolis force. The wall velocity strongly affects the flow field. When the wall velocity is equal to the free stream velocity, the skin friction coefficients for the primary and secondary flows vanish, but this does not imply separation. (C) 2002 Published by Elsevier Science Ltd.

Unsteady free convection flow over an infinite vertical porous plate due to the combined effects of thermal and mass diffusion, magnetic field and Hall currents

The unsteady free convection flow over an infinite vertical porous plate, which moves with time-dependent velocity in an ambient fluid, has been studied. The effects of the magnetic field and Hall current are included in the analysis. The buoyancy forces arise due to both the thermal and mass diffusion. The partial differential equations governing the flow have been solved numerically using both the implicit finite difference scheme and the difference-differential method. For the steady case, analytical solutions have also been obtained. The effect of time variation on the skin friction, heat transfer and mass transfer is very significant. Suction increases the skin friction coefficient in the primary flow, and also the Nusselt and Sherwood numbers, but the skin friction coefficient in the secondary flow is reduced. The effect of injection is opposite to that of suction. The buoyancy force, injection and the Hall parameter induce an overshoot in the velocity profiles in the primary flow which changes the velocity gradient from a negative to a positive value, but the magnetic field and suction reduce this velocity overshoot.

New nine-node Lagrangian quadrilateral plate element based on Mindlin-Reissner theory using IFM

This paper presents a new nine-node Lagrangian quadrilateral plate bending element (MQP9) using the Integrated Force Method (IFM) for the analysis of thin and moderately thick plate bending problems. Three degrees of freedom: transverse displacement w and two rotations theta(x) and theta(y) are considered at each node of the element. The Mindlin-Reissner theory has been employed in the formulation which accounts the effect of shear deformation. Many standard plate bending benchmark problems have been analyzed using the new element MQP9 for various grid sizes via Integrated Force Method to estimate defections and bending moments. These results of the new element MQP9 are compared with those of similar displacement-based plate bending elements available in the literature. The results are also compared with exact solutions. It is observed that the presented new element MQP9 is free from shear locking and produced, in general, excellent results in all plate bending benchmark problems considered.

A novel free floating accelerometer force balance system for shock tunnel applications

In order to overcome the interference of the model mounting system with the external aerodynamics of the body during shock tunnel testing, a new free floating internally mountable balance system that ensures unrestrained model motion during testing has been designed, fabricated and tested. Minimal friction ball bearings are used for ensuring the free floating condition of the model during tunnel testing. The drag force acting on a blunt leading edge flat plate at hypersonic Mach number has been measured using the new balance system. Finite element model (FEM) and CFD are exhaustively used in the design as well as for calibrating the new balance system. The experimentally measured drag force on the blunt leading edge flat plate at stagnation enthalpy of 0.7 and 1.2 MJ/kg and nominal Mach number of 5.75 matches well with FEM results. The concept can also be extended for measuring all the three fundamental aerodynamic forces in short duration test facilities like free piston driven shock tunnels.

Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics

The way in which basal tractions, associated with mantle convection, couples with the lithosphere is a fundamental problem in geodynamics. A successful lithosphere-mantle coupling model for the Earth will satisfy observations of plate motions, intraplate stresses, and the plate boundary zone deformation. We solve the depth integrated three-dimensional force balance equations in a global finite element model that takes into account effects of both topography and shallow lithosphere structure as well as tractions originating from deeper mantle convection. The contribution from topography and lithosphere structure is estimated by calculating gravitational potential energy differences. The basal tractions are derived from a fully dynamic flow model with both radial and lateral viscosity variations. We simultaneously fit stresses and plate motions in order to delineate a best-fit lithosphere-mantle coupling model. We use both the World Stress Map and the Global Strain Rate Model to constrain the models. We find that a strongly coupled model with a stiff lithosphere and 3-4 orders of lateral viscosity variations in the lithosphere are best able to match the observational constraints. Our predicted deviatoric stresses, which are dominated by contribution from mantle tractions, range between 20-70 MPa. The best-fitting coupled models predict strain rates that are consistent with observations. That is, the intraplate areas are nearly rigid whereas plate boundaries and some other continental deformation zones display high strain rates. Comparison of mantle tractions and surface velocities indicate that in most areas tractions are driving, although in a few regions, including western North America, tractions are resistive. Citation: Ghosh, A., W. E. Holt, and L. M. Wen (2013), Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics.

Molecular Origin of the Intrinsic Bending Force for Helical Morphology Observed in Chiral Amphiphilic Assemblies: Concentration and Size Dependence

The formation of the helical morphology in monolayers and bilayers of chiral amphiphilic assemblies is believed to be driven at least partly by the interactions at the chiral centers of the amphiphiles. However, a detailed microscopic understanding of these interactions and their relation with the helix formation is still not clear. In this article a study of the molecular origin of the chirality-driven helix formation is presented by calculating, for the first time, the effective pair potential between a pair of chiral molecules. This effective potential depends on the relative sizes of the groups attached to the two chiral centers, on the orientation of the amphiphile molecules, and also on the distance between them. We find that for the mirror-image isomers (in the racemic modification) the minimum energy conformation is a nearly parallel alignment of the molecules. On the other hand, the same for a pair of molecules of one kind of enantiomer favors a tilt angle between them, thus leading to the formation of a helical morphology of the aggregate. The tilt angle is determined by the size of the groups attached to the chiral centers of the pair of molecules considered and in many cases predicted it to be close to 45 degrees. The present study, therefore, provides a molecular origin of the intrinsic bending force, suggested by Helfrich (J. Chem. Phys. 1986, 85, 1085-1087), to be responsible for the formation of helical structure. This effective potential may explain many of the existing experimental results, such as the size and the concentration dependence of the formation of helical morphology. It is further found that the elastic forces can significantly modify the pitch predicted by the chiral interactions alone and that the modified real pitch is close to the experimentally observed value. The present study is expected to provide a starting point for future microscopic studies.