Micro/macro-crack growth due to creep–fatigue dependency on time–temperature material behavior

The implicit character of micro-structural degradation is determined by specifying the time history of crack growth caused by creep–fatigue interaction at high temperature. A dual scale micro/macro-equivalent crack growth model is used to illustrate the underlying principle of multiscaling which can be applied equally well to nano/micro. A series of dual scale models can be connected to formulate triple or quadruple scale models. Temperature and time-dependent thermo-mechanical material properties are developed to dictate the design time history of creep–fatigue cracking that can serve as the master curve for health monitoring.In contrast to the conventional procedure of problem/solution approach by specifying the time- and temperature-dependent material properties as a priori, the desired solution is then defined for a class of anticipated loadings. A scheme for matching the loading history with the damage evolution is then obtained. The results depend on the initial crack size and the extent of creep in proportion to fatigue damage. The path dependent nature of damage is demonstrated by showing the range of the pertinent parameters that control the final destruction of the material. A possible scenario of 20 yr of life span for the 38Cr2Mo2VA ultra-high strength steel is used to develop the evolution of the micro-structural degradation. Three micro/macro-parameters μ^*, d^* and σ^* are used to exhibit the time-dependent variation of the material, geometry and load effects. They are necessary to reflect the scale transitory behavior of creep–fatigue damage. Once the algorithm is developed, the material can be tailor made to match the behavior. That is a different life span of the same material would alter the time behavior of μ^*, d^* and σ^* and hence the micro-structural degradation history. The one-to-one correspondence of the material micro-structure degradation history with that of damage by cracking is the essence of path dependency. Numerical results and graphs are obtained to demonstrate how the inherently implicit material micro-structure parameters can be evaluated from the uniaxial bulk material properties at the macroscopic scale.The combined behavior of creep and fatigue can be exhibited by specifying the parameter ξ with reference to the initial defect size a₀. Large ξ (0.90 and 0.85) gives critical crack size a_cr = 11–14 mm (at t < 20 yr) for a₀ about 1.3 mm. For small ξ (0.05 and 0.15), there results critical a_cr = 6–7 mm (at t < 20 yr) for a₀ about 0.7–0.8 mm. The initial crack is estimated to increase its length by an order of magnitude before triggering global to the instability. This also applies ξ ≈ 0.5 where creep interacts severely with fatigue. Fine tuning of a_cr and a₀ can be made to meet the condition oft = 20 yr.Trade off among load, material and geometric parameters are quantified such that the optimum conditions can be determined for the desired life qualified by the initial–final defect sizes. The scenario assumed in this work is indicative of the capability of the methodology. The initial–final defect sizes can be varied by re-designing the time–temperature material specifications. To reiterate, the uniqueness of solution requires the end result to match with the initial conditions for a given problem. This basic requirement has been accomplished by the dual scale micro/macro-crack growth model for creep and fatigue.     

Double diffusive natural convection in a partially heated enclosure with Soret and Dufour effects

The effect of double-diffusive natural convection of water in a partially heated enclosure with Soret and Dufour coefficients around the density maximum is studied numerically. The right vertical wall has constant temperature θ_c, while left vertical wall is partially heated θ_h, with θ_h > θ_c. The concentration in right wall is maintained higher than left wall (C_c < C_h) for case I, and concentration is lower in right wall than left wall (C_h > C_c) for case II. The remaining left vertical wall and the two horizontal walls are considered adiabatic. Water is considered as the working fluid. The governing equations are solved by control volume method using SIMPLE algorithm with QUICK scheme. The effect of the various parameters (thermal Rayleigh number, center of the heating location, density inversion parameter, Buoyancy ratio number, Schmidt number, and Soret and Dufour coefficients) on the flow pattern and heat and mass transfer has been depicted. Comprehensive Nusselt and Sherwood numbers data are presented as functions of the governing parameters mentioned above.     

High-frequency dielectric properties of BSTO ceramic prepared with hydrothermal synthesized SrTiO3 and BaTiO3 powders

BSTO dielectric ceramic was prepared from SrTi03 and BaTi03 powders synthesized by hydrothermal method, as well as from Bao.sSro.4TiO3 powder synthesized by conventional solid-state reaction. The former can be sintered at a relatively low temperature of 1120 ℃. Characterization by SEM showed that the grain shapes of both ceramics are cubical, though the grain size of the former is much smaller. Dielectric constants measured at 20℃ were shown to vary with frequency in the range from I kHz to 2 MHz and dc bias field, and further that the dielectric loss of the former to be less than 2 × 10^- 3 in the frequency range of 20 kHz to 1 MHz, much smaller than that of the latter sample. For the former, temperature dependence of dielectric constant is much flatter and there exists an extended phase transition diffusion covering a wide temperature range of Curie temperature To. The smaller grain size of the former depresses the dc bias electrical field dependence of dielectric constant. The tunability is 7% under a bias field of 0.6 kV/mm dc.     

