Transient thermopiezoelectric response of a one-dimensional functionally graded piezoelectric medium to a moving heat source

The thermopiezoelectricity problem of a one-dimensional (1-D), finite length, functionally graded medium excited by a moving heat source is investigated in this paper. The Lord and Shulman theory of generalized coupled thermoelasticity is employed to account for both the finite speed of thermal waves and coupling of temperature field with displacement and electric fields. Except thermal relaxation time and specific heat, which are taken to be constant for simplicity, all other properties are assumed to vary exponentially along the length through an arbitrary non-homogeneity index. Laplace transform has been used to eliminate the time effect, and three coupled fields, namely, displacement, temperature, and electric fields are obtained analytically in the Laplace domain. The solutions are then inverted to time domain using a numerical Laplace inversion method. Numerical examples are displayed to illustrate the effects of non-homogeneity index, length and thermal relaxation time on the results. When the medium is homogeneous, the results of the current paper are reduced to exactly the same results available in the literature.     

Fundamental solutions of the thermoelastic equations for transversely isotropic plates

The problem of thermoelasticity for transversely isotropic plates acted upon by concentrated heat sources is solved. The {1, 2}-order equations of thermoelasticity that incorporate the transverse shear and normal stresses are used. A bending heat source with symmetric heat transfer is considered. The dependence of thermal stress components on the thermal and thermomechanical parameters of transversely isotropic plates is studied     

A numerical simulation of combined radiation and natural convection heat transfer in a square enclosure heated by a centric circular cylinder

A numerical simulation of combined natural convection and radiation in a square enclosure heated by a centric circular cylinder and filled with absorbing-emitting medium is presented. The ideal gas law and the discrete ordinates method are used to model the density changes due to temperature differences and the radiation heat transfer correspondingly. The influence of Rayleigh number, optical thickness and temperature difference on flow and temperature fields along with the natural convection, radiation and total Nusselt number at the source surfaces is studied. The results reveal that the radiation heat transfer as well as the optical thickness of the fluid has a distinct effect on the fluid flow phenomena, especially at high Rayleigh number. The heat transfer and so the Nusselt number decreases with increase in optical thickness, while increases greatly with increase in temperature difference. The variation in radiation heat transfer with optical thickness and temperature difference is much more obvious as comparison with convection heat transfer.     

Numerical modelling of combined conductive and radiative heat transfer within square enclosure using discrete ordinate method

The current study addresses the mathematical modeling aspects of coupled conductive and radiative heat transfer in presence of absorbing, emitting and isotropic scattering gray medium within two-dimensional enclosure. The walls of the enclosure are considered to be opaque, diffuse and gray. The enclosure comprises isothermal vertical walls and insulated horizontal walls. The discrete ordinate method has been employed for modeling the radiative transport equation and the finite volume method has been adopted as the numerical technique. The effect of various parameters, i.e., radiation-conduction parameter, surface emissivity single scattering albedo and optical thickness has been illustrated.     

Thermal gradient effect on the free vibration of a cardioid cross-section cylindrical shell

For the first time, a new cross-section profile and efficient method are developed for the vibration analysis of isotropic and orthotropic cylindrical shells having circumferentially varying profile of a cardioid cross-section expressed as an arbitrary function, under thermal gradient effect. The governing equations of orthotropic cylindrical shells with varying thermal gradient around its circumference are derived as a boundary-value problem and solved numerically as an initial-value problem, based on the framework of Flügge's shell theory, transfer matrix approach and Romberg integration method. As a semi-analytical procedure, the trigonometric functions are used with Fourier's approach to approximate the solution in the longitudinal direction and also to reduce the two-dimensional problem to one-dimensional one. The thermal gradient is assumed to arise due to the variation of Young's moduli and shear modulus, along the circumferential direction of the shell. The results are obtained to indicate the effects of cardioid cross-section on the natural frequencies and corresponding mode shapes in the thermal environment as well as the sensitivity of the vibration behavior to the thermal gradient ratio and the orthotropy of the shell is also investigated for different types of vibration modes. In general, close agreement between the obtained results and those of other researchers has been found.     

温度梯度对热致型形状记忆聚合物板力学性能的影响

