Thermal buckling and post-buckling response of imperfect temperature-dependent sandwich FGM plates resting on elastic foundation

Post-buckling behaviour of sandwich plates with functionally graded material (FGM) face sheets under uniform temperature rise loading is considered. It is assumed that the plate is in contact with a Pasternak-type elastic foundation during deformation, which acts in both compression and tension. The derivation of equations is based on the first-order shear deformation plate theory. Thermomechanical non-homogeneous properties of FGM layers vary smoothly by the distribution of power law across the thickness, and temperature dependency of material constituents is taken into account. Using the non-linear von-Karman strain-displacement relations, the equilibrium and compatibility equations of imperfect sandwich plates with FGM face sheets are derived. The boundary conditions for the plate are assumed to be simply supported in all edges. The governing equations are reduced to two coupled equation in terms of stress function and lateral deflection. Employing the single mode approach combined with Galerkin technique, an approximate closed-form solution is presented to calculate the critical buckling temperature and post-buckling equilibrium path of the plate. Presented numerical examples contain the influences of power law index, sandwich plate geometry, geometrical imperfection, temperature dependency, and the elastic foundation coefficients.     

On the layered structure in a stable solute gradient heated from below

When a linear stable solute gradient, contained in a tall tank, is heated uniformly from below, double diffusive layers may be formed successively from the bottom to the top of the fluid. A simple model is developed which considers the evolution of the bottom mixed layer, and a set of experiments was performed when the temperature at the solid boundary is specified. The model and the experiments show that the instantaneous thickness of the bottom mixed layer is directly proportional to the difference between the instantaneous temperature of the solid boundary and that of the far fluid, and inversely proportional to the ambient solute gradient. The experiments reveal that during its formation the second layer persists above the bottom one notwithstanding the advancing interface below. Similar behavior is observed during the formation of subsequent secondary layers.     

粘贴压电层功能梯度材料Timoshenko梁的热过屈曲分析

研究了上下表面粘贴压电层的功能梯度材料Timoshenko梁在升温及电场作用下的过屈曲行为。在精确考虑轴线伸长和一阶横向剪切变形的基础上,建立了压电功能梯度Timoshenko层合梁在热-电-机械载荷作用下的几何非线性控制方程。其中,假设功能梯度的材料性质沿厚度方向按照幂函数连续变化,压电层为各向同性均匀材料。采用打靶法数值求解所得强非线性边值问题,获得了在均匀电场和横向非均匀升温场内两端固定Timoshenko梁的静态非线性屈曲和过屈曲数值解。并给出了梁的变形随热、电载荷及材料梯度参数变化的特性曲线。结果表明,通过施加电压在压电层产生拉应力可以有效地提高梁的热屈曲临界载荷,延缓热过屈曲发生。由于材料在横向的非均匀性,即使在均匀升温和均匀电场作用下,也会产生拉-弯耦合效应。但是对于两端固定的压电-功能梯度材料梁,在横向非均匀升温下过屈曲变形仍然是分叉形的。     

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

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

INVESTIGATIONS ON THE MECHANICAL CHARACTERISTICS OF NANOBEAMS BASED ON THE NONLOCAL STRAIN GRADIENT THEORY1)

基于非局部应变梯度理论,建立了一种具有尺度效应的高阶剪切变形纳米梁的力学模型. 其中,考虑了应变场和一阶应变梯度场下的非局部效应. 采用哈密顿原理推导了纳米梁的控制方程和边界条件,并给出了简支边界条件下静弯曲、自由振动和线性屈曲问题的纳维级数解. 数值结果表明,非局部效应对梁的刚度产生软化作用,应变梯度效应对纳米梁的刚度产生硬化作用,梁的刚度整体呈现软化还是硬化效应依赖于非局部参数与材料特征尺度的比值. 梁的厚度与材料特征尺度越接近,非局部应变梯度理论与经典弹性理论所预测结果之间的差异越显著.     

Correlation between interfacial interactions and internal stresses in filled high density polyethylene

The stress relaxation and the creep behaviour of high density polyethylene (HDPE) filled with glass fibres, clay (plate-like particles) or CaCO₃ (particles with irregular shape) were measured in uniaxial extension at room temperature. It was observed that the addition of filler increased the internal stress level, as evaluated from stress relaxation data. This increase was larger than the corresponding increase in the (short-term) elastic modulus. This behaviour may be attributed to a reduced macromolecular mobility in the matrix material close to the filler surface, i.e. to formation of an interphase region in the HDPE-matrix. From the internal stress values, the thickness of this interphase region around each filler particle was estimated, assuming a uniform coverage of the particles. It was suggested that the amount of matrix material with reduced mobility (or the thickness of the interphase region) reflected the degree of adhesion between the filler and HDPE. The change in the internal stress level due to the incorporation of different fillers, which were surface treated in some cases, was also consistent with the observed creep behaviour.     

