蜂窝夹层复合材料结构非线性传热分析.pdf

《蜂窝夹层复合材料结构非线性传热分析.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

复合材料学报第22卷第5期10月2005年ActaMateriaeCompositaeSinicaVol122No15October2005文章编号:100023851(2005)0520147208蜂窝夹层复合材料结构非线性传热分析3邱志平,林强(北京航空航天大学固体力学研究所,北京100083)摘要:考虑材料热物性随温度变化的影响,针对蜂窝夹层复合材料结构的特点,假定其内部的瞬态温度场的分布状况,并在此假定的基础上基于三维各向异性热传导理论及复合材料层合结构传热分析的热叠层方法,提出一种针对蜂窝夹层复合材料结构的改进的计算方法,并推导出蜂窝夹层复合材料结构瞬态温度场分析的有限元方程。应用这种方法对热物性非线性蜂窝夹层复合材料板的瞬态温度场进行分析,数值算例结果显示了材料特性的非线性对蜂窝夹层复合材料板热传导的影响,只有非线性程度很小时忽略非线性因素导致的误差才可被忽略;当非线性程度较强烈时,考虑非线性因素的计算结果与忽略非线性因素的计算结果就会产生很大偏差。关键词:蜂窝夹层复合材料;非线性传热;瞬态温度场中图分类号:V21414;TB33011文献标识码:ANONLINEARHEATTRANSFEROFHONEYCOMBSANDWICHCOMPOSITEMATERIALSTRUCTURE3QIUZhiping,LINQiang(InstituteofSolidMechanics,BeijingUniversityofAeronauticsandAstronautics,Beijing100083,China)Abstract:Basedonthecharacteristicsofthehoneycombsandwichcompositematerialstructures,thetransienttemperaturefieldofthestructuresisassumed.Dependingonthisassumptionandwiththeconcernoftheeffectoftemperatureontheheatcharacteristicofthematerial,amethodforhoneycombsandwichcompositematerialstruc2tureheattransferisobtained.Itisanimprovedapproachbasedonthe32Danisotropyheattransfertheoryandthecompositematerialthermallaminatedtheory.ThecorrespondingFEMequationsarealsoconstructed.Anexamplesolvedwiththismethodisgivenattheend.Therearethreedifferentlevelsfornonlinearheatcharacteristicsofthematerialdiscussedinthisexample.Theresultsofthenumericalexampleshowthatthenonlinearfactorscanbeneg2lectedonlywhenthedegreeofthenonlinearityisveryweak,andinmostinstancesthenonlinearfactorscannotbeneglectedatall.Keywords:honeycombsandwichcompositematerial;nonlinearheattransfer;transienttemperaturefield[3,4]复合材料结构由于具有强的热各向异性、非线等。最近RRolfes提出了与有限元结合的热叠层性、热耦合性等特点,与各向同性均质结构相比,理论,这种理论是解决复合材料叠层结构瞬态热分其热传导分析十分复杂。传统的解析方法只能解决析的一种有效方法,与三维热计算相比具有计算效[5][6,7,8]复合材料的一些简单热传导问题,无法求解复杂边率高的优点。李松年、杨自春等学者也应用界条件及复杂结构的传热问题。许多研究人员采用这种方法对复合材料叠层结构非线性瞬态热分析进有限元法对复合材料夹层结构的热传导进行了研究行了研究。但是以上各种方法都是对一般的复合材[1]计算,比较突出的结果有Padovan考虑各向异性料层合结构的分析,没有特别针对复合材料蜂窝夹的影响,研究了正交各向异性无限大层合板和层合层结构进行研究,因此这些方法并不十分适用于对[2]圆柱壳的热传导问题,Tamma等人对三维热分复合材料蜂窝夹层结构非线性瞬态热传导的计算。析模型进行了简化,提出了准三维温度场分析模型为了更好地分析和计算蜂窝夹层复合材料结构的热收稿日期:2004210213;收修改稿日期:2005201204基金项目:航天支撑技术基金项目(20032.2BH03)资助通讯作者:邱志平,教授,研究方向:结构可靠性E2mail:zpqiu@buaa1edu1cn©1995-2005TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved. ·148·复合材料学报传导问题,在前人的基础上还有一些工作要做。的,并近似认为面板中的温度分布规律为事先假定在进行复合材料传热分析的计算时通常假定复的蜂窝结构部分温度分布的延伸,即忽略由材料性合材料结构热物性是恒定的,一般主要是指热物性质变化而导致的温度分布规律的变化。有关面板中不随温度变化,即只限于线性热传导问题的计算。的温度分布规律的假定显然会给计算结果带来误但是实际上,上述假定对于每种复合材料都只在特差,但是考虑到面板处在于环境直接接触的边界定的温度范围内才是正确的,虽然有一些复合材料上,其温度分布受环境温度影响很大,而物性的影的热物性在较小的温度范围内不随温度变化,但一响处于次要地位,并且相比于芯层厚度面板的厚度旦超过这个温度界限,则表现出强烈的热非线性行一般都很小,即使温度梯度误差较大,面板上的温为。复合材料的热物性参数主要是指材料的热导率度误差也不会很大,所以忽略物性变化而直接认为kij、对流换热系数α和比热容c。一般情况下,大温度梯度是连续的而导致的误差在一般情况下是可多数复合材料的比热容随着温度的升高而增大,而以接受的。由传热时蜂窝夹层结构中蜂窝结构是主热导率则随着温度的升高而减小。而实践证明,很体成分的假设,简化面板部分对传热的影响,将整多情况下假定材料热物性是恒定的将导致不能忽视个结构的传热看作是蜂窝部分的性质起决定作用,的误差,必须考虑材料的热物性随温度变化造成的面板部分性质起次要作用的传热就是本文作者的中影响。心思想。本文作者在前人工作的基础上,针对复合材料对于结构中材料的各种热物性的变化规律,无蜂窝夹层结构的特性,并考虑热物性随温度变化的法作出通用的假定,实际应用中可根据实验或者经影响,作出合理假定,并基于三维各向异性热传导验作出适当的假定。另外,一般的蜂窝夹层结构可理论及复合材料的一般理论,采用与有限元方法结看作以中面为对称面的对称结构,当外界环境或边合的热叠层理论方法,得出分析复合材料蜂窝夹层界条件已知时,其内部温度分布规律的大致趋势是结构非线性热传导的一般方法,并且利用数值算例可以大致判断的,如只有一侧面板被加热时,沿厚显示了热物性的变化对结构瞬态温度场的影响,以度方向温度分布一般不可能是对称的,所以假设蜂及非线性的强弱对分析结果的影响。窝结构内部温度分布规律时需要根据实际情况具体分析,只要温度及温度梯度变化连续即可。本文作1基本假定者在以下的分析中,以一种具体情况为例,为了推和一般普通的复合材料层合结构相比,蜂窝夹导结论的过程简便起见,采取以下假定:层复合材料结构有其独有的特性。首先,蜂窝夹层(1)假定结构在下表面承受热载,上表面为自结构层数较少,并且两侧面板厚度与芯层厚度相差由边界,温度沿蜂窝夹层结构芯层厚度方向为奇函很大,整个结构中芯层占有很大比重;而一般的复数多项式合材料层合结构层数较多,每层厚度相差不大或者T(x,y,z,t)=T0(x,y,t)+T0z(x,y,t)·z3基本相同,且每层的材料特性都不相同。其次,虽+T03z(x,y,t)·z(1)然蜂窝夹层板壳结构不是连续体,严格说来有关材其中:T0(x,y,t)表示夹层板中面温度;T0(x,y,t)、料连续的假设在这里不再适应,但就一般常见的夹T0z(x,y,t)、T03z(x,y,t)为待定系数,它们是x、y层板壳而言,夹芯的格子都很小,制造夹芯的金属和时间t的函数,均与z无关。波或非金属波也很薄。因此近似认为是连续的,不(2)假定芯层材料的热导率kxij、对流换热系致引入较大误差。所以宏观上看蜂窝夹层结构中占数αx和比热容cx均随温度T作线性变化,在这种很大比重的蜂窝部分可看作连续的热物性各向异性假定下有材料,可用折合后的连续体的热物性系数代表原来kxij=kx0ij(1+λxij·T)的蜂窝结构的传热性质。αx=αx0(1+βx·T)(2)针对以上分析,本文作者对蜂窝夹层复合材料cx=cx0(1+γx·T)结构中的瞬态温度场作出以下假设:根据蜂窝芯层式中:kx0ij(i,j=x,y,z)表示当前坐标系下的热传材料本身的性质,宏观上假设其内部(即芯层)的温导系数,它可由材料主向的热传导率通过转轴变换度的变化是连续的,并且温度的变化率也是连续加以确定。©1995-2005TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved. 邱志平,等:蜂窝夹层复合材料结构非线性传热分析·149·(3)假定面板材料的热导率kmij、对流换热系数αm和比热容cm均随温度Tm作线性变化,在这种假定下有kmij=km0ij(1+λmij·Tm)αm=αm0(1+βm·Tm)(3)cm=cm0(1+γm·Tm)式中:km0ij(i,j=x,y,z)表示当前坐标系下的热传导系数,它可由材料主向的热传导率通过转轴变换加以确定。这里面板材料的热导率kmij、对流换热系数αm和比热容cm是指复合材料面板层合结构的折合数值,可按照一般复合材料层合结构的折合处图1示意图与坐标系Fig11Sketchandcoordinate理方法得到。以上关于温度分布及热物性随温度变化的规律式中:qx、qy、qz表示沿坐标轴x、y、z方向的热的假定是一种非常简单的假定,只是为了下面分析流密度。线性热传导系数矩阵[k]、与温度相关的和推导方便。在实际中对于不同的具体问题可以有非线性热导率矩阵[k(T)]及{T′}分别为不同的假定,可根据具体情况来处理,但是最后解k0xxk0xy0决问题的基本方法是相同的,均可按照下面的过程[k]=k0xyk0yy0(5)来处理。00k0zzk0xxλxxTk0xyλxyT02基础理论[k(T)]=k0xyλxyTk0yyλyyT0(6)根据Fourier热传导定律,由式(1)、式(2)、00k0zzλzzT式(3)三式可得T9T9T9Tqx{T′}=9x9y9z(7){q}=qy=([k]+[k(T)]){T′}(4)将式(1)代入式(7)可得qz9T9T09T0z9z9x9x9x9x9T9T09T0z9z{T′}==+·z+·T0z9y9y9y9y9T9T09T0z9z9z9z9z9zT0zT03z9T09x9T039T03z9z9y9x9x10100z0z309T039T03z39z239z+·z+·T03z=03·z0100z0z(8)9y9y9T0z9T30000100009x03z9z9z9z9T0z9y9T03z9x9T03z9y©1995-2005TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved. ·150·复合材料学报22令9Tx19Tm+kxzzdv1+kmxx9T09T09T09z∫v229x{T0′}=T0zT03z9x9y9z9T22m9TmT+kmyy+kmzzdv29T0z9T0z9T03z9T03z9y9z(9)9x9y9x9y9T221m9Tm3+kmxx+kmyy10100z0z0∫v329x9y[SŠ]=03·z20100z0z329Tm+kmzzdv3-h(T·Tw0000100009z∫Sc(10)12-Tw)dSc-T·qSdS2-T·qBdv2∫2可得S2∫v9T{T′}=[SŠ]{T0′}(11)+ρ·cx·T·dv1∫v19t9T3有限元方程的推导+ρ·cm·T·dv2∫v29t芯层导热系数主轴的方向与面板导热系数主轴9T+ρ·cm·T·dv3(17)的方向可能是不重合的,并且这两者与分析所用的∫v39t坐标系的方向也可能是不重合的,这时芯层与面板上式中v1、v2、v3分别表示单元芯层、单元上面的导热系数矩阵就都包含有耦合项导热系数。为了板、单元下面板的体积。推导的简便和结果表达的简单,在这里假设芯层导令热系数主轴的方向、面板导热系数主轴的方向、分9T09T09T0析所用的坐标系的方向三者是重合的,这对问题的{T}=T0T0zT03z9x9y9z本质没有影响。9TT0z9T0z9T03z9T03z(18)热传导问题的控制方程及边界条件分别为:9x9y9x9y99T99T99T110100z0z30控制方程kxx+kyy+kzz+qB9x9x9y9y9z9z23[SŠ]=003·z0100z0z9T=ρ·c9t(12)0000010000边界条件:(19)3给定温度边界S1T=T(x,y,z,t)(13){R}={1zz0000000}9T(20)热流边界S2kn=hr1(Tr-Tw)(14)9n取单元上的T0(x,y,t)、T0z(x,y,t)、T03z(x,9T对流边界Sckn=α(Tf-Tw)(15)y,t)的插值函数为相同的平面插值函数N,将式9n(18)至式(20)代入泛函表达式可得初始条件Tt=0=T(x,y,z)(16)1TT将结构在中面上划分为有限个平面单元,在单πe=[T][S][kx][T][S]dv1∫v12元中的控制方程和边界条件仍和式(12)至式(16)一1TT致。根据这个原则,单元上的控制方程及其边界条+[T][S][kx(T)][T][S]dv1∫v12件对应的泛函是1TT9T29T29T2+[T][S][km][T][S]dv2π1∫v22e=kxx+kyy+kzzdv∫v29x9y9z1TT12+[T][S][km(T)][T][S]dv2-h(T·Tw-Tw)dSc-T·qSdS2∫v22∫22Sc∫S21TT9T+[T][S][km][T][S]dv3-T·qBdv+ρ·c·T·dv∫v32∫v∫v9t9T221TT=1kx9Tx+[T][S][km(T)][T][S]dv3xxx+kxyy∫v32∫v129x9y©1995-2005TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved. 邱志平,等:蜂窝夹层复合材料结构非线性传热分析·151·TTTTT-∫α0([T][R]-[Tw])[Tw]dSc+[N][S][km][S][N]dv3Sc∫v3TTTTT3-∫α0β[T][R]([T][R]-[Tw])[Tw]dSc[K′1]=(T0+z·T0z+z·T03z)Sc∫v1TTTTTT[N][S][kx(T)][S][N]dv1-[T][R]qSdS2-[T][R]qBdv∫2S2∫v3+(T0+z·T0z+z·T03z)TT∫9[T][R]v2+ρx·cx0·dv1∫v19tTT[N][S][km(T)][S][N]dv2TT9[T][R]+ρx·cx0γx[R]·[T]·dv1(T3∫9t+0+z·T0z+z·T03z)v1∫v3TT9[T][R][N]T[S]T[k+ρm·cm0·dv2m(T)][S][N]dv3∫v29tTTTT[K2]=α0[N]{R}{R}[N]dSc+ρm·cm0γm9[T][R][R]·[T]·dv2∫Sc∫v29t3TT[K′2]=α0β(T0+z·T0z+z·T03z)+ρm·cm09[T][R]·dv3∫Sc∫v39tTT[N]{R}{R}[N]dScTT9[T][R]+ρm·cm0γm[R]·[T]·dv3(21){P}={P1}+{P2}+{P′2}∫v39tTT取δπe=0,并且在单元i上对δπe进行变分{P1}=qS[N]{R}dS2∫2S2时,认为另外单元j上的温度T为常数。经过整TT理可得有限元方程{P2}=∫α0·Tw[N]{R}dScSc[C]{TÛ}+[K]{T}={P}(22)3{P′2}=∫α0β(T0+z·T0z+z·T03z)·Sc可知{T}中未知量只有3个:T0、T0z、T03z。TTTw[N]{R}dSc式中各矩阵的具体表达式为得到单元的热刚度矩阵、热载向量后,组装即[C]=[C1]+[C2]可得到总体热刚度矩阵和总体热载列向量,最后可TT[C1]=ρx·cx0[N]{R}{R}[N]dv1∫v1得到一个非线性方程组,可用数值解法来进行求解,如牛顿法及其各种改进形式等。TT+∫ρm·cm0[N]{R}{R}[N]·dv2v2这种方法的实质是在一般复合材料结构传热的TT热叠层理论的基础上,根据蜂窝夹层复合材料结构+∫ρm·cm0[N]{R}{R}[N]·dv3v3整体来看分层较少的特点,略去热导率kij、对流ρ3换热系数α和比热容c的折合计算而分别严格计入[C2]=∫x·cx0·γx(T0+z·T0z+z·T03z)v1芯层与面板热导率kij、对流换热系数α和比热容TT[N]{R}{R}[N]dv1cm的一种改进的计算方法,并且根据蜂窝夹层复3+∫ρm·cm0·γm(T0+z·T0z+z·T03z)合材料本身的特点,作出了一些不同于一般复合材v2料层合结构的假设,如结构温度分布及温度梯度分TT[N]{R}{R}[N]dv2布连续的假设和减弱面板芯层之间物性变化的影响3+∫ρm·cm0·γm(T0+z·T0z+z·T03z)的假设,而根据这些假设来改进复合材料结构传热v3TT热叠层理论得到的本文中的方法更接近蜂窝夹层结[N]{R}{R}[N]dv3构的实际,比一般的复合材料结构传热的热叠层理[K]=[K1]+[K′1]+[K2]+[K′2]论更适合以蜂窝夹层复合材料结构传热的计算。TT[K1]=∫[N][S][kx][S][N]dv1v14数值算例TT+∫[N][S][km][S][N]dv2v2利用本文作者提供的方法,对以下的算例进行©1995-2005TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved. ·152·复合材料学报研究。如图2所示,夹层为石墨环氧45°角铺层方板,芯层为蜂窝夹层,结构的几何参数为a=b=1m,夹层厚度h0=01005m,芯层厚度h=0104m,2四周绝热,外加热流垂直作用在板下表面0109m上,热流作用中心在板中心偏左0115m处。热流2q=2000W/m,环境温度T∞=300K,结构具有均匀的初始温度T0=400K。材料的热物性参数分别为:夹层材料主轴的热导率K11=35W/(m·K),K12=K21=30W/(m·K),K22=30W/(m·K),K33=25W/(m·K);芯层蜂窝材料主轴的热导率K11=40W/(m·K),K12=K21=30W/(m·K),K22=30W/(m·K),K33=25W/(m·K)。2夹层材料对流换热系数α0jc=10W/(m·K);芯2层材料对流换热系数α0x=10W/(m·K)。夹层材料比热容c0jc=1000J/(kg·K);芯层材料比热容c0x=1000J/(kg·K)。对以下3种情况进行计算:(a)非线性材料热系数为:λij=β=γ=01001;(b)非线性材料热系数为:λij=β=γ=01005;(c)非线性材料热系数为:λij=β=γ=0101。图2蜂窝夹层复合材料板结构示意图Fig12Sketchforhoneycombsandwichcompositematerialplate图3和图4分别是3种条件下10s内下表面中点处温度变化曲线和10s后中面45°对角线上温度分布曲线。图中实线表示的是考虑材料热物性随温度变化时的计算结果,而短直线虚线表示的是不考虑材料热物性随温度变化、假定物性值为常数时的图310s内下表面中点处温度变化计算结果。由图3(a)和图4(a)可见,当材料的非Fig13Temperaturevariationforthemiddlepoint线性特性并不十分显著时,虽然非线性特性对结构ofbottomsurfacein10s瞬态温度场有影响,但是影响很小,且没有改变其态温度场的变化趋势及分布,并且两者之间已经有原有的变化趋势。所以当材料的非线性特性并不显明显的差值出现。我们可以看到这时的材料的非线著时,可以忽略其影响。性特性其实并不是十分得显著,但是形成的影响却当材料的非线性特性逐渐增大时,如图3(b)已经是不能忽略的了,如果这时继续忽略其影响,与图4(b)所示,这时非线性特性已经开始影响瞬将对结果造成很大的误差。而且通过前两组图片,©1995-2005TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved. 邱志平,等:蜂窝夹层复合材料结构非线性传热分析·153·时的情况。通过图3(c)与图4(c),我们可以清楚地看到由于材料的非线性特性导致的结果,这时如果忽略材料非线性特性的影响,得到的结果根本不能相信。综上,当材料物性的非线性系数较小时,与假定材料物性为常值时相比,结构温度变化所受影响较小,但是随着时间的增加,其累积结果也相当显著;而当非线性系数较大时,结构温度所受的影响十分明显,即便加倍计入计算误差所造成的影响,由较大的非线性系数所造成的影响也是令人吃惊的,已经到了不能不考虑的地步。由此可见,考虑传热过程中的非线性因素的影响是十分必要的。在此指出,本文作者采用的具体的有关温度分布规律及物性变化规律的假定是相当简陋的,计算结果不可避免地带有误差,只能作为示例使用。5结论本文作者研究了蜂窝夹层复合材料结构热物性参数(热传导率、对流换热系数和比热容)随温度变化的非线性热传导问题,基于对一般复合材料结构传热热叠层理论,针对蜂窝夹层复合材料结构的特点作出改进,给出了复合材料蜂窝夹层板瞬态温度场分析的有限元方程。数值算例表明:材料非线性对复合材料蜂窝夹层板的热传导有显著影响。最后再次指出,按照本文作者的思想并不难导出基于其他各种假定的复合材料蜂窝夹层结构热传导理论及其有限元方程。参考文献:[1]PadovanJ.Steadyconductionofheatinlinearandnonlinearfullyanisotropicmediabyfiniteelements[J].JournalofHestTransfer,1974,96:313-318.[2]TammaKK,YurkoAA.Aunifiedfiniteelementmodelinganalysisapproachforthermal2structuralresponseinlayeredcomposites[J].ComputStruct,1988,29:743-754.[3]RolfesR.Higherordertheoryandfiniteelementforheatconductionincomposites[A].In:LewisRW,ChinJH,HomsyGM,eds.StanfordUSA7thIntConfonNumer图410s后中面45°对角线上温度分布MethinThermalProbl[C].Swansea:PineridgePress,Fig14Temperaturedistributionforthediagonal1991.880-889.ofmiddlesurfaceafter10s[4]RolfesR.Nonlineartransientthermalanalysisofcompositesstructures[A].In:HirschCh,ZienkiewiczOC,OnateE,我们可以看到即使材料的非线性特性十分微弱,但eds.NumericalMethodsinEngineering92.Brussels,Bel2是随着时间历程的推移,其造成的影响也将会变得gium:ProceedingsoftheFirstEuropeanConferenceonNu2不可忽略。mericalMethodsinEngineering[C].Brussels,Belgium,最后来看一下材料的非线性特性已经较为明显1992.653-659.©1995-2005TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved. ·154·复合材料学报[5]李松年,邓忠民.复合材料叠层板的非线性热分析[J].北compositestructure[J].JournalofNavalAcademyofEn2京航空航天大学学报,1996,22(3):290-295.gineering,1998(1):1-7.LiSongnian,DengZhongmin.Nonlineartransientthermala2[8]杨自春.复合材料层合结构的非线性传热分析[J].海军工nalysisofcompositelaminate[J].JournalofBeijingUni2程大学学报,2000(5):9-13.versityofAeronauticsandAstronautics,1996,22(3):YangZichun.Nonlinearheattransferofcompositematerial290-295.laminatedstructure[J].JournalofNavalAcademyofEn2[6]杨自春,黄玉盈.复合材料新的热层合理论及其有限元方法gineering,2000(5):9-13.研究[J].海军工程学院学报,1995(3):29-34.[9]俞昌铭.热传导及其数值分析[M].北京:清华大学出版YangZichun,HuangYuying.Anewthermallaminatetheo2社,1982.ryandfiniteelementapproachforthetransienttemperatureYuChangming.HeatTransferandItsNumericalAnalysisfieldanalysisofcompositelaminates[J].JournalofNaval[M].Beijing:TsinghuaUniversityPress,1982.AcademyofEngineering,1995(3):29-34.[10]张志民,张开达,杨乃宾.复合材料结构力学[M].北京:[7]杨自春,曹跃云.复合材料结构传热分析新的高阶热层合理北京航空航天大学出版社,1993.论研究[J].海军工程大学学报,1998(1):1-7.ZhangZhimin,ZhangKaida,YangNaibin.CompositeMate2YangZichun,CaoYueyun.Anewhigh2orderthermallami2rialStructuralMechanics[M].Beijing:BeijingUniversityofnatetheoryforthetransienttemperaturefieldanalysisofAeronauticsandAstronauticsPress,1993.©1995-2005TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved.