Acoustic radiation and scattering from cylindrical bodies and analysis of transducer arrays

New mathematical methods to analytically investigate linear acoustic radiation and scattering from cylindrical bodies and transducer arrays are presented. Three problems of interest involving cylinders in an infinite fluid are studied. In all the three problems, the Helmholtz equation is used to model propagation through the fluid and the beam patterns of arrays of transducers are studied. In the first problem, a method is presented to determine the omni-directional and directional far-field pressures radiated by a cylindrical transducer array in an infinite rigid cylindrical baffle. The solution to the Helmholtz equation and the displacement continuity condition at the interface between the array and the surrounding water are used to determine the pressure. The displacement of the surface of each transducer is in the direction of the normal to the array and is assumed to be uniform. Expressions are derived for the pressure radiated by a sector of the array vibrating in-phase, the entire array vibrating in-phase, and a sector of the array phase-shaded to simulate radiation from a rectangular piston. It is shown that the uniform displacement required for generating a source level of 220 dB ref. μPa @ 1m that is omni directional in the azimuthal plane is in the order of 1 micron for typical arrays. Numerical results are presented to show that there is only a small difference between the on-axis pressures radiated by phased cylindrical arrays and planar arrays. The problem is of interest because cylindrical arrays of projectors are often used to search for underwater objects. In the second problem, the errors, when using data-independent, classical, energy and split beam correlation methods, in finding the direction of arrival (DOA) of a plane acoustic wave, caused by the presence of a solid circular elastic cylindrical stiffener near a linear array of hydrophones, are investigated. Scattering from the effectively infinite cylinder is modeled using the exact axisymmetric equations of motion and the total pressures at the hydrophone locations are computed. The effect of the radius of the cylinder, a, the distance between the cylinder and the array, b, the number of hydrophones in the array, 2H, and the angle of incidence of the wave, α, on the error in finding the DOA are illustrated using numerical results. For an array that is about 30 times the wavelength and for small angles of incidence (α<10), the error in finding the DOA using the energy method is less than that using the split beam correlation method with beam steered to α; and in some cases, the error increases when b increases; and the errors in finding the DOA using the energy method and the split beam correlation method with beam steered to α vary approximately as a7 / 4 . The problem is of interest because elastic stiffeners – in nearly acoustically transparent sonar domes that are used to protect arrays of transducers – scatter waves that are incident on it and cause an error in the estimated direction of arrival of the wave. In the third problem, a high-frequency ray-acoustics method is presented and used to determine the interior pressure field when a plane wave is normally incident on a fluid cylinder embedded in another infinite fluid. The pressure field is determined by using geometrical and physical acoustics. The interior pressure is expressed as the sum of the pressures due to all rays that pass through a point. Numerical results are presented for ka = 20 to 100 where k is the acoustic wavenumber of the exterior fluid and a is the radius of the cylinder. The results are in good agreement with those obtained using field theory. The directional responses, to the plane wave, of sectors of a circular array of uniformly distributed hydrophones in the embedded cylinder are then computed. The sectors are used to simulate linear arrays with uniformly distributed normals by using delays. The directional responses are compared with the output from an array in an infinite homogenous fluid. These outputs are of interest as they are used to determine the direction of arrival of the plane wave. Numerical results are presented for a circular array with 32 hydrophones and 12 hydrophones in each sector. The problem is of interest because arrays of hydrophones are housed inside sonar domes and acoustic plane waves from distant sources are scattered by the dome filled with fresh water and cause deterioration in the performance of the array.

An exact series solution for the vibration analysis of cylindrical shells with arbitrary boundary conditions

In this paper, an exact series solution for the vibration analysis of circular cylindrical shells with arbitrary boundary conditions is obtained, using the elastic equations based on Flügge's theory. Each of the three displacements is represented by a Fourier series and auxiliary functions and sought in a strong form by letting the solution exactly satisfy both the governing differential equations and the boundary conditions on a point-wise basis. Since the series solution has to be truncated for numerical implementation, the term exactly satisfying should be understood as a satisfaction with arbitrary precision. One of the important advantages of this approach is that it can be universally applied to shells with a variety of different boundary conditions, without the need of making any corresponding modifications to the solution algorithms and implementation procedures as typically required in other techniques. Furthermore, the current method can be easily used to deal with more complicated boundary conditions such as point supports, partial supports, and non-uniform elastic restraints. Numerical examples are presented regarding the modal parameters of shells with various boundary conditions. The capacity and reliability of this solution method are demonstrated through these examples. © 2012 Elsevier Ltd. All rights reserved.

