Erratum to: Rock Dynamic Fracture Toughness Tested with Holed-Cracked Flattened Brazilian Discs Diametrically Impacted by SHPB and its Size Effect

Dynamic fracture initiation toughness of marble was tested using two types of the holed-cracked flattened Brazilian disc (HCFBD) specimens, which were diametrically impacted at the flat end of the disc by the split Hopkinson pressure bar (SHPB) of 100 mm diameter. One type of the discs is geometrically similar with different outside diameter of 42 mm, 80 mm, 122 mm and 155 mm respectively, and with crack length being half the diameter; another type of the discs has identical 80 mm diameter and different crack length. Issues associated with determination of the stress wave loading by the SHPB system and the crack initiation time in the disc specimen were resolved using strain gage technique. The stress waves recorded on the bars and the disc failure patterns are shown and explained. The tested dynamic fracture toughness increases obviously with increasing diameter for the geometrically similar HCFBD specimens. It changes moderately for the one-size specimens of identical diameter and different crack length. The size effect of rock dynamic fracture toughness is mainly caused by the fracture process zone length l and fracture incubation time τ, the latter being an additional influencing factor for the dynamic loading as compared with the counterpart static situation. Hence a method is proposed to determine a unique value for the dynamic fracture initiation toughness, the approach takes average of the local distribution and time history for dynamic stress intensity factor in the spatial-temporal domain, which is defined by l and τ jointly. In this way the dynamic size effect is minimized.     

Modeling of macroscopic response of phase transforming materials under quasi-static loading

Based on a continuum model of solid-solid phase transformations, the macroscopic response of a bar of a thermoelastic phase transforming material loaded quasistatically is investigated. A critical loading rate is identified for the evolution of a single phase boundary in the bar during an isothermal process. It is shown that, when the loading rate is larger than this critical loading rate, nucleation occurs either continuously or at multiple sites; when the loading rate is lower than this critical loading rate, the size of the hysteresis loop increases with increasing loading rate, and decreases with an increase in the mobility of the phase boundary. The heat conduction due to the heat generated by the latent heat of the phase transformation is considered for a special case.     

Assessment of Non-adiabatic Behaviour in Thermoelastic Stress Analysis of Composite Sandwich Panels

Thermoelastic stress analysis (TSA) is used to derive the surface stresses in large sandwich structure panels with honeycomb core and carbon fibre face sheets. The sandwich panels are representative of those used for secondary aircraft structure. The panels were subjected to a pressure load, similar to that experienced in-service, using a custom designed test rig. To achieve the necessary adiabatic conditions for TSA, cyclic loading is regarded as an essential feature. As the panels were full-scale, the maximum loading frequency that could be imparted to the panels by the rig was 1 Hz, which is below the usual range recommended to achieve adiabatic behaviour. To assess the effectiveness of TSA at low frequencies two approaches to calibration are investigated and compared with the stress distribution obtained from independently validated FE models. The thermoelastic response was calibrated into stress data using thermoelastic constants derived experimentally from tensile strips of the sandwich panel face sheet material. It is shown that by using thermoelastic constants obtained from the tensile strips manufactured with the same lay-up as the sandwich panel face sheets, and at the same cyclic load frequency used in the full-scale tests, quantitative stress metrics can be derived from the TSA data. More significantly, a deeper insight into the importance of the thermal characteristics in TSA of laminated materials is provided. It is demonstrated that, for the material used in this work, it is possible to use the global material behaviour to obtain quantitative results when adiabatic conditions do not prevail.     

Transient thermal cracking associated with non-classical heat conduction in cylindrical coordinate system

This paper studies the fracture behavior of a thermoelastic cylinder subjected to a sudden temperature change on its outer surface within the framework of non-classical heat conduction.The heat conduction equation is solved by separation of variable technique.Closed form solution for the temperature field and the associated thermal stress are established.The critical parameter governing the level of the transient thermal stress is identified.Exact expression for the transient stress intensity factor is obtained for a crack in the cylinder.The difference between the non-classical solutions and the classical solution are discussed.It is found that the traditional classical heat conduction considerably underestimates the transient thermal stress and thermal stress intensity factor.     

On the Stability of Thermoelastic Materials

A notion of material stability is introduced and discussed in the setting of nonlinear thermoelasticity. Necessary and sufficient conditions are established for the stability of a general thermoelastic material. The adiabatic theory and the theory that accounts for heat conduction are considered separately. This revised version was published online in August 2006 with corrections to the Cover Date.     