形状记忆聚合物作为一种新型的智能材料,由于质量轻、成本低以及变形回复率大等优势,已经在航空、医学等领域得到广泛应用。当前对于热致型形状记忆聚合物力学行为的研究,大都集中在整体变温的情况下,随着应用环境的越来越广泛,温度梯度对材料力学性能的影响效果越发重要。本文在均匀应力的假设下,给出材料在不同初始和传热条件下的温度梯度分布情况,结合传热学和热致型形状记忆聚合物相变理论本构模型,分别讨论了不同温度梯度对存储应变、弹性模量等力学性能的影响,通过理论分析和实验数据对比验证了模型正确性。本文研究可为不同导热情况下,对热致型形状记忆聚合物力学行为监测提供思路,也为形状记忆聚合物的进一步工程应用提供理论依据。     

Hopkinson拉杆试验的优化与高导无氧铜拉伸本构关系的确定

为有效测量试件中的应力、应变及应变率,Hopkinson拉伸试验(TSHB)必须作优化分析,所进行的数值模拟涉及试件与杆件等连接对于实验结果的影响.为减小上升前沿、惯性效应且使试件处于一维应力及均匀应力与应变状态,优化的试件具有一定的长度与形状要求.对于高导无氧铜,由准静态试验及不同应变率与温度的优化TSHB试验得到一...     

一种散体材料SHPB被动围压试验体应力修正方法

本文利用有限元仿真给出了一种修正方法，并用数值仿真和试验验证了该方法的可靠性。研究表明:散体材料SHPB被动围压试验中，试样厚度远小于厚壁圆筒长度时，端部效应会导致厚壁圆筒不均匀凸出变形，计算材料的体应力-应变关系不能将厚壁圆筒应力状态简化为平面应力问题；厚壁圆筒处于弹性状态下，通过厚壁圆筒理论计算出的径向力与真实径向力存在一定比例关系，在一定范围内，折算系数与试样实时厚度呈二次函数关系。     

Heat transfer by conduction,natural convection and radiation across a rectangular cellular structure

This paper describes results on the effects of wall conduction and radiation heat exchange among surfaces on laminar natural convection heat transfer in a two-dimensional rectangular cavity modelling a cellular structure. Parametric heat transfer calculations have been performed, and numerical results are presented in graphical and tabular form. Local and average Nusselt numbers along the cavity walls are reported for a range of parameters of physical interest. The findings suggest that the local or the average Nusselt number is one of many parameters that control conjugate heat transfer problems. The results indicate that natural convection heat transfer in the cavity is reduced by heat conduction in the walls and radiation exchange among surfaces. The results obtaibed for the total heat transfer rate through the system using the two-dimensional model are compared with those based on a one-dimensional model.     

Solution of the elastic boundary value problems for a layer with tunnel stress raisers

A novel procedure for solving three-dimensional problems for elastic layer weakened by through-thickness tunnel cracks has been developed and is presented in this paper. This procedure reduces the given boundary value problem to an infinite system of one-dimensional singular integral equations and is based on a system of homogeneous solutions for a layer. Integral representations of single- and double-layer potentials are used for metaharmonic and harmonic functions entering in the singular integral equations. These representations provide a continuous extendibility of the stress vector while allowing a jump in the displacement vector in the transition through the cut.Expanding the potential and biharmonic solutions in the Fourier series over the thickness coordinate yields the integral representations of the displacement vector and stress tensor. The problem of reducing a denumerable set of the integral equations of the given boundary value problem to one-to-one correspondence with the set of unknown densities appearing in the Fourier’s coefficient representations has been settled efficiently. Numerical investigations show a rapid convergence of the proposed reduction procedure as applied to the solution of the infinite system of one-dimensional integral equations. Numerical examples illustrate the proposed method and demonstrate its advantages.     

Effects of specimen thickness, hardening and crack closure for the plastic strip model

The triaxial stress constraint T_z and the effective yield stress distribution in the plastic zone for strain hardening materials are considered for the strip where the thickness effect is investigated by introducing a plastic constraint factor α. This factor depends on the specimen thickness, crack length, load level and hardening exponent. A simple expression of the plastic zone length and an expression involving α are obtained. Application of the strip model to crack closure shows that the specimen thickness has an influence and the results are compared with those found finite element.     

Thermal Analysis of Natural Convection Porous Fins