Temperature fields induced by direct contact condensation of steam in a cross-flow in a channel

The temperature fields in the center plane of a channel with a square cross-section have been measured. Steam injected at relatively low mass fluxes through a small hole in one of the walls of the channel condensed intermittently in a small area close to the inlet. The upstream temperature of the liquid cross-flow, T _L , the momentum ratio, J, and the Prandtl number proved to be important for the single-phase temperature field induced in the jet further away from the steam inlet. Jet centerlines of velocity and temperature are measured and positions are compared. Different locations for J < 100 and low T _L are explained from dependencies on Reynolds and Prandtl numbers. Next to the jet centerline a second high-temperature zone was found to occur, close to the wall and downstream of the steam inlet. The importance of capillary forces is investigated with the aid of 3D CFD computations.     

Oblique stagnation-point flow of a viscoelastic fluid with heat transfer

The two-dimensional forced convection stagnation-point flow and heat transfer of a viscoelastic second grade fluid obliquely impinging on an infinite plane wall is considered as an exact solution of the full partial differential equations. This oblique flow consists of an orthogonal stagnation-point flow to which a shear flow whose vorticity is fixed at infinity is added. The relative importance of these flows is measured by a parameter γ. The viscoelastic problem is reduced to two ordinary differential equations governed by the Weissenberg number W_e, two parameters α and β, the later being a free parameter β, introduced by Tooke and Blyth [A note on oblique stagnation-point flow, Physics of Fluids 20 (2008) 033101-1–3], and the Prandtl number Pr. The two cases when α=β and α≠β are, respectively, considered. Physically the free parameter may be viewed as altering the structure of the shear flow component by varying the magnitude of the pressure gradient. It is found that the location of the separation point x_s of the boundary layer moves continuously from the left to the right of the origin of the axes (x_s<0).     

Three-dimensional flow of Newtonian and Boger fluids in square–square contractions

The flow of a Newtonian fluid and a Boger fluid through sudden square–square contractions was investigated experimentally aiming to characterize the flow and provide quantitative data for benchmarking in a complex three-dimensional flow. Visualizations of the flow patterns were undertaken using streak-line photography, detailed velocity field measurements were conducted using particle image velocimetry (PIV) and pressure drop measurements were performed in various geometries with different contraction ratios. For the Newtonian fluid, the experimental results are compared with numerical simulations performed using a finite volume method, and excellent agreement is found for the range of Reynolds number tested (Re₂ ≤ 23). For the viscoelastic case, recirculations are still present upstream of the contraction but we also observe other complex flow patterns that are dependent on contraction ratio (CR) and Deborah number (De₂) for the range of conditions studied: CR = 2.4, 4, 8, 12 and De₂ ≤ 150. For low contraction ratios strong divergent flow is observed upstream of the contraction, whereas for high contraction ratios there is no upstream divergent flow, except in the vicinity of the re-entrant corner where a localized atypical divergent flow is observed. For all contraction ratios studied, at sufficiently high Deborah numbers, strong elastic vortex enhancement upstream of the contraction is observed, which leads to the onset of a periodic complex flow at higher flow rates. The vortices observed under steady flow are not closed, and fluid elasticity was found to modify the flow direction within the recirculations as compared to that found for Newtonian fluids. The entry pressure drop, quantified using a Couette correction, was found to increase with the Deborah number for the higher contraction ratios.     

CFB cyclones at high temperature: Operational results and design assessment