轴对称横观各向同性热弹性圆柱的分解

为了得到精确的应力场、位移场、温度场,将扭转圆轴的精化理论研究方法推广到轴对称横观各向同性热弹性圆柱。利用Bessel函数以及轴对称横观各向同性热弹性圆柱的通解,给出了轴对称横观各向同性热弹性圆柱的分解定理。根据柱面齐次边界条件获得了精确的精化方程,精化方程可以分解为一阶方程、超越方程、温度方程,从而将横观各向同性热弹性圆柱的轴对称问题分解为轴向拉压问题、超越问题、热-应力耦合问题。超越部分对应端部自平衡情况,可以清晰地了解到端部应力分布对内部应力场的影响,热-应力耦合部分对应无外加应力场时圆柱内部因温度变化引起的热应力。     

超短激光脉冲加热薄板的广义热弹扩散问题

由于超短激光脉冲具有功率密度高、持续时间短、加工精度高等优势,近年来被广泛应用于超精细加工、光学储存和微电子器件制造等领域.本文基于L-S型广义热弹扩散理论,建立了考虑材料记忆依赖效应和空间非局部效应的记忆依赖型非局部广义热弹扩散耦合理论,它能够准确预测几何尺寸与内部特征尺寸相近结构的热弹扩散瞬态响应.推导了所建理论的控制方程,并基于拉普拉斯积分变换获得了控制方程的解.作为算例,利用所建理论和求解方法研究了半无限大薄板受非高斯激光脉冲加热和化学冲击联合作用下的热弹扩散瞬态响应问题,得到了薄板的温度、化学势、位移、应力和浓度等随非局部参数、热时间迟滞因子和扩散时间迟滞因子等参数变化的分布规律.结果表明:传热对传质影响显著,传质对传热影响甚微;非局部参数对位移、应力影响显著,对温度、化学势和浓度几乎没有影响.该理论及求解方法的建立,旨在实现材料在机械、热、化学势等冲击作用下传热传质瞬态响应的准确预测.     

A GENERALIZED THERMOELASTIC DIFFUSION PROBLEM OF THIN PLATE HEATED BY THE ULTRASHORT LASER PULSES 1)

由于超短激光脉冲具有功率密度高、持续时间短、加工精度高等优势, 近年来被广泛应用于超精细加工、光学储存和微电子器件制造等领域. 本文基于L-S型广义热弹扩散理论, 建立了考虑材料记忆依赖效应和空间非局部效应的记忆依赖型非局部广义热弹扩散耦合理论, 它能够准确预测几何尺寸与内部特征尺寸相近结构的热弹扩散瞬态响应. 推导了所建理论的控制方程, 并基于拉普拉斯积分变换获得了控制方程的解. 作为算例, 利用所建理论和求解方法研究了半无限大薄板受非高斯激光脉冲加热和化学冲击联合作用下的热弹扩散瞬态响应问题, 得到了薄板的温度、化学势、位移、应力和浓度等随非局部参数、热时间迟滞因子和扩散时间迟滞因子等参数变化的分布规律. 结果表明: 传热对传质影响显著, 传质对传热影响甚微; 非局部参数对位移、应力影响显著, 对温度、化学势和浓度几乎没有影响. 该理论及求解方法的建立, 旨在实现材料在机械、热、化学势等冲击作用下传热传质瞬态响应的准确预测.     

Rational design of thermoelastic damping in microresonators with phase-lagging heat conduction law

The design of thermoelastic damping (TED) affected by the phase-lagging non-Fourier heat conduction effects becomes significant but challenging for enlarging the quality factor of widely-used microresonators operating in extreme situations, including ultra-high excitation frequency and ultra-low working temperature. However, there does not exist a rational method for designing the TED in the framework of non-Fourier heat conduction law. This work, therefore, proposes a design framework to achieve low thermoelastic dissipation of microresonators governed by the phase-lagging heat conduction law. The equation of motion and the heat conduction equation for phase-lagging TED microresonators are derived first, and then the non-Fourier TED design problem is proposed. A topology optimization-based rational design method is used to resolve the design problem. What is more, a two-dimensional (2D) plain-strain-based finite element method (FEM) is developed as a solver for the topology optimization process. Based on the suggested rational design technique, numerical instances with various phase lags are investigated. The results show that the proposed design method can remarkably reduce the dissipation of microresonators by tailoring their substructures.     