One and two degrees-of-freedom Vortex-Induced Vibration experiments with yawed cylinders

Vortex-Induced Vibration (VIV) experiments were carried out with yawed cylinders. The purpose was to investigate the validity of the Independence Principle (IP) for properly describing the flow characteristics and the dynamics of structures subjected to oblique flow. Five yaw angles in relation to the direction perpendicular to the free stream velocity were tested, namely View the MathML sourceθ=0°,10°,20°,30° and 45°. Both the upstream and downstream orientations were tested. The models were mounted on a leaf spring apparatus that allows experiments with one or two degrees of freedom. The Reynolds numbers based on the component normal to the cylinder axis fell in the interval 3×103

On the flexural vibration of cylinders under axial loads:Numerical and experimental study

The flexural vibration of a homogeneous isotropic linearly elastic cylinder of any aspect ratio is analysed in this paper. Natural frequencies of a cylinder under uniformly distributed axial loads acting on its bases are calculated numerically by the Ritz method with terms of power series in the coordinate directions as approximating functions. The effect of axial loads on the flexural vibration cannot be described by applying infinitesimal strain theory, therefore, geometrically nonlinear strain–displacement relations with second-order terms are considered here. The natural frequencies of free–free, clamped–clamped, and sliding–sliding cylinders subjected to axial loads are calculated using the proposed three-dimensional Ritz approach and are compared with those obtained with the finite element method and the Bernoulli–Euler theory. Different experiments with cylinders axially compressed by a hydraulic press are carried out and the experimental results for the lowest flexural frequency are compared with the numerical results. An approach based on the Ritz formulation is proposed for the flexural vibration of a cylinder between the platens of the press with constraints varying with the intensity of the compression. The results show that for low compressions the cylinder behaves similarly to a sliding–sliding cylinder, whereas for high compressions the cylinder vibrates as a clamped–clamped one.

A canonical UTD solution for electromagnetic scattering by an electrically large impedance circular cylinder illuminated by an obliquely incident plane wave

A uniform geometrical theory of diffraction (UTD) solution is developed for the canonical problem of the electromagnetic (EM) scattering by an electrically large circular cylinder with a uniform impedance boundary condition (IBC), when it is illuminated by an obliquely incident high frequency plane wave. A solution to this canonical problem is first constructed in terms of an exact formulation involving a radially propagating eigenfunction expansion. The latter is converted into a circumferentially propagating eigenfunction expansion suited for large cylinders, via the Watson transform, which is expressed as an integral that is subsequently evaluated asymptotically, for high frequencies, in a uniform manner. The resulting solution is then expressed in the desired UTD ray form. This solution is uniform in the sense that it has the important property that it remains continuous across the transition region on either side of the surface shadow boundary. Outside the shadow boundary transition region it recovers the purely ray optical incident and reflected ray fields on the deep lit side of the shadow boundary and to the modal surface diffracted ray fields on the deep shadow side. The scattered field is seen to have a cross-polarized component due to the coupling between the TEz and TMz waves (where z is the cylinder axis) resulting from the IBC. Such cross-polarization vanishes for normal incidence on the cylinder, and also in the deep lit region for oblique incidence where it properly reduces to the geometrical optics (GO) or ray optical solution. This UTD solution is shown to be very accurate by a numerical comparison with an exact reference solution.

Discontinuous solutions and boundary value problems for non-linearly elastic fibre-reinforced materials

