饱和多孔弹性Timoshenko梁的大挠度分析

基于微观不可压饱和多孔介质理论和弹性梁的大挠度变形假设,考虑梁剪切变形效应,在梁轴线不可伸长和孔隙流体仅沿轴向扩散的限定下,建立了饱和多孔弹性Timoshenko梁大挠度弯曲变形的非线性数学模型.在此基础上,利用Galerkin截断法,研究了两端可渗透简支饱和多孔Timoshenko梁在突加均布横向载荷作用下的拟静态弯曲,给出了饱和多孔 Timoshenko梁弯曲变形时固相挠度、弯矩和孔隙流体压力等效力偶等随时间的响应.比较了饱和多孔Timoshenko梁非线性大挠度和线性小挠度理论以及饱和多孔 Euler-Bernoulli梁非线性大挠度理论的结果,揭示了他们间的差异,指出当无量纲载荷参数q＞l0时,应采用饱和多孔Timoshenko梁或Euler-Bernoulli梁的大挠度数学模型进行分析,特别的,当梁长细比λ＜30时,应采用饱和多孔Timoshenko梁大挠度数学模型进行分析.     

饱和多孔弹性Timoshenko悬臂梁的拟静力弯曲

在经典单相Timoshenko梁变形和孔隙流体仅沿多孔梁轴向运动的假定下,基于不可压饱和多孔介质的三维理论,本文首先建立了横观各向同性饱和多孔弹性Timoshenko悬臂梁拟静力弯曲的一维数学模型,并给出了相应的边界条件。其次,利用Laplace变换及其数值逆变换,分析了端部不同渗透条件下,饱和多孔弹性Timoshenko悬臂梁在端部梯载荷作用下的拟静力响应,给出了饱和多孔Timoshenko悬臂梁弯曲时挠度、弯矩以及孔隙流体压力等效力偶等随时间的响应曲线,并与饱和多孔Euler-Bernoulli悬臂梁的响应进行了比较,考察了梁长细比对弯曲的影响。数值结果表明：固相骨架与孔隙流体的相互作用具有粘性效应,梁弯曲的拟静态挠度具有蠕变行为,端部渗透条件对梁的弯曲响应有显著的影响,并且,饱和多孔弹性Timoshenko悬臂梁的拟静态响应亦存在Mandel-Cryer现象。     

不可压饱和多孔Timoshenko梁动力响应的数学模型

基于饱和多孔介质理论,假定孔隙流体仅沿梁的轴向运动,本文建立了横观各向同性饱和多孔弹性Timoshenko梁动力响应的一维数学模型,通过不同的简化,该模型可分别退化为饱和多孔梁的Euler-Bernoulli模型、Rayleigh模型和Shear模型等。研究了两端可渗透Timoshenko简支梁自由振动的固有频率、衰减率和阶梯载荷作用下的动力响应特征,给出了梁弯曲时挠度、弯矩以及孔隙流体压力等效力偶等随时间的响应曲线,并与饱和多孔Euler-Bernoulli简支梁响应进行了比较,考察了固相与流相相互作用系数、梁长细比等的影响。可见,固相骨架与孔隙流体的相互作用具有粘性效应,随着作用系数的增加,梁挠度振动幅值衰减加快,并最终趋于静态响应,Euler-Bernoulli梁的挠度幅值和振动周期小于Timoshenko梁的挠度幅值和周期,而Euler-Bernoulli梁的弯矩极限值等于Timoshenko梁的弯矩极限值。     

不可压饱和多孔弹性梁、杆动力响应的数学模型

基于多孔介质理论,首先建立了饱和多孔弹性杆件弯曲与轴向变形时动力响应的数学模型.其次,基于多孔弹性梁弯曲变形的数学模型,利用Laplace变换,分析了两端可渗透的饱和多孔弹性悬臂梁在自由端受阶梯载荷作用下的动静力响应,给出了梁弯曲时挠度、弯矩以及孔隙流体压力等效力偶等物理量随时间的响应曲线.发现不可压多孔弹性梁的拟静态响应亦存在Mandel-Cryer现象,多孔弹性梁的挠度具有与粘弹性梁挠度类似的蠕变特征,然而,其应力响应不同于粘弹性梁,随着时间的增加,梁拟静态响应的弯矩逐渐增加,并达到一个稳态值.这些结果有助于揭示植物根茎等力学行为的机理.     

