Regularized numerical integration of multibody dynamics with the generalized alpha method*

This paper discusses the consistent regularization property of the generalized α method when applied as an integrator to an initial value high index and singular differential-algebraic equation model of a multibody system. The regularization comes from within the discretization itself and the discretization remains consistent over the range of values the regularization parameter may take. The regularization involves increase of the smallest singular values of the ill-conditioned Jacobian of the discretization and is different from Baumgarte and similar techniques which tend to be inconsistent for poor choice of regularization parameter. This regularization also helps where pre-conditioning the Jacobian by scaling is of limited effect, for example, when the scleronomic constraints contain multiple closed loops or singular configuration or when high index path constraints are present. The feed-forward control in Kane's equation models is additionally considered in the numerical examples to illustrate the effect of regularization. The discretization presented in this work is adopted to the first order DAE system (unlike the original method which is intended for second order systems) for its A-stability and same order of accuracy for positions and velocities.

L2-stability of nonstationary feedback systems: Frequency-domain criteria

A frequency-domain positivity condition is derived for linear time-varying operators in2and is used to develop2stability criteria for linear and nonlinear feedback systems. These criteria permit the use of a very general class of operators in2with nonstationary kernels, as multipliers. More specific results are obtained by using a first-order differential operator with a time-varying coefficient as multiplier. Finally, by employing periodic multipliers, improved stability criteria are derived for the nonlinear damped Mathieu equation with a forcing function.

Multipliers of Sobolev spaces

Let Wm,p denote the Sobolev space of functions on Image n whose distributional derivatives of order up to m lie in Lp(Image n) for 1 less-than-or-equals, slant p less-than-or-equals, slant ∞. When 1 < p < ∞, it is known that the multipliers on Wm,p are the same as those on Lp. This result is true for p = 1 only if n = 1. For, we prove that the integrable distributions of order less-than-or-equals, slant1 whose first order derivatives are also integrable of order less-than-or-equals, slant1, belong to the class of multipliers on Wm,1 and there are such distributions which are not bounded measures. These distributions are also multipliers on Lp, for 1 < p < ∞. Moreover, they form exactly the multiplier space of a certain Segal algebra. We have also proved that the multipliers on Wm,l are necessarily integrable distributions of order less-than-or-equals, slant1 or less-than-or-equals, slant2 accordingly as m is odd or even. We have obtained the multipliers from L1(Image n) into Wm,p, 1 less-than-or-equals, slant p less-than-or-equals, slant ∞, and the multiplier space of Wm,1 is realised as a dual space of certain continuous functions on Image n which vanish at infinity.

Finite element analysis of curved anisotropic laminated hollow bars

Curved hollow bars of laminated anisotropic construction are used as structural members in many industries. They are used in order to save weight without loss of stiffness in comparison with solid sections. In this paper are presented the details of the development of the stiffness matrices of laminated anisotropic curved hollow bars under line member assumptions for two typical sections, circular and square. They are 16dof elements which make use of one-dimensional first-order Hermite interpolation polynomials for the description of assumed displacement state. Problems for which analytical or other solutions are available are first solved using these elements. Good agreement was found between the results. In order to show the capability of the element, application is made to carbon fibre reinforced plastic layered anisotropic curved hollow bars.

Ultrasonic velocities and thermal expansion coefficients of amorphous Se80Te20 and Se90Te10 alloys near glass transitions

Precise measurements of the ultrasonic velocities and thermal expansivities of amorphous Se80Te20 and Se90Te10 alloys are reported near the glass transition. The samples are produced by liquid quenching. The longitudinal and transverse velocities are measured at 10 MHz frequency using the McSkimin pulse superposition technique. The thermal expansivities,agr, are measured using a three-terminal capacitance bridge. Theagr-values show a sharp maximum near the glass transition temperature,T g. The ultrasonic velocities also show a large temperature derivative, dV/dT nearT g. The data are discussed in terms of existing theories of the glass transition. The continuous change inagr shows that the glass transition is not a first-order transition, as suggested by some theories. The samples are found to be deformed by small loads nearT g. The ultrasonic velocities and dV/dT have contributions arising from this deformation.

Finite element analysis of laminated anisotropic beams of bimodulus materials

This paper presents finite element analysis of laminated anisotropic beams of bimodulus materials. The finite element has 16 d.o.f. and uses the displacement field in terms of first order Hermite interpolation polynomials. As the neutral axis position may change from point to point along the length of the beam, an iterative procedure is employed to determine the location of zero strain points along the length. Using this element some problems of laminated beams of bimodulus materials are solved for concentrated loads/moments perpendicular and parallel to the layering planes as well as combined loads.

Intrinsic two-phonon-exchange mechanism of superconductivity and enhancement of Tc

It is shown that the intrinsic two-phonon terms occurring in first order in the electron-phonon interaction Hamiltonian can give rise to (i) an essential doubling of the interaction phase space (BCS cutoff) and (ii) an attractive pairing interaction proportional to the phonon occupation numbers. This suggests a possible enhancement of the superconductive transition temperature in the presence of high-frequency acoustic field.

Thermal decomposition of titanyl oxalates II : Kinetics of decomposition of barium titanyl oxalate