Nanoindentation of hyperelastic polymer layers at finite deformation and parameter re-identification

Thin polymer layers on substrates have a wide range of application in important areas. However, it is impossible to measure the mechanical properties with the traditional testing methods. Recently, nanoindentation became a new but primary testing technique of thin layers. In the present work, based on a finite element model of contact mechanics and hyperelastic materials, nanoindentation of polymer layers is simulated with the finite element code ABAQUS^?. Three often used hyperelastic models, that is, the neo-Hookean, Mooney–Rivlin and Yeoh models are investigated. The behaviour of these three models is compared to each other in different boundary value problems of nanoindentation in order to get some feeling of the different behaviour of various hyperelastic models under nanoindentation. In contrast to the traditional analytical method, the penetration depth is not restrained to avoid the influence of the substrate. A parameter re-identification strategy is employed to extract the parameters of the material models at small and finite deformation based on the principle of biological evolution. Furthermore, it is investigated how large the penetration depth has to be chosen in order to distinguish different models in reference to the load–displacement curves. Finally, the possibility is discussed of describing the data obtained by a non-linear complex model using the relatively simple approach based on the neo-Hookean model.     

DNS,LES and RANS of turbulent heat transfer in boundary layer with suddenly changing wall thermal conditions

The objectives of this study are to investigate a thermal field in a turbulent boundary layer with suddenly changing wall thermal conditions by means of direct numerical simulation (DNS), and to evaluate predictions of a turbulence model in such a thermal field, in which DNS of spatially developing boundary layers with heat transfer can be conducted using the generation of turbulent inflow data as a method. In this study, two types of wall thermal condition are investigated using DNS and predicted by large eddy simulation (LES) and Reynolds-averaged Navier–Stokes equation simulation (RANS). In the first case, the velocity boundary layer only develops in the entrance of simulation, and the flat plate is heated from the halfway point, i.e., the adiabatic wall condition is adopted in the entrance, and the entrance region of thermal field in turbulence is simulated. Then, the thermal boundary layer develops along a constant temperature wall followed by adiabatic wall. In the second case, velocity and thermal boundary layers simultaneously develop, and the wall thermal condition is changed from a constant temperature to an adiabatic wall in the downstream region. DNS results clearly show the statistics and structure of turbulent heat transfer in a constant temperature wall followed by an adiabatic wall. In the first case, the entrance region of thermal field in turbulence can be also observed. Thus, both the development and the entrance regions in thermal fields can be explored, and the effects upstream of the thermal field on the adiabatic region are investigated. On the other hand, evaluations of predictions by LES and RANS are conducted using DNS results. The predictions of both LES and RANS almost agree with the DNS results in both cases, but the predicted temperature variances near the wall by RANS give different results as compared with DNS. This is because the dissipation rate of temperature variance is difficult to predict by the present RANS, which is found by the evaluation using DNS results.     

功能梯度矩形厚板的三维热弹性分析

直接从三维热弹性力学基本方程出发，通过引入两个位移函数和两个应力函数，导出了一个二阶的齐次状态方程和一个四阶的非齐次状态方程。分析中采用了层合近似模型，即将板划分成厚度足够小的若干薄层，从面可将每一层内的材料常数近似为常数。给出了任意厚度的四边简支横观各向同性功能梯度矩形板的热弹性分析，特别当板较薄时，与薄板理论进行了数值对比，发现两者结果吻合很好。最后研究了材料梯度指标对热弹性场的影响，结果显示梯度指标对热应力和位移都有着显著的作用：在不同的区间，梯度指标对它们有不同的影响；并且在同一区间，梯度指标对两者的影响程度也有所不同。     

功能梯度材料微梁的热弹性阻尼研究

基于Euler-Bernoulli梁理论和单向耦合的热传导理论,研究了功能梯度材料(functionally graded material,FGM)微梁的热弹性阻尼(thermoelastic damping,TED).假设矩形截面微梁的材料性质沿厚度方向按幂函数连续变化,忽略了温度梯度在轴向的变化,建立了单向耦合的变系数一维热传导方程.热力耦合的横向自由振动微分方程由经典梁理论获得.采用分层均匀化方法将变系数的热传导方程简化为一系列在各分层内定义的常系数微分方程,利用上下表面的绝热边界条件和界面处的连续性条件获得了微梁温度场的分层解析解.将温度场代入微梁的运动方程,获得了包含热弹性阻尼的复频率,进而求得了代表热弹性阻尼的逆品质因子.在给定金属-陶瓷功能梯度材料后,通过数值计算结果定量分析了材料梯度指数、频率阶数、几何尺寸以及边界条件对TED的影响.结果表明:(1)若梁长固定不变,梁厚度小于某个数值时,改变陶瓷材料体积分数可以使得TED取得最小值;(2)固有频率阶数对TED的最大值没有影响,但是频率阶数越高对应的临界厚度越小;(3)不同的边界条件对应的TED的最大值相同,但是随着支座约束刚度增大对应的临界厚度减小;(4)TED的最大值和对应的临界厚度随着金属组分的增大而增大.     

Free vibration of functionally graded material beams with surface-bonded piezoelectric layers in thermal environment

Free vibration of statically thermal postbuckled functionally graded material （FGM） beams with surface-bonded piezoelectric layers subject to both temperature rise and voltage is studied. By accurately considering the axial extension and based on the Euler-Bernoulli beam theory, geometrically nonlinear dynamic governing equations for FGM beams with surface-bonded piezoelectric layers subject to thermo-electro- mechanical loadings are formulated. It is assumed that the material properties of the middle FGM layer vary continuously as a power law function of the thickness coordinate, and the piezoelectric layers are isotropic and homogenous. By assuming that the amplitude of the beam vibration is small and its response is harmonic, the above mentioned non-linear partial differential equations are reduced to two sets of coupled ordinary differential equations. One is for the postbuckling, and the other is for the linear vibration of the beam superimposed upon the postbuckled configuration. Using a shooting method to solve the two sets of ordinary differential equations with fixed-fixed boundary conditions numerically, the response of postbuckling and free vibration in the vicinity of the postbuckled configuration of the beam with fixed-fixed ends and subject to transversely nonuniform heating and uniform electric field is obtained. Thermo-electric postbuckling equilibrium paths and characteristic curves of the first three natural frequencies versus the temperature, the electricity, and the material gradient parameters are plotted. It is found that the three lowest frequencies of the prebuckled beam decrease with the increase of the temperature, but those of a buckled beam increase monotonically with the temperature rise. The results also show that the tensional force produced in the piezoelectric layers by the voltage can efficiently increase the critical buckling temperature and the natural frequency.     

轴对称自由振动功能梯度材料微圆板中的热弹性阻尼

基于经典弹性薄板理论和单向耦合热传导理论,研究了材料性质沿厚度连续变化的功能梯度微圆板的热弹性阻尼特性．首先,考虑热力耦合效应,建立了功能梯度微圆板轴对称横向自由振动微分方程．然后,忽略温度梯度在面内的变化,建立了单向耦合变系数一维热传导方程．采用分层均匀化近似方法,将变系数热传导方程转化为一系列常系数的微分方程,利用上下表面的热边界条件和层间连续性条件获得了微圆板温度场解析解．将所得温度场代入微圆板的自由振动微分方程,得到了包含热弹性阻尼的复频率,从而获得了反映热弹性阻尼水平的逆品质因子．最后,针对材料性质沿板厚按幂函数变化的陶瓷-金属功能梯度微圆板,定量地分析材料梯度指数、几何尺寸、边界条件、温度环境等对微圆板热弹性阻尼的影响．     

Plasticity of formable all-metal sandwich sheets: Virtual experiments and constitutive modeling

The mechanical behavior of a metallic sandwich sheet material composed of two flat face sheets and two bi-directionally corrugated core layers is analyzed in detail. The manufacturing of the sandwich material is simulated to obtain a detailed unit cell model which accounts for the non-uniform thickness distribution and residual stresses associated with the stamping of the core layers. Virtual experiments are performed by subjecting the unit cell model to various combinations of bi-axial in-plane loading including the special cases of uniaxial tension, uniaxial compression, equi-biaxial tension and shear. The results demonstrate that the core structure’s contribution to the in-plane load carrying capacity of the sandwich sheet material is similar to that of the face sheets. The numerical results are also used to identify the effective yield surface and hardening response of both the core layer and the face sheets. An anisotropic yield function with linear pressure dependency is proposed to approximate the equal-plastic work surfaces for the core structure and face sheets. Furthermore, a new two-surface model with non-linear interpolation based on plastic work density is presented to describe the observed combined isotropic-distortional hardening of the core structure.     

Solutions for axial and transverse boundary layers in the case of steady laminar flow past a horizontal plate embedded in a saturated porous medium

The present paper is concerned with the interesting situation of the existence of transverse, velocity and thermal, boundary layers in the case of flow past a horizontal plate embedded in a saturated porous medium. The analysis shows that these layers are thicker than the axial boundary layers. However, the increase in the porous parameterB reduces the thickness. The expressions for the boundary layer thickness, the shear-stress at the wall and the specific heat flux are obtained and the behaviour of the above quantities whenB takes different values are discussed in detail.     

On Maxwell fluids with relaxation time and viscosity depending on the pressure