轴向扩散不可压饱和多孔弹性梁的拟静态弯曲

建立了横观各向同性不可压饱和多孔弹性梁拟静态弯曲的数学模型,并给出了一般的求解方法。作为例子,研究了端部不同渗透条件对梁中点承受突加常集中载荷作用的饱和多孔悬臂梁拟静态弯曲的影响,给出了挠度和孔隙流体压力等效力偶沿梁轴线的分布以及随时间的响应曲线。结果表面:端部渗透条件对饱和多孔弹性梁的弯曲行为有显著的影响,梁的弯曲挠度既可随时间单调递增、亦可单调递减,其性态依赖于梁端部的渗透条件。同时发现不同于经典单相弹性梁,由于孔隙流体压力的作用,不承受载荷作用的梁段亦发生弯曲,并且Man-del-Cryer效应亦存在于不可压饱和多孔弹性梁的拟静态响应中,这些结果有助于揭示传热管道、植物根茎等力学行为的机理。     

A nonlinear mathematical model for large deflection of incompressible saturated poroelastic beams

Nonlinear governing equations are established for large deflection of incom- pressible fluid saturated poroelastic beams under constraint that diffusion of the pore fluid is only in the axial direction of the deformed beams.Then,the nonlinear bend- ing of a saturated poroelastic cantilever beam with fixed end impermeable and free end permeable,subjected to a suddenly applied constant concentrated transverse load at its free end,is examined with the Gaierkin truncation method.The curves of deflections and bending moments of the beam skeleton and the equivalent couples of the pore fluid pressure are shown in figures.The results of the large deflection and the small deflection theories of the cantilever poroelastic beam are compared,and the differences between them are revealed.It is shown that the results of the large deflection theory are less than those of the corresponding small deflection theory,and the times needed to approach its stationary states for the large deflection theory are much less than those of the small deflection theory.     

简支饱和多孔弹性梁的非线性弯曲

基于饱和多孔介质理论和弹性梁的大挠度弯曲假设,在多孔弹性梁轴线不可伸长,孔隙流体仅沿轴向方向扩散的限制下,建立了微观不可压饱和多孔弹性梁大挠度拟静态响应的一维非线性数学模型.在此基础上,利用Galerkin截断法,分析了两端可渗透的简支多孔弹性梁在突加横向均布载荷作用下的非线性弯曲,给出了梁弯曲时挠度、弯矩以及孔隙流体压力等效力偶随时间的响应曲线.数值结果表明:当载荷较小时,大挠度非线性与小挠度线性理论的结果相差很小,而当载荷较大时,非线性大挠度理论的结果小于相应线性小挠度理论的结果,并且这种差异随着载荷的增大而增大.同时,在载荷突加于梁上时,多孔弹性梁骨架起初不变形,孔隙流体压力等效力偶由零突增为非零,其值与外载荷保持平衡.随着时间的增加,梁的挠度增加,等效力偶逐渐减小为零,最终多孔梁骨架承担全部的外载荷.     

轴向扩散下简支饱和多孔弹性梁的大挠度分析

基于多孔介质理论和弹性梁的大挠度理论,并考虑轴向变形,在孔隙流体仅沿轴向扩散的假设下,建立了微观不可压饱和多孔弹性梁大挠度弯曲变形的一维非线性数学模型.在此基础上,忽略饱和多孔弹性梁的轴向应变,并利用Galerkin截断法,研究了两端可渗透的简支饱和多孔弹性梁在突加横向均布载荷作用下的拟静态弯曲,给出了饱和多孔梁弯曲时挠度、弯矩和轴力以及孔隙流体压力等效力偶等沿轴线的分布曲线.揭示了大挠度非线性和小挠度线性模型的结果差异,指出大挠度非线性模型的结果小于相应小挠度线性模型的结果,并且这种差异随着载荷的增大而增大.计算表明:当无量纲载荷参数q>5时,应该采用大挠度非线性数学模型进行研究.     