Pressure drop and cut size measurements are reported for a full scale cyclone operating within a 58 MWth CFB-combustor unit at 775 ℃.
The paper reviews the vast number of equations to calculate the pressure drop and separation efficiency of cyclones, generally for operation at ambient temperature and at low Cs [〈0.5]. None of the literature correlations predicts the pressure drop with a fair accuracy within the range of experimental operating conditions. The cut size d50 can be estimated using direct empirical methods or using the Stokes number, Stk50. Both methods were used to compare measured and predicted values of d50. With the exception of Muschelknautz and Krambrock, none of the equations made accurate predictions.
Finally, an alternative method to determine the friction factor of the pressure drop equation （Euler number, Eu） and of the cut size is proposed. The Eu number is determined from the geometry of common cyclones, and the derived value of Stk50 defines more accurate cut sizes. The remaining discrepancy of less than 5%, when compared with the measured values, is tentatively explained in terms of a reduced cyclone diameter due to the solids layer formed near its wall. Further measurements, mostly using positron emission particle tracking, elucidate the particle motion in the cyclone and both tracking results and the influence of the particle movement on Eu and Stk50 will be discussed in a follow-up paper.     

Oscillating flow in a 2-D diffuser

Separating oscillating flows in an internal, adverse pressure gradient geometry are studied experimentally. Simultaneous velocity and pressure measurements demonstrate that the minor losses associated with oscillating flow in an adverse pressure gradient geometry can be smaller or larger than those for steady flow. Separation is found to begin high in the diffuser and propagate downward. The flow is able to remain attached further into the diffuser with larger Reynolds numbers, small displacement amplitudes, and smaller diffuser angles. The extent of separation grows with L ₀/h. The minor losses grow with increasing displacement amplitude in the measured range 10 < L ₀/h < 40. Losses decrease with increasing Re _δ in the measured range of 380 < Re _δ < 740. It is found that the losses increase with increasing diffuser angle over the measured range of 12° < θ < 30°. The nondimensional acoustic power dissipation increases with Reynolds number in the measured range and decreases with displacement amplitude.     

The ultra-low Reynolds number airfoil wake

Lift force and the near wake of an NACA 0012 airfoil were measured over the angle (α) of attack of 0°–90° and the chord Reynolds number (Re _c ), 5.3 × 10³–5.1 × 10⁴, with a view to understand thoroughly the near wake of the airfoil at low- to ultra-low Re _c . While the lift force is measured using a load cell, the detailed flow structure is captured using laser-Doppler anemometry, particle image velocimetry, and laser-induced fluorescence flow visualization. It has been found that the stall of an airfoil, characterized by a drop in the lift force, occurs at Re _c  ≥ 1.05 × 10⁴ but is absent at Re _c  = 5.3 × 10³. The observation is connected to the presence of the separation bubble at high Re _c but absence of the bubble at ultra-low Re _c , as evidenced in our wake measurements. The near-wake characteristics are examined and discussed in detail, including the vortex formation length, wake width, spanwise vorticity, wake bubble size, wavelength of K–H vortices, Strouhal numbers, and their dependence on α and Re _c .     

Fractal model for virgin compression of pure clays

The well-known linear relation between the void ratio and the logarithm of pressure for the one-dimensional consolidation has widely incorporated into elastoplastic constitute equation of soils, however, it is not at all straight line for more compressible clays and has several physically unaccepted properties. The compression of clays is explained by fractal theory, and is seen as a balance between the vertical pressure and repulsive force between clay plates. The distance between clay plates depends on the structure formation of the clay flocculation and clay surface. The structure formation of the clay flocculation and clay surface can be modeled by fractal approach. A new formulation of one-dimensional virgin compression of pure clays, expressed by the linear function of ln ν–ln p or ln e–ln p, is proposed based on the fractal model for clay structure. The linear function of ln ν–ln p is suitable to express the virgin compression of the extremely high plasticity clays up to the high pressure of 5000 kPa. For low plasticity clays, the virgin compression satisfies with the linear relationship of ln e–ln p at the pressure higher than 100 kPa. These too are show to be in satisfactory agreement with experimental data. The gradient of the plots of ln ν–ln p and ln e–ln p is related to the fractal dimensions of clay structure, the fundamental material parameters for the first time.     

Study of vibrational deactivation of molecules of carbon dioxide gas during cooling of stream in a supersonic nozzle

The vibrational temperature of the antisymmetrical type of vibrations (v ₃) of the CO₂ molecule at the exit of a supersonic nozzle is measured in the present work using the method of recording the infrared emission. Freezing in of thev ₃-type vibrations was observed during the flow of undiluted carbon dioxide in a nozzle. In this case the vibrational temperature T₃ considerably exceeded the translational temperature. On the basis of a comparison of the experimental results with calculation it can be concluded that vibrational deactivation of CO₂ molecules occurs three to five times faster than the excitation of the vibrations during heating in a shock wave. All the experiments were conducted under the following conditions: maximum expansion of gas in nozzle A/A* = 115, temperature range 1900–2400 °K, pressure range 1–17.5 atm.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 6, pp. 32–40, November–December, 1973.The authors are grateful to U. G. Pirumov and É. A. Ashratov for the calculation of the nozzle profile and the distribution of streamlines as well as for a discussion of the results.     