Evaporation flow pattern and heat transfer of R-22 and R-134a in small diameter tubes

The flow patterns and heat transfer coefficients of R-22 and R-134a during evaporation in small diameter tubes were investigated experimentally. The evaporation flow patterns of R-22 and R-134a were observed in Pyrex sight glass tubes with 2 and 8 mm diameter tube, and heat transfer coefficients were measured in smooth and horizontal copper tubes with 1.77, 3.36 and 5.35 mm diameter tube, respectively. In the flow patterns during evaporation process, the annular flows in 2 mm glass tube occurred at a relatively lower vapor quality compared to 8 mm glass tube. The flow patterns in 2 mm glass tube did not agree with the Mandhane’s flow pattern maps. The evaporation heat transfer coefficients in the small diameter tubes (d _i  < 6 mm) were observed to be strongly affected by tube diameters, and to differ from those in the large diameter tubes. The heat transfer coefficients of 1.77 mm tube were higher than those of 3.36 mm and 5.35 mm tube. Most of the existing correlations failed to predict the evaporation heat transfer coefficient in small diameter tubes. Therefore, based on the experimental data, the new correlation is proposed to predict the evaporation heat transfer coefficients of R-22 and R-134a in small diameter tubes.     

Free convection heat transfer from a horizontal fin attached cylinder between confined nearly adiabatic walls

Steady state two-dimensional free convection heat transfer from a horizontal, isothermal fin attached cylinder, located between nearly two adiabatic walls is studied experimentally using a Mach–Zehnder interferometer. Effects of the walls inclination angel (θ) on heat transfer from the cylinder is investigated for Rayleigh number ranging from 1000 to 15,500. Two cylinders with different diameters of D = 10 and 20 mm are used to cover wide Rayleigh range. Results indicate that, heat transfer phenomena differ for different Rayleigh number. For Rayleigh numbers lower than 5500, heat transfer rate from cylinder surface is lower than the heat transfer from a single cylinder. In this range by the use of walls, heat transfer from the cylinder decreases slightly and walls’ inclination does not change heat transfer rate from the cylinder surface. For Rayleigh number ranging from 5500 to 15,500, amount of heat transfer from the cylinder surface is less than that of a single cylinder. However, by adding nearly adiabatic walls to experimental model heat transfer mechanism differs and chimney effect between fin and walls increases the heat transfer rate from the cylinder surface. By increasing the walls inclination angel from 0° to 20°, the chimney effect between walls and fin diminishes and heat transfer rate from the cylinder surface is approaching to the heat transfer rate of fin attached cylinder without adiabatic walls.     

Image enhancement and stress separation of thermoelastically measured data under random loading

Traditional thermoelastic stress analysis (TSA) presupposes that the structure being analyzed is cyclically loaded at a constant amplitude and frequency. This approach typically has been used to satisfy the adiabatic reversible assumptions. The authors employ an alternative signal analysis technique that enables one to evaluate the magnitude of the individual components of stress in a component subjected to a loading that is random in both frequency and magnitude. However, the nature of the measured information does not change; i.e., data are inherently noisy, and edge information is unreliable. The latter two aspects have caused many thermoelastic stress analyses to be more qualitative than quantitative. The present paper emphasizes developing the TSA technique into a practical, noncontacting quantitative method for stress analyzing actual engineering structures that are randomly loaded. In particular, ability to determine the individual stresses thermoelastically under random loading is demonstrated.     

Adiabatic shear banding in plane strain tensile deformations of 11 thermoelastoviscoplastic materials with finite thermal wave speed