This work introduces a simple method of analysis to study the performance of porous fins in a natural convection environment. The method is based on using energy balance and Darcy’s model to formulate the heat transfer equation. The thermal performance of porous fins is then studied for three types of fins: long fin, finite-length fin with insulated tip and a finite-length fin with tip exposed to a known convection coefficient. It is found from the analysis that the effect of different design and operating parameters such as: Ra number, Da number, thermal conductivity ratio, Kr and length thickness ratio on the temperature distribution along the fin is grouped into one newly defined parameter called S_H. The effect of the variation of S_H on the porous fin thermal performance is established. The effect of varying the fin length and thermal conductivity ratio on the heat transfer rate from the fin is investigated and compared with that for a solid fin at certain conditions. It is found that the heat transfer rate from porous fin could exceed that of a solid fin. It is also found that increasing the fin length and effective thermal conductivity enhances the heat transfer from the fin up certain limit, where a further increase in these parameters adds no improvement to the fin performance. On Leave from Jordan University of Science and Technology, Irbid-Jordan     

Review and evaluation of the double-torsion technique for fracture toughness and fatigue testing of brittle materials

This paper reviews the double-torsion (DT) test as an experimental technique for the measurement of fracture toughness and slow-crack-growth behavior in brittle materials based on the authors' experiences and an evaluation of current literature. The DT technique has numerous advantages due primarily to the fact that the stress intensity is independent of crack length, at least for the central half of the specimen. Although the technique was first proposed about 20 years ago, and has been used extensively since then, there are a number of important unresolved questions concerning the methodology. To date there has been no standardization of test procedure or specimen geometry. A review of specimen geometries in use indicates that the proportions that are most commonly employed (based on the literature and experience) are width:W, length: 3W, thickness:W/6–W/15. Grooves on both the top and the bottom surfaces have been used to guide the crack, however it has been found that ungrooved, but very accurately aligned, specimens give the best results. Theoretical thickness-correction factors which account for both relatively thick specimens (with respect to width) and the effect of the size of the loading points, have been considered. The effect of crack-front profile on measured values of crack velocity and stress intensity is contentious. Althogh the stress intensity,K, varies along the crack front, the front merely translates axially. For the presentation of crackvelocity stress-intensity (V-K) data, the consensus seems to be that the only rational velocity to use is that based on the crack's intersection with the tensile surface. Despite some of the shortcomings mentioned above, the DT technique is widely accepted and gaining in popularity. It is particularly useful under cyclic fatigue conditions for investigating the effect of a change in a single parameter on crack-growth rate, using the very effective and elegant ‘changeover’ method.     

Stress intensity factor for an edge crack in two bonded dissimilar materials under convective cooling

The transient thermal stress crack problem for two bonded dissimilar materials subjected to a convective cooling on the surface containing an edge crack perpendicular to the interface is considered. The problem is solved using the principle of superposition and the uncoupled quasi-static thermoelasticity. The crack problem is formulated by applying the transient thermal stresses obtained from the uncracked medium with opposite sign on the crack surfaces to be the only external loads. Fourier integral transform is used to solve the perturbation problem resulting in a singular integral equation of Cauchy type in which the derivative of the crack surface displacement is the unknown function. The numerical results of the stress intensity factors are calculated for both the edge crack and the crack terminating at the interface using two different composite materials and illustrated as a function of time, crack length, coefficient of heat transfer, and the thickness ratio.     

有限厚度板穿透裂纹前缘附近三维弹性应力场分析

通过三维有限元计算来研究有限宽度、有限厚度含有穿透裂纹板的裂纹前缘应力场，从中找出应力强度因子与板的厚度、裂纹长度之间的关系，同时还分析了裂尖的三维约束程度和三维约束区的大小。分析结果表明：应力强度因子沿厚度的分布是不均匀的，应力强度因子的最大值及其位置与厚度有关；有限厚度板中面应力强度因子(KI)m-p及最大应力强度因子(KI)max均大于平面应力或平面应变的应力强度因子。对有限厚度裂纹问题，按平面应力或平面应变来考虑是不安全的；板中面的应力强度因子(KI)m-p及最大应力强度因子(KI)max是厚度B／a的函数；板的中面离面约束系数Tx最大，自由面(z=B)Tx=0。沿厚度方向裂尖附近的离面约束系数Tx也是z／B和B／a的函数，随着厚度的增加离面约束系数Tx增大，离中面越近离面约束系数Tx越大。Tx随着x的增大急剧减小，三维约束影响区域大小大约为板厚的一半，且裂纹长度a／W对应力强度因子沿厚度变化规律及Tx影响区域大小影响较小。     

Conjugate heat transfer within a concentric annulus filled with micropolar fluid

