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

在经典单相Timoshenko梁变形和孔隙流体仅沿饱和多孔弹性梁轴向运动的假定下,基于不可压饱和多孔介质的三维Gurtin型变分原理,首先建立了饱和多孔弹性Timoshenko悬臂梁动力响应的一维数学模型．在若干特殊情形下,该模型可分别退化为饱和多孔弹性梁的Euler-Bernoulli模型、Rayleigh模型和Shear模型等．其次,利用Laplace变换,分析了固定端不可渗透、自由端可渗透的饱和多孔弹性Timoshenko悬臂梁在自由端阶梯载荷作用下的动静力响应,给出了梁自由端处挠度随时间的响应曲线,考察了固相与流相相互作用系数、梁长细比等参数对悬臂梁动静力行为的影响．结果表明:饱和多孔弹性梁的拟静态挠度具有与粘弹性梁挠度类似的蠕变特征．在动力响应中,随着梁长细比的增大,自由端挠度的振动周期和幅值增大,且趋于稳态值的时间增长,而随着两相相互作用系数的增大,梁挠度振动衰减加快,并最终趋于经典单相弹性Timoshenko梁的静态挠度．     

不可压饱和多孔弹性杆动力响应的多辛方法

研究了不可压饱和多孔弹性杆的一维动力响应问题.基于多孔介质理论,在流相和固相微观不可压、固相骨架小变形的假定下,建立了不可压流体饱和多孔弹性杆一维轴向动力响应的数学模型.利用Hamilton空间体系的多辛理论,构造了不可压饱和多孔弹性杆轴向振动方程的多辛形式及其多种局部守恒律.采用中点Box离散方法得到轴向振动方程的多辛离散格式和局部能量守恒律以及局部动量守恒律的离散格式;数值模拟了不可压饱和多孔弹性杆的轴向振动过程,记录了每一时间步的局部能量数值误差和局部动量数值误差.结果表明,已构造的多辛离散格式具有很高的精确性和较长时间的数值稳定性,这为解决饱和多孔介质的动力响应问题提供了新的途径.     

饱和多孔弹性杆流固耦合动力响应的保结构算法

研究了不可压饱和多孔弹性杆的流固耦合动力响应问题.基于多孔介质理论,根据多孔介质流固混合物动量方程、孔隙流体动量方程及体积分数方程,建立流固耦合不可压饱和多孔弹性杆的轴向振动方程;引入正则变量,构造饱和多孔弹性杆轴向振动方程的广义多辛保结构形式、广义多辛守恒律及广义多辛局部动量误差;采用中点Box离散方法得到轴向振动方程的广义多辛离散格式、广义多辛守恒律数值误差及局部动量数值误差;数值模拟不可压饱和多孔弹性杆的轴向振动过程及流相渗流速度分布,考察了流固两相耦合系数对轴向振动过程及广义多辛守恒律误差和局部动量误差的影响.结果表明,已构造的广义多辛保结构算法具有很高的精确性和长时间的数值稳定性.     

多孔固体充满黏性流体时的边界条件

基于物理学基本原理和能量守恒定律的精确检查,导出充满黏性流体多孔固体边界呈连续性要求的边界条件.当孔隙流体具有黏性时,多孔弹性固体就是一个耗散的充满黏性流体的多空固体.孔隙流体的黏性造成的耗散应力准确地表达了边界条件.边界上两种固体连接的不完全,导致孔隙流体的流出,多孔骨料两边微粒运动的不平衡.导出多孔.多孔固体界面孔隙局部连接时的数学模型.在该界面上,滑移的松.紧,以及孔隙开.合,能造成一部分应变能的耗散.数值结果表明,在水和饱和油砂岩之间的界面上,修正的边界条件将影响各向同性多孔介质中折射波的能量.     

饱和多孔介质粘弹性理论

本文讨论了饱和多孔介质的粘弹体本构关系。把Lee的弹性和粘弹性比拟理论,拓展到饱和多孔介质领域内。用比拟理论解出了饱和粘性土的常载和变载的一维固结问题,并利用比拟理论的属性,可以由简单的粘弹性本构关系来推求复杂的本构关系,得到了粘弹性地基上的粘弹性梁板问题的解答。由此初步形成了饱和多孔介质粘弹性比拟法的完整理论。     

饱和粘土非线性渗流规律与径向固结

