Large deformation finite element analysis of laminated circular composite plates

A 48 d.o.f., four-noded quadrilateral laminated composite shell finite element is particularised to a sector finite element and is used for the large deformation analysis of circular composite laminated plates. The strain-displacement relationships for the sector element are obtained by reducing those of the quadrilateral shell finite element by substituting proper values for the geometric parameters. Subsequently, the linear and tangent stiffness matrices are formulated using conventional methods. The Newton-Raphson method is employed as the nonlinear solution technique. The computer code developed is validated by solving an isotropic case for which results are available in the literature. The method is then applied to solve problems of cylindrically orthotropic circular plates. Some of the results of cylindrically orthotropic case are compared with those available in the literature. Subsequently, application is made to the case of laminated composite circular plates having different lay-up schemes. The computer code can handle symmetric/unsymmetric lay-up schemes. The large displacement analysis is useful in estimating the damage in composite plates caused by low-velocity impact.

Crystal structure of Boc-LAla-Delta Phe-Delta Phe-Delta Phe-Delta Phe-NHMe: a left-handed helical peptide

alpha,beta-Dehydrophenylalanine residues constrain the peptide backbone to beta-bend conformation. A pentapeptide containing four consecutive (Delta Phe) residues has been synthesised and crystallised. The peptide Boc-LAla-Delta Phe-Delta Phe-Delta Phe-Delta Phe-NHMe (C45H46N6O7, MW = 782.86) was crystallised from an acetonitrile/methanol mixture. The crystal belongs to the orthorhombic space group P2(1)2(1)2(1) With a = 19.455(6), b = 20.912(9), c = 11.455(4) Angstrom and Z = 4. The X-ray (MoKalpha, lambda = 0.7107 Angstrom) intensity data were collected using the Rigaku-AFC7 diffractrometer. The crystal structure was determined by direct methods and refined using the least-squares technique, R = 8.41% for 1827 reflections with \F-o\ > 4 sigma\F-o\. The molecule contains the largest stretch of consecutive dehydrophenylalanine residues whose crystal structure has been determined so far. The peptide adopts left-handed 3(10)-helical conformation despite the presence of LAla at the N-terminus. The mean phi, psi values, averaged across the last four residues are 56.8 degrees and 17.5 degrees, respectively. There are four 4-->1 intramolecular hydrogen bonds, characteristic of the 3(10)-helix. In the crystal each molecule interacts with four crystallographically symmetric molecules with one hydrogen bond each.

Lopsidedness in WHISP galaxies I. Rotation curves and kinematic lopsidedness

The frequently observed lopsidedness of the distribution of stars and gas in disc galaxies is still considered as a major problem in galaxy dynamics. It is even discussed as an imprint of the formation history of discs and the evolution of baryons in dark matter haloes. Here, we analyse a selected sample of 70 galaxies from the Westerbork Hi Survey of Spiral and Irregular Galaxies. The Hi data allow us to follow the morphology and the kinematics out to very large radii. In the present paper, we present the rotation curves and study the kinematic asymmetry. We extract the rotation curves of the receding and approaching sides separately and show that the kinematic behaviour of disc galaxies can be classified into five different types: symmetric velocity fields where the rotation curves of the receding and approaching sides are almost identical; global distortions where the rotation velocities of the receding and approaching sides have an offset that is constant with radius; local distortions leading to large deviations in the inner and negligible deviations in the outer parts (and vice versa); and distortions that divide the galaxies into two kinematic systems that are visible in terms of the different behaviour of the rotation curves of the receding and approaching sides, which leads to a crossing and a change in side. The kinematic lopsidedness is measured from the maximum rotation velocities, averaged over the plateau of the rotation curves. This gives a good estimate of the global lopsidedness in the outer parts of the sample galaxies. We find that the mean value of the perturbation parameter denoting the lopsided potential as obtained from the kinematic data is 0.056. Altogether, 36% of the sample galaxies are globally lopsided, which can be interpreted as the disc responding to a halo that was distorted by a tidal encounter. In Paper II, we study the morphological lopsidedness of the same sample of galaxies.