不可压饱和多孔弹性简支梁的动力响应

在杆件弯曲小变形的假定下,考虑杆件的侧向变形因素,根据多孔介质理论,本文首先建立了不可压饱和多孔弹性梁弯曲变形时动力响应的控制方程。其次,基于所建立的控制微分方程,利用变量分离法,研究了两端可渗透的饱和多孔弹性简支梁在梁中间集中载荷作用下的动力响应,得到了不同物性参数下简支梁动态弯曲时挠度和孔隙流体压力等效力偶等随时间的响应曲线。研究发现由于孔隙流体和固相骨架的相互作用,不可压饱和多孔弹性梁挠度的动力响应具有粘性特征,同时,随着时间的增加,饱和多孔弹性梁的挠度、弯矩等最终趋于经典弹性梁的静挠度、弯矩,此时,孔隙流体压力为零,梁的固相骨架承担所有的外载荷。     

不可压流体饱和多孔弹性梁的变分原理及有限元方法

基于不可压饱和多孔弹性梁动力弯曲的数学模型,建立了以多孔弹性梁挠度和孔隙流体压力等效力偶为宗量的Gurtin型变分原理,并给出了特殊边界条件下解耦时的仅以挠度为宗量的变分原理.同时,作为动力响应的退化情形,讨论了拟静态情形下的相应变分原理.根据所建立的变分原理,导出了一个有限元离散公式.由于Gurtin型变分原理是关于时间的卷积型的泛函,空间的有限元离散导致一个关于时间的对称微分一积分方程组,此方程组可进一步转化为常微分方程组.利用隐式Euler法,给出了时间区域的计算格式.作为一个数值例子,分析了饱和多孔弹性悬臂梁在自由端简谐载荷作用下的动力响应,分析了流相与固相相互作用对饱和多孔弹性悬臂梁动力响应的影响.     

Acoustical characterization of fluid-saturated porous media with local heterogeneities: Theory and application

A linear dynamic model of fully saturated porous media with local (either microscopic or mesoscopic) heterogeneities is developed within the context of Biot’s theory of poroelasticity. Viscoporoelastic behavior associated with local fluid flow is characterized by the notion of the dynamic compatibility condition on the interface between the solid and the fluid. Complex, frequency-dependent material parameters characterizing the viscoporoelasticity are derived. The complex properties can be obtained through determining the quasi-static poroelastic parameters, the properties of individual constituents, and the relaxation time of the dynamic compatibility condition on the interface. Relationships among various quasi-static poroelastic parameters are developed. It is shown that local fluid flow mechanism is significant only in the porous media with local heterogeneities. The relaxation time of the compatibility condition on the interface depends upon the details of local structure of porous media that control local fluid pressure diffusion. The new model is used to describe the velocity dispersion and attenuation in fully saturated porous media. The proposed model provides a theoretical framework to simulate the acoustical behavior of fully saturated porous media over a wide range of frequencies without making any explicit assumption about the structure of local heterogeneities.     

DYNAMIC STABILITY ANALYSIS OF EMBEDDED MULTI-WALLED CARBON NANOTUBES IN THERMAL ENVIRONMENT

In the present paper, the dynamic stability of multi-walled carbon nanotubes(MWCNTs) embedded in an elastic medium is investigated including thermal environment effects. To this end, a nonlocal Timoshenko beam model is developed which captures small scale effects.Dynamic governing equations of the carbon nanotubes are formulated based on the Timoshenko beam theory including the effects of axial compressive force. Then a parametric study is conducted to investigate the influences of static load factor, temperature change, nonlocal parameter, slenderness ratio and spring constant of the elastic medium on the dynamic stability characteristics of MWCNTs with simply-supported end supports.     

Thermo-Osmosis Effect in Saturated Porous Medium