Experimental study of sedimentation characteristics of spheroidal particles

The drag of non-spherical particles is a basic, important parameter for multi-phase flow. As the first step in research in this area, the terminal velocities, Ut, of hemispherical and spherical segment particles with maximal diameters of 6-21 mm were measured in static fluids by using a high-speed video camera. The drag coefficient, CD, measured for Reynolds number, Re of 10^1-10^5, has been obtained and compared with those for a sphere. The Re based on the terminal velocity has a logarithmic linear relationship with Ar number for both the facet facing upwards or downwards for the two experimental spheroidal particles, and their Co values are greater than those of spheres. A shape function that depends on the initial orientation of the particle facet is presented to correct for the shape effects.     

Absorption spectroscopy of mid and neighboring Z plasmas: Iron, nickel,copper and germanium

Opacities of four medium Z element plasmas (iron, nickel, copper and germanium) have been measured at the LULI-2000 facility in similar conditions: temperatures between 15 and 25 eV and densities between 2 and 10 mg/cm³, in a wavelength range (8–18 Å) including the strong 2p–3d structures.Two laser beams from the LULI facility were used in the nanosecond-picosecond configuration. The NANO-2000 beam (at λ = 0.53 μm) heated a gold hohlraum with an energy between 30 and 150 J with a duration of 0.6 ns. Samples covering half a hohlraum hole were thus radiatively heated. The picosecond pulse PICO-2000 beam (at λ = 1.053 μm) has been used to produce a short (about 10 ps) X-ray backlighter in order to reduce time variations of temperatures and densities during the measurement. A crystal high-resolution spectrometer was used as the main diagnostic to record at the same time the non-absorbed and the absorbed backlighter spectra. Radiation temperatures were measured using a broadband spectrometer. 1D and 2D simulations have been performed in order to estimate hydrodynamic plasmas parameters.The measured spectra have been compared with theoretical ones obtained using either the superconfiguration code SCO or the detailed term accounting code HULLAC. These comparisons allow us to check the modeling of the statistical broadening and of the spin-orbit splitting of the 2p–3d transitions and related effects such as the interaction between relativistic subconfigurations belonging to the same non-relativistic configuration.     

Comparisons of different procedures of pre-compaction stress determination on weakly structured soils

Compaction is an important component of soil degradation. In this regard, the pre-compaction stress (σ_pc) concept is considered useful in mechanized agriculture nowadays. When the external forces exceed the internal strength (σ_pc) of soil, soil structure and soil physical quality will deteriorate. This concept was introduced at first for confined consolidation of non-structured, homogenized and saturated subsoils in civil engineering, though it is also suitable for agricultural conditions where the topsoil and subsoil are considered and both are often structured, heterogeneous and unsaturated. The best method for predicting σ_pc is by the plate sinkage test (PST) in the field, but it is expensive and time-consuming. This study was conducted to find an alternative laboratory method besides the confined compaction test (CCT) for predicting σ_pc. The CCT may not be a good method, especially at higher water contents, and for soils with low organic matter content, because of low sharpness of the critical region on the stress–strain curves. The study was performed on five soil types with a range of soil textures and organic matter content from central Iran using three loading types and three pF (i.e. Log [soil matric suction in cm]) values of 2.3, 2.7 and 2.9 with two replicates. Loading types consisted of CCT, the semi-confined compaction test (SCCT) and the kneading compaction test (KCT) at three maximum (or pre-compaction) stresses of 200, 400 and 600 kPa. The experiment was a completely randomized factorial design. The aim was to determine how accurately each loading type test could predict σ_pc of soil pre-compacted by one of the other methods. The applied combinations of CCT–SCCT, SCCT–CCT and KCT–CCT mean that the soil was pre-compacted by the first loading type and evaluated by the second one. The results showed that σ_pc and the sharpness of the σ_pc region were significantly affected by loading types as well as the soil conditions. The sharpest σ_pc region was observed in SCCT and the least sharp in CCT with the overall order being CCT–SCCT > SCCT–CCT > KCT–CCT. The sharpness of the σ_pc region was reduced for the soil samples with higher water content and coarser texture. Regardless of the soil and loading conditions, the prediction by SCCT was consistently more accurate than by CCT. The prediction of σ_pc by SCCT was more precise in comparison with CCT especially at higher stresses and soil water contents. However, the prediction of σ_pc by SCCT was very accurate at pF values of 2.7 and 2.9, and with low σ_pc values, when compared with the actual values of the σ_pc. For the clay soil at a pF value of 2.3, no pre-compaction region (i.e. zero σ_pc) could be determined by CCT at a maximum axial stress of 600 kPa. This was because of the incompressibility of soil water at this near-saturated soil condition at high stress. However, the sharpness of the critical region in SCCT was high enough to predict σ_pc satisfactorily. There was no significant difference between the combinations of SCCT–CCT and KCT–CCT in predicting σ_pc. The SCCT is a compromise method that lies between CCT and PST. SCCT has the advantage of using a limited and definite soil volume that can be modeled as a soil element. Marginal effects of disturbance caused by coring/sampling as well as pre-test sample preparation seem to have minor effects on the stress–strain curve determined by SCCT in comparison with CCT. Moreover, the soil volume needed for this test is the same as for CCT.     

