牛顿-拉弗逊法求解非线性问题迭代过程的简单程序例子.docVIP

牛顿－拉弗逊法求解非线性问题迭代过程的简单程序例子C 一个牛顿－拉弗逊法求解非线性问题迭代过程的简单程序例子 可以进行荷载步长自动调整 读者可参考该模块加入相应的子程序编制自己的非线性计算程序module Main_Iteration use TypeDef ! 数据结构定义模块 use Data_Input ! 数据输入模块 use Data_Output ! 数据输出模块 use Elem_Prop ! 单元属性模块 use MatSolve ! 矩阵求解模块 implicit none contains subroutine RC2D_Main() type(typ_GValue) :: GValue !总体控制变量数据结构数组 type(typ_Node),pointer :: Node(:) !节点数据结构数组 type(typ_Elem),pointer :: Elem(:) !混凝土单元数据结构数组 type(typ_Rebar),pointer :: Rebar(:) !钢筋单元数据结构数组 type(typ_Load),pointer :: Load(:) !荷载数据结构数组 type(typ_Load),pointer :: Initial_Load(:) !初始荷载 type(typ_Material),pointer :: Material(:) !材料数据结构数组 type(typ_Support),pointer :: Support(:) !支座数据结构数组 call DataInput(GValue,Node,Elem,Rebar,Material,Load,Support, Initial_Load) !读入数据 call Iteration(GValue,Node,Elem,Rebar,Load,Support,Initial_Load) !迭代核心程序 call DataOutput(GValue,Node,Elem,Rebar,Material,Load,Support) !输出数据 write(*,*) 计算成功结束 read(*,*) deallocate(node) deallocate(Elem) deallocate(Rebar) deallocate(Load) deallocate(Support) return end subroutine subroutine Iteration(GValue0,Node0,Elem0,Rebar0,Load0,Support0,Initial_Load0) type(typ_GValue) :: GValue0 type(typ_Node) :: Node0(:) type(typ_Elem) :: Elem0(:) type(typ_Rebar) :: Rebar0(:) type(typ_Load):: Load0(:) type(typ_Load):: Initial_Load0(:) type(typ_Support) :: Support0(:) type(typ_GValue) :: GValue type(typ_Node) :: Node(GValue0%NNode) type(typ_Elem) :: Elem(GValue0%NElem) type(typ_Rebar) :: Rebar(GValue0%NRebar) type(typ_Load):: Load(GValue0%NLoad) type(typ_Load):: Initial_Load(GValue0%NInitial_Load) type(typ_Support) :: Support(GValue0%NSupport) Integer(ikind) :: I,J,II,III GValue=GValue0 GValue%FinishedStep=0 II=1 !开始计算加载步数 LOOP_A: do write(*,*) ======================================================= write(*,*) Step,II,总加载步数,GValue%NStep write(*,*) 完成,(GValue%FinishedStep)/real(GValue%NStep)*100.0,% !对数据进行备份，在荷载步调整操作时恢复上一步的参数 do I=1,GValue%NNode; Node(I)=Node0(I); end do do I=1,GValue%NElem; Elem(I)=Elem0(I); end do do I=1,GValue%NRebar; Rebar(I)=Rebar