In this paper, the thermo-poroelasticity theory is used to investigate the quasi-static response of temperatures, pore pressure, stress, displacement, and fluid flux around a cylindrical borehole subjected to impact thermal and mechanical loadings in an infinite saturated poroelastic medium. It has been reported in literatures that coupled flow known as thermo-osmosis by which flux is driven by temperature gradient, can significantly change the fluid flux in clay, argillaceous and many other porous materials whose permeability coefficients are very small. This study presents a mathematical model to investigate the coupled effect of thermo-osmosis in saturated porous medium. The energy balance equations presented here fulfill local thermal non-equilibrium condition (LTNE) which is different from the local thermal equilibrium transfer theory, accounting for that temperatures of solid and fluid phases are not the same and governed by different heat transfer equations. Analytical solutions of temperatures, pore pressure, stress, displacement, and fluid flux are obtained in Laplace transform space. Numerical results for a typical clay are used to investigate the effect of thermo-osmosis. The effects of LTNE on temperatures, pore pressure, and stress are also studied in this paper.     

Dynamic Model of a Periodic Medium with Double Porosity

The paper presents a unified mathematical approach for describing the dynamic stressstrain state of mechanical structures from heterogeneous materials possessing a double coupled system of pore channels filled with fluid. New dynamic equations describing the oscillations of poroelastic systems based on the developed model of a continuous medium with additional degrees of freedom in the form of various pressures of the components constituting the liquid phase of the material are obtained. The equations and the method of obtaining them have a greater degree of generalization than those encountered in the literature. Theoretical results can be used to study the propagation of vibrations in fractured geological rocks saturated with liquid, to develop technical systems of new structural materials with a porous structure, for the analysis of micro streams of fluid in the hierarchical system of microporous bone tissue.     

Homogenized modeling for vascularized poroelastic materials

A new mathematical model for the macroscopic behavior of a material composed of a poroelastic solid embedding a Newtonian fluid network phase (also referred to as vascularized poroelastic material), with fluid transport between them, is derived via asymptotic homogenization. The typical distance between the vessels/channels (microscale) is much smaller than the average size of a whole domain (macroscale). The homogeneous and isotropic Biot’s equation (in the quasi-static case and in absence of volume forces) for the poroelastic phase and the Stokes’ problem for the fluid network are coupled through a fluid-structure interaction problem which accounts for fluid transport between the two phases; the latter is driven by the pressure difference between the two compartments. The averaging process results in a new system of partial differential equations that formally reads as a double poroelastic, globally mass conserving, model, together with a new constitutive relationship for the whole material which encodes the role of both pore and fluid network pressures. The mathematical model describes the mutual interplay among fluid filling the pores, flow in the network, transport between compartments, and linear elastic deformation of the (potentially compressible) elastic matrix comprising the poroelastic phase. Assuming periodicity at the microscale level, the model is computationally feasible, as it holds on the macroscale only (where the microstructure is smoothed out), and encodes geometrical information on the microvessels in its coefficients, which are to be computed solving classical periodic cell problems. Recently developed double porosity models are recovered when deformations of the elastic matrix are neglected. The new model is relevant to a wide range of applications, such as fluid in porous, fractured rocks, blood transport in vascularized, deformable tumors, and interactions across different hierarchical levels of porosity in the bone.     

Large thermal deflections of Timoshenko beams under transversely non-uniform temperature rise

Geometrically non-linear deformation of axially extensional Timoshenko beams subjected mechanical as well thermal loadings were characterized by a system of 7 coupled and highly non-linear ordinary differential equations, which results in a complicated two-point boundary-value problem. By using shooting method this kind of problem can be numerically solved efficiently. Based on the above-mentioned mathematical formulation and numerical procedure, analysis of large thermal deflections for Timoshenko beams, subjected transversely non-uniform temperature rise and with immovably pinned–pinned as well as fixed–fixed ends, is presented. Characteristic curves showing the relationships between the beam deformation and temperature rise are illustrated. Especially, the effects of shear deformation on the bending and buckling response are quantitatively investigated. The numerical results show, as we know, that shear deformation effects become significant with the decrease of the slenderness and with the increase of the shear flexibility.     

Quasi-static response of linear viscoelastic cantilever beams subject to a concentrated harmonic end load