Velocity distribution for a dilute vibrated granular material

The velocity distribution for a vibrated granular material is determined in the dilute limit where the frequency of particle collisions with the vibrating surface is large compared to the frequency of binary collisions. The particle motion is driven by the source of energy due to particle collisions with the vibrating surface, and two dissipation mechanisms-inelastic collisions and air drag-are considered. In the latter case, a general form for the drag force is assumed. First, the distribution function for the vertical velocity for a single particle colliding with a vibrating surface is determined in the limit where the dissipation during a collision due to inelasticity or between successive collisions due to drag is small compared to the energy of a particle. In addition, two types of amplitude functions for the velocity of the surface, symmetric and asymmetric about zero velocity, are considered. In all cases, differential equations for the distribution of velocities at the vibrating surface are obtained using a flux balance condition in velocity space, and these are solved to determine the distribution function. It is found that the distribution function is a Gaussian distribution when the dissipation is due to inelastic collisions and the amplitude function is symmetric, and the mean square velocity scales as [[U-2](s)/(1 - e(2))], where [U-2](s) is the mean square velocity of the vibrating surface and e is the coefficient of restitution. The distribution function is very different from a Gaussian when the dissipation is due to air drag and the amplitude function is symmetric, and the mean square velocity scales as ([U-2](s)g/mu(m))(1/(m+2)) when the acceleration due to the fluid drag is -mu(m)u(y)\u(y)\(m-1), where g is the acceleration due to gravity. For an asymmetric amplitude function, the distribution function at the vibrating surface is found to be sharply peaked around [+/-2[U](s)/(1-e)] when the dissipation is due to inelastic collisions, and around +/-[(m +2)[U](s)g/mu(m)](1/(m+1)) when the dissipation is due to fluid drag, where [U](s) is the mean velocity of the surface. The distribution functions are compared with numerical simulations of a particle colliding with a vibrating surface, and excellent agreement is found with no adjustable parameters. The distribution function for a two-dimensional vibrated granular material that includes the first effect of binary collisions is determined for the system with dissipation due to inelastic collisions and the amplitude function for the velocity of the vibrating surface is symmetric in the limit delta(I)=(2nr)/(1 - e)much less than 1. Here, n is the number of particles per unit width and r is the particle radius. In this Limit, an asymptotic analysis is used about the Limit where there are no binary collisions. It is found that the distribution function has a power-law divergence proportional to \u(x)\((c delta l-1)) in the limit u(x)-->0, where u(x) is the horizontal velocity. The constant c and the moments of the distribution function are evaluated from the conservation equation in velocity space. It is found that the mean square velocity in the horizontal direction scales as O(delta(I)T), and the nontrivial third moments of the velocity distribution scale as O(delta(I)epsilon(I)T(3/2)) where epsilon(I) = (1 - e)(1/2). Here, T = [2[U2](s)/(1 - e)] is the mean square velocity of the particles.

A differential geometric method for kinematic analysis of two- and three-degree-of-freedom rigid body motions

In this paper, we present a novel differential geometric characterization of two- and three-degree-of-freedom rigid body kinematics, using a metric defined on dual vectors. The instantaneous angular and linear velocities of a rigid body are expressed as a dual velocity vector, and dual inner product is defined on this dual vector, resulting in a positive semi-definite and symmetric dual matrix. We show that the maximum and minimum magnitude of the dual velocity vector, for a unit speed motion, can be obtained as eigenvalues of this dual matrix. Furthermore, we show that the tip of the dual velocity vector lies on a dual ellipse for a two-degree-of-freedom motion and on a dual ellipsoid for a three-degree-of-freedom motion. In this manner, the velocity distribution of a rigid body can be studied algebraically in terms of the eigenvalues of a dual matrix or geometrically with the dual ellipse and ellipsoid. The second-order properties of the two- and three-degree-of-freedom motions of a rigid body are also obtained from the derivatives of the elements of the dual matrix. This results in a definition of the geodesic motion of a rigid body. The theoretical results are illustrated with the help of a spatial 2R and a parallel three-degree-of-freedom manipulator.