This paper investigates numerically the conjugate heat transfer in an annulus between two concentric cylinders. The annulus contains micropolar fluid and is heated isothermally from its inner wall. The effect of Rayleigh number, thickness of inner wall, inner wall-fluid thermal conductivity ratio, and material parameters of micropolar fluid on heat transfer rate within the annulus has been investigated. The study has shown that for low Rayleigh number regimes and for thermal conductivity of the inner wall greater than that of the fluid, the increase of inner wall thickness increases the heat transfer rate through the annulus and vice versa. While for convection dominating regimes Ra ≥ 10⁴ the increase of inner wall thickness decreases the heat transfer rate. Moreover, the study has shown that for fixed geometrical and flow parameters the heat transfer decreases in case of micropolar fluids in comparison with that of Newtonian fluids.     

Analysis of coupled conduction and radiation heat transfer in presence of participating medium- using a hybrid method

The current study addresses the mathematical modeling aspects of coupled conductive and radiative heat transfer in presence of absorbing, emitting and isotropic scattering gray medium within two-dimensional square enclosure. The walls of the enclosure are considered to be opaque, diffuse and gray. The enclosure comprised of isothermal vertical walls and insulated horizontal walls. A new hybrid method where the concepts of modified differential approximation employed by blending discrete ordinate method and spherical harmonics method, has been developed for modeling the radiative transport equation. The finite volume method has been adopted as the numerical technique. The effect of various influencing parameters i.e., radiation-conduction parameter, surface emissivity, single scattering albedo and optical thickness has been illustrated. The compatibility of the method with regard to solving coupled conduction and radiation has also been addressed.     

Necking evolution in dynamically stretched bars: New experimental and computational insights

This paper presents new results on dynamic neck evolution in steel bars of varying diameters. Dynamic tensile tests were carried out in a Kolsky apparatus using cylindrical steel specimens with various cross-section diameters ranging from 1.5 mm to 4 mm. A high speed digital camera was used to record the deformation of the specimen during the loading process. Video recording of the tests enabled accurate experimental measurements of the necking evolution, specifically its growth rate as a function of the diameter. The experiments show that increasing the specimen cross-section slows down the neck development. This behavior has been further investigated using two different kinds of numerical calculations: (1) axisymmetric finite element simulations and (2) one-dimensional finite difference computations. While the finite difference model only considers the normal stress along the longitudinal direction of the bar, the finite element model does not entail any simplification on the stress state of the specimen during the loading process. In agreement with the experiments, the finite element calculations show a decrease of the necking growth rate with the increase in the cross-section of the sample. On the contrary, the damping effect of the specimen cross-section on the necking evolution is not captured by the finite difference computations. We postulate that this result comes from the one-dimensional nature of the finite difference model. This work uncovers, by means of combined experiments and modelling, the key role played by stress multiaxiality in the growth rate of dynamic necks.     

Experiments and homogeneous turbulence model of boiling flowin narrow channels

The boiling heat transfer experiments have been carried out in vertical narrow annular channels with pure water. A two-dimensional homogeneous turbulence model of boiling flow has been developed and solved numerically to yield pressure gradient, and velocity, thermal and turbulence fields, together with local heat transfer coefficient along the length of the tube. Predictions are compared with the data of experiments and agreed well with it. The model results show that the heat transfer coefficient increases as the gap size decreases in annular channels. This model can be used to predict heat transfer of boiling flow in narrow channels.     

The reduction of spatial aliasing by long hot-wire anemometer probes

 Experiment and numerical analysis are presented to demonstrate that a hot-wire anemometer probe reduces spatial aliasing of turbulent velocity fluctuations because of the filtering property of the probe sensing element. The experiment focuses on the one-dimensional turbulent velocity spectrum and utilizes a long sensing length hot-wire probe to exaggerate the effect of the sensing element on the turbulent field. The numerical analysis utilizes a model of the hot-wire probe from Wyngaard (1968) along with isotropic turbulence relations to obtain an equation for the hot-wire response in a turbulent velocity field. The model can be used to determine the effect of hot-wire length on the one and three-dimensional turbulent spectra. The experimental study demonstrates that the finite length, hot-wire probe filters out energy in the high wave number region of the one-dimensional spectrum thereby verifying its ability to reduce spatial aliasing. Interestingly, the study also shows that energy in the low wave numbers of the one-dimensional spectrum is attenuated. The numerical study of the hot-wire probe demonstrates that this low wave-number attenuation is purely an artifact of the one-dimensional spectrum and not an effect of the hot-wire probe. Received: 20 May 996 / Accepted: 14 November 1996