ANSYS-第六章热边界条件.ppt

第6章,复杂的，时间和/或空间变化的边界条件,6-2,表格化热边界条件(载荷),以前，施加在几何模型和/或有限元模型上的载荷都是常数或随温度变化的数值。当要施加复杂的随时间和空间变化的载荷时，表格化的热边界条件就更加方便。定义-表格化边界条件就是通过表格施加的边界条件，其数值经过计算并在求解过程中施加。,6-3,表格化边界条件-基本变量,基本变量时间全局笛卡尔坐标系中坐标温度速度,用于Fluid116单元，1-D热流单元,相应ANSYS标记TIMEX,Y,ZTEMPVELOCITY,表格由下列基本变量组成:,6-4,其它表格定义,独立变量直接或间接随基本变量变化并形成表格的变量，如，hf=f(time)换热系数(独立变量)是基本变量时间的函数hf=f(x,rpm)rpm=g(time)hf和rpm是独立变量,6-5,其它表格定义,多元表格表格有多于一个的变量，如hf=f(x,y,time)循环表格表格要引用其它带有独立变量的表格，如hf=f(x,y,rpm)rpm=g(time),6-6,表格输入的应用,表格化边界条件可以使用于所有实体模型和/或有限元模型热边界条件:使用D,DK,DL,DA命令施加的温度使用F,FK命令施加的热流使用SF,SFE,SFL,SFA命令施加的对流使用SF,SFE,SFL,SFA命令施加的热流量使用BF,BFE,BFL,BFA,BFV命令施加的内部热生成,使用HELP命令查看上述命令相应的菜单项。,6-7,表格化边界条件的重要说明,如果表格化边界条件不是时间的函数，他们可以很好的随时间阶跃施加，而不受KBC命令的约束(见第4章)。,6-8,基本变量-边界条件应用,热边界条件施加温度热流换热系数介质温度热流量内部热生成均匀热生成,基本变量TIME,X,Y,ZTIME,X,Y,Z,TEMPTIME,X,Y,Z,TEMP,VELOCITYTIME,X,Y,ZTIME,X,Y,Z,TEMPTIME,X,Y,Z,TEMPTIME,6-9,单元-特殊表格支持,Fluid1161-D热流单元载荷热流是基本变量TIME的函数压力是基本变量TIME,X,Y,Z的函数SURF151,152平面效果单元实常数旋转速度和滑动因子(主要用于汽轮机应用如转子)是基本变量TIME,X,Y,Z的函数,6-10,定义ANSYS表格,表格可以使用下面的方法定义:ANSYS命令(*DIM)在GUI,使用UtilityMenuParametersArrayParametersDefine/Edit在GUI,在施加边界条件时，软件允许使用现有的表格或定义新的表格。,用外部编辑器如EXCEL生成的文件。数据必须是tab-delimited格式才能用ANSYS的*TREAD命令读入。,6-11,ANSYSAPDL数组复习,ANSYS有两种数组参数:数组类型定义离散的函数表格类型使用线性插值方法定义连续的函数。数组可以是一、二、三维的。两维数组有m行n列;行用行号I标志,从1到m列用列号J标志，从1到n。组成数组的数称为数组元素。每个数组元素用下标(I,J)标志，I是行号，J是列号。,6-12,ANSYSAPDL数组复习,对于ARRAY类型，所有数都按照下标存储和引用。所有下标为整数。一个典型的5x3数组如下。,6-13,ANSYSAPDL数组复习,例如，在ANSYS命令中使用数组元素NSTR(2,4)会将数字19.1代入，如果数组定义如下:,6-14,ANSYSAPDL表格复习,TABLE同样由一系列排列为矩阵的数组成，但TABLE允许使用数组元素的中间值，这是通过插值方法实现的。而且,J=0列,通常包括行号数字(I=1到m),必须变化为单向增加的数字系列(不一定是整数)。同样,I=0行,通常包括列号数字(J=1到n),必须变化为单向增加的数字系列(不一定是整数)。缺省情况下，所有标志数字从接近0的数字开始。,6-15,1-D表格例子,例如,1-D表格A可以解释如下:A(1.5)表示20.0(12.0和28.0的中间数值)A(2.5)表示87.2(28.0和146.4的中间数值),6-16,2-D表格例子,2-D表格PQ可以解释如下:PQ(1,1.5)表示3.5(2.8和4.2的中间数值)PQ(3.5,1.3)表示14.88(42.0和-4.5平均,9.7和2.0平均,然后乘上30%加到42.0和-4.5的平均值上)。,6-17,载荷数组的通常习惯,尽管参数名可以达到8个字符长，表格参数由于应用于命令的域中，只能限制在7个字符以下。,6-18,例:使用表格定义A“载荷”,在载荷中使用表格允许用户在表格元素之间自动进行线性，双线性和三次插值。考虑下面的载荷(FORCE可以是节点热流):,6-19,使用命令定义APDL表格,表格(FORCE)元素为载荷数值，行标号为时间数值，使用下面的ANSYS命令定义:*DIM,FORCE,TABLE,5,1,1,TIMEFORCE(1,1)=0,560,560,238.5,0FORCE(1,0)=1.E-6,0.8,7.2,8.5,9.3FORCE(0,1)=1.0(不必须)见右图所示。,6-20,使用GUI表格编辑器定义APDL表格,使用GUI,intheUtilityMenuParametersArrayParametersDefine/EditAdd,1,2,6-21,使用外部模板定义APDL表格(如,EXCEL),在EXCEL中定义表格元素，包括“0”号行,“0”号列。用户可以使用空行将数据平面分割开，数据会可读性更强。一个3-D表格的例子见下，注意y平面变量有两个元素:在EXCEL中用tab-delimited文本格式写出，使用菜单FileSaveAschangeSaveAsTypetotext(tabdelimited),6-22,使用外部模板定义APDL表格(续),用户可以给EXCEL数据加上题目，而且空白行可以使得数据更加容易识别。在ANSYS中,定义数组使用UtilityMenuParametersArrayParametersDefine/EditAdd,thisisa3-DTable,6-23,使用外部模板定义APDL表格(续),要读入EXCEL数据，使用UtilityMenuParametersArrayParametersReadfromFileTablearray并将对话框中所需数据给出。,注-跳过EXCEL数据的前四行。,6-24,在GUI中检查表格,数据可以通过UtilityMenuParametersDefine/Edit来检查,数据平面,6-25,表格可以用UtilityMenuPlotArrays绘制。使用前面的1-D例子，将对话框中填充，生成的图见后面。,绘制APDL表格,6-26,绘制APDL表格(续),填充完对话框并单击OK,下面的图形出现,6-27,使用命令操作表格,表格可以使用标准的ANSYS加载命令施加。使用表格输入的方式是在命令的数值域将表格名填进去。不用使用其他表格参数。下面有几个例子:SF,all,HFLUX,%fluxtab%D,all,TEMP,%temptab%SFL,all,CONV,%hftab%,1000,6-28,通过GUI施加操作表格,在GUI,当用户选择加载并单击APPLY,对话框可以让用户指出施加的方式:常数或现存的表格或新表格,6-29,检查通过表格输入的边界条件,当表格边界条件施加到实体模型或有限元模型上后，在前处理器中列出边界条件时表格名将出现。如果使用SBCTRAN(SolutionLoadsOperateTransfertoFE)命令将载荷传递到有限元模型上时,在求解器中可以列出或绘制出结点和单元的表格边界条件数值。如果表格边界条件是时间的函数，显示的是TIME=0的数值.。如果载荷步文件(jobname.SXX)是在传递以后写的，结点和单元的载荷步文件边界条件(实体模型载荷不存储在LS文件中)还是由表格引用，而不是数值。,6-30,检查通过表格输入的边界条件(续),求解结束后，列出或绘制单元和结点的边界条件(在求解器中)会显示出时间为最后载荷子步的数值。在POST1(通用后处理器)中,边界条件可以在结点和单元上列出或绘制，时间是存储在数据库中的数值。,6-31,其他选项-表格输入,在许多情况下，载荷类型可以得到但载荷的表格数值在建模时并不知道。在这种情况下，用户可能只想确定不同的载荷区域。要做到这一点:使用新表格施加边界条件，表格大小和数值都取缺省数值(空白)。如果后来知道了边界条件，只要重新定义表格大小和数值即可。不用重新施加载荷。,6-32,例子：换热系数是长度的函数HF=f(x),假设面上的对流换热系数是非线性函数，并只与面的全局X坐标有关。假设有6个插值点，可以在GUI中定义如下表格:,6-33,例子：换热系数是长度的函数HF=f(x)(续),假设全局笛卡儿坐标从1到6变化，换热系数表格可以定义为:,6-34,例子：换热系数是长度的函数HF=f(x)(续),使用GUI在线上施加对流，出现下面的对话框。选择“existingtable”为换热系数，“constantvalue”为介质温度(因为没有用表格定义)。输入TBULK数值并单击OK。,6-35,例子：换热系数是长度的函数HF=f(x)(续),下面的对话框允许用户选择合适的表格(在本题中,只有一个表格,HF,可以选择)。,6-36,例子：换热系数是长度的函数HF=f(x)(续),要显示对流边界条件，我们打开右图选项UtilityMenuunderPlotCtlsSymbols,6-37,例子：换热系数是长度的函数HF=f(x)(续),要在绘图时显示表格名，在UtilityMenu,underPlotCtlsNumbering,打开TableNames,6-38,例子：换热系数是长度的函数HF=f(x)(续),画线显示表格化换热系数,HF,使用表格名。,6-39,例子：换热系数是长度的函数HF=f(x)(续),在通用后处理器(POST1)中,数值可以通过将表格名关闭并将数值打开显示(UtilityMenuPlotCtlsNumbering),Turnoff,Turnon,6-40,例子：换热系数是长度的函数HF=f(x)(续),数值不能在几何模型显示;只能在有限元模型显示:,这是单元绘图。,6-41,时间步长和输出控制,表格边界条件允许用户在一个载荷步中定义和求解复杂，瞬态热传递问题。多载荷步同样可以使用表格化边界条件。但是，当自动时间步(ATS)打开，求解器可能会跳过时间的关键部分(如,载荷突然变化)。要避免这一点，用户可以定义关键时间(时间步中的某一时间将设置为初始的NSUB/DELTIM数值)。这些设置可以通过ANSYS命令TSRES来设置。结果写入结果文件(jobname.rth)的频率可以使用OUTRES命令控制。,6-42,关键时间和输出时间,RPM,Time,输出时间,关键时间,6-43,时间步重置(TSRES命令),在特定时间点重新设置时间步大小为初始数值(使用NSUB/DELTIM命令设置)可以参照数组(不是表格)，使用TSRES命令:TSRES,数组名其中的数组必须是ANSYS一维数组Nx1,N为关键时间的数目，数组元素就是相应的要重新设置的时间步的时间数值。数组指定的时间点之间的时间增量必须比NSUB或DELTIM指定的时间初始数值大。,6-44,时间步重置-使用GUI,3,1,2,4,6-45,结果文件输出控制(OUTRES),OUTRES命令允许用户控制计算结果输出到热结果文件(jobname.rth)中的频率。命令格式为,OUTRES,item,FREQ,Cname要使用数组参数(不是表格)控制输出频率，在命令的