Aerodynamics of delta wings with flaps at hypersonic speeds

Control surface effectiveness is an important parameter for any aeroplane. For a hypersonic aircraft, though the power required to operate the flaps is determined by low speed flying conditions, it is imperative to know the effect of flaps at hypersonic speeds. Hence, studies have been done on this topic by aerodynamicists for over 40 years. In spite of this, only a limited data is available in the literature on this subject. This paper discusses the experimental study of the effect of sweep on the aerodynamic characteristics of thin slab delta wings with flaps at hypersonic speeds. For the purpose of this investigation, a novel special thin six-component balance, which has a thickness of 4mm and can be housed inside wings with 8mm thickness, has been designed. The wings had a sweep of 76degrees, 70degrees and 65degrees, t/c of 0.053 and flaps with 12% of wing area and 12% of wing chord. Testing were done at Mach 8.2, Re number of 2.13 x 10(6) (based on chord), from alpha = -12degrees to 12degrees and flap angle of 20degrees, 30degrees and 40degrees. Separation lengths, measured from Schlieren pictures, clearly show that there is 'no appreciable' effect of sweep on them. Also, using a simple local flow field calculation, the separation has been identified to be transitional in nature. These features of separation reflect in the force data. Because of the small separation length, the flaps (inspite of their small size) were very effective in generating additional C-N, C-M and C-l, which increased with increase in flap angle. In general, the C-N, C-M and X-CP were unaffected by sweep for symmetric flap deflection at positive incidences and asymmetric flap case, For symmetric flap case at negative incidences, only C-N was not influenced by the sweep but C-M decreased and X-CP moved upstream as the sweep is decreased, The wing with lower sweep produces higher CA and lower (L/D)(max) for both symmetric and asymmetric flaps. The rolling moment and adverse yaw increased with decrease in sweep for asymmetric flap deflection. Newtonian theory is shown to be incapable of predicting the effect of sweep on C-l, C-n and on the incremental values of C-N, C-M and C-A. In conclusion, it can be said that a small flap is generally adequate for hypersonic aeroplanes provided they operate at altitudes where transitional and turbulent separation can be expected to occur. This would make the flaps effective and thus enable ample control authority.

Receiver-only optimized Vector Quantization for fading channels

This paper considers the design and analysis of a filter at the receiver of a source coding system to mitigate the excess distortion caused due to channel errors. The index output by the source encoder is sent over a fading discrete binary symmetric channel and the possibly incorrect received index is mapped to the corresponding codeword by a Vector Quantization (VQ) decoder at the receiver. The output of the VQ decoder is then processed by a receive filter to obtain an estimate of the source instantiation. The distortion performance is analyzed for weighted mean square error (WMSE) and the optimum receive filter that minimizes the expected distortion is derived for two different cases of fading. It is shown that the performance of the system with the receive filter is strictly better than that of a conventional VQ and the difference becomes more significant as the number of bits transmitted increases. Theoretical expressions for an upper and lower bound on the WMSE performance of the system with the receive filter and a Rayleigh flat fading channel are derived. The design of a receive filter in the presence of channel mismatch is also studied and it is shown that a minimax solution is the one obtained by designing the receive filter for the worst possible channel. Simulation results are presented to validate the theoretical expressions and illustrate the benefits of receive filtering.

Effect of convection on domain growth during demixing transitions in fluids

