功能梯度梁静力与动力分析的格林函数法.docx

该【功能梯度梁静力与动力分析的格林函数法 】是由【科技星球】上传分享，beplayapp体育下载一共【27】页，该beplayapp体育下载可以免费在线阅读，需要了解更多关于【功能梯度梁静力与动力分析的格林函数法 】的内容，可以使用beplayapp体育下载的站内搜索功能，选择自己适合的beplayapp体育下载，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此beplayapp体育下载到您的设备，方便您编辑和打印。功能梯度梁静力与动力分析的格林函数法??邓先琪,苏成,2,马海涛(,广州510640;,广州510640;,广州510405)功能梯度材料(Functionallygradedmaterial)[1]指由两种或两种以上不同材料复合而成,且各组成成分的结构和性能在空间上沿着某一方向呈连续梯度变化的一种新型非匀质复合材料。由于其新颖的材料可设计性思想和优良的结构性能,近年来引起学者的广泛关注和研究[2?5]。文献[6?8]对功能梯度材料的力学性能与研究进展进行了综述。目前,关于材料性能沿梁高度方向梯度变化的功能梯度梁静动态响应分析问题,解析方法和数值方法已有大量的研究。Sankar[9]假定材料弹性模量沿梁高度方向按指数函数变化,采用Euler-Bernoulli梁理论给出了功能梯度简支梁在正弦分布荷载作用下的弹性力学解;而Zhong和Yu[10]则基于Airy应力函数给出了功能梯度悬臂梁的二维弹性力学解。Li[11]提出了功能梯度Timoshenko梁和Euler-Bernoulli梁的静力与动力分析统一列式。Simsek[12]利用不同的高阶梁理论分析了功能梯度梁的基本频率。Chakraborty等[13]基于一阶剪切变形理论提出了一种新的梁单元来研究功能梯度梁的静力、自由振动和波动力学行为;而Kadoli等[14]则基于高阶剪切变形理论,提出了高阶梁单元来进行功能梯度梁的静力分析。Filippi等[15]基于一类统一列式提出了一种功能梯度梁单元模型,并对功能梯度梁静力问题做了相关的研究。Pradhan和Murmu[16]采用微分求积法研究了材料参数沿高度方向呈幂函数形式分布的弹性地基功能梯度梁和层合梁的振动问题。Jiao等[17]基于一阶剪切变形与经典梁理论,采用微分求积法分析了功能梯度梁的线性弯曲问题。李世荣和刘平[18]基于Euler-Bernoulli梁理论,比较了功能梯度梁和均匀梁的控制微分方程,得到了功能梯度梁与均匀梁静动态解之间的相似转换关系。王瑄和李世荣[19]基于Levinson高阶剪切变形理论,研究了两端简支功能梯度梁自由振动频率与经典梁理论下参考匀质梁自由振动频率之间的解析转换关系。以上文献中仅考虑了材料沿梁高度梯度变化的情况。相比于材料性能沿梁高度方向梯度变化的功能梯度梁,关于轴向功能梯度梁静力与动力分析的研究相对较少。轴向功能梯度梁的控制微分方程本质上是一个变系数常微分方程,通常情况下很难获得该类方程的解析解。在数值方法方面,Murín和Kuti?[20]采用传递系数描述功能梯度梁在空间上变化的弹性模量,建立了功能梯度Euler-Bernoulli梁的等效刚度矩阵。Alshorbagy等[21]提出了材料沿轴向或高度方向梯度变化的功能梯度梁动力特性分析的有限元法。汪亚运等[22]基于切比雪夫多项式展开技术研究了轴向功能梯度变截面梁的自由振动问题。本文主要研究弹性模量沿轴向梯度变化的功能梯度Euler-Bernoulli等截面梁的静力与动力分析问题。由于这种梁的弹性模量沿轴向呈函数变化,因此采用传统有限元法进行结构分析时,不可避免地要进行有限单元网格细分,随之带来的问题是计算规模的增大与计算效率的降低。为了提高功能梯度梁静动力问题的计算精度与效率,本文首先将功能梯度梁静力分析的控制微分方程转化为与匀质材料等截面梁静力分析控制微分方程相一致的形式,然后利用无限长匀质材料等截面梁静力问题的格林函数[23],提出了功能梯度梁静力分析的格林函数法。在此基础上,进一步推导获得功能梯度梁的柔度矩阵,据此建立功能梯度梁的运动方程,开展功能梯度梁的动力特性分析和动力响应分析。:式中:x为梁的轴向坐标;u(x)和v(x)分别为梁的轴向位移和横向位移;E为参考匀质梁的弹性模量;η(x)为反映弹性模量沿轴向变化的函数;A和I分别为梁的截面面积和截面惯性矩;p(x)和q(x)分别为梁的轴向和横向分布荷载。根据微分运算关系,式(1)可以转化为:式中,F1(x)和F2(x)分别为梁的等效分布荷载,表示为:式中:由式(2)可以看出,功能梯度梁静力分析的控制微分方程可以写成与匀质梁静力分析的控制微分方程相似的形式,唯一的差别在于式(2)中的等效分布荷载与方程左端梁的位移相关联。尽管如此,仍可以利用匀质梁静力问题的格林函数[23]进行求解。利用匀质梁静力问题的格林函数进行求解时,在内力计算上存在一些不同。考虑到材料的非匀质性,功能梯度梁内力与位移之间的微分关系为:式中,b(x)见式(4)。由式(5)可见,在功能梯度梁的静力分析中,内力与位移微分关系在形式上与匀质梁问题的内力与位移微分关系并不完全一致。因此,在采用匀质梁静力问题的格林函数进行求解时,需要考虑函数η(x)及修正项b(x)M(x)对内力的影响。利用式(5),可将式(3)给出的等效荷载表示为:由式(6)可见,等效荷载可以写成与内力N(x)、Q(x)和M(x)相关的形式,因此在采用格林函数法分析时,需要建立关于内力响应的补充方程,方可求解。,可以将功能梯度梁转化为如图1所示的等效匀质梁ab进行求解。其中,x?y为梁的坐标系,梁端点轴向坐标为xa和xb,弹性模量取为参考匀质梁的弹性模量E,截面面积和截面惯性矩分别取为功能梯度梁的截面面积A和截面惯性矩I,轴向和横向分布荷载分别取为式(6)所示的等效分布荷载F1(x)和F2(x)。[24],将等效匀质梁ab嵌入到具有相同材料性质与截面性质的无限长梁中,并在梁ab域外ξj(j=1,2)处布设6个大小未知的集中虚荷载Xij(i=1,2,3;j=1,2)。根据叠加原理,在域内等效荷载Fi(x)(i=1,2)与域外虚荷载Xij(i=1,2,3;j=1,2)的共同作用下,梁ab上任意一点的位移与内力响应为理论上,此式与式(10)在同一点处应给出相同的内力值。基于此,在梁域内选择(l+1)个点,利用轴力N(x)相等可以得到(l+1)个方程;类似的,在梁域内选择(m+1)个点,利用弯矩M(x)相等可以得到(m+1)个方程。最终可得:式中,H的各行取决于式(13)系数矩阵在所选择点处的取值,而A?、B?和C?的各行取决于式(10)中对应影响矩阵在所选择点处的取值。联立式(12)和式(14),可以解得:式中:式中,I为6×6的单位矩阵。从以上推导可见,由式(15)求得的基本未知向量X和Z分别包含6个虚荷载和(l+m+2)个内力系数,基本未知量数目与梁的子段划分无关,即与子段数n无关。将X和Z代入式(10),即可得到功能梯度梁内任意一点的位移和内力响应为:式中: