【文件属性】：
文件名称：用matlab生成谐波代码-fem_1:fem_1
文件大小：287KB
文件格式：ZIP
更新时间：2021-05-21 17:47:25
系统开源
用matlab生成谐波代码fem_1固定码（正在进行中）
在这里，我将我的上一份工作（于2017年12月1日发布）发布在有限元方法上。
我的目标是观察亚音速下一对（机翼，气流）的行为。
出于仿真目的，我将考虑到机翼在气流绕动时是固定的。
这与自然现象完全相反，但它是模型的合理首次尝试。
为了能够继续前进，我们必须深入一点物理方面的考虑，我不会详细介绍，但我们必须牢记以下几点：
该模拟基于不可压缩的Navier-Stokes方程，但解决方案简化为谐波函数（位移）
我们认为机翼位于方盒中。
上边界和下边界被认为在物理上远离机翼，左边界是气流的来源，其轨迹是水平的并且与高度成反比。
右边界迫使解与左边界相同
说到解决方案，我们在这里计算流线函数
关于机翼，我正在研究可提供的NACA0012机翼
我使用了一系列的fortran工具，这些工具最初是在我的培训中提供的，并且基于一些fortran库和GMSH
meshtool。
代码中提到了所有积分
我修改了代码的某些部分并整理了目录。
它仍然不像我想要的那样模块化，因为一方面模块管理（库+其他一些..）对我来说有些混乱，另一方面，代码仍需修改以处
【文件预览】：
fem_1-master
----libraries()
--------tunesolver.mod(658B)
--------riters.f90(5KB)
--------matcsr2.o(74KB)
--------gestionerr.mod(415B)
--------precond.mod(3KB)
--------f90_kind.mod(446B)
--------riter2.f90(3KB)
--------syslincsr.mod(3KB)
--------gestionerr.f90(638B)
--------resol.o(11KB)
--------iters.o(116KB)
--------initsolver.f90(5KB)
--------tri.f90(7KB)
--------blas1.f(16KB)
--------iters_utils.mod(1KB)
--------iters.f(102KB)
--------libsparsekit.a(335KB)
--------gestionerr.o(3KB)
--------makefile(2KB)
--------blas1.mod(1KB)
--------blassm.o(5KB)
--------blas1.o(13KB)
--------iters.mod(4KB)
--------f90_kind.o(952B)
--------f90_kind.f90(2KB)
--------initsolver.o(14KB)
--------riters_ver3.f90(2KB)
--------matcsr.f90(39KB)
--------resol.f90(4KB)
--------matcsr2.f90(39KB)
--------tri_mod.mod(600B)
--------riters_ver2.f90(3KB)
--------initsolver.mod(992B)
--------riters_ver1.f90(5KB)
--------blassm.f(6KB)
--------tri.o(8KB)
--------resolrun.mod(939B)
--------blass.mod(729B)
----mod_lec_fic.f90(13KB)
----README.md(3KB)
----makefile(572B)
----pw2mts.f90(3KB)
----EFLib.f90(11KB)
----streamlines.PNG(65KB)
----wing_naca.geo(967B)
----load_gmsh.m(7KB)
----trace_naca.m(506B)