This study evaluates the response of a uniform cantilever beam with a symmetric cross-section fixed at one end, and submitted to a lateral concentrated sinusoidal load at the free extremity. The beam material is assumed to be homogeneous, isotropic and linear viscoelastic. Due to the nature of the loading and the beam slenderness, large displacements are developed but the strains are considered small. Consequently, the mathematical formulation only involves geometrical non-linearity. It is also assumed that the beam is inextensible (neutral axis length is constant) and that inertial forces are negligible, i.e., dynamic effects are insignificant and the system can thus be modeled quasi-statically. The beam is therefore subject to oscillations caused by the sinusoidal time-dependent load, leading to a transient response until the material stabilizes and the system exhibits a periodic response, which can be conveniently described in the frequency domain. The time domain solution of this problem is elaborated by considering the quasi-static response for each time interval. The mathematical equations are presented in dimensional and dimensionless forms, and for the latter case, a numerical solution is generated and several case studies are presented. The problem is governed by a set of non-linear ordinary differential equations encompassing functions of space and time that relate the curvature, rotation angle, bending moment and geometrical coordinates. In this study, an elegant solution is deduced using perturbation theory, yielding a precise steady-state solution in the frequency domain with considerable computational economy. The solutions for both time and frequency domain methods are developed and compared using a case study for a series of dimensionless parameters that influence the response of the system.     

饱和土中端承桩纵向振动特性研究

基于饱和多孔介质理论研究了三维轴对称条件下端承桩在饱和土中的纵向耦合振动. 首先通 过引入势函数对Biot动力固结方程解耦，采用算子分解方法及分离变量法求得饱和土层振动 解，进而利用桩土完全耦合条件得到桩土系统定解. 然后对桩顶的频率和时域响应进行参数 研究，结果表明，桩的长径比和桩土模量比对桩顶动力响应有较大影响，而渗透力对其影响 较小. 最后将幅频及时域反射波理论拟合结果与桩基实测结果加以对比，结果表明理论曲线 与实测曲线规律一致.     

Transient wave analysis of a cantilever Timoshenko beam subjected to impact loading by Laplace transform and normal mode methods

This study applies two analytical approaches, Laplace transform and normal mode methods, to investigate the dynamic transient response of a cantilever Timoshenko beam subjected to impact forces. Explicit solutions for the normal mode method and the Laplace transform method are presented. The Durbin method is used to perform the Laplace inverse transformation, and numerical results based on these two approaches are compared. The comparison indicates that the normal mode method is more efficient than the Laplace transform method in the transient response analysis of a cantilever Timoshenko beam, whereas the Laplace transform method is more appropriate than the normal mode method when analyzing the complicated multi-span Timoshenko beam. Furthermore, a three-dimensional finite element cantilever beam model is implemented. The results are compared with the transient responses for displacement, normal stress, shear stress, and the resonant frequencies of a Timoshenko beam and Bernoulli–Euler beam theories. The transient displacement response for a cantilever beam can be appropriately evaluated using the Timoshenko beam theory if the slender ratio is greater than 10 or using the Bernoulli–Euler beam theory if the slender ratio is greater than 100. Moreover, the resonant frequency of a cantilever beam can be accurately determined by the Timoshenko beam theory if the slender ratio is greater than 100 or by the Bernoulli–Euler beam theory if the slender ratio is greater than 400.     

Upscaling of the Coupling of Hydromechanical and Thermal Processes in a Quasi-static Poroelastic Medium

We undertake a formal derivation of a linear poro-thermo-elastic system within the framework of quasi-static deformation. This work is based upon the well-known derivation of the quasi-static poroelastic equations (also known as the Biot consolidation model) by homogenization of the fluid-structure interaction at the microscale. We now include energy, which is coupled to the fluid-structure model by using linear thermoelasticity, with the full system transformed to a Lagrangian coordinate system. The resulting upscaled system is similar to the linear poroelastic equations, but with an added conservation of energy equation, fully coupled to the momentum and mass conservation equations. In the end, we obtain a system of equations on the macroscale accounting for the effects of mechanical deformation, heat transfer, and fluid flow within a fully saturated porous material, wherein the coefficients can be explicitly defined in terms of the microstructure of the material. For the heat transfer we consider two different scaling regimes, one where the Péclet number is small, and another where it is unity. We also establish the symmetry and positivity for the homogenized coefficients.