An improved iterative finite element solution for pin joints

Accurate, reliable and economical methods of determining stress distributions are important for fastener joints. In the past the contact stress problems in these mechanically fastened joints using interference or push or clearance fit pins were solved using both inverse and iterative techniques. Inverse techniques were found to be most efficient, but at times inadequate in the presence of asymmetries. Iterative techniques based on the finite element method of analysis have wider applications, but they have the major drawbacks of being expensive and time-consuming. In this paper an improved finite element technique for iteration is presented to overcome these drawbacks. The improved iterative technique employs a frontal solver for elimination of variables not requiring iteration, by creation of a dummy element. This automatically results in a large reduction in computer time and in the size of the problem to be handled during iteration. Numerical results are compared with those available in the literature. The method is used to study an eccentrically located pin in a quasi-isotropic laminated plate under uniform tension.

Analysis of reflection characteristics of a normal incidence plane wave on resonant sound absorbers: A finite element approach

Resonant sound absorbers are used widely as anechoic coatings in underwater applications. In this paper a finite element scheme based on the Galerkin technique is used to analyze the reflection characteristics of the resonant absorber when insonified by a normal incidence plane wave. A waveguide theory coupled with an impedance matching condition in the fluid is used to model the problem. It is shown in this paper that the fluid medium encompassing the absorber can be modeled as an elastic medium with equivalent Lamé constants. Quarter symmetry conditions within the periodic unit cell are exploited. The finite element results are compared with analytical results, and with results published elsewhere in the literature. It is shown in the process that meshing of the fluid domain can be obviated if the transmission coefficients or reflection coefficients only are desired as is often the case. Finally, some design curves for thin resonant absorbers with water closure are presented in this paper.

Analysis of the evanescent field coupling between a fiber and a planar waveguide using the explicit finite difference method with a nonuniform mesh (Journal Paper)

We investigate an optical waveguide system consisting of an unclad fiber core suspended at a constant distance parallel to the surface of a planar waveguide. The coupling and propagation of light in the combined system is studied using the three-dimensional explicit finite difference beam propagation method with a nonuniform mesh configuration. The power loss in the fiber and the field distribution in the waveguide are studied as a function of various parameters, such as index changes, index profile, and propagation distance, for the combined system.

Dynamic response of a beam on elastic foundation of finite depth under a moving force

In this paper, dynamic response of an infinitely long beam resting on a foundation of finite depth, under a moving force is studied. The effect of foundation inertia is included in the analysis by modelling the foundation as a series of closely spaced axially vibrating rods of finite depth, fixed at the bottom and connected to the beam at the top. Viscous damping in the beam and foundation is included in the analysis. Steady state response of the beam-foundation system is obtained. Detailed numerical results are presented to study the effect of various parameters such as foundation mass, velocity of the moving load, damping and axial force on the beam. It is shown that foundation inertia can considerably reduce the critical velocity and can also amplify the beam response.

Finite element estimation of elastic interlaminar stresses in laminates

Low interlaminar strength and the consequent possibility of interlaminar failures in composite laminates demand an examination of interlaminar stresses and/or strains to ensure their satisfactory performance. As a first approximation, these stresses can be obtained from thickness-wise integration of ply equilibrium equations using in-plane stresses from the classical laminated plate theory. Implementation of this approach in the finite element form requires evaluation of third and fourth order derivatives of the displacement functions in an element. Hence, a high precision element developed by Jayachandrabose and Kirkhope (1985) is used here and the required derivatives are obtained in two ways. (i) from direct differentiation of element shape functions; and (ii) by adapting a finite difference technique applied to the nodal strains and curvatures obtained from the finite element analysis. Numerical results obtained for a three-layered symmetric and a two-layered asymmetric laminate show that the second scheme is quite effective compared to the first scheme particularly for the case of asymmetric laminates.

Flap‐Lag Damping in Hover and Forward Flight with a Three‐Dimensional Wake

Prediction of lag damping is difficult owing to the delicate balance of drag, induced drag and Coriolis forces in the in‐plane direction. Moreover, induced drag” is sensitive to dynamic wake, bath shed and trailing components, and thus its prediction requires adequate unsteady‐wake representation. Accordingly, rigid‐blade flap‐lag equations are coupled with a three‐dimensional finite‐state wake model; three isolatcd rotor canfigurations with three, four and five blades are treated over a range of thrust levels, tack numbers, lag frequencies and advance ratios. The investigation includes convergence characteristics of damping with respect to the number of radial shape functions and harmonics of the wake model for multiblade modes of low frequency (< 1/ rev.) to high frequency (> 1/rev.). Predicted flap and lag damping levels are then compared with similar predictions with 1) rigid wake (no unsteady induced now), 2) Loewy lift deficiency and 3) dynamic inflow. The coverage also includes correlations with the measured lag regressive‐mode damping in hover and forward flight and comparisons with similar correlations with dynamic inflow. Lag‐damping predictions with the dynamic wake model are consistently higher than the predictions with the dynamic inflow model; even for the low frequency lag regressive mode, the number of wake harmonics should at least be equal to twice the number of blades.

Nonlinear finite element analysis of reinforced concrete corbel

Reinforced concrete corbels have been analysed using the nonlinear finite element method. An elasto-plastic-cracking constitutive formulation using Huber-Hencky-Mises yield surface augmented with a tension cut-off is employed. Smeared-fixed cracking with mesh-dependent strain softening is employed to obtain objective results. Multiple non-orthogonal cracking and opening and closing of cracks are permitted. The model and the formulation are verified with respect to available numerical solution for an RC corbel. Results of analyses of nine reinforced concrete corbels are presented and compared with experimental results. Nonlinear finite element analysis of reinforced concrete structures is shown to be a complement and also a feasible alternative to laboratory testing.

