A general procedure for modified crack closure integral in 3D problems with cracks

In linear elastic fracture mechanics (LEFM), Irwin's crack closure integral (CCI) is one of the signficant concepts for the estimation of strain energy release rates (SERR) G, in individual as well as mixed-mode configurations. For effective utilization of this concept in conjunction with the finite element method (FEM), Rybicki and Kanninen [Engng Fracture Mech. 9, 931 938 (1977)] have proposed simple and direct estimations of the CCI in terms of nodal forces and displacements in the elements forming the crack tip from a single finite element analysis instead of the conventional two configuration analyses. These modified CCI (MCCI) expressions are basically element dependent. A systematic derivation of these expressions using element stress and displacement distributions is required. In the present work, a general procedure is given for the derivation of MCCI expressions in 3D problems with cracks. Further, a concept of sub-area integration is proposed which facilitates evaluation of SERR at a large number of points along the crack front without refining the finite element mesh. Numerical data are presented for two standard problems, a thick centre-cracked tension specimen and a semi-elliptical surface crack in a thick slab. Estimates for the stress intensity factor based on MCCI expressions corresponding to eight-noded brick elements are obtained and compared with available results in the literature.

Unsteady free convection flow under the influence of a magnetic field

The unsteady free convection boundary layer at the stagnation point of a two-dimensional body and an axisymmetric body with prescribed surface heat flux or temperature has been studied. The magnetic field is applied parallel to the surface and the effect of induced magnetic field has been considered. It is found that for certain powerlaw distribution of surface heat flux or temperature and magnetic field with time, the governing boundary layer equations admit a self-similar solution locally. The resulting nonlinear ordinary differential equations have been solved using a finite element method and a shooting method with Newton's corrections for missing initial conditions. The results show that the skin friction and heat transfer coefficients, and x-component of the induced magnetic field on the surface increase with the applied magnetic field. In general, the skin friction, heat transfer and x-component of the induced magnetic field for axisymmetric case are more than those of the two-dimensional case. Also they change more when the surface heat flux or temperature decreases with time than when it increases with time. The skin friction, heat transfer and x-component of the induced magnetic field are significantly affected by the magnetic Prandtl number and they increase as the magnetic Prandtl number decreases. The skin friction and x-component of the magnetic field increase with the dissipation parameter, but heat transfer decreases.

A numerical study of division of flow in open channels

A two-dimensional numerical model which employs the depth-averaged forms of continuity and momentum equations along with k-e turbulence closure scheme is used to simulate the flow at the open channel divisions. The model is generalised to flows of arbitrary geometries and MacCormack finite volume method is used for solving governing equations. Application of cartesian version of the model to analyse the flow at right-angled junction is presented. The numerical predictions are compared with experimental data of earlier investigators and measurements made as part of the present study. Performance of the model in predicting discharge distribution, surface profiles, separation zone parameters and energy losses is evaluated and discussed in detail. To illustrate the application of the numerical model to analyse the flow in acute angled offtakes and streamlined branch entries, a few computational results are presented.

Plates on Elastic Foundation

Analysis of rectangular plates resting on a Winkler-type, one-parameter foundation is studied. The finite element method is applied and a 12-degree-of-freedom, nonconforming rectangular plate element is adopted. Based on shape functions of the plate element, an energy approach is used to derive a closed-form, 12-by-12, consistent foundation stiffness matrix for a rectangular plate on an elastic subgrade. A commonly used method of modeling structural elements on an elastic foundation is the application of discrete springs at the element nodes. The model developed in this paper is compared with the discrete spring model and the convergence of both models is discussed. The convergence of the models is compared with the well-known classical solution of plates on elastic foundation developed in the 1950s. Both models show good convergence to the classical solution. The continuous subgrade response model converges in a manner more consistent with the flexibility of the plate element.

Stability of laminated composite plates with cut-outs

Critical buckling loads of laminated fibre-reinforced plastic square panels have been obtained using the finite element method. Various boundary conditions, lay-up details, fibre orientations, cut-out sizes are considered. A 36 degrees of freedom triangular element, based on the classical lamination theory (CLT) has been used for the analysis. The performance of this element is validated by comparing results with some of those available in literature. New results have been given for several cases of boundary conditions for [0°/ ± 45°/90°]s laminates. The effect of fibre-orientation in the ply on the buckling loads has been investigated by considering [±?]6s laminates.

Stiffened composite cylindrical panels�Optimum lay-up for buckling by ranking

