[转帖]FEPG的二次开发，在瞬态传热中的应用

jiangfan2008

FEPG是一种优秀的有限元软件，由中科院研究员的伟大创举。她的优势就在于微分方程组、解法自定义，使用者有很大的灵活性。综合而言，和Ansys等大型有限元比较，毫不逊色。基于FEPG，完全可以开发大型的专用有效元程序。
然而，和当前流行的大多数有限元软件一样，FEPG采用fortran开发。fortran是一种老式的结构化程序设计语言，当前已经是过时的产品。要编写大程序，要采用面向对象的语言，如c++，vc++,c++builder,vb等.采用C++语言的另一个好处是可以使用专业公司如微软\autodesk,borland等公司提供的库和控件,简化程序的编制.
基于以上原因,本软件为用户提供了一个FEPG c++的转化程序.熟悉C语言的用户可以轻松的开发Fepg,快速开发大型的专用有限元系统.
//
//
使用法：在用FEPG生成了fortran源程序后,打开FEPG_C++_Translator.exe, 在其中输入FEPG工程目录,运行,软件将自动把所有的for程序转化为c,c++程序,并在FEPG工程中生成VC++工程文件main.dsw.
假如安装了VC++6.0以上,双击工程文件,就将打开该工程.
最新更新：05.09.03
修改了生成程序的结构
05.9.4
可模拟弹塑性变形
在user_functions.cpp中智能化增加函数。
http://www.fegensoft.com/phpwind/read.php?tid=1500&keyword=
程序编写者：华工的Eric博士。

jiangfan2008

1 打开FEPG---samples------heat---2dxy----nonliner_dynamic_heat, 拷贝你们的文件到工作目录，这里是d:\temp
2 用写字板打开里面的heat.pde文件，删除里面的子程序。子程序将用c语言写。
3 在dos或者fepg的run窗口中运行：gio np,或者test np。dos窗口显示程序正在生成，并提示缺少两个子程序eku 和ecu。
4打开本系统的FEPG_VC++Transfer.exe，在窗口中输入工程目录，这里是：d:\temp, 运行程序。假如一切正常，将出现dos窗口，提示程序的生成过程，并显示：vc++ 程序生成，其工程文件为：d:\temp\main.dsw. 关闭dos.
5 假如系统中已经安转有vc++,双击d:\temp\main.dsw,打开vc工程。其中显示了所有的c,c++文件。
6 双击dll项目中的user_function.cpp文件，在其中增加函数：
extern "C" __declspec( dllexport ) void eku_(double& tn,double& ek){
    AFX_MANAGE_STATE(AfxGetStaticModuleState());
    //Add your function code here...
    //Notice: in pde file, the function is eku, this is eku_
    ek=1.44+1.236-3*tn-4.63-7*tn*tn;
}
假如比热的计算采用样条函数插值，定义一样条类，定义见附件。CSpline。 将温度和比热的关系写入一文本文件，如附件中的ec.txt. 增加函数：
CSpline sp_ec("ec.txt");
extern "C" __declspec( dllexport ) void ecu_(double& tn,double& ec){
    AFX_MANAGE_STATE(AfxGetStaticModuleState());   
    //Add your function code here...
    //Notice: in pde file, the function is ecu,here is ecu_
    if(tn>1000) {
        ec=-0.44-5-0.95-5*tn+0.11-7*tn*tn;
        return;
    }
    if(tn<100) tn=100;
    ec=sp_ec.get_y(tn);
}
7 编译，运行程序。run
8 在fepg中显示处理结果

jiangfan2008

程序的运行过程

hp270388026

请大家注意软件的更新。
http://www.fegensoft.com/phpwind/read.php?tid=1500