Active flow control for reduction of fluctuating aerodynamic forces of a blunt trailing edge profiled body

Vortex shedding in the wake of two-dimensional bluff bodies is usually accompanied by three dimensional instabilities. These instabilities result in streamwise and vertical vorticity components which occur at a certain spanwise wavelength. The spanwise wavelength of the instabilities (λ_Z) depends on several parameters, including profile geometry and Reynolds number. The objective of the present work is to study the three dimensional wake instabilities for a blunt trailing edge profiled body, comprised of an elliptical leading edge and a rectangular trailing edge, and to manipulate these instabilities to control the aerodynamic forces. Results of numerical simulations of flow around the body at Re(d) = 400, 600, and 1000, as well as planar Laser Induced Fluorescence (LIF) flow visualizations at Re(d) = 600 and 1000 are analyzed to determine the wake vorticity structure and λ_Z. Based on the findings of these analyses, an active flow control mechanism for attenuation of the fluctuating aerodynamic forces on the body is proposed. The flow control mechanism is comprised of a series of trailing edge injection ports distributed across the span, with a spacing equal to λ_Z. Injection of a secondary flow leads to amplification of the three dimensional instabilities and disorganization of the von Kármán vortex street. Numerical simulations indicate that the flow control mechanism can attenuate the fluctuating aerodynamic forces at lower Reynolds numbers (Re(d) = 400 and 600) where λ_Z is constant in time. However, the control mechanism loses its effectiveness at Re(d) = 1000, due to the temporal variations of λ_Z.     

Experimental determination of the free-stream disturbance field in a short-duration supersonic wind tunnel

The free-stream disturbance field in a short-duration supersonic wind tunnel is investigated at a nominal Mach number of Ma=2.54. A specially designed constant-temperature anemometer is used to be able to draw a complete fluctuation diagram within one wind tunnel run (testing time: 120 ms). It is shown that the disturbance field is dominated by acoustic waves radiated from the turbulent boundary layer on the nozzle and the sidewalls, like in conventional supersonic wind tunnels. The acoustic field appears to be composed of highly localized shivering Mach waves superimposed on a background of eddy Mach waves.Abbreviations  a constant in the thermal conductivity/temperature power law of air: k/k_r=(T/T_r)^a (dimensionless) - b constant in the viscosity/temperature power law of air: /_r=(T/T_r)^b (dimensionless) - Be bandwidth (Hz) - A, B constants in the wire heat transfer relation (Eq. (7), dimensionless) -   (dimensionless) - c_p specific heat at constant temperature (kJ/kg K) - c_v specific heat at constant volume (kJ/kg K) -  boundary layer thickness (m) - D function of the overheat ratio (dimensionless) - e anemometer output voltage (V) - _F end-loss attenuation factor for mass flow sensitivity (dimensionless) - _G end-loss attenuation factor for total temperature sensitivity (dimensionless) -  recovery factor (dimensionless) - f frequency (Hz) - f₁ normalized frequency (dimensionless) - F anemometer nondimensional sensitivity to mass flow fluctuations (dimensionless) - G anemometer nondimensional sensitivity to total temperature fluctuations (dimensionless) - F_AC F×F (dimensionless) - G_AC G×G (dimensionless) - f,g functions in the wire heat transfer relation (Eq. (7), dimensionless) -  c_p/c_v (dimensionless) - k thermal conductivity of air (W/m K) - k_r thermal conductivity of air at temperature T_r (W/m K) - k anemometer sensitivity to total temperature fluctuations (V/K) - l Mach rhombus half-length (Fig. 1, m) - Ma Mach number (dimensionless) -  viscosity of air (kg/m·s) - _r viscosity of air at temperature T_r (kg/m·s) - n constant in the wire heat transfer relation (Eq. (7), dimensionless) - Nu Nusselt number (dimensionless) - p pressure (Pa) - p₀ stagnation pressure (Pa) - r –F/G (dimensionless) - R unit Reynolds number (1/m) - Re Reynolds number (dimensionless) -  correlation coefficient between mass flow and total temperature fluctuations (dimensionless) -  density (kg/m³) - T time span (s) - T₀ total temperature (K) - T_r reference temperature (K) - T_w hot wire temperature (K) -  overheat ratio: =(T_w–T₀)/T₀ (dimensionless) -  –<e>/G (%) - u x-component of the flow velocity (m/s) - u_s source velocity at acoustic origin (m/s) - u inviscid velocity at acoustic origin (m/s) - x wind tunnel axis (Fig. 1, m)Symbols  x̄ temporal mean value of a fluctuating quantity x - x fluctuating part of x: x=x–x̄ - x_RMS^' root mean square of x - <x> x_RMS^'/x̄ - (X) relative uncertainty of a random variable X