En este trabajo se han analizado varios problemas en el contexto de la elasticidad no lineal basándose en modelos constitutivos representativos. En particular, se han analizado problemas relacionados con el fenómeno de perdida de estabilidad asociada con condiciones de contorno en el caso de material reforzados con fibras. Cada problema se ha formulado y se ha analizado por separado en diferentes capítulos. En primer lugar se ha mostrado el análisis del gradiente de deformación discontinuo para un material transversalmente isótropo, en particular, el modelo del material considerado consiste de una base neo-Hookeana isótropa incrustada con fibras de refuerzo direccional caracterizadas con un solo parámetro. La solución de este problema se vincula con instabilidades que dan lugar al mecanismo de fallo conocido como banda de cortante. La perdida de elipticidad de las ecuaciones diferenciales de equilibrio es una condición necesaria para que aparezca este tipo de soluciones y por tanto las inestabilidades asociadas. En segundo lugar se ha analizado una deformación combinada de extensión, inación y torsión de un tubo cilíndrico grueso donde se ha encontrado que la deformación citada anteriormente puede ser controlada solo para determinadas direcciones de las fibras refuerzo. Para entender el comportamiento elástico del tubo considerado se ha ilustrado numéricamente los resultados obtenidos para las direcciones admisibles de las fibras de refuerzo bajo la deformación considerada. En tercer lugar se ha estudiado el caso de un tubo cilíndrico grueso reforzado con dos familias de fibras sometido a cortante en la dirección azimutal para un modelo de refuerzo especial. En este problema se ha encontrado que las inestabilidades que aparecen en el material considerado están asociadas con lo que se llama soluciones múltiples de la ecuación diferencial de equilibrio. Se ha encontrado que el fenómeno de instabilidad ocurre en un estado de deformación previo al estado de deformación donde se pierde la elipticidad de la ecuación diferencial de equilibrio. También se ha demostrado que la condición de perdida de elipticidad y ^W=2 = 0 (la segunda derivada de la función de energía con respecto a la deformación) son dos condiciones necesarias para la existencia de soluciones múltiples. Finalmente, se ha analizado detalladamente en el contexto de elipticidad un problema de un tubo cilíndrico grueso sometido a una deformación combinada en las direcciones helicoidal, axial y radial para distintas geotermias de las fibras de refuerzo . In the present work four main problems have been addressed within the framework of non-linear elasticity based on representative constitutive models. Namely, problems related to the loss of stability phenomena associated with boundary value problems for fibre-reinforced materials. Each of the considered problems is formulated and analysed separately in different chapters. We first start with the analysis of discontinuous deformation gradients for a transversely isotropic material under plane deformation. In particular, the material model is an augmented neo-Hookean base with a simple unidirectional reinforcement characterised by a single parameter. The solution of this problem is related to material instabilities and it is associated with a shear band-type failure mode. The loss of ellipticity of the governing differential equations is a necessary condition for the existence of these material instabilities. The second problem involves a detailed analysis of the combined non-linear extension, inflation and torsion of a thick-walled circular cylindrical tube where it has been found that the aforementioned deformation is controllable only for certain preferred directions of transverse isotropy. Numerical results have been illustrated to understand the elastic behaviour of the tube for the admissible preferred directions under the considered deformation. The third problem deals with the analysis of a doubly fibre-reinforced thickwalled circular cylindrical tube undergoing pure azimuthal shear for a special class of the reinforcing model where multiple non-smooth solutions emerge. The associated instability phenomena are found to occur prior to the point where the nominal stress tensor changes monotonicity in a particular direction. It has been also shown that the loss of ellipticity condition that arises from the equilibrium equation and ^W=2 = 0 (the second derivative of the strain-energy function with respect to the deformation) are equivalent necessary conditions for the emergence of multiple solutions for the considered material. Finally, a detailed analysis in the basis of the loss of ellipticity of the governing differential equations for a combined helical, axial and radial elastic deformations of a fibre-reinforced circular cylindrical tube is carried out.

Thermal buckling of circular cylindrical and conical thin-walled shells

"Aeronautical Research Laboratory. Contract AF 33(616)-3885, Project no. 7(8-1367), Task no. 70524."

On the buckling of circular cylindrical shells under external pressure,

"Project no. 9782. Contract AF 49(638)-430... CU-23-62-AF-430-CE."

Modeling of Cylindrical Couette Flow of Rarefied Gas Between Rotating Cylinders

Петър Господинов, Добри Данков, Владимир Русинов, Стефан Стефанов - Изследвано е цилиндрично течение на Кует за разреден газ между два въртящи се цилиндъра. Получени са профилите на налягането, скоростта и температурата по метода на прякото статистическо моделиране (DSMC) и чрез числено решаване на уравненията на Навие-Стокс за свиваем флуид. Резултатите сочат много добро съвпадение за малки числа на Кнудсен Kn = 0.02. Показано е, че при различни кинематични гранични условия, газът изостава или избързва спрямо скоростта на стената, или има поведение на твърдо еластично тяло. Получените резултати са важни при решаването на неравнинни, задачи от микрофлуидиката с отчитане на ефектите на кривината.

