userb_cij.for

来自「两维弹性边界元程序」· FOR 代码 · 共 99 行

!For User
! to prescribe MATERIALS CONSTANTS
!  of the Nth BODY (NthBd)
!ES0: SHEAR MODULUS
!PR0: POISON RATIO
***************************************************
	Subroutine UserElasConst(ES0,PR0,NthBd)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
	Common/Example/NthExample
*=======================================================
	If(NthExample.eq.1) then
	 Call Ex1_UserElasConst(ES0,PR0,NthBd)
	 return
	elseif(NthExample.eq.2) then
	 Call Ex2_UserElasConst(ES0,PR0,NthBd)
	 return
	elseif(NthExample.eq.3) then
	 Call Ex3_UserElasConst(ES0,PR0,NthBd)
	 return
	elseif(NthExample.eq.4) then
	 Call Ex4_UserElasConst(ES0,PR0,NthBd)
	 return
	endif
*===================================================
	If(NthBd.eq.1) then
	 ES0=1.d0
	 PR0=0.3d0
	endif
*========================================================
	return
	end
***************************************************
	Subroutine Ex1_UserElasConst(ES0,PR0,NthBd)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
*===================================================
	If(NthBd.eq.1) then
	 ES0=1.d0
	 PR0=0.3d0
	endif
*========================================================
	return
	end
***************************************************
	Subroutine Ex2_UserElasConst(ES0,PR0,NthBd)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
*===================================================
	If(NthBd.eq.1) then
	 ES0=1.d0
	 PR0=0.3d0
	endif
*========================================================
	return
	end
***************************************************
	Subroutine Ex3_UserElasConst(ES0,PR0,NthBd)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
*===================================================
	If(NthBd.eq.1) then
	 ES0=1.5d0
	 PR0=0.3d0
	elseif(NthBd.eq.2) then
	 ES0=1.d0
	 PR0=0.3d0
	elseif(NthBd.eq.3) then
	 ES0=1.5d0
	 PR0=0.3d0
	elseif(NthBd.eq.4) then
	 ES0=1.d0
	 PR0=0.3d0
	endif
*========================================================
	return
	end
***************************************************
	Subroutine Ex4_UserElasConst(ES0,PR0,NthBd)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
*===================================================
	If(NthBd.eq.1) then
	 ES0=1.d0
	 PR0=0.3d0
	elseif(NthBd.eq.2) then
	 ES0=1.5d0
	 PR0=0.3d0
	elseif(NthBd.eq.3) then
	 ES0=1.5d0
	 PR0=0.3d0
	elseif(NthBd.eq.4) then
	 ES0=1.5d0
	 PR0=0.3d0
	endif
*========================================================
	return
	end