---

Validation of a discrete element model using magnetic resonance measurements

The discrete element model （DEM） is a very promising modelling strategy for two-phase granular systems. However, owing to a lack of experimental measurements, validation of numerical simulations of two-phase granular systems is still an important issue. In this study, a small two-dimensional gas- fluidized bed was simulated using a discrete element model. The dimensions of the simulated bed were 44mm × 10mm × 120 mm and the fluidized particles had a diameter dp = 1.2 mm and density ρp = 1000 kg/m^3. The comparison between DEM simulations and experiments are performed on the basis of time-averaged voidage maps. The drag-law of Beetstra et al. [Beetstra, R., van der Hoef, M.A., ＆ Kuipers,J. A. M. （2007b）. Drag force of intermediate Reynolds number flow past mono- and bidispersed arrays of spheres. AIChE Journal, 53,489-501 ] seems to give the best results. The simulations are fairly insensitive to the coefficient of restitution and the coefficient of friction as long as some route of energy dissipation during particle-particle and particle-wall contact is provided. Changing the boundary condition of the gas phase at the side-walls from zero-slip to full-slip does not affect the simulation results. Care is to be taken that the cell sizes are chosen so that a reasonable number of particles can be found in a fluid cell.     

2-D Simulation of Non-Isothermal Fate and Transport of a Drip-Applied Fumigant in Plastic-Mulched Soil Beds. II. Sensitivity Analysis and Model Application

Temperature dependence of selected parameters was defined and illustrated. The order of sensitivity among tested model parameters (sorption coefficient, K _D; diffusion coefficient, D _g; vapor pressure, P; degradation rate, k; Henry’s constant, H) was k > P > K _D > H > D _g for 1-day simulations and k > P and K _D > D _g > H for 2-day simulations. The most sensitive parameter for both days was k. Mass balance errors were calculated using NRMSE (normalized root mean square error) to evaluate the level of agreement. The numerical model was applied to predict fumigant movement in the soil bed.     

Analytical and numerical analysis for the FGM thick sphere under combined pressure and temperature loading

Thermo-mechanical analysis of functionally graded hollow sphere subjected to mechanical loads and one-dimensional steady-state thermal stresses is carried out in this study. The material properties are assumed to vary non-linearly in the radial direction, and the Poisson’s ratio is assumed constant. The temperature distribution is assumed to be a function of radius, with general thermal and mechanical boundary conditions on the inside and outside surfaces of the sphere. In the analysis presented here, the effect of non-homogeneity in FGM thick sphere was implemented by choosing a dimensionless parameter, named β _i (i = 1, . . . , 3), which could be assigned an arbitrary value affecting the stresses in the sphere. It is observed that the solution process for β _i (i = 3) = −1 are different from those obtained for other values of β _i (i = 1, . . . , 3). Cases of pressure, temperature, and combined loadings were considered separately. It is concluded that by changing the value of β _i (i = 1 . . . 3), the properties of FGM can be so modified that the lowest stress levels are reached. The present results agree well with existing results. Using FEM simulations, the analytical findings for FGM spheres under the influence of internal pressure and temperature gradient were compared to finite element results.