Thermomechanical generalization of pin joints in finite plates

The understanding of thermoelastic behaviour of joints is significant in order to ensure the integrity of large and complex structures exposed to a thermal environment, particularly in fields such as aerospace and nuclear engineering. Thermomechanical generalization of partial contact behaviour of a pin joint under combined in-plane mechanical loading and on-axis unidirectional heat flow has already been established by the authors for the analytically simpler domains of large plates. This paper successfully extends the on-going investigation to a single pin in a finite rectangular isotropic plate as a two-dimensional abstraction from a practical situation of a multipin fastener joint. The finite element method is used to analyse the joint problem under on-axis thermomechanical loading and unified load-contact relationships are established for a class of loading conditions.

Time-domain-finite-wave analysis of the engine exhaust system by means of the stationary-frame method of characteristics. Part I. Theory

Time-domain-finite-wave analysis of the engine exhaust system is usually done using the method of characteristics. This makes use of either the moving frame method, or the stationary frame method. The stationary frame method is more convenient than its counterpart inasmuch as it avoids the tedium of graphical computations. In this paper (part I), the stationary-frame computational scheme along with the boundary conditions has been implemented. The analysis of a uniform tube, cavity-pipe junction including the engine and the radiation ends, and also the simple area discontinuities has been presented. The analysis has been done accounting for wall friction and heat-transfer for a one-dimensional unsteady flow. In the process, a few inconsistencies in the formulations reported in the literature have been pointed out and corrected. In the accompanying paper (part II) results obtained from the simulation are shown to be in good agreement with the experimental observations.

Time-domain-finite-wave analysis of the engine exhaust system by means of the stationary-frame method of characteristics. Part II. Computed results and experimental corroboration thereof

Time-domain-finite-wave analysis of engine exhaust systems is usually carried out by means of the method of characteristics. The theory and the computational details of the stationary-frame method have been worked out in the accompanying paper (part I). In this paper (part II), typical computed results are given and discussed. A setup designed for experimental corroboration is described. The results obtained from the simulation are found to be in good agreement with experimental observations.

Bearing capacity factors of circular foundations for a general c-phi soil using lower bound finite elements limit analysis

By using the lower bound limit analysis in conjunction with finite elements and linear programming, the bearing capacity factors due to cohesion, surcharge and unit weight, respectively, have been computed for a circular footing with different values of phi. The recent axisymmetric formulation proposed by the authors under phi = 0 condition, which is based on the concept that the magnitude of the hoop stress (sigma(theta)) remains closer to the least compressive normal stress (sigma(3)), is extended for a general c-phi soil. The computational results are found to compare quite well with the available numerical results from literature. It is expected that the study will be useful for solving various axisymmetric geotechnical stability problems. Copyright (C) 2010 John Wiley & Sons, Ltd.

On Two-User Gaussian Multiple Access Channels With Finite Input Constellations

Constellation Constrained (CC) capacity regions of two-user Single-Input Single-Output (SISO) Gaussian Multiple Access Channels (GMAC) are computed for several Non-Orthogonal Multiple Access schemes (NO-MA) and Orthogonal Multiple Access schemes (O-MA). For NO-MA schemes, a metric is proposed to compute the angle(s) of rotation between the input constellations such that the CC capacity regions are maximally enlarged. Further, code pairs based on Trellis Coded Modulation (TCM) are designed with PSK constellation pairs and PAM constellation pairs such that any rate pair within the CC capacity region can be approached. Such a NO-MA scheme which employs CC capacity approaching trellis codes is referred to as Trellis Coded Multiple Access (TCMA). Then, CC capacity regions of O-MA schemes such as Frequency Division Multiple Access (FDMA) and Time Division Multiple Access (TDMA) are also computed and it is shown that, unlike the Gaussian distributed continuous constellations case, the CC capacity regions with FDMA are strictly contained inside the CC capacity regions with TCMA. Hence, for finite constellations, a NO-MA scheme such as TCMA is better than FDMA and TDMA which makes NO-MA schemes worth pursuing in practice for two-user GMAC. Then, the idea of introducing rotations between the input constellations is used to construct Space-Time Block Code (STBC) pairs for two-user Multiple-Input Single-Output (MISO) fading MAC. The proposed STBCs are shown to have reduced Maximum Likelihood (ML) decoding complexity and information-losslessness property. Finally, STBC pairs with reduced sphere decoding complexity are proposed for two-user Multiple-Input Multiple-Output (MIMO) fading MAC.

Decoupled three-dimensional finite element computation of thermoelastic damping using Zener's approximation

We consider three dimensional finite element computations of thermoelastic damping ratios of arbitrary bodies using Zener's approach. In our small-damping formulation, unlike existing fully coupled formulations, the calculation is split into three smaller parts. Of these, the first sub-calculation involves routine undamped modal analysis using ANSYS. The second sub-calculation takes the mode shape, and solves on the same mesh a periodic heat conduction problem. Finally, the damping coefficient is a volume integral, evaluated elementwise. In the only other decoupled three dimensional computation of thermoelastic damping reported in the literature, the heat conduction problem is solved much less efficiently, using a modal expansion. We provide numerical examples using some beam-like geometries, for which Zener's and similar formulas are valid. Among these we examine tapered beams, including the limiting case of a sharp tip. The latter's higher-mode damping ratios dramatically exceed those of a comparable uniform beam.