以多孔介质等效渗流概念得出粘土平均孔隙微尺度数量级范围在0.01 μm至0.1μm,与测试结果一致．实验结果表明饱和粘土微尺度孔隙渗流为非线性流．理论推导表明固液界面作用与渗透率平方根或孔隙半径成反比,固液界面相互作用是导致饱和粘土非线性渗流的重要原因．提出了精确描述饱和粘土微尺度孔隙非线性渗流基本规律的数学模型,其参数量纲明确,物理意义清楚．建立了饱和粘土非线性渗流新定律,统一描述了从较低到较高水力梯度全过程渗流曲线特征,达西定律是其特例．基于新定律,建立了饱和粘土非线性渗流定流量径向固结数学模型．以粘性边界层思想与稳态依序替换法及积分方程法,导出了粘土非线性渗流平均质量守恒方程及活动边界运动方程,给出了饱和粘土非线性渗流超孔隙水压力分布公式与平均固结度计算公式,获得了粘土层压力分布规律和平均固结度随时间变化规律．结果表明:饱和粘土非线性渗流使活动边界运动速度减小．研究结果为粘土地质工程与岩土工程应用提供了新的科学依据．达西渗流径向固结计算是新的非线性渗流固结计算的特例．     

多孔介质平板通道发展传热中非局部热平衡时的温度分布特征

研究了多孔介质平板通道中,Darcy流体发展传热强迫对流非局部热平衡下,固相骨架和孔隙流体的温度分布特征．考虑流体流动方向的热传导以及固相和流相相互作用的粘性耗散,根据非局部热平衡的两能量方程模型,得到了常壁温度时多孔介质固相骨架温度和孔隙流体温度的解析解．证明了当两相间的热交换系数趋于无穷大时,两能量方程的温度解趋于局部热平衡时一能量方程的温度解．针对不同的无量纲参数,给出了固相和流相的温度分布状态,通过参数研究,揭示了非局部热平衡强迫对流时温度对无量纲参数的依赖关系．     

轴向复合载荷作用下变截面悬臂弹性立柱后屈曲分析

基于轴线可伸长弹性杆的几何非线性理论,建立了同时作用端部轴向集中荷载和沿轴线作用分布轴向载荷的变截面弹性悬臂柱的后屈曲控制方程。采用打靶法直接求解了所得强非线性边值问题,给出了截面线性变化的圆截面柱的二次平衡路径及其过屈曲位形曲线。     

流体饱和多孔隙介质波动方程小波有限差分法

研究流体饱和多孔隙介质中波动方程的数值模拟．针对求解二维弹性波方程问题,提出小波有限差分法．该方法综合了小波多分辨分析计算灵活、计算效率高特性和有限差分易于实现的优点．数值模拟的结果显示,此方法对于求解流体饱和多孔隙介质方程的数值模拟是有效稳定的．     

横观各向同性含液饱和多孔介质中应力波传播的特征分析

根据广义特征理论,对横观各向同性含液饱和多孔介质中应力波传播特性进行了特征分析．给出了特征曲面的微分方程以及沿次特征线的相容条件,得到了波阵面的解析表达式．详细地讨论了应力波在横观各向同性含液饱和多孔介质中传播时,其速度曲面和波阵面的形状及性质．分析结果亦表明,纯固体中应力波传播的特征方程,是含液饱和多孔介质中应力波特征方程的特例．     

Two-material optimal design for nonlinear elastica

The work considers an optimal design problem in the context of nonlinear elastica. More specifically, we deal with finding the best way of mixing fixed amounts of two different elastic materials, so as to minimize the tip deflection of a cantilever beam loaded on its free extreme under the assumption of large deflections. Applying an optimality criteria method to the relaxed problem, simulations give us numerical evidence that the original design problem admits classical solutions (i.e. there is no microstructure) and those are the same as the respective ones for the case of small deflections.     

A homotopy analysis solution to large deformation of beams under static arbitrary distributed load

In this article, homotopy analysis method (HAM) is employed to investigate non-linear large deformation of Euler–Bernoulli beams subjected to an arbitrary distributed load. Constitutive equations of the problem are obtained. It is assumed that the length of the beam remains constant after applying external loads. Different auxiliary parameters and functions of the HAM and the extra auxiliary parameter, which is applied to initial guess of the solution, are employed to procure better convergence rate of the solution. The results of the solution are obtained for two different examples including constant cross sectional beam subjected to constant distributed load and periodic distributed load. Special base functions, orthogonal polynomials e.g. Chebyshev expansion, are employed as a tool to improve the convergence of the solution. The general solution, presented in this paper, can be used to attain the solution of the beam under arbitrary distributed load and flexural stiffness. Ultimately, it is shown that small deformation theory overestimates different quantities such as bending moment, shear force, etc. for large deflection of the beams in comparison with large deformation theory. Finally, it is concluded that solution of small deformation theory is far from reality for large deflection of straight Euler–Bernoulli beams.     

压电梁的多项式解（Ⅱ）——典型问题解析解

从正交各向异性压电介质平面问题,对于材料3个特征根互不相等情况下,以3个拟调和位移函数表达位移、电势、应力和电位移的通解出发,利用调和多项式的显式表达式,结合试凑法,给出了平面压电梁的若干典型问题的解析解,包括悬臂压电梁自由端作用横向集中力和点电荷,悬臂压电梁表面作用线性电势和均布载荷,以及两端简支压电梁作用均布载荷等的解析解．     

