内容概要：本文针对黄平《润滑数值计算方法》一书中随机粗糙线接触弹流问题的Fortran代码存在的语法与数值计算错误进行修正，重点解决了数组越界、迭代收敛条件不合理等问题，并通过Matlab实现数据可视化。修正后压力分布与油膜厚度更符合物理实际，揭示了经典教材代码在实际运行中可能存在的严重偏差。
适合人群：从事润滑理论、弹流润滑数值模拟、Fortran编程或机械工程相关研究的研究生及科研人员。
使用场景及目标：①学习弹流润滑仿真中Fortran代码的常见错误及调试方法；②掌握Fortran与Matlab联合进行数值计算与可视化的技术路径；③理解数值稳定性与收敛条件对仿真结果的影响。
阅读建议：在复现代码时应重点关注数组维度设置与文件读写格式，建议结合hexdump等工具验证二进制输出的正确性，同时采用全局误差判断提升迭代稳定性。
2025-09-22 21:12:47
467KB
1
利用迭代的思想求解线性方程组，适用于齐次和非齐次的情况，并附上了运行测试代码，能够支持三种情况的方程组求解，并与正确结果进行比对，基本一致，在小数点后位有些许差异。
2022-11-27 19:24:57
2KB
1
Fortran代码分析工具-understand for fortran 1.4 (含注册机)
2021-12-23 16:05:30
16.98MB
1
拉普拉斯2D
用于在矩形网格上以二维方式求解Laplace（或Poisson）方程的fortran代码。
该代码用在域边界上定义的u（x，y）的值求解方程u_ {xx} + u_ {yy} = f（x，y）。
2021-12-11 18:23:48
3.7MB
1
该程序按照http://blog.csdn.net/zc02051126 有限元算法1-3实现
2021-11-21 10:05:27
1.52MB
1
C-----------------------------------------------------------------------
program runocc
C-----------------------------------------------------------------------
c
C OCCAM 2.0: Steven Constable IGPP/SIO La Jolla CA 92093-0225
c Program Revision 2.01, 13 Jan 1993
c
c runocc is the calling program for OCCAM v2.0.
C
C REFERENCES: CONSTABLE, PARKER & CONSTABLE, 1987: GEOPHYSICS 52, 289-300.
C DEGROOT-HEDLIN & CONSTABLE, 1990: GEOPHYSICS 55, 1613-1624.
C CONSTABLE, 1991: GEOPHYS. J. INT. 106, 387-388.
C CONSTABLE, 1992: OCCAM DISTRIBUTION NOTES.
C
C REVISION HISTORY:
C MARCH 1986, (VERSION 1.2)
C OCTOBER 1987, (VERSION 1.3)
C AUGUST 1988, (VERSION 1.4)
C JANUARY 1989 (WITH THANKS TO C.DEGROOT-HEDLIN), (VERSION 1.5)
C SEPTEMBER 1989, (VERSION 1.5.2)
C APRIL 1992 (MORE THANKS TO CATHERINE AND ALSO WSE) (VERSION 2.0)
C
C IF YOU OBTAIN THIS CODE FROM A THIRD PERSON, PLEASE SEND YOUR NAME AND ADDRESS
C TO S. CONSTABLE. YOU WILL THEN RECEIVE UPDATES, NEWS ON BUGS, ETC.
C
C VERSION 2.0 IS A RE-WRITE TO MAKE OCCAM INDEPENDENT OF THE DIMENSIONALITY
C OF THE FORWARD PROBLEM, AND TO MAKE MAXIMUM USE OF DYNAMIC MEMORY ALLOCATION,
C AS THE ARRAY SIZES FOR THE 2D PROBLEM ARE NOW GETTING TOO LARGE FOR COMFORT.
C INCLUDE FILES ARE NOW USED TO STREAMLINE CHANGES OF DIMENSION.
C IMPLICIT DOUBLE PRECISION HAS BEEN REMOVED, AGAIN FOR STORAGE REASONS.
C CALLING PROGRAM IS NOW AN INTEGRAL PART OF THE PACKAGE AS ALL THE MODEL-
C DEPENDENT STUFF HAS BEEN SPUN OFF.
C
C IF YOU OBTAIN THIS CODE FROM A THIRD PERSON, PLEASE SEND YOUR NAME AND ADDRESS
C TO S. CONSTABLE. YOU WILL THEN RECEIVE UPDATES, NEWS ON BUGS, ETC.
C
C SUBROUTINES WHICH MUST BE SUPPLIED BY THE USER:
C
C FORMOD(NP,ND,PM,DP,DM), COMPUTES THE FORWARD FUNCTION FOR MODEL PM() AT
C THE DATA PARAMETERS DP() AN
2021-11-10 11:19:30
55KB
1