Kinetics of the thermal decomposition of barium titanyl oxalate have been studied. Decomposition of the anhydrous oxalate is complex and deceleratory throughout. Kinetics of decomposition of the intermediate carbonate Ba2Ti2O5CO3 is greatly influenced by the thermal effects during its formation. The sigmoidal (α, t) curves obey a power law equation followed by first order decay. Presence of carbon in the vacuum prepared carbonate has a strong deactivating effect. Decomposition of the carbonate is accompanied by growth in particle size of the product, barium titanate.

Dynamics of a class of social interaction systems

A mathematical model of social interaction in the form of two coupler! first-order non-linear differential equations, forms the topic of this study. This non-conservative model io representative of such varied social interaction problems as coexisting sub-populations of two different species, arms race between two rival countries and the like. Differential transformation techniques developed elsewhere in the literature are seen to be effective tools of dynamic analysis of this non-linear non-conservative mode! of social interaction process.

Liquid phase hydrogenation of aniline in a trickle bed reactor

A rate equation is developed for the liquid phase hydrogenation of aniline over cylindrical catalyst pellets of 30% nickel deposited on clay in a trickle bed reactor. The equation takes into account external and internal diffusional limitations, and describes the experimental data adequately. The hydrogenation reaction is first order with respect to hydrogen and zero order with respect to aniline. Effectiveness factors are in the range 0.003-0.03. Apparent activation energy of the reaction is 12.7 kcal/mol and true activation energy is 39.6 kcal/mol.

Approximation of the perturbation equations of a quasi-linear hyperbolic system in the neighborhood of a bicharacteristic

In 1956 Whitham gave a nonlinear theory for computing the intensity of an acoustic pulse of an arbitrary shape. The theory has been used very successfully in computing the intensity of the sonic bang produced by a supersonic plane. [4.] derived an approximate quasi-linear equation for the propagation of a short wave in a compressible medium. These two methods are essentially nonlinear approximations of the perturbation equations of the system of gas-dynamic equations in the neighborhood of a bicharacteristic curve (or rays) for weak unsteady disturbances superimposed on a given steady solution. In this paper we have derived an approximate quasi-linear equation which is an approximation of perturbation equations in the neighborhood of a bicharacteristic curve for a weak pulse governed by a general system of first order quasi-linear partial differential equations in m + 1 independent variables (t, x1,…, xm) and derived Gubkin's result as a particular case when the system of equations consists of the equations of an unsteady motion of a compressible gas. We have also discussed the form of the approximate equation describing the waves propagating upsteam in an arbitrary multidimensional transonic flow.

On the computation of two-dimensional compressible flow fields by the modified local stability scheme

The modified local stability scheme is applied to several two-dimensional problems—blunt body flow, regular reflection of a shock and lambda shock. The resolution of the flow features obtained by the modified local stability scheme is found to be better than that achieved by the other first order schemes and almost identical to that achieved by the second order schemes incorporating artificial viscosity. The scheme is easy for coding, consumes moderate amount of computer storage and time. The scheme can be advantageously used in place of second order schemes.

Kinetics of formalization of poly(vinyl acetate)

Kinetic information on the formation of poly(vinyl formal) by the reaction of poly(vinyl acetate) and formaldehyde in presence of aqueous acid has been derived from the spectroscopic analysis of polymer samples after different periods of reaction. The hydroxyl content of poly(vinyl formal) is found to be nearly independent of reaction time and only slightly affected by temperature while the fall of acetate content and the increase in formal content are most rapid in the initial period and are largely influenced by temperature. The rate expression formulated on the assumption that the formalization reaction is of first order with respect to both poly(vinyl acetate) and formaldehyde explains the observed variation of polymer composition with reaction time. The activation energy for the reaction is found to be 17.3 kcal/mol.

Phase diagram of vortex matter in layered superconductors with tilted columnar pinning centers

We study the vortex matter phase diagram of a layered superconductor in the presence of columnar pinning defects, tilted with respect to the normal to the layers. We use numerical minimization of the free energy written as a functional of the time-averaged vortex density of the Ramakrishnan-Yussouff form, supplemented by the appropriate pinning potential. We study the case where the pin density is smaller than the areal vortex density. At lower pin concentrations, we find, for temperatures of the order of the melting temperature of the unpinned lattice, a Bose glass type phase which at lower temperatures converts, via a first-order transition, to a Bragg glass, while, at higher temperatures, it crosses over to an interstitial liquid. At somewhat higher concentrations, no transition to a Bragg glass is found even at the lowest temperatures studied. While qualitatively the behavior we find is similar to that obtained using the same procedures for columnar pins normal to the layers, there are important and observable quantitative differences, which we discuss.

Automata and logics over finitely varying functions

We extend some of the classical connections between automata and logic due to Büchi (1960) [5] and McNaughton and Papert (1971) [12] to languages of finitely varying functions or “signals”. In particular, we introduce a natural class of automata for generating finitely varying functions called View the MathML source’s, and show that it coincides in terms of language definability with a natural monadic second-order logic interpreted over finitely varying functions Rabinovich (2002) [15]. We also identify a “counter-free” subclass of View the MathML source’s which characterise the first-order definable languages of finitely varying functions. Our proofs mainly factor through the classical results for word languages. These results have applications in automata characterisations for continuously interpreted real-time logics like Metric Temporal Logic (MTL) Chevalier et al. (2006, 2007) [6] and [7].