We study the initiation and propagation of adiabatic shear bands (ASBs) in 11 homogeneous materials each modeled as microporous, isotropic and thermoelastoviscoplastic, and deformed in plane strain tension. The heat conduction in each material is assumed to be governed by a hyperbolic heat equation; thus thermal and mechanical waves propagate with finite speeds. The decrease in the thermophysical parameters due to the increase in porosity is considered. An ASB is assumed to initiate at a material point when the maximum shear stress there has dropped to 80% of its peak value for that material point and it is deforming plastically. An approximate solution of the coupled nonlinear partial differential equations subject to suitable initial and boundary conditions is found by the finite element method (FEM). In contrast to the Considerè and the Hart criterion, it is found that an ASB initiates when the axial load drops rapidly and not when it peaks. The refinement of the 40 × 40 uniform FE mesh to 120 × 120 uniform elements decreased the ASB initiation time by 2.1% while increasing the CPU time by a factor of ∼26. By locating points where the ASB has initiated we find its current length, width and speed. The 11 materials are ranked according to the time of initiation of an ASB under otherwise identical geometric and loading conditions with the same initial nonuniform porosity distribution. This ranking of materials is found to differ somewhat from that ascertained by Batra and Kim (1992) who studied simple shearing deformations, and by Batra et al. (1995) who analyzed three-dimensional torsional deformations of thin-walled tubular specimens. The average axial strain determined from the maximum axial load condition differs noticeably from that when an ASB initiates.     

Thermoelastic state of a half-space having a thermally active crack parallel to its boundary

Using the boundary integral equation method, the problem of stationary heat conduction and thermoelasticity for a semi-infinite body with a crack parallel to its boundary is solved. Temperature or heat flow on the crack is prescribed. The body boundary is heat-insulated or is at zero temperature. The dependence of the stress intensity factor on the depth of occurrence of a circular crack at a constant temperature or under a constant heat flow is studied. In contrast to mechanical loading, thermal loading shows less SIF values than in an infinite body __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 4, pp. 46–54, April 2007.     

Rock Dynamic Fracture Toughness Tested with Holed-Cracked Flattened Brazilian Discs Diametrically Impacted by SHPB and its Size Effect

Dynamic fracture initiation toughness of marble was tested using two types of the holed-cracked flattened Brazilian disc (HCFBD) specimens, which were diametrically impacted at the flat end of the disc by the split Hopkinson pressure bar (SHPB) of 100 mm diameter. One type of the discs is geometrically similar with different outside diameter of 42 mm, 80 mm, 122 mm and 155 mm respectively, and with crack length being half the diameter; another type of the discs has identical 80 mm diameter and different crack length. Issues associated with determination of the stress wave loading by the SHPB system and the crack initiation time in the disc specimen were resolved using strain gage technique. The stress waves recorded on the bars and the disc failure patterns are shown and explained. The tested dynamic fracture toughness increases obviously with increasing diameter for the geometrically similar HCFBD specimens. It changes moderately for the one-size specimens of identical diameter and different crack length. The size effect of rock dynamic fracture toughness is mainly caused by the fracture process zone length l and fracture incubation time τ, the latter being an additional influencing factor for the dynamic loading as compared with the counterpart static situation. Hence a method is proposed to determine a unique value for the dynamic fracture initiation toughness, the approach takes average of the local distribution and time history for dynamic stress intensity factor in the spatial-temporal domain, which is defined by l and τ jointly. In this way the dynamic size effect is minimized.     

Effects of three-phase-lag on two-temperature generalized thermoelasticity for infinite medium with spherical cavity

The thermoelastic interaction for the three-phase-lag (TPL) heat equation in an isotropic infinite elastic body with a spherical cavity is studied by two-temperature generalized thermoelasticity theory (2TT). The heat conduction equation in the theory of TPL is a hyperbolic partial differential equation with a fourth-order derivative with respect to time. The medium is assumed to be initially quiescent. By the Laplace transformation, the fundamental equations are expressed in the form of a vector-matrix differential equation, which is solved by a state-space approach. The general solution obtained is applied to a specific problem, when the boundary of the cavity is subjected to the thermal loading (the thermal shock and the ramp-type heating) and the mechanical loading. The inversion of the Laplace transform is carried out by the Fourier series expansion techniques. The numerical values of the physical quantity are computed for the copper like material. Significant dissimilarities between two models (the two-temperature Green-Naghdi theory with energy dissipation (2TGN-III) and two-temperature TPL model (2T3phase)) are shown graphically. The effects of two-temperature and ramping parameters are also studied.     

On a general treatment for coupled thermal stress problems in a cylindrical coordinate system

Summary In this paper, we study the three-dimensional coupled thermoelastic problem governed by the cylindrical coordinate system under the consideration of the thermo-mechanical coupling effect. The basic relations for the three-dimensional displacement field and the corresponding stress components are derived from the field equations of motion. These expressions enable us to determine both the temperature and the stress distributions simultaneously. As an illustration, numerical calculations are carried out for an axisymmetrical coupled thermoelastic problem. From the numerical results it is shown that the thermomechanical coupling term appearing in the corrected heat conduction equation gives a fairly large difference between the coupled theory and the uncoupled one.