A reduced model to locate low contrast inclusions in poroelastic media

This paper introduces a numerical method to localize inclusions having slightly different elastic coefficients than those of a fully saturated poroelastic matrix, whose detection is often difficult. This method can be used to find weakly stiffer or softer objects in saturated soils or diseased biological tissues at early stages. To this end, we propose a reduced model from the Biot’s equations by splitting the fluid pressure into two parts: one embedded into an elasticity model and the other one used as a corrector term. By applying the small amplitude homogenization method, we can successfully retrieve the position and extension of inclusions in poroelastic media employing this simplified model. Numerical results show a good agreement for the location of inclusions when the contrast is below 30% stiffer or softer than the matrix, and for a noise level up to 5% for frequencies below 50 Hz.     

Transient wave propagation in a one dimensional partially saturated poroelastic column

Based on the theory of mixtures, a dynamic three phase model for partially saturated poroelasticity is established. This model is applied to a one-dimensional unsaturated poroelastic column and an analytical solution in the Laplace domain is deduced. By using the Convolution Quadrature Method, the solution in the time domain is obtained. Using some widely distributed porous materials as rock, soil and sediment, the wave propagation behavior in terms of displacement and pore pressure is examined. By neglecting the viscosity of the fluid, assuming very large fluid permeabilities, the second and third compressional waves are identified. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analytical solution for large deflections of a cantilever beam under nonconservative load based on homotopy analysis method

In this article, large deflection and rotation of a nonlinear beam subjected to a coplanar follower static loading is studied. It is assumed that the angle of inclination of the force with respect to the deformed axis of the beam remains unchanged during deformation. The governing equation of this problem is solved analytically for the first time using a new kind of analytic technique for nonlinear problems, namely, the homotopy analysis method (HAM). The present solution can be used in wide range of load and length for beams under large deformations. The results obtained from HAM are compared with those results obtained by fourth order Range Kutta method. Finally, the load‐displacement characteristics of a uniform cantilever under a follower force normal to the deformed beam axis are presented. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27:541–553, 2011     

Elastic solutions of a functionally graded cantilever beam with different modulus in tension and compression under bending loads

This paper considers a functionally graded cantilever beam with different modulus in tension and compression. The beam is subjected to bending loads, including pure bending, shear force at the free end and uniform pressure on the upper lateral, respectively. Its modulus values in tension and compression both change with the thickness coordinate as arbitrary functions, which could bring the beam a broader range of applications in engineering. The problem is treated as a plane stress case and described by Airy stress function. By using semi-inverse method, the elastic solutions for the beam are obtained, which can be easily degenerated into the ones for homogeneous beams. An example is finally presented to show the effect of nonhomogeneous materials with different modulus on the elastic field in a cantilever beam.     

考虑非局部剪切效应的碳纳米管弯曲特性研究

基于Hamilton(哈密顿)变分原理和非局部连续介质弹性理论,建立了新型非局部Timoshenko(铁木辛柯)梁模型（ANT）,推导了碳纳米管（CNT）的ANT弯曲平衡方程以及两端简支梁、悬臂梁和简支 固定梁的边界条件表达式,分析了剪切变形效应和非局部微观尺度效应对碳纳米管弯曲特性的影响.数值计算结果显示,碳纳米管的弯曲刚度随着小尺度效应的增强而升高.其次,这种小尺度效应对自由端受集中力的悬臂梁碳纳米管有明显作用,其刚度变化规律和其它约束条件的碳纳米管一样,这一点是ANT模型区别于普通非局部纳米梁模型的主要特点.经分子动力学模拟验证,ANT模型是合理分析碳纳米管力学特性的有效方法.     

Stability and vibration analysis of axially-loaded shear beam-columns carrying elastically restrained mass

Classical shear beams only consider the deflection resulting from sliding of parallel cross-sections, and do not consider the effect of rotation of cross-sections. Adopting the Kausel beam theory where cross-sectional rotation is considered, this article studies stability and free vibration of axially-loaded shear beams using Engesser’s and Haringx’s approaches. For attached mass at elastically supported ends, we present a unified analytical approach for obtaining a characteristic equation. By setting natural frequencies to be zero in this equation, critical buckling load can be determined. The resulting frequency equation reduces to the classical one when cross-sections do not rotate. The mode shapes at free vibration and buckling are given. The frequency equations for shear beam-columns with special free/pinned/clamped ends and carrying concentrated mass at the end can be obtained from the present. The influences of elastic restraint coefficients, axial loads and moment of inertia on the natural frequencies and buckling loads are expounded. It is found that the Engesser theory is superior to the Haringx theory.