与坐标或时间有关的边值问题(适用于FEPG单机版和pFEPG版)20140410.docVIP

与坐标或时间有关的边值问题 一、适用于pPFEG版本（元计算刘钊老师总结） 修改gcn文件 在待求解的场之前，添加bft元件程序，后面跟场名，表示要修改该场的边 值。例如，对于一个瞬态温度场问题，场名为a，gcn文件如下 defi a par startsin a if exist stop del stop :1 bft a solvsin a post if not exist stop goto 1 准备boundbft.for文件 其中定义边值关于坐标和时间的函数，名为bound，例如 real*8 function bound r,t,j,nf implicit real*8 a-h,o-z dimension r 2 x r 1 y r 2 if abs x .lt.1.0e-8 then bound t*t+x*x+y*y endif return end 其中nf为场的编号，j表示该场的第j个未知量，r为坐标，t为时间 生成程序 前处理 在GiD中将约束的边界规格数赋为-1，约束的值随便写，后面计算的时候bft会进行修改 二、适用于FEPG6.0单机版（在元计算刘钊老师的基础上和指导下，由中科院南海所杨小秋总结） 具体操作过程： 1. 在gcn文件中的bft命令行中添加a场，也即运行bft.exe程序时，自动根据bft中的相关bound函数对a场的边界进行约束。 defi a par b str a START%1 a if exist stop del stop :1 bft a SOLV%1 a stress b post a if not exist stop goto 1 2. 运行gcn.bat lps, 重新生成所有程序 3. 进入bft.for源程序，根据实际情况修改相关bound函数，并在其相应位置调用，例如在E:\FEPGsy\Tt-bft测试算例中对bound函数修改成： real*8 function bound r,t,j,nf implicit real*8 a-h,o-z dimension r 2 c bound t*t+r 1 *r 1 +r 2 *r 2 bound sin t*2*3.1415926/36000.0 c write *,* bound ,bound return end 4. 进入gid前处理，根据实际情况，对温度随时间变化的边界进行设置，约束的边界规格数赋为-1，约束的值随便写，后面计算的时候bft会进行修改。 5. 运行实例见E:\FEPGsy\Tt-bft，计算结果符合预想情况，具体如下： 计算100个小时后的温度场分布 1-50 h（间隔5 h），温度-深度剖面，与解析解得到的一致。 0-5 m深度处（间隔0.5 m） 6. 将实测数据加载到温度边界E:\FEPGsy\PLHF，成功！ 主要在上面的基础上，修改了bft.for中的bound函数 real*8 function bound r,t,j,nf implicit real*8 a-h,o-z dimension r 2 c bound t*t+r 1 *r 1 +r 2 *r 2 c bound sin t*2*3.1415926/36000.0 c write *,* bound ,bound OPEN 1,FILE PLHF.txt,STATUS old tnum 0.0 temp1 0.0 1000 read 1,* tnum,temp1 if tnum*20.0*60.0 .eq.t then bound temp1 goto 1001 else goto 1000 1001 endif CLOSE 1 return end 需要换成直接读取文件的语句，效率会提高很多。 但暂时还受999步的限制，有待修改源程序，实现实际采样步数。