---

Zur allgemeinen Behandlung gekoppelter Wärmespannungsprobleme in Zylinderkoordinaten

Übersicht Es wird das dreidimensionale Wärmespannungsproblem unter Berücksichtigung der thermomechanischen Kopplung untersucht. Ausgehend von den Bewegungsgleichungen werden Temperatur und Spannungen angegeben. Zur Illustration der Methode wird ein axialsymmetrisches Problem numerisch behandelt und der Einfluß der thermo-meehanischen Kopplung auf die Wärmeleitgleichung gezeigt.

     

Flow boiling heat transfer of R134a and R404A in a microfin tube at low mass fluxes and low heat fluxes

An experimental investigation of flow boiling heat transfer in a commercially available microfin tube with 9.52 mm outer diameter has been carried out. The microfin tube is made of copper with a total fin number of 55 and a helix angle of 15°. The fin height is 0.24 mm and the inner tube diameter at fin root is 8.95 mm. The test tube is 1 m long and is electrically heated. The experiments have been performed at saturation temperatures between 0 and −20°C. The mass flux was varied between 25 and 150 kg/m²s, the heat flux from 15,000 W/m² down to 1,000 W/m². All measurements have been performed at constant inlet vapour quality ranging from 0.1 to 0.7. The measured heat transfer coefficients range from 1,300 to 15,700 W/m²K for R134a and from 912 to 11,451 W/m²K for R404A. The mean heat transfer coefficient of R134a is in average 1.5 times higher than for R404A. The mean heat transfer coefficient has been compared with the correlations by Koyama et al. and by Kandlikar. The deviations are within ±30% and ±15%, respectively. The influence of the mass flux on the heat transfer is most significant between 25 and 62.5 kg/m²s, where the flow pattern changes from stratified wavy flow to almost annular flow. This flow pattern transition is shifted to lower mass fluxes for the microfin tube compared to the smooth tube.     

On the nature of the relaxation time,the Maxwell–Cattaneo and Fourier law in the thermodynamics of a continuous medium,and the scale effects in thermal conductivity

We present the theory of space–time elasticity and demonstrate that it is the extended reversible thermodynamics and gives the coupled model of thermoelasticity and heat conductivity and involves traditional thermoelasticity. We formulate the generally covariant variational model’s dynamic thermoelasticity and heat conductivity in which the basic kinematic and static variables are unified tensor objects (subject, matter). Variation statement defines the whole set of the initial-boundary problems for the 4D vector governing equation (Euler equation), the spatial projections of which define motion equations and the time projection gives the heat conductivity equation. We show that space–time elasticity directly implies the Fourier and the Maxwell–Cattaneo laws of heat conduction. However, space–time elasticity is richer than classical thermoelasticity, and it advocates its own equations of motion for coupled thermoelasticity. Moreover, we establish that the Maxwell–Cattaneo law and Fourier law can be defined for the reversible processes as compatibility equations without introducing dissipation. We argue that the present framework of space–time elasticity should prove adequate to describe the thermoelastic phenomena at low temperatures for interpreting the results of molecular simulations of heat conduction in solids and for the optimal heat and stress management in the microelectronic components and the thermoelectric devices.     

Performance characteristics of HFC-134a and HFC-410A refrigeration system using a short-tube orifice as an expansion device

In the present article, the effect of heat source temperature, heat sink temperature, short-tube orifice diameter and short-tube orifice length on the performance characteristics of HFC-140A and HFC-134a refrigeration system using a short-tube orifice as expansion device, i.e., mass flow rate, cooling capacity, compressor pressure ratio, power consumption, and second law efficiency are experimentally studied. The short-tube orifices diameters ranging from 0.849 to 1.085 mm with length ranging from 10 to 20 mm are used in this examination. The test run are done at heat source temperature ranging between 16.5 and 18.5°C, and heat sink temperature ranging between 30 and 35°C. The results show that the tendency of second law efficiency is increased as the short-tube orifice diameter and heat source temperature are enhanced, but it is decreased by increasing the short-tube orifice length and heat sink temperature. Under the similar conditions, the mass flow rate, cooling capacity, and compressor power consumption obtained from HFC-410A are higher than those obtained from HFC-134a.