Buckling of discretely stiffened composite cylindrical panels made of repeated sublaminate construction is studied using a finite element method. In repeated sublaminate construction, a full laminate is obtained by repeating a basic sublaminate, which has a smaller number of plies. This paper deals with the determination of the optimum lay-up for buckling by ranking of such stiffened (longitudinal and hoop) composite cylindrical panels. For this purpose we use the particularized form of a four-noded, 48 degrees of freedom doubly curved quadrilateral thin shell finite element together with a fully compatible two-noded, 16 degrees of freedom composite stiffener element. The computer program developed has been used, after extensive checking for correctness, to obtain an optimum orientation scheme of the plies in the sublaminate so as to achieve maximum buckling load for a specified thickness of typical stiffened composite cylindrical panels.

Role of surface tension driven convection in pulsed - laser melting and solidification

Surface melting by a stationary, pulsed laser has been modelled by the finite element method. The role of the surface tension driven convection is investigated in detail. Numerical results are presented for a triangular laser pulse of durations 10, 50 and 200 ms. Though the magnitude of the velocity is high due to the surface tension forces, the present results indicate that a finite time is required for convection to affect the temperature distribution within the melt pool. The effect of convection is very significant for pulse durations longer than 10 ms.

Synthesis of path generating flexible-link mechanisms

Flexible-link mechanisms are those linkage mechanisms (or structures) which are capable of motion by virtue of elastic deformation of one or more;links. In such mechanisms a single flexible link; can replace several rigid links and joints resulting in fewer links, fewer pin joints, reduced overall weight and reduced mechanical error. In spite of such clear advantages, contributions toward flexible-link mechanisms remain very scarce. The area of flexible-link mechanisms offers much scope for further exploration. This paper attempts to show the potential of flexible-link mechanisms in accomplishing a kinematic task like path generation. Synthesis of a four-bar mechanism with a flexible rocker for circular and straight line path generation is carried out. Displacement analysis of the structure is carried out using finite element method (FEM) and synthesis is formulated and solved as an optimization problem. Several numerical examples are presented for illustration. Based on the results obtained with these examples, the flexible-link mechanism considered shows good promise for-path generation.

Microwave heating in multiphase systems: evaluation of series solutions

The 1D electric field and heat-conduction equations are solved for a slab where the dielectric properties vary spatially in the sample. Series solutions to the electric field are obtained for systems where the spatial variation in the dielectric properties can be expressed as polynomials. The series solution is used to obtain electric-field distributions for a binary oil-water system where the dielectric properties are assumed to vary linearly within the sample. Using the finite-element method temperature distributions are computed in a three-phase oil, water and rock system where the dielectric properties vary due to the changing oil saturation in the rock. Temperature distributions predicted using a linear variation in the dielectric properties are compared with those obtained using the exact nonlinear variation.

Modelling, simulation, and graphical user interface for industrial gas carburising process

A mathematical model has been developed for the gas carburising (diffusion) process using finite volume method. The computer simulation has been carried out for an industrial gas carburising process. The model's predictions are in good agreement with industrial experimental data and with data collected from the literature. A study of various mass transfer and diffusion coefficients has been carried out in order to suggest which correlations should be used for the gas carburising process. The model has been interfaced in a Windows environment using a graphical user interface. In this way, the model is extremely user friendly. The sensitivity analysis of various parameters such as initial carbon concentration in the specimen, carbon potential of the atmosphere, temperature of the process, etc. has been carried out using the model.

Rayleigh-Benard convection during solidification of an eutectic solution cooled from the top

The interaction between laminar Rayleigh-Benard convection and directional solidification is studied for the case of an eutectic solution kept in a rectangular cavity cooled from the top. Experiments and numerical simulations are carried out using an NH4Cl-H2O solution as the model fluid. The flow is visualized using a sheet of laser light scattered by neutrally buoyant, hollow-glass spheres seeded in the fluid. The numerical modeling is performed using a pressure-based finite-volume method according to the SIMPLER algorithm. The present configuration enables us to visualize flow vortices in the presence of a continuously evolving solid/liquid interface. Clear visualization of the Rayleigh-Benard convective cells and their interaction with the solidification front are obtained. It is observed that the convective cells are characterized by zones of up-flow and down-flow, resulting in the development of a nonplanar interface. Because of the continuous advancement of the solid/liquid interface, the effective liquid height of the cavity keeps decreasing. Once the height of the fluid layer falls below a critical value, the convective cells become weaker and eventually die out, leading to the growth of a planar solidification front. Results of flow visualization and temperature measurement are compared with those from the numerical simulation, and a good agreement is found.

Study of gas-fines in decreasing gas velocity in a packed bed

