专利名称：定向凝固叶片生产的自适应变速抽拉仿真方法
技术领域：
本发明涉及一种定向凝固叶片生产的自适应变速抽拉的方法，特别涉及一种采用数值模拟技术制定定向凝固叶片变速抽拉工艺的速度曲线的方法。
背景技术：
定向凝固叶片(也成定向叶片)在航空和民用领域都有重要应用，常被用来作为飞机发动机的导向叶片、涡轮叶片或者民用燃气机的叶片。定向叶片的制造技术也备受各国材料加工方向专家所重视。目前，仅从速度工艺角度来分，定向凝固叶片生产分为定速工艺和变速工艺两种，以下是当前两种工艺的特点对比
权利要求
1.定向凝固叶片生产用的自适应变速抽拉仿真方法，其特征在于，依次含有以下步骤步骤(1)，构建一个基于数值模拟方法的计算机-实验装置系统 所述实验装置，包括定向凝固炉及钨/铼热电偶，其中定向凝固炉，炉体简化为用挡板隔离的位于上部的加热区和位于下部的冷却区，在所述的冷却区的抽拉机构由圆形结晶器和抽拉推杆连接组成，所述圆形结晶器上装有浇注用的叶片型壳，在所述抽拉机构作用下，所述叶片型壳穿过所述挡板上的孔，能够在所述抽拉推杆作用下，在所述加热区和冷却区之间做上下往复运动，在所述叶片型壳的空腔内对应叶片的引晶段、叶身、上部缘板、下部缘板以及榫头处设定测温点，至少5个点，所述挡板水平地连接在定向凝固炉内炉壁两侧的中部位置，钨/铼热电偶，在测定所述叶片型壳测温点的温度时，该钨/铼热电偶的输出端向所述计算机输入对应测温点的温度，所述计算机，预置有FT-Star软件；步骤0)，依次按以下步骤进行定向凝固叶片生产用的自适应变速抽拉仿真过程 步骤(2. 1)，操作员向所述计算机输入叶片生产所用材料的模拟用的基本参数，其中包括热传导系数，辐射换热系数，所用合金的物理参数比热、潜热、密度和固/液相线温度， 定向凝固炉体参数所述加热区直径和高度、所述冷却区直径和高度、挡板厚度和圆盘结晶器直径，同时，输入固定抽拉速度值3mm/min，步骤(2.幻，操作员向所述计算机输入所述叶片型壳的三维简化模型，所述计算机用所述FTItar软件对该三维简化模型用正六面体单元进行三维离散化，步骤(2. 3)，操作员向所述计算机下达用所述FT-Star软件对叶片凝固过程的温度场进行测试模拟的指令，并输出步骤⑴所述至少5个测温点的温度值随时间变化的数据，步骤(2. 4)，操作员以所述叶片型壳为样片，以3mm/min为固定抽拉速度值，在相同于步骤(2. 1)中的模拟用的基本参数条件下，在所述实验装置中进行浇注实验，实时测定各测温点实际的温度随时间变化的数据，步骤(2.幻，操作员把步骤(2. 3)数值模拟得到在同一时刻条件下的所述各测温点的数据与步骤(2.4)实测的所述各测温点的温度数据进行对比分析，在任意相同时刻，只要有一个测温点的相对温度误差δ > 10%，便需要增加或减小上一次确定的热传导系数和辐射换热系数值10%，并作为新的所述基本参数输入到所述FT-Star软件，直到在不同时刻下所有所述测温点的相对温度误差都不大于10%为止，从而确定最终用于模拟计算的基本参数，所述相对温度误差是指温度模拟值与温度实测值之差的绝对值与所述温度实测值之比，步骤(2. 6)，按照以下步骤对所述定向凝固叶片生产用的自适应变速抽拉过程进行预测模拟步骤(2. 6. 0)，建立实际生产用叶片型壳的三维模型，并采用FTItar软件进行三维离散化，步骤(2. 6. 1)，确定抽拉速度曲线ν (t)，v(t)表示t时刻的抽拉速度值 步骤(2. 6. 1. 1)，确定变速判据所述变速判据用来决定何时对抽拉速度值进行改变， 其包括两方面判据减速判据和加速判据，定义叶片的离散单元C(i，j，k)的热流方向与x，y，ζ轴的夹角分别为a(i,j,k)，日⑴」』，￥(^10，其温度梯度大小为6^10，1，」，1^是对应坐标轴下的坐标值，定义糊状区为凝固过程中固相温度线和液相温度线之间的区域，1)减速判据，对糊状区离散单元单元C、(i，j，k)，γ 、(i,j,k) > α 、(i,j,k)或者Y 、(i,j,k) > β ^i.j.k)，同时G ^^k) <2°C/mm，对于凝固过程的糊状区，定义当满足上述条件的单元数η比糊状区的总离散单元数m的值为λ，当λ > 0. 0509时满足判据；2)加速判据，糊状区所有单元Y 、(i,j,k) < 0.785且G、(ijJ,k) >3°C/mm，为满足加速判据。 步骤(2. 6. 1. 2)，确定截面积&和最大抽拉速度Vmax的函数关系， 步骤(2. 6. 1. 2. 1)，建立第一组5个不同截面的棒状试样模型，所有棒状试样模型长度300mm，作为第一组标准算例，记为标准模型系列1，S0 = 78. 5，314. 2，1256. 6，1963. 5， 2827. 4，单位mm2，其中，S0为截面积值系列，步骤(2. 6. 1. 2. 2)，计算所述标准试样的最大抽拉速度Vmax 对所述的5个不同截面的棒状试样模型采用所述FT-Star软件进行三维离散化，并以所述确定的模拟用基本参数为输入量，对所述的5个不同截面的棒状试样模型分别进行模拟计算，并以步骤(2. 6. 1. 1)所述的变速判据为依据，当计算数据满足所述变速判据时，对抽拉速度进行改变，每次减小或增加的速度增量绝对值为lmm/min，不断进行判断，直至确定的抽拉速度不满足所述变速判据，且凝固完全结束，此时确定的抽拉速度值为所述最大抽拉速度值Vmax，从而建立所述标准模型系列1的5个所述试样模型对应的截面积值系列&和最大抽拉速度Vmax的关系对 tS0(n)~Vmax(n)], η = 1,2,3,4,5,步骤(2. 6. 1.2. 3)，根据步骤(2. 6. 1.2. 2)所述的关系对[Stlfc)-Vmax(n)]，采用BP神经网络算法，以&为输入层变量，以Vmax为输出层变量，根据所述关系对，对所述第一组标准算例进行系统训练，得到输入量与输出量之间内联关系，即 F(S0) = [VfflaJ其中，F(S0)表示根据第一组标准算例的5个棒状试样模型通过BP神经网络算法确定的统一函数关系，步骤O. 6. 1.3)，确定突变截面[S1 S2]和最大抽拉速度Vmax以及变速时间提前量At 之间的函数关系步骤(2. 6. 1. 3. 1)，建立第二组5个截面突变的棒状试样模型，作为第二组标准算例， 记为标准模型系列 2，[S1 S2] =314.2 706. 9,314. 2 1256. 6,314. 2 1963.5， 706.9 1963.5,1256.6 2827. 4，单位mm2 mm2，所有棒状试样模型的长度为300mm，截面开始突变位置ζ = 150mm，其中，S1为截面变化前的截面积值，其ζ方向坐标为Z1, S2为截面变化后的截面积值，其ζ方向坐标为^， 且满足，Z1 < ζ < z2, (Z2-Z1)A < 5%,步骤(2. 6. 1. 3.幻，对应于所述第二组标准算例中的每一个棒状试样模型，确定最大抽拉速度随时间变化曲线，办法如下步骤(2.6. 1.3. 2. 1)，采用步骤(2.6. 1.2. 3)所述的最大速度公式F (Stl) = [Vmax]，当截面积变化时，抽拉速度根据所述函数关系，相应变化，步骤(2. 6. 1. 3. 2. 2)，确定最大抽拉速度与时间关系，根据截面开始突变位置ζ = 150mm，确定变速点时间为、=z/vi; t2 = z/v2，其中V1为截面突变前面积为S1所对应的最大抽拉速度值，V2为截面突变后面积为&所对应的最大抽拉速度值，从而确定的抽拉速度曲线为，以 V1速度抽拉，持续时间为、，再改变速度值，开始以V2速度进行抽拉，持续时间为t2，步骤(2. 6. 1.3.幻，根据计算模拟，修正第二组标准算例的最大抽拉速度随时间变化关系以步骤(2. 6. 1. 3. 2. 2)所确定的最大抽拉速度随时间变化曲线为输入，采用所述 FT-Star软件进行模拟计算，并结合所述变速判据进行判断，当满足变速判据时，并不改变速度值，而是减小变速点时间、的大小，每次减小量为0. 5min，并再次计算，直至凝固结束，累加所有时间减小值，记为At，最终确定第二组标准算例中每个模型所对应的关系对 tSi(n) · S2(n)-v1(n) . v2(n) . Δ t(n)] ,η = 1,2,3,4,5,步骤(2. 6. 1. 3. 4)，根据步骤(2. 6. 1. 3. 3)所述的关系对，采用BP神经网络算法，以S1, &为输入层变量，以Vl，v2, At为输出层变量，对所述第二标准算例进行系统训练，得到输入量与输出量之间内联关系，即 F (S1, S2) = [V1, V2, Δ t]其中，F(SnS2)表示根据第二组标准算例的5个棒状试样模型通过BP神经网络算法确定的统一函数关系，步骤(2. 6. 1.4)，根据所述实际叶片型壳的三维模型，分析不同截面的截面积值，并根据步骤(2. 6. 1. 3. 4)确定的截面积、速度、时间函数关系，计算所需要的抽拉速度随时间变化曲线，步骤(2. 6. 2)，优化凝固抽拉速度曲线ν (t)以步骤(2. 5)所确定的模拟用基本参数为输入，以步骤(2. 6. 0)所建立的三维离散化模型为输入，以步骤(2. 6. 1)所建立的抽拉速度随时间变化曲线为输入，采用FT-Mar软件进行模拟计算，并采用变速判据进行变速判断， 当满足判据时变速方法如下对抽拉速度进行自动优化，调整其速度改变量和速度改变时间点，变速规则用来确定每次变速的速度改变量，变速规则分为两部分，减速规则和增速规则 A减速规则当满足减速判据时，减速量Δ、= (I-LJvm, Lm= (λ-0.0509)/0.0509减速时间位置程序回退时间量At = tprt-tbf, 其中，λ，当满足减速判据的步骤(2. 6. 1. 1)所述糊状区的单元数η与所述糊状区的总单元数m的比值t_，当前凝固时刻，tbf，当前糊状区处于固相线的单元的温度值从开始浇注温度降至液相线温度时所对应的时刻，满足tbf < tPrt，vffl为tbf时刻对应的速度值， Lm，速度减小调整系数， B增速规则当满足加速判据时，增速量 Δνρ = (1+Lp)vp, Lp = (0. 785-γ 、(ijJ,k))/0. 785，增速时间位置当前时刻增速，其中，Y 、(i,j,k)，在步骤(2.6. 1. 1)中定义，Vp，为当前时刻速度， Lp，速度增大调整系数，步骤(2. 6. 3)，当根据变速判据对计算结果验证通过，满足全叶片凝固结束时，根据速度对时间绘制曲线，得到v(t)凝固抽拉速度曲线。
全文摘要
定向凝固叶片生产的自适应变速抽拉仿真方法，属于材料加工技术领域，其特征在于，通过对Bridgeman法定向凝固生产涡轮发动机叶片的过程进行模拟仿真，集中模拟金属液浇铸后，抽拉机构运转过程中，铸件凝固过程的温度场演变过程，同时本方法以凝固过程的温度场模拟结果为指导，结合工业对定向凝固叶片生产的标准要求，采用神经网络算法和PID控制原理确定凝固过程的抽拉速度的改变时间和改变大小，对全叶片凝固过程的抽拉速度工艺进行优化。本方法提高了铸件生产合格率，提高了生产效率，降低了工艺调试周期，降低了生产成本，将定向凝固生产的单级变速抽拉工艺升级为多级变速抽拉工艺，具有广阔的应用前景。
文档编号G06F17/50GK102508965SQ20111034369
公开日2012年6月20日 申请日期2011年11月3日 优先权日2011年11月3日
发明者唐宁, 张航, 柳百成, 许庆彦 申请人:清华大学