Current analytical work on the effect of convection on the late stages of spinodal decomposition in liquids is briefly described. The morphology formed during the spinodal decomposition process depends on the relative composition of the two species. Droplet spinodal decomposition occurs when the concentration of one of the species is small. Convective transport has a significant effect on the scaling laws in the late-stage coarsening of droplets in translational or shear flows. In addition, convective transport could result in an attractive interaction between non-Brownian droplets which could lead to coalescence. The effect of convective transport for the growth of random interfaces in a near-symmetric quench was analysed using an area distribution function, which gives the distribution of surface area of the interface in curvature space. It was found that the curvature of the interface decreases proportional to time t in the late stages of spinodal decomposition, and the surface area also decreases proportional to t.

Investigation of short-time isomerization dynamics in p-nitroazobenzene from resonance Raman intensity analysis

Resonance Raman (RR) spectra are presented for p-nitroazobenzene dissolved in chloroform using 18 excitation Wavelengths, covering the region of (1)(n --> pi*) electronic transition. Raman intensities are observed for various totally symmetric fundamentals, namely, C-C, C-N, N=N, and N-O stretching vibrations, indicating that upon photoexcitation the excited-state evolution occurs along all of these vibrational coordinates. For a few fundamentals, interestingly, in p-nitroazobenzene, it is observed that the RR intensities decrease near the maxima of the resonant electronic (1)(n --> pi*) transition. This is attributed to the interference from preresonant scattering due to the strongly allowed (1)(pi --> pi*) electronic transition. The electronic absorption spectrum and the absolute Raman cross section for the nine Franck-Condon active fundamentals of p-nitroazobenzene have been successfully modeled using Heller's time-dependent formalism for Raman scattering. This employs harmonic description of the lowest energy (1)(n --> pi*) potential energy surface. The short-time isomerization dynamics is then examined from a priori knowledge of the ground-state normal mode descriptions of p-nitroazobenzene to convert the wave packet motion in dimensionless normal coordinates to internal coordinates. It is observed that within 20 fs after photoexcitation in p-nitroazobenzene, the N=N and C-N stretching vibrations undergo significant changes and the unsubstituted phenyl ring and the nitro stretching vibrations are also distorted considerably.

Persistent passive hopping and juggling is possible even with plastic collisions

We describe simple one-dimensional models of passive (no energy input, no control), generally dissipative, vertical hopping and one-ball juggling. The central observation is that internal passive system motions can conspire to eliminate collisions in these systems. For hopping, two point masses are connected by a spring and the lower mass has inelastic collisions with the ground. For juggling, a lower point-mass hand is connected by a spring to the ground and an upper point-mass ball is caught with an inelastic collision and then re-thrown into gravitational free flight. The two systems have identical dynamics. Despite inelastic collisions between non-zero masses, these systems have special symmetric energy-conserving periodic motions where the collision is at zero relative velocity. Additionally, these special periodic motions have a non-zero sized, one-sided region of attraction on the higher-energy side. For either very large or very small mass ratios, the one-sided region of attraction is large. These results persist for mildly non-linear springs and non-constant gravity. Although non-collisional damping destroys the periodic motions, small energy injection makes the periodic motions stable, with a two-sided region of attraction. The existence of such special energy conserving solutions for hopping and juggling points to possibly useful strategies for both animals and robots. The lossless motions are demonstrated with a table-top experiment.

Photochemical behaviour of 4-styrylcoumarin polymorphs in the crystalline state

Addition of a small amount of coumarin during crystallization produces new polymorphic modifications in 4-styrylcoumarin and 4-(3-fluorostyryl)coumarin, which are photolabile. Interestingly, upon irradiation polymorphic modification of 4-(3-fluorostyryl)coumarin produces a mirror-symmetric photodimer in contrast to the centrosymmetric photodimer obtained without addition of coumarin during crystallization.

Turbulent near-wake studies behind an infinitely swept wing

