Axisymmetric Thermoelastic Interactions without Energy Dissipation in an Unbounded Body with Cylindrical Cavity

The linear theory of thermoelasticity without energy dissipation is employed to study thermoelastic interactions in a homogeneous and isotropic unbounded body containing a cylindrical cavity. The interactions are supposed to be due to a constant step in radial stress or temperature applied to the boundary of the cavity, which is maintained at a constant temperature or zero radial stress (as the case may be). By using the Laplace transform technique, it is found that the interactions consist of two coupled waves both of which propagate with a finite speed but with no attenuation. The discontinuities that occur at the wavefronts are computed. Numerical results applicable to a copper-like material are presented

Fully coupled peridynamic thermomechanics

This study concerns the derivation of the coupled peridynamic (PD) thermomechanics equations based on thermodynamic considerations. The generalized peridynamic model for fully coupled thermomechanics is derived using the conservation of energy and the free-energy function. Subsequently, the bond-based coupled PD thermomechanics equations are obtained by reducing the generalized formulation. These equations are also cast into their nondimensional forms. After describing the numerical solution scheme, solutions to certain coupled thermomechanical problems with known previous solutions are presented

Solution of Coupled Thermoelasticity Problem in Rotating Disks

The main purpose of this dissertation is to study coupled thermoelastic behaviors in disks subjected to thermal shock loads based on the generalized and classic theories of coupled thermoelasticity. To this end, this research has been carried out in two stages. In the first stage, thermoelasticity problems in an axisymmetric rotating disk with constant thickness made of a homogeneous isotropic material are analytically solved and closed-form formulations are presented for temperature and displacement fields. Since, the analytical solution is not always feasible, the finite element (FE) method can be employed for more sophisticated coupled thermoelasticity problems. Accordingly, in the second stage of the research, a novel refined 1D finite element approach with 3D-like accuracies are developed for theories of coupled thermoelasticity. Then, the developed FE models are applied for a 3D solution of the dynamic generalized coupled thermoelasticity problem in disks. Use of the reduced models with low computational costs may be of interest in a laborious time history analysis of the dynamic problems. The obtained analytical and numerical solutions are in good agreement with the results available in the literature. It is further shown that the proposed analytical and FE methods are quite efficient with very high rate of convergence

Laser thermoelastic generation in metals above the melt threshold

An approach is presented for calculating thermoelastic generation of ultrasound in a metal plate exposed to nanosecond pulsed laser heating, sufficient to cause melting but not ablation. Detailed consideration is given to the spatial and temporal profiles of the laser pulse, penetration of the laser beam into the sample, the appearance and subsequent growth and then contraction of the melt pool, and the time dependent thermal conduction in the melt and surrounding solid throughout. The excitation of the ultrasound takes place during and shortly after the laser pulse and occurs predominantly within the thermal diffusion length of a micron or so beneath the surface. It is shown how, because of this, the output of the thermal simulations can be expressed as axially symmetric transient radial and normal surface force distributions. The epicentral displacement response to these force distributions is obtained by two methods, the one based on the elastodynamic Green’s functions for plate geometry determined by the Cagniard generalized ray method and the other using a finite element numerical method. The two approaches are in very close agreement. Numerical simulations are reported on the epicentral displacement response of a 3.12mm thick tungsten plate irradiated with a 4 ns pulsed laser beam with Gaussian spatial profile, at intensities below and above the melt threshol

Development of an integrated BEM approach for hot fluid structure interaction

A comprehensive boundary element method is presented for transient thermoelastic analysis of hot section Earth-to-Orbit engine components. This time-domain formulation requires discretization of only the surface of the component, and thus provides an attractive alternative to finite element analysis for this class of problems. In addition, steep thermal gradients, which often occur near the surface, can be captured more readily since with a boundary element approach there are no shape functions to constrain the solution in the direction normal to the surface. For example, the circular disc analysis indicates the high level of accuracy that can be obtained. In fact, on the basis of reduced modeling effort and improved accuracy, it appears that the present boundary element method should be the preferred approach for general problems of transient thermoelasticity

Fully-coupled 3D Simulation of Electromagnetic Forming

Electromagnetic metal forming is a contact-free high-speed forming process in which strain rates of more than 103 s^(-1) are achieved. The deformation of the workpiece is driven by a material body force, the Lorentz force, that results from the interaction of a pulsed magnetic field with eddy currents induced in the workpiece by the magnetic field itself. The purpose of this work is to present a fully-coupled 3D simulation of the process. For the mechanical structure a thermoelastic, viscoplastic, electromagnetic material model is relevant, which is incorporated in a large-deformation dynamic formulation. The evolution of the electromagnetic fields is governed by Maxwell s equations under quasistatic conditions. Their numerical solution in 3D requires particular arrangements due to a reduced regularity at material interfaces. Hence, Nédélec elements are employed. Coupling between the thermomechanical and electromagnetic subsystems takes the form of the Lorentz force, the electromotive intensity, and the current geometry of the workpiece. A staggered scheme based on a Lagrangian mesh for the workpiece and an ALE formulation for the electromagnetic field is utilized to solve the coupled system, guaranteeing the efficiency and accuracy of the data transfer between the two meshes

Development of an integrated BEM approach for hot fluid structure interaction

The progress made toward the development of a boundary element formulation for the study of hot fluid-structure interaction in Earth-to-Orbit engine hot section components is reported. The convective viscous integral formulation was derived and implemented in the general purpose computer program GP-BEST. The new convective kernel functions, in turn, necessitated the development of refined integration techniques. As a result, however, since the physics of the problem is embedded in these kernels, boundary element solutions can now be obtained at very high Reynolds number. Flow around obstacles can be solved approximately with an efficient linearized boundary-only analysis or, more exactly, by including all of the nonlinearities present in the neighborhood of the obstacle. The other major accomplishment was the development of a comprehensive fluid-structure interaction capability within GP-BEST. This new facility is implemented in a completely general manner, so that quite arbitrary geometry, material properties and boundary conditions may be specified. Thus, a single analysis code (GP-BEST) can be used to run structures-only problems, fluids-only problems, or the combined fluid-structure problem. In all three cases, steady or transient conditions can be selected, with or without thermal effects. Nonlinear analyses can be solved via direct iteration or by employing a modified Newton-Raphson approach