It has been well recognized now that the blast furnace aerodynamics can be represented more accurately under the decreasing gas velocity condition. Therefore, gas-fines study has been carried out in a packed bed under the decreasing gas velocity condition. Gas and fines flow equations have been developed and solved, for two-dimensional case using finite volume method. To take into account the turbulence, k-e turbulent flow model has also been developed in two-dimension. The model's predictions have been validated against the published experimental data for the increasing gas velocity case, as no experimental data are available in open literature for the decreasing gas velocity. This study shows the difference in the results for increasing and decreasing gas velocity cases under various conditions which have been reported here. Implication of the results to the blast furnace condition has also been discussed.

Effects of phase inhomogeneity and boundary conditions on the dynamic response of SMA wire actuators

This paper reports the simulation results from the dynamic analysis of a Shape Memory Alloy (SMA) actuator. The emphasis is on understanding the dynamic behavior under various loading rates and boundary conditions, resulting in complex scenarios such as thermal and stress gradients. Also, due to the polycrystalline nature of SMA wires, presence of microstructural inhomogeneity is inevitable. Probing the effect of inhomogeneity on the dynamic behavior can facilitate the prediction of life and characteristics of SMA wire actuator under varieties of boundary and loading conditions. To study the effect of these factors, an initial boundary value problem of SMA wire is formulated. This is subsequently solved using finite element method. The dynamic response of the SMA wire actuator is analyzed under mechanical loading and results are reported. Effect of loading rate, micro-structural inhomogeneity and thermal boundary conditions on the dynamic response of SMA wire actuator is investigated and the simulation results are reported.

Modified lattice model for mode-I fracture analysis of notched plain concrete beam using probabilistic approach

A modified lattice model using finite element method has been developed to study the mode-I fracture analysis of heterogeneous materials like concrete. In this model, the truss members always join at points where aggregates are located which are modeled as plane stress triangular elements. The truss members are given the properties of cement mortar matrix randomly, so as to represent the randomness of strength in concrete. It is widely accepted that the fracture of concrete structures should not be based on strength criterion alone, but should be coupled with energy criterion. Here, by incorporating the strain softening through a parameter ‘α’, the energy concept is introduced. The softening branch of load-displacement curves was successfully obtained. From the sensitivity study, it was observed that the maximum load of a beam is most sensitive to the tensile strength of mortar. It is seen that by varying the values of properties of mortar according to a normal random distribution, better results can be obtained for load-displacement diagram.

Implementation of a phase array diffuse optical tomographic imager

Diffuse optical tomography (DOT) using near-infrared (NIR) light is a promising tool for noninvasive imaging of deep tissue. This technique is capable of quantitative reconstructions of absorption coefficient inhomogeneities of tissue. The motivation for reconstructing the optical property variation is that it, and, in particular, the absorption coefficient variation, can be used to diagnose different metabolic and disease states of tissue. In DOT, like any other medical imaging modality, the aim is to produce a reconstruction with good spatial resolution and accuracy from noisy measurements. We study the performance of a phase array system for detection of optical inhomogeneities in tissue. The light transport through a tissue is diffusive in nature and can be modeled using diffusion equation if the optical parameters of the inhomogeneity are close to the optical properties of the background. The amplitude cancellation method that uses dual out-of-phase sources (phase array) can detect and locate small objects in turbid medium. The inverse problem is solved using model based iterative image reconstruction. Diffusion equation is solved using finite element method for providing the forward model for photon transport. The solution of the forward problem is used for computing the Jacobian and the simultaneous equation is solved using conjugate gradient search. The simulation studies have been carried out and the results show that a phase array system can resolve inhomogeneities with sizes of 5 mm when the absorption coefficient of the inhomogeneity is twice that of the background tissue. To validate this result, a prototype model for performing a dual-source system has been developed. Experiments are carried out by inserting an inhomogeneity of high optical absorption coefficient in an otherwise homogeneous phantom while keeping the scattering coefficient same. The high frequency (100 MHz) modulated dual out-of-phase laser source light is propagated through the phantom. The interference of these sources creates an amplitude null and a phase shift of 180° along a plane between the two sources with a homogeneous object. A solid resin phantom with inhomogeneities simulating the tumor is used in our experiment. The amplitude and phase changes are found to be disturbed by the presence of the inhomogeneity in the object. The experimental data (amplitude and the phase measured at the detector) are used for reconstruction. The results show that the method is able to detect multiple inhomogeneities with sizes of 4 mm. The localization error for a 5 mm inhomogeneity is found to be approximately 1 mm.