边界元法程序
姚振汉 王海涛 编著
高等教育出版社 高等教育电子音像出版社
《边界元法》程序光盘目录
第一部分：弹性力学平面问题边界元程序BEM2 1. FORTRAN 程序BIEBE2（简称BEM2）简介 2. F90源程序
BEM2.F90 BEM2PRP.F90 BEMPN.F90 BEM2POP.F90 BEM2PLT.F90
3. 执行文件及输入数据文件
BEM2.EXE CGSDAT.DAT TEQPN.DAT Input.dat 4. 使用说明
BEM2主要变量说明 输入文件模板 5. 考题与算例
BEM2计算例题汇总
(1) 第一组（方形域简单考题） 第一组算例
159MB
132KB 约1530行 约880行 约880行 约390行 约590行 533KB 529KB 2KB 2KB （简例） 580KB
10KB 94.8MB
4.38MB
输出文件 TecPlot 输出文件 图形文件 中间文件
第二组（薄板梁纯弯）
第二组算例 输入文件 输出文件 TecPlot 输出文件 图形文件 中间文件
第三组（悬臂薄板梁弯曲）
第三组算例 输入文件 输出文件 TecPlot 输出文件 图形文件
第四组（直杆与刚性基础垂直接触）第四组算例 输入文件 输出文件 TecPlot 输出文件
9个 36个 16个 3个 8.29MB
10个 10个 40个 36个 1个 7.93MB
11个 11个 44个 43个 10.8MB
11个 11个 44个
(2)(3)(4)
中间文件
(5) 第五组（方板中心圆孔应力集中）
第五组算例 输入文件 输出文件 TecPlot 输出文件 图形文件
(6) 第六组（圆形边界内部和外部问题）
第六组算例 输入文件 输出文件 TecPlot 输出文件 图形文件
(7) 第七组（链状子域法及其在无滑动摩擦接触中的应用）第七组算例 输入文件 输出文件 TecPlot 输出文件 图形文件
(8) 第八组（相同材料和不同材料链状子域计算）
第八组算例 输入文件 输出文件
1个 12.4MB
10个 10个 40个 50个 8.65MB
10个 10个 40个 40个 9.60MB
10个 10个 40个 54个 13.2MB
12个 12个
TecPlot 输出文件 图形文件
(9) 第九组（较大规模算例和二次等参边界段算例等）
第九组算例 输入文件 输出文件 TecPlot 输出文件 图形文件 中间文件
(10)第十组（内点变量计算）
第十组算例 输入文件 输出文件 TecPlot 输出文件 中间文件
6. 边界应力图和域内应力云图比较
边界应力图和域内应力分布云图比较 (1) 算例1（带一中心圆孔方板单向拉伸）
算例1 输入文件 输出文件 TecPlot 输出文件 图形文件
48个 51个 17.8MB
12个 11个 47个 49个 1个 1.21MB
10个 10个 20个 5个 62.1MB
8.31MB
1个 1个 4个 5个
(2) 算例2（带5个圆孔方板单向拉伸）
算例2 输入文件 输出文件 TecPlot 输出文件 图形文件
(3) 算例3（薄板梁三点弯曲）
算例3 输入文件 输出文件 TecPlot 输出文件 图形文件
(4) 算例4（薄板梁四点弯曲）
算例4 输入文件 输出文件 TecPlot 输出文件 图形文件
(5) 算例5（圆盘对径压缩）
算例5 输入文件 输出文件
8.78MB
1个 1个 4个 5个 8.28MB
1个 1个 4个 5个 3.66MB
1个 1个 4个 5个 6.57MB
1个 1个
TecPlot 输出文件 图形文件
(6) 算例6（含100个随机分布圆孔方板单向拉伸）
算例6 输入文件 输出文件 中间文件 TecPlot 输出文件 图形文件
第二部分：二维弹性常值边界元程序BEM2C++ 1. C++源程序和头文件
main.CPP mat_vec.CPP solver.CPP complex.H mat_vec.H solver.H
2. 执行文件和输入输出文件
BEM2C++.EXE BEM_INP BEM_OUT 3. 数值算例
4个 5个 30.4MB
3个 2个 2个 11个 10个 2.38MB 29.3KB 231行 114行 271行 220行 96行 110行 48.5KB 48.0KB (简例) (简例) 2.30MB
BEM2C++计算例题汇总
(1) TestE001（方板四边给定均匀法向位移）
(2) TestE002（方形域平面应变四边给定均匀法向位移） (3) TestE003（含中心圆孔方板四边给定均匀法向位移） (4) TestE004（含中心圆孔方板单向拉伸取四分之一结构计算）
532KB 532KB 580KB 520KB
每个算例有输入、输出文件。此程序虽然只有常值直线单元，但结构简明，可供学习边界元编程的读者参考。
第三部分：三维位势问题快速多极边界元程序FMBEM_LAPLACE3D 1. 三维Laplace 方程快速多极边界元软件FMBEM_LAPLACE3D用户手册 2. 执行文件
FMBEM_LAPLACE3D.EXE PAT_TRANSLATOR.EXE 3. 数值算例
FMBEM_LAPLACE3D计算例题汇总
(1) 第一组（长方体区域稳态温度场计算，使用常规求解器）
TestE001 TestE002 TestE003
(2) 第二组（长方体区域稳态温度场计算，使用快速求解器）
TestE004 TestE005 TestE006 TestE007
121MB
432KB 236KB 196KB 120MB
2.38MB 803KB 817KB 821KB 6.39MB 806KB 821KB 819KB 815KB
TestE008 TestE009 TestE010 TestE011
(3) 第三组（核电站某结构件稳态温度场计算，使用快速求解器）
TestE012 TestE013 TestE014 TestE015 TestE016
816KB 823KB 822KB 824KB 89.2MB 1.15MB 3.05MB 16.7MB 34.1MB 34.1MB
(4) 第四组（理想不可压无旋流场圆球扰流速度势计算，用快速求解器）13.6MB
TestE017 TestE018
(5) 第五组（阴极保护静电场计算，用快速求解器）
TestE019 TestE020
1.37MB 12.3MB 8.31MB 1.38MB 7.19MB
以上20个算例均含参数卡片文件BEM_DATACARD.DAT，前处理文件BEM_INP.DAT，输出数据文件BEM_OUT.DAT，图形输出文件TEC_OUT.DAT，以及计算日志文件BEM_LOG.DAT，为便于没有TecPlot 软件的读者查看输出图形，还给出了边界温度（或其它未知函数）分布云图OUT_T和相应的法向导数分布云图OUT_DTDN。
另外，在“执行文件备份”子目录有三个程序的执行文件和固定输入文件的备份，所
有文字文档也同时给出MSWord 文件和PDF 文件，以便读者的使用。在“彩色插图”子目录下，提供了书中部分插图的彩图电子文档，以弥补黑白印刷的不足，使读者能将这些插图看得更清楚些。
边界元法程序
姚振汉 王海涛 编著
高等教育出版社 高等教育电子音像出版社
《边界元法》程序光盘目录
第一部分：弹性力学平面问题边界元程序BEM2 1. FORTRAN 程序BIEBE2（简称BEM2）简介 2. F90源程序
BEM2.F90 BEM2PRP.F90 BEMPN.F90 BEM2POP.F90 BEM2PLT.F90
3. 执行文件及输入数据文件
BEM2.EXE CGSDAT.DAT TEQPN.DAT Input.dat 4. 使用说明
BEM2主要变量说明 输入文件模板 5. 考题与算例
BEM2计算例题汇总
(1) 第一组（方形域简单考题） 第一组算例
159MB
132KB 约1530行 约880行 约880行 约390行 约590行 533KB 529KB 2KB 2KB （简例） 580KB
10KB 94.8MB
4.38MB
输出文件 TecPlot 输出文件 图形文件 中间文件
第二组（薄板梁纯弯）
第二组算例 输入文件 输出文件 TecPlot 输出文件 图形文件 中间文件
第三组（悬臂薄板梁弯曲）
第三组算例 输入文件 输出文件 TecPlot 输出文件 图形文件
第四组（直杆与刚性基础垂直接触）第四组算例 输入文件 输出文件 TecPlot 输出文件
9个 36个 16个 3个 8.29MB
10个 10个 40个 36个 1个 7.93MB
11个 11个 44个 43个 10.8MB
11个 11个 44个
(2)(3)(4)
中间文件
(5) 第五组（方板中心圆孔应力集中）
第五组算例 输入文件 输出文件 TecPlot 输出文件 图形文件
(6) 第六组（圆形边界内部和外部问题）
第六组算例 输入文件 输出文件 TecPlot 输出文件 图形文件
(7) 第七组（链状子域法及其在无滑动摩擦接触中的应用）第七组算例 输入文件 输出文件 TecPlot 输出文件 图形文件
(8) 第八组（相同材料和不同材料链状子域计算）
第八组算例 输入文件 输出文件
1个 12.4MB
10个 10个 40个 50个 8.65MB
10个 10个 40个 40个 9.60MB
10个 10个 40个 54个 13.2MB
12个 12个
TecPlot 输出文件 图形文件
(9) 第九组（较大规模算例和二次等参边界段算例等）
第九组算例 输入文件 输出文件 TecPlot 输出文件 图形文件 中间文件
(10)第十组（内点变量计算）
第十组算例 输入文件 输出文件 TecPlot 输出文件 中间文件
6. 边界应力图和域内应力云图比较
边界应力图和域内应力分布云图比较 (1) 算例1（带一中心圆孔方板单向拉伸）
算例1 输入文件 输出文件 TecPlot 输出文件 图形文件
48个 51个 17.8MB
12个 11个 47个 49个 1个 1.21MB
10个 10个 20个 5个 62.1MB
8.31MB
1个 1个 4个 5个
(2) 算例2（带5个圆孔方板单向拉伸）
算例2 输入文件 输出文件 TecPlot 输出文件 图形文件
(3) 算例3（薄板梁三点弯曲）
算例3 输入文件 输出文件 TecPlot 输出文件 图形文件
(4) 算例4（薄板梁四点弯曲）
算例4 输入文件 输出文件 TecPlot 输出文件 图形文件
(5) 算例5（圆盘对径压缩）
算例5 输入文件 输出文件
8.78MB
1个 1个 4个 5个 8.28MB
1个 1个 4个 5个 3.66MB
1个 1个 4个 5个 6.57MB
1个 1个
TecPlot 输出文件 图形文件
(6) 算例6（含100个随机分布圆孔方板单向拉伸）
算例6 输入文件 输出文件 中间文件 TecPlot 输出文件 图形文件
第二部分：二维弹性常值边界元程序BEM2C++ 1. C++源程序和头文件
main.CPP mat_vec.CPP solver.CPP complex.H mat_vec.H solver.H
2. 执行文件和输入输出文件
BEM2C++.EXE BEM_INP BEM_OUT 3. 数值算例
4个 5个 30.4MB
3个 2个 2个 11个 10个 2.38MB 29.3KB 231行 114行 271行 220行 96行 110行 48.5KB 48.0KB (简例) (简例) 2.30MB
BEM2C++计算例题汇总
(1) TestE001（方板四边给定均匀法向位移）
(2) TestE002（方形域平面应变四边给定均匀法向位移） (3) TestE003（含中心圆孔方板四边给定均匀法向位移） (4) TestE004（含中心圆孔方板单向拉伸取四分之一结构计算）
532KB 532KB 580KB 520KB
每个算例有输入、输出文件。此程序虽然只有常值直线单元，但结构简明，可供学习边界元编程的读者参考。
第三部分：三维位势问题快速多极边界元程序FMBEM_LAPLACE3D 1. 三维Laplace 方程快速多极边界元软件FMBEM_LAPLACE3D用户手册 2. 执行文件
FMBEM_LAPLACE3D.EXE PAT_TRANSLATOR.EXE 3. 数值算例
FMBEM_LAPLACE3D计算例题汇总
(1) 第一组（长方体区域稳态温度场计算，使用常规求解器）
TestE001 TestE002 TestE003
(2) 第二组（长方体区域稳态温度场计算，使用快速求解器）
TestE004 TestE005 TestE006 TestE007
121MB
432KB 236KB 196KB 120MB
2.38MB 803KB 817KB 821KB 6.39MB 806KB 821KB 819KB 815KB
TestE008 TestE009 TestE010 TestE011
(3) 第三组（核电站某结构件稳态温度场计算，使用快速求解器）
TestE012 TestE013 TestE014 TestE015 TestE016
816KB 823KB 822KB 824KB 89.2MB 1.15MB 3.05MB 16.7MB 34.1MB 34.1MB
(4) 第四组（理想不可压无旋流场圆球扰流速度势计算，用快速求解器）13.6MB
TestE017 TestE018
(5) 第五组（阴极保护静电场计算，用快速求解器）
TestE019 TestE020
1.37MB 12.3MB 8.31MB 1.38MB 7.19MB
以上20个算例均含参数卡片文件BEM_DATACARD.DAT，前处理文件BEM_INP.DAT，输出数据文件BEM_OUT.DAT，图形输出文件TEC_OUT.DAT，以及计算日志文件BEM_LOG.DAT，为便于没有TecPlot 软件的读者查看输出图形，还给出了边界温度（或其它未知函数）分布云图OUT_T和相应的法向导数分布云图OUT_DTDN。
另外，在“执行文件备份”子目录有三个程序的执行文件和固定输入文件的备份，所
有文字文档也同时给出MSWord 文件和PDF 文件，以便读者的使用。在“彩色插图”子目录下，提供了书中部分插图的彩图电子文档，以弥补黑白印刷的不足，使读者能将这些插图看得更清楚些。