复合材料层合的结构的非线性传热分析

发布 2021-06-02 21:16:28 阅读 4327

复合材料层合结构朔阴溯眼撇汾析。

杨自春。侮军工程大学.武汉,43

摘要本文基于复合材料传热分析的一阶热层合理论,考虑材料非线性的影。

响,建立了复合材料层合结构非线性传热分析的热层合理论及其有限元方程。数值算例。

研究了材料非线性对反对称复合材料层合板瞬态温度场的影响。

关健词友合材权非线性传热。

引言。复合材料热分析的主要内容是考虑强热各向异性,对复合材料结构的热传递进行稳态和瞬态计算。复合材料由于具有强的热各向异性、非均质性、非线性、热藕合性等特点,与各向同性结构相比,复合材料结构的热传导分析是一项十分复杂的工作。

作者从微元热传导的热流平衡方程和fou热传导定律出发,导出了复合材料层合结构传热分析的一。

阶和高阶热层合理论及有限元算法[[1但是该理论算法是建立在复合材料热物性不随温。

度变化的基本假定之上,即只限于研究线性热传导问题。实际上,上述假定对于每种复合材料都只在一特定的温度使用范围内才是正确的,一但超过这个温度界限,则表现出强烈的热非线性行为。目前,对复合材料结构的非线性热传导问题研究较少[31因此,进一。

步考虑材料非线性,深人研究复合材料结构非线性传热分析的理论和计算方法是十分必要的,而且具有重要的实际意义。

复合材料的热物性参数主要是指材料的热导率、对流换热系数和比热容。一般情况。

下,大多数复合材料的比热容随着温度的升高而增大,而热导率则随着温度的升高而减小[31本文将在文[[1文[[2研究的基础上,重点研究材料的热导率、对流换热系数及比热容。

均随温度变化的复合材料层合板的非线性热传导问题。数值算例将分析材料非线性对复合材料层合板非线性热传导瞬态温度场的影响。

层合板非线性热传导分析的基本理论及有限元方程。

为了叙述和计算的简便,假定复合材料的热导率kv、对流换热系数a和比热容c

均随温度作线性变化,这与多数复合材料热物性实验结论相一致[ku十、。约。

41。这种情况下有。

式中,入。,p非线性参数可由实验确定。这里采用一阶热层合理论[[1即假定温度沿厚度方向线性分布。

式中表示层合板中面温度表示沿层合板厚度方向的温度梯度,它们是x,y和时间t的函数,均与z无关;下标”表示对后一个坐标的求导。

借助((2试,并将((1式代人fou澈传导定律,得。

护l仓l一rrl一j

式中,4,表示沿坐标轴x,y方向的热流密度。线性热传导矩阵【们的元素祥见文。一一。肛。

十。1l,为节省篇幅,以下凡在文[[1中有详细推导、说明和列式的相关内容,本文不在。

赞述;而医(t)则是与温度相关的非线性热导率矩阵:

式中表示当前坐标系下的热传导系数,它可由材料主向的热传导。

率通过转轴变换加以确定[[1

.ls破|||

将((3试代人下式,并沿厚度方向对z积分,可得广义热力。

仁卜(5)一叭一。

y龙中熟流力”,熟流力弯矩”,两者统称为广义热式,qx鸟,q,称为m,表示。

5)式右边的最后一项即为材料热导率随温度变化而引起的非线力。

之z么乌么从叽。

k,』求解祥见文[[1

性项,它可表示为。一。公卜‘

一十-.工...

卜。中式,如株示温度场梯度列向量:

赚[几。}洁佃w=y二二,t"

么。乌-么一峡一叽。

了。..l记矛,r*十。

乳马。一么1么。

久。一峨一叭。

将式(4)和上式代人(6)式,并沿z向积分得。

一。佃。

们|列川」式中为层合板的非线性热传导矩阵,它们分别为。弋。厉。

一-刁。.洲。

0弋。y阴。

此成凡成0。d呱巩。

knl的各元素分别为。

弋=叉寸,.t冲h(k成d一。

一香n}k干,y.

刀。kql的各5m分别为:月。伙。

一(z(川。

叮一((z一一二,y.一。

仿照文[1]中传热有限元方程的推导步骤,将以上各式代入三维各向异性热传导的泛。么。于。

一一。1.,口二,,.

}0f芬0f;

叱g二。gam函表达式[[5并进行离散,由驻值条件可得到考虑材*[4线性后层合板瞬态温度场分析。

的有限元方程(推导过程略)为[[1

k}卜[k,

式中。各矩阵符号意义分别为:

co[了{ry枷pv知圈一一妙c刚一〕十因风]区-同。一 ̄

闻。刚圈。

韧。了。

nf一-升j

了rnj一弄。

ao-仇十:to知无〕

何}=又aup十:to二[nt

线性热传导问题,原则上并无任何困难。

上述各式中,凡未给出表达式的矩阵,在文[[1中均有说明而且有相应的计算公式。

应该指出,虽然上述算式是基于一阶热层合理论,但若用高阶热层合理论【2]分析非。

求解方法。非线性热传导平衡方程(12在求解方法上要比线性情况复杂,这是因为在每一个迭代。

时间步上,都需要形成新的热传导矩阵、非线性热容矩阵、非线性对流矩阵以及新的有效热流向量,从而使得计算比较费时。本节采用逐步积分方法求解方程(12具体步骤如。

下:)非线性热流平衡方程的线性化。

在每一个微元时间步上,利用增量原理有下列关系:

o十a_-十△at扩+ac相应有。

kc]一。砂护=砂i+o沪}

将上式代人到方程(12即可将非线性热传导平衡方程逐步线性化。

)初始步计算。

-0,此时二。,因此非线性项消失,故而形成线性热容矩阵。

k小线性热传导矩阵[k,及线性热流载荷识}仇}.

)用前一步算得的和t.,代)1.线性项,计算t时刻的非线性热传导。

矩阵{刘‘、非线性热流an.矩阵护}及非线性热viv阵【引‘,对时间采用差5}*式。

)计算‘时刻的温度及温度梯度。

)计算‘十△‘时刻的非线性热传导矩阵【k'非线性热流载荷矩阵lr*及非线性热容矩阵【kt对时间采用差f*式。

)计算‘十at时刻的温度及温度梯度to.艰+oa

)在每一个时间步上,需要进行多次的热流平衡迭代。若前后两次迭代的结果误差小于给定的误差容限,则终止该时间步上的迭代,然后进人下一时间步的计算。

数值算例。利用本文非线性热层合理论及有限元方法,研究材料非线性对复合材料层合板传热的。

影响。分析对象为四层石墨环氧反对称角铺层方板[10一其几何参。

数是二0.0外加热流垂直作用在层合板上表面中心区域,热流作用。

区面积为热流温度为500材料参数分别为:材料主轴的热导率对流换热系数比热容。=1刀kgk非线性材料热系数为。

层合板几何形状及加载热流如图1所示:个z口口。

图t所研究层合板的几何形状及加载热流。

层合板共划分64个平面矩形单元,节点总数81个。瞬态温度场分析时时间步长为。

.01秒。利用本文给出的非线性热层合理论和有限元方法,对其进行非线性瞬态热传导数值分析,并将计算结果与线性分析结果相比较。图2给出的是点1位置上,温度随时间变化的曲线,图3给出的是在时间t=1秒时,沿层合板上表面线c的温度分布。

由图2、图3可以明显看出材料的热物性〔复合材料的热导率、对流换热系数及比热容)对复合材料层合板非线性热传导有显著的影响,因而,在进行重要的受热复合材料结构设。

计时,材料热物性的非线性不容忽略。

比。图2层合板非线性热传导瞬态温度场分析。

刁3-0刁。

la(卜10s

图3层合板上表面的温度分布(线c)

结束语。本文研究了复合材料热物性参数(热传导率、对流换热系数和比热容)随温度变化的非。

线性热传导问题,基于一阶热层合理论,给出了复合材料层合板瞬态温度场分析的有限元。

方程及逐步积分算法。数值算例表明材料非线性对复合材料层合板的温度分布有显著影。

响。参考文献。

1(杨自春,黄玉盈,复合材料热层合理论及其有限元方法研究,海军工程学院学报。

杨自春,曹跃云,复合材料高阶热层合理论及其有限元方法研究,海军工程学院学报。

杨自春,复合材料层合板壳中一些非线性热相关问题的理论和实验研究,华中理工大学博士论。

文,19孔详谦.有限单元法在传热学中的应用.科学出版社。北京。19

80一。

复合材料的结构及作用

二 举例说明。聚乳酸 纳米碳管防静电复合材料。此材料是以纳米碳管为导电料通过球磨和密炼2种方法添加到聚乳酸基体中制备的防静电复合材料。具体工艺流程如下 纳米碳管的纯化处理 p cnt 纳米碳管功能化 f cnt 球磨法或密炼法混合 热压 成型。聚乳酸可以看做复合材料的基层,是复合材料的基材框架。pl...

rPETR MFMBGF复合材料的结构及力学性能研究

工程塑料应用 年,第 卷,第 期。复合材料的结构及力学性能研究。郑丽娜吕亚男魏金刚。代佳丽。王经武。郑州大学材料科学与 程学院,郑州。摘要采用自制的反应性多功能母料 及玻璃纤维 与废旧聚对苯二甲酸乙二酯 热机械反。应性共混,制得复合材料,对其化学结构 形态结构及力学性能进行研究。结果表明,中的环氧官...

线性结构在非线性结构中的应用

摘要 数据结构课程中数据的逻辑结构分为线性结构和非线性结构。数据结构中线性结构指的是数据元素之间存在着 一对一 的线性关系的数据结构。相对应于线性结构,非线性结构的逻辑特征是一个结点元素可能对应多个直接前驱和多个后驱。关键字 线性表,树,图。1.引言。数据元素相互之间的关系称为结构。有四类基本结构 ...