数值模拟程序.docx

《数值模拟程序.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、油藏数值模拟程序设计报告班级：油工12-3班姓名：孙雷学号：2程序一：1、程序内容求解三对角方程组开始输入a,b,c,d值输入主对角线个数令计算令输出值结束图1程序一程序框图i=2,3…ni=n-1,n-2,…,12、程序框图3、程序代码PrivateSubCommand1_Click()n=InputBox("输入主对角线个数")Dima(),b(),c(),d()Diml(),y(),u(),x()ReDima(n),d(n),b(n-1),c(n)Fori=1Tona(i)=Text1(i-1).Textd(i)=Text4(i-1).TextNextiFori=1T

2、on-1b(i)=Text2(i-1).TextNextiFori=2Tonc(i)=Text3(i-2).TextNextiReDiml(n),y(n),u(n)l(1)=a(1)y(1)=d(1)/a(1)u(1)=b(1)/a(1)Fori=2Ton-1l(i)=a(i)-c(i)*u(i-1)u(i)=b(i)/l(i)y(i)=(d(i)-c(i)*y(i-1))/l(i)Nextil(n)=a(n)-c(n)*u(n-1)y(n)=(d(n)-c(n)*y(n-1))/l(n)ReDimx(n)x(n)=y(n)Fori=n-1To1Step-1x(i)=y(i

3、)-u(i)*x(i+1)NextiFori=1TonPicture1.Print"x(";i;")=";Format(x(i),0)NextiEndSub4、计算结果图2程序一运行界面程序运行结果为：程序二：1、程序内容已知一维均质等厚单相渗流，已知边界条件定压、初始压力分布，求任意时刻的压力分布，渗流微分方程为：开始输入已知的三个边界条件输入值输入和系数矩阵和右端项利用追赶法求解压力方程组打印结束是图3程序二程序框图2、程序框图3、程序代码Dima(),b(),c(),d()Diml(),y(),u(),x()DimnSubjs()ReDima(n),d(n),b(n-

4、1),c(n)Fori=1Tona(i)=Text1(i-1).Textd(i)=Text4(i-1).TextNextiFori=1Ton-1b(i)=Text2(i-1).TextNextiFori=2Tonc(i)=Text3(i-2).TextNextiReDiml(n),y(n),u(n)l(1)=a(1)y(1)=d(1)/a(1)u(1)=b(1)/a(1)Fori=2Ton-1l(i)=a(i)-c(i)*u(i-1)u(i)=b(i)/l(i)y(i)=(d(i)-c(i)*y(i-1))/l(i)Nextil(n)=a(n)-c(n)*u(n-1)y(n

5、)=(d(n)-c(n)*y(n-1))/l(n)ReDimx(n)x(n)=y(n)Fori=n-1To1Step-1x(i)=y(i)-u(i)*x(i+1)NextiEndSubPrivateSubCommand1_Click()Text5=""DimdetDimchushi()n=Text9/Text10det=Text7/Text10^2Fori=0Ton-2Text3(i)=detText2(i)=detNextiFori=0Ton-1Text1(i)=-(1+2*det)NextiReDimchushi(n)Fori=1Tonchushi(i)=Text11N

6、extitmax=Val(Text8)dt=Val(Text7)t=0t=t+dtWhile(t

7、4、计算结果图4程序二运行界面本次计算得出十五个时刻的压力分布，具体结果如下表表一程序二计算结果时间p1p2p3p418.629.499.337.8327.658.348.466.5336.937.417.435.6246.396.726.634.9756.006.206.034.4965.705.825.594.1475.495.545.273.8985.335.335.033.7095.215.184.853.56105.125.064.723.46115.054.984.623.39125.014.924.5