We study a variant of the well-known Maxwell model for viscoelastic fluids, namely we consider a Maxwell fluid with viscosity and relaxation time depending on the pressure. Such a model is relevant for example in modelling the behaviour of some polymers and geomaterials. Although it is experimentally known that the material moduli of some viscoelastic fluids can depend on the pressure, most of the studies concerning the motion of viscoelastic fluids do not take such effects into account despite their possible practical significance in technological applications. Using a generalized Maxwell model with pressure dependent material moduli we solve a simple boundary value problem and we demonstrate interesting non-classical features exhibited by the model.     

Experiments on localized disturbances in a flat plate boundary layer. Part 1. The receptivity and evolution of a localized free stream disturbance

The receptivity of a laminar boundary layer to free stream disturbances has been experimentally investigated through the introduction of deterministic localized disturbances upstream of a flat plate mounted in a wind tunnel. Hot-wire measurements indicate that the spanwise gradient of the normal velocity component (and hence the streamwise vorticity) plays an essential role in the transfer of disturbance energy into the boundary layer. Inside the laminar boundary layer the disturbances were found to give rise to the formation of longitudinal structures of alternating high and low streamwise velocity. Similar streaky structures exist in laminar boundary layers exposed to free stream turbulence, in which the disturbance amplitude increases in linear proportion to the displacement thickness. In the present study the perturbation amplitude of the streaks was always decaying for the initial amplitudes used, in contrast to the growing fluctuations that are observed in the presence of free stream turbulence. This points out the importance of the continuous influence from the free stream turbulence along the boundary layer edge.     

考虑尺度效应的微梁静态弯曲数值分析

基于修正偶应力理论及表面弹性理论,本文提出了一种新的双曲线剪切变形梁模型,用于均匀微尺度梁的静态弯曲分析。该理论可以直接利用本构关系获得横向剪切应力,满足梁顶部和底部的无应力边界条件,避免了引入剪切修正因子。根据广义Young-Laplace方程建立了梁的内部与表面层的应力连续性条件,单一的变量场可以描述梁的位移模式。通过在位移场中考虑表面层厚度以及表面层的应力连续条件,可以使新模型能够更准确地预测微尺寸和表面能相关的尺度效应。通过Hamilton原理推导出了梁的控制方程和边界条件。应变能除了考虑经典弹性理论,还要考虑微结构效应和表面能。Navier-type的解析解适用于简支边界条件,而基于拉格朗日插值的微分求积法(DQEM)可以研究在不同边界条件下的力学响应。把该数值解与Navier方法得出的解析解作了对比,得出:微尺度梁在考虑表面能或微尺寸效应、不同载荷和梁高变化下的响应一致;当不考虑微结构相关性和表面能效应时,该模型退化为经典的欧拉梁模型。     

热冲击下物性与温度相关性对热弹性响应的渐进分析

计及材料物性与温度的相关性,基于Green-Naghdi能量无耗散广义热弹性理论(G-N II理论),对热冲击下具有变物性特征材料的热弹性响应进行了求解分析。借助Laplace正、反变换技术以及Krichhoff变换,在热物性参数随真实温度呈线性规律的前提下,推导了半无限大体受热冲击作用时热弹性响应的解析表达式,通过求解分析,得到了热冲击下热波、热弹性波的传播规律,位移场、温度场以及应力场的分布情况,以及物性随温度相关性对热弹性响应的影响效果。结果表明：当考虑材料物性随温度的变化时,热波、热弹性波的传播以及各物理场的分布均受到不同程度的影响,且物性随温度相关性对热弹性响应的作用效果将受到材料热-力耦合特性的影响。     

A mathematical model for the waste region of a long-wall mineworking

The mechanical behaviour of the broken rock material occupying the waste region of a long-wall mineworking is investigated. A mathematical model is proposed, developed from physical considerations, which is based upon the continuum, plane strain, rigid-plastic theory for double-shearing flows of compressible granular materials. The model comprises the stress equilibrium equations, the Coulomb yield criterion, kinematic equations due to Mehrabadi and Cowin (1978) and the continuity equation, together with the waste region geometry and the boundary conditions imposed upon the field variables. Statical indeterminancy is a property of the formulation and this results in a connection between the stress and velocity fields through the boundary conditions. A strategy is presented for the iterative construction of the stress and velocity fields and approximations to solutions of the equations of the model are obtained numerically.     

The dissipative mechanism in flowing polymers: theory and experiments

The application of modern non-linear continuum thermodynamics to the analysis of energy dissipation in both isothermal and non-isothermal flow of polymers is discussed. A fundamental simplification of the field equation for the balance of energy arises if the assumption is made that the flowing polymer is a material with entropic elasticity. This assumption is partly justified by structural arguments, and its validity can be checked experimentally. The results of the experiments performed support the validity of the assumption.