The three-dimensional asymmetric turbulent near-Rake behind an infinitely swept wing with GAW(2) airfoil has been investigated at low speeds. The near-wake in the present study is asymmetric because the boundary layers on the top and bottom surfaces of the model develop under different streamwise pressure gradients. Distributions of mean velocity, three turbulent normal stresses, and two important Reynolds shear stresses have been measured using hot-wire anemometry. The profiles of mean velocity and Reynolds shear stress exhibit asymmetry near the trailing edge and seem to have become symmetric within a short distance of 6 trailing edge momentum thicknesses. Results of computation using K-epsilon turbulence model with a simple scheme to predict the near-wake behind the swept wing have also been presented and compared with the experimental data. The agreement of the predicted mean How development with the experiment is fair considering the simplicity of the scheme.

Motion of a charged particle in superposed Heliotron and biconical cusp fields

We have studied the behaviour of a charged particle in an axially symmetric magnetic field having a neutral point, so as to find a possibility of confining a charged particle in a thermonuclear device. In order to study the motion we have reduced a three-dimensional motion to a two-dimensional one by introducing a fictitious potential. Following Schmidt we have classified the motion, as an ‘off-axis motion’ and ‘encircling motion’ depending on the behaviour of this potential. We see that the particle performs a hybrid type of motion in the negative z-axis, i.e. at some instant it is in ‘off-axis motion’ while at another instant it is in ‘encircling motion’. We have also solved the equation of motion numerically and the graphs of the particle trajectory verify our analysis. We find that in most of the cases the particle is contained. The magnetic moment is found to be moderately adiabatic.

Engineered Biferroic 0.7Pb(Mg1/3Nb2/3)O3-0.3PbTiO3/La0.6Sr0.4MnO3 Epitaxial Superlattices

Symmetric and asymmetric superlattices (SLs) composed of ferromagnetic La0.6Sr0.4MnO3 (LSMO) and ferroelectric 0.7Pb(Mg1/3Nb2/3)O3 – 0.3PbTiO3 (PMN-PT) with different periodicities have been fabricated on LaNiO3 (LNO) coated LaAlO3 (100) (LAO) substrates by pulsed laser ablation deposition. Structural, ferromagnetic and ferroelectric properties have been studied for all the SLs. All the heterostructures exhibited good ferromagnetic response over a wide range of temperatures (10K – 300K), whereas only the asymmetric SLs exhibited reasonably good ferroelectric behaviour. Ferromagnetic and ferroelectric hysteresis loops observed in the asymmetric SLs confirmed their biferroic nature. Studies were conducted towards understanding the influence of LSMO layers on the electrical responses of the heterostructures. Absence of ferroelectricity in the symmetric SL structures has been attributed to their high leakage characteristics. Strong influence of an applied magnetic field of 1.2T was observed on the ferroelectric properties of the asymmetric SLs. The effect of magnetic field on the ferroelectric properties of the SLs indicated possibility of strong interfacial effect.

Dynamics and transport properties of Kondo insulators

A many-body theory of paramagnetic Kondo insulators is described, focusing specifically on single-particle dynamics, scattering rates, dc transport and optical conductivities. This is achieved by development of a non-perturbative local moment approach to the symmetric periodic Anderson model within the framework of dynamical mean-field theory. Our natural focus is the strong-coupling, Kondo lattice regime, in particular the resultant 'universal' scaling behaviour in terms of the single, exponentially small low-energy scale characteristic of the problem. Dynamics/transport on all relevant (ω, T)-scales are considered, from the gapped/activated behaviour characteristic of the low-temperature insulator through to explicit connection to single-impurity physics at high ω and/or T; and for optical conductivities emphasis is given to the nature of the optical gap, the temperature scale responsible for its destruction and the consequent clear distinction between indirect and direct gap scales. Using scaling, explicit comparison is also made to experimental results for dc transport and optical conductivities of Ce3Bi4Pt3, SmB6 and YbB12. Good agreement is found, even quantitatively; and a mutually consistent picture of transport and optics results.