均布载荷作用下功能梯度悬臂梁的弹性半逆解.docxVIP

均布载荷作用下功能梯度悬臂梁的弹性半逆解 0 功能梯度材料的弹性力学分析方法 功能梯度材料作为材料设计的概念，日本科学家在上个世纪80年代中期首次提出了该问题。为了解决航天技术行业在高温环境下的材料热应力问题，第一个在宇宙高温环境下使用的功能梯度材料通常由耐高温工程陶瓷和高承载能力的金属组成。其设计理念是在材料制备过程中连续控制陶瓷和金属体积的分布，并在空间位置上显示宏观材料的特性。没有明显的材料性能界面。为了缓解高温和高使用环境下材料的内在热应力，现在日本科学家发现了电子、化学、光学、声音、医学和其他技术领域的功能梯度材料。 梁、板、壳是功能梯度材料部件的主要结构形式,主要的理论分析方法有简化模型法、层合模型法和精确解法等.简化模型法将经典的或考虑剪切变形的一维梁理论或二维板、壳理论直接套用到功能梯度材料结构的分析上,但是这些理论的基本假设是基于均匀或层合材料做出的,并不能完全适用于功能梯度材料.层合模型法将功能梯度材料梁、板、壳等结构沿厚度方向分成若干层,每一层等效为均匀材料,当分层足够多时,分析结果将趋近于弹性力学精确解.精确解法则是严格求解问题的控制方程和边界条件,目前只能获得特殊边界条件(如简单支撑)下的结果.Sankar给出了功能梯度简支梁受横向载荷作用问题的二维弹性解,并采用Euler-Bernoulli梁理论的基本假设建立了功能梯度梁的一维简化模型.但是Sankar的弹性解是在简支边界条件下并假设材料弹性模量沿梁厚度方向呈指数函数变化的前提下获得的. 本文从正交各向异性材料基本方程出发,考虑材料常数沿梁厚度方向梯度变化,采用弹性力学半逆解法,求取了按任意梯度函数分布的悬臂梁在均布载荷作用下的通解.并给出了指数函数梯度变化时的算例. 1 材料本构模型 考虑如图1所示的正交各向异性功能梯度悬臂梁(长度为l,宽度为b,厚度为h).当三维弹性介质在y轴方向的尺寸较其它两轴方向尺寸小得多,且全部载荷均作用于x-z平面内且不沿y轴方向变化时,问题可简化为x-z平面内的平面应力问题,其平衡方程为 ?σx?x+?τzx?z=0,?τzx?x+?σz?z=0(1) 物理方程为 εx=s11σx+s13σz,εz=s13σx+s33σz,γzx=s44τzx(2) 几何方程为 εx=?u?x,εz=?w?z,γzx=?u?z+?w?x(3) 应变协调方程为 ?2εx?z2+?2εz?x2-?2γzx?z?x=0(4) 上述各式中u和w为位移,σx、σz和τzx为应力,εx、εz和γzx为应变.s11、s33、s13、和s44为弹性材料的柔度系数,对于功能梯度材料,其材料常数沿厚度方向连续变化,因此是坐标z的连续函数.在本文的研究中,我们假设上述材料常数沿厚度方向按同一函数规律变化 sij=s0ijF(z)(5) sij代表s11、s33、s13和s44,s0ij为相应的材料常数在z0(z0满足F(z0)=1)处的值. 根据圣维南原理,自由端(x=0)的边界条件可写为 ∫h/2-h/2σxdz=0,∫h/2-h/2σxzdz=0,τzx=0(6) 下表面(z=h/2)的自由边界条件为 σz=τzx=0(7) 上表面(z=-h/2)作用均布竖向荷载,边界条件为 σz=-q,τzx=0(8) 固支端的边界条件可近似为:在z=0,x=1处 u=w=?w?x=0(9) 2 应力、应变和位移的公式 引入应力函数U,使满足 σx=?2U?z2,σz=?2U?x2,τzx=-?2U?z?x(10) 将式(2)和(10)代入应变协调方程(4)得 ?2?z2(s11?2U?z2+s13?2U?x2)+??z(s44?3U?x2?z)+s13?4U?x2?z2+s33?4U?x4=0(11) 令 U=x22f(z)+xf1(z)+f2(z)(12) 则有 {σx=x22d2f(z)dz2+xd2f1(z)dz2+d2f2(z)dz2τzx=-xdf(z)dz-df1(z)dz,σz=f(z)(13) 由力边界条件:x=0时,τzx=0,得到df1(z)/dz=0,则f1(z)为常数项,对应力无影响,可略去. 为了表达方便,引入如下符号 {Ηi(z)=∫z0zidzF(z)?Ηi(z)=∫z0Ηi(z)dz?Η1i(z)=∫z0?Η1i(z)dz?Η11i(z)=∫z0?Ηi(z)zdz?i=0?1,2,3,4{G0(z)=∫z0F(z)dz,G1(z)=∫z0F(z)zdz?G0(z)=∫z0G0(z)dz,G10(z)=∫z0G0(z)F(z)dzG110(z)=∫z0G0(z)F(z)zdz,G11(z)=∫z0G1(z)F(z)dz{Τ0(z)=∫z0Η0(z)F(z)dzΤ1(z)=∫z0Η1(z)F(z)dz?Τ0(z)=∫z0Τ0(z)d