Analysis of co-current hydraircooling of food products in bulk

A mathematical model is developed to describe the hydraircooling process when the water and air are flowing in the same direction. The governing equations for the simultaneous heat and mass transfer are solved using finite-difference numerical methods. The half cooling time of the food products is correlated as a function of the dimensionless process parameters. It is observed that a process time of approximately double the half cooling time will result in the food products attaining almost a steady state. The process times of the bulk hydraircooling process and the bulk air precooling process are compared.

Electrostatic Field Analysis of Electrooptic Devices

We describe a Finite Difference Method for the determination of the electrostatic field in a multilayered electrooptic device. The Laplace equation is solved, assuming a suitable closed area, by taking into account the different permittivities of the various layers. The effect of a higher permittivity in the guiding layer has been explicitly considered. As a practical example, we calculate the phase shift of a guided optical wave within an electrooptic modulator. A review of the various methods in use for the field analysis is given. Some criteria for the selection of the appropriate method are also mentioned.

Analysis of heat and mass transfer during bulk hydraircooling of spherical food products

Hydraircooling is a technique used for precooling food products. In this technique chilled water is sprayed over the food products while cold unsaturated air is blown over them. Hydraircooling combines the advantages of both air- and hydrocooling. The present study is concerned with the analysis of bulk hydraircooling as it occurs in a package filled with several layers of spherical food products with chilled water sprayed from the top and cold unsaturated air blown from the bottom. A mathematical model is developed to describe the hydrodynamics and simultaneous heat and mass transfer occurring inside the package. The non-dimensional governing equations are solved using the finite difference numerical methods. The results are presented in the form of time-temperature charts. A correlation is obtained to calculate the process time in terms of the process parameters.

Method of Analysis if Uniaxial Deformation Data

A comprehensive scheme for analysing uniaxial deformation data, taking into account the finite stiffness of the testing machine is presented. Equations relevant to tension and stress relaxation tests carried out under cross head speed control, and to creep testing under constant load, are described. For the first two cases, the implications of not using gauge length extensometry but relying upon cross head displacement for inferring specimen extension, and the role of uncertainty in machine stiffness are also examined. The final section touches upon the extension of the present scheme to account for specimen anelasticity.

Coupled electrostatic-elastic analysis for topology optimization using material interpolation

In this paper, we present a novel analytical formulation for the coupled partial differential equations governing electrostatically actuated constrained elastic structures of inhomogeneous material composition. We also present a computationally efficient numerical framework for solving the coupled equations over a reference domain with a fixed finite-element mesh. This serves two purposes: (i) a series of problems with varying geometries and piece-wise homogeneous and/or inhomogeneous material distribution can be solved with a single pre-processing step, (ii) topology optimization methods can be easily implemented by interpolating the material at each point in the reference domain from a void to a dielectric or a conductor. This is attained by considering the steady-state electrical current conduction equation with a `leaky capacitor' model instead of the usual electrostatic equation. This formulation is amenable for both static and transient problems in the elastic domain coupled with the quasi-electrostatic electric field. The procedure is numerically implemented on the COMSOL Multiphysics (R) platform using the weak variational form of the governing equations. Examples have been presented to show the accuracy and versatility of the scheme. The accuracy of the scheme is validated for the special case of piece-wise homogeneous material in the limit of the leaky-capacitor model approaching the ideal case.

Fracture Analysis of Concrete Structural Components Accounting for Tension Softening Effect

This paper presents methodologies for fracture analysis of concrete structural components with and without considering tension softening effect. Stress intensity factor (SIF) is computed by using analytical approach and finite element analysis. In the analytical approach, SW accounting for tension softening effect has been obtained as the difference of SIP obtained using linear elastic fracture mechanics (LEFM) principles and SIP due to closing pressure. Superposition principle has been used by accounting for non-linearity in incremental form. SW due to crack closing force applied on the effective crack face inside the process zone has been computed using Green's function approach. In finite element analysis, the domain integral method has been used for computation of SIR The domain integral method is used to calculate the strain energy release rate and SIF when a crack grows. Numerical studies have been conducted on notched 3-point bending concrete specimen with and without considering the cohesive stresses. It is observed from the studies that SW obtained from the finite element analysis with and without considering the cohesive stresses is in good agreement with the corresponding analytical value. The effect of cohesive stress on SW decreases with increase of crack length. Further, studies have been conducted on geometrically similar structures and observed that (i) the effect of cohesive stress on SW is significant with increase of load for a particular crack length and (iii) SW values decreases with increase of tensile strength for a particular crack length and load.

Surface-ray tracing on hybrid surfaces of revolution for UTD mutual coupling analysis

An angle invariance property based on Hertz's principle of particle dynamics is employed to facilitate the surface-ray tracing on nondevelopable hybrid quadric surfaces of revolution (h-QUASOR's). This property, when used in conjunction with a Geodesic Constant Method, yields analytical expressions for all the ray-parameters required in the UTD formulation. Differential geometrical considerations require that some of the ray-parameters (defined heuristically in the UTD for the canonical convex surfaces) be modified before the UTD can be applied to such hybrid surfaces. Mutual coupling results for finite-dimensional slots have been presented as an example on a satellite launch vehicle modeled by general paraboloid of revolution and right circular cylinder.

Analysis of a cracked lug loaded by an interference fit pin

In the present study, a lug joint fitted with an interference fit (oversized) pin is considered with radial through cracks situated at diametrically opposite points perpendicular to the loading direction. A finite element contact stress algorithm is developed with linear elastic assumptions to deal with varying partial contact/separation at the pin-plate interface using a marching solution. Stress Intensity Factor (SIF) at the crack tips is evaluated using the Modified Crack Closure Integral (MCCI) method. The effect of change in crack length and edge distance on the load-contact relation, SIFs and stress distributions are studied. A rigorous plane stress elasticity solution of the pin-plate interface at the crack mouth confirmed the existence of the stress concentration leading to a local peak in the radial stress at the crack mouth and provided a method of estimating it quantitatively. Copyright (C) 1996 Elsevier Science Ltd.

Elasto-plastic contact stress analysis of joints subjected to cyclic loading

Rigorous elastic-plastic finite element analysis of joints subjected to cyclic loading is carried out. An incremental-iterative algorithm is developed in a modular form combining elasto-plastic material behaviour and contact stress analysis. For the case of the interference fit, the analysis sequentially carries out insertion of the pin and application of the load on the joint, covering possible initiation of separation (and/or yielding) and progressively the receding/advancing contact at the pin-plate interface. Deformations of both the plate and the pin are considered in the analysis. Numerical examples are presented for the case of an interference fit pin in a large plate under remote cyclic tension, and for an interference fit pin lug joint subjected to cyclic loading. A detailed study is carried out for the latter problem considering the effect of change in contact/separation at the pin-plate interface on local stresses, strains and redistribution of these stresses with the spread of a plastic zone. The results of the study are a useful input for the estimation of the fatigue life of joints. Copyright (C) 1996 Elsevier Science Ltd

Parallel computing concepts and methods for floquet analysis of helicopter trim and stability

Floquet analysis is widely used for small-order systems (say, order M < 100) to find trim results of control inputs and periodic responses, and stability results of damping levels and frequencies, Presently, however, it is practical neither for design applications nor for comprehensive analysis models that lead to large systems (M > 100); the run time on a sequential computer is simply prohibitive, Accordingly, a massively parallel Floquet analysis is developed with emphasis on large systems, and it is implemented on two SIMD or single-instruction, multiple-data computers with 4096 and 8192 processors, The focus of this development is a parallel shooting method with damped Newton iteration to generate trim results; the Floquet transition matrix (FTM) comes out as a byproduct, The eigenvalues and eigenvectors of the FTM are computed by a parallel QR method, and thereby stability results are generated, For illustration, flap and flap-lag stability of isolated rotors are treated by the parallel analysis and by a corresponding sequential analysis with the conventional shooting and QR methods; linear quasisteady airfoil aerodynamics and a finite-state three-dimensional wake model are used, Computational reliability is quantified by the condition numbers of the Jacobian matrices in Newton iteration, the condition numbers of the eigenvalues and the residual errors of the eigenpairs, and reliability figures are comparable in both the parallel and sequential analyses, Compared to the sequential analysis, the parallel analysis reduces the run time of large systems dramatically, and the reduction increases with increasing system order; this finding offers considerable promise for design and comprehensive-analysis applications.

Role of elasto-plastic analysis under cyclic loading in fatigue crack growth studies

Linear Elastic Fracture Mechanics (LEFM) has been widely used in the past for fatigue crack growth studies, but this is acceptable only in situations which are within small scale yielding (SSY). In many practical structural components, conditions of SSY could be violated and one has to look for fracture criteria based on elasto-plastic analysis. Crack closure phenomenon, one of the most striking discoveries based on inelastic deformations during crack growth, has significant effect on fatigue crack growth rate. Numerical simulation of this phenomenon is computationally intensive and involved but has been successfully implemented. Stress intensity factors and strain energy release rates lose their meaning, J-integral (or its incremental) values are applicable only in specific situations, whereas alternate path independent integrals have been proposed in the literature for use with elasto-plastic fracture mechanics (EPFM) based criteria. This paper presents certain salient features of two independent finite element (numerical) studies of relevance to fatigue crack growth, where elasto-plastic analysis becomes significant. These problems can only be handled in the current day computational environment, and would have been only a dream just a few years ago.

Vibration Signature Analysis Using Artificial Neural Networks

Damage detection by measuring and analyzing vibration signals in a machine component is an established procedure in mechanical and aerospace engineering. This paper presents vibration signature analysis of steel bridge structures in a nonconventional way using artificial neural networks (ANN). Multilayer perceptrons have been adopted using the back-propagation algorithm for network training. The training patterns in terms of vibration signature are generated analytically for a moving load traveling on a trussed bridge structure at a constant speed to simulate the inspection vehicle. Using the finite-element technique, the moving forces are converted into stationary time-dependent force functions in order to generate vibration signals in the structure and the same is used to train the network. The performance of the trained networks is examined for their capability to detect damage from unknown signatures taken independently at one, three, and five nodes. It has been observed that the prediction using the trained network with single-node signature measurement at a suitability chosen location is even better than that of three-node and five-node measurement data.

Stochastic Aeroelastic Analysis of Composite Helicopter Rotor

The effect of uncertainty in composite material properties on the aeroelastic response, vibratory loads, and stability of a hingeless helicopter rotor is investigated. The uncertainty impact on rotating natural frequencies of the blade is studied with Monte Carlo simulations and first-order reliability methods. The stochastic aeroelastic analyses in hover and forward flight are carried out with Monte Carlo simulations. The flap, lag, and torsion responses show considerable scatter from their baseline values, and the uncertainty impact varies with the azimuth angle. Furthermore, the blade response shows finite probability of resonance-type conditions caused by modal frequencies approaching multiples of the rotor speed. The 4/rev vibratory forces show large deviations from their baseline values. The lag mode damping shows considerable scatter due to uncertain material properties with an almost 40% probability of instability in hover.

Nonlinear analysis of adhesively bonded lap joints considering viscoplasticity in adhesives

This paper presents nonlinear finite element analysis of adhesively bonded joints considering the elastoviscoplastic constitutive model of the adhesive material and the finite rotation of the joint. Though the adherends have been assumed to be linearly elastic, the yielding of the adhesive is represented by a pressure sensitive modified von Mises yield function. The stress-strain relation of the adhesive is represented by the Ramberg-Osgood relation. Geometric nonlinearity due to finite rotation in the joint is accounted for using the Green-Lagrange strain tensor and the second Piola-Kirchhoff stress tensor in a total Lagrangian formulation. Critical time steps have been calculated based on the eigenvalues of the transition matrices of the viscoplastic model of the adhesive. Stability of the viscoplastic solution and time dependent behaviour of the joints are examined. A parametric study has been carried out with particular reference to peel and shear stress along the interface. Critical zones for failure of joints have been identified. The study is of significance in the design of lap joints as well as on the characterization of adhesive strength. (C) 1999 Elsevier Science Ltd. All rights reserved.

Performance analysis of a modified spatial smoothing technique for direction estimation

The statistical performance analysis of ESPRIT, root-MUSIC, minimum-norm methods for direction estimation, due to finite data perturbations, using the modified spatially smoothed covariance matrix, is developed. Expressions for the mean-squared error in the direction estimates are derived based on a common framework. Based on the analysis, the use of the modified smoothed covariance matrix improves the performance of the methods when the sources are fully correlated. Also, the performance is better even when the number of subarrays is large unlike in the case of the conventionally smoothed covariance matrix. However, the performance for uncorrelated sources deteriorates due to an artificial correlation introduced by the modified smoothing. The theoretical expressions are validated using extensive simulations. (C) 1999 Elsevier Science B.V. All rights reserved.