Impact response and parametric studies of reinforced concrete circular columns

This paper presents a detailed description of the influence of critical parameters that govern the vulnerability of columns under lateral impact loads. Numerical simulations are conducted by using the Finite Element program LS-DYNA, incorporating steel reinforcement, material models and strain rate effects. A simplified method based on impact pulse generated from full scale impact tests is used for impact reconstruction and effects of the various pulse loading parameters are investigated under low to medium velocity impacts. A constitutive material model which can simulate failures under tri-axial state of stresses is used for concrete. Confinement effects are also introduced to the numerical simulation and columns of Grade 30 to 50 concrete under pure axial loading are analysed in detail. This research confirmed that the vulnerability of the axially loaded columns can be mitigated by reducing the slenderness ratio and concrete grade, and by choosing the design option with a minimal amount of longitudinal steel. Additionally, it is evident that approximately a 50% increase in impact capacity can be gained for columns in medium rise buildings by enhancing the confinement effects alone. Results also indicated that the ductility as well as the mode of failure under impact can be changed with the volumetric ratio of lateral steel. Moreover, to increase the impact capacity of the vulnerable columns, a higher confining stress is required. The general provisions of current design codes do not sufficiently cover this aspect and hence this research will provide additional guidelines to overcome the inadequacies of code provisions.

Determination of elastic modulus of thin coating materials using nanoindentation and finite element analysis

A deconvolution method that combines nanoindentation and finite element analysis was developed to determine elastic modulus of thin coating layer in a coating-substrate bilayer system. In this method, the nanoindentation experiments were conducted to obtain the modulus of both the bilayer system and the substrate. The finite element analysis was then applied to deconvolve the elastic modulus of the coating. The results demonstrated that the elastic modulus obtained using the developed method was in good agreement with that reported in literature.

Elastic lateral buckling of simply supported LiteSteel beams subject to transverse loading

The buckling strength of a new cold-formed hollow flange channel section known as LiteSteel beam (LSB) is governed by lateral distortional buckling characterised by simultaneous lateral deflection, twist and web distortion for its intermediate spans. Recent research has developed a modified elastic lateral buckling moment equation to allow for lateral distortional buckling effects. However, it is limited to a uniform moment distribution condition that rarely exists in practice. Transverse loading introduces a non-uniform bending moment distribution, which is also often applied above or below the shear centre (load height). These loading conditions are known to have significant effects on the lateral buckling strength of beams. Many steel design codes have adopted equivalent uniform moment distribution and load height factors to allow for these effects. But they were derived mostly based on data for conventional hot-rolled, doubly symmetric I-beams subject to lateral torsional buckling. The moment distribution and load height effects of transverse loading for LSBs, and the suitability of the current design modification factors to accommodate these effects for LSBs is not known. This paper presents the details of a research study based on finite element analyses on the elastic lateral buckling strength of simply supported LSBs subject to transverse loading. It discusses the suitability of the current steel design code modification factors, and provides suitable recommendations for simply supported LSBs subject to transverse loading.

Control issues for velocity sourced Series Elastic Actuators

The Series Elasic Actuator has been proposed as a method for providing safe force or torque based acutation for robots that interact with humans. In this paper we look at some outstanding issues in the implementation and control of Series Elastic Actuators. The study addresses issues in making the Series Elastic Actuator respond effectively in the presence of physical difficulties such as restriction, using a computation efficient controller. The improvement over previous implementations is achieved by treating the motor as a velocity source to the elastic element, rather than as a torque source.

Assessing the safety of a velocity sourced series elastic actuator using the head injury criterion

The Velocity Sourced Series Elastic Actuator has been proposed as a method for providing safe force or torque based actuation for robots without compromising the actuator performance. In this paper we assess the safety of Velocity Sourced Series Elastic Actuators by measuring the Head Injury Criterion scores for collisions with a model head. The study makes a comparative analysis against stiff, high impedance actuation using the same motor without the series elastic component, showing that the series elastic component brings about a massive reduction in the chance of head injury. The benefits of a collision detection and safe reaction system are shown to be limited to collisions at low speeds, providing greater interaction comfort but not necessarily contributing to safety from injury.