usere_eigensn.for

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

!For User
! to prescribe hydrostatic EIGENSTRAIN for each body
! eij=e0*dij
!It may be a function of loading step.
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!NthBd: Nth Boby
!NLoad: Nth loading step
!UserIsoEigenSn: value of eigenstrain 
***************************************************
	Function UserIsoEigenSn(NLoad,NthBd)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
	Common/Example/NthExample
*================================================================
	If(NthExample.eq.1) then
	 UserIsoEigenSn=Ex1_UserIsoEigenSn(NLoad,NthBd)
	 return
	elseif(NthExample.eq.2) then
	 UserIsoEigenSn=Ex2_UserIsoEigenSn(NLoad,NthBd)
	 return
	elseif(NthExample.eq.3) then
	 UserIsoEigenSn=Ex3_UserIsoEigenSn(NLoad,NthBd)
	 return
	elseif(NthExample.eq.4) then
	 UserIsoEigenSn=Ex4_UserIsoEigenSn(NLoad,NthBd)
	 return
	endif
*===================================================
	UserIsoEigenSn=0.d0
*===================================================
	return
	end
***************************************************
	Function Ex1_UserIsoEigenSn(NLoad,NthBd)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
*===================================================
	Ex1_UserIsoEigenSn=0.d0
*===================================================
	return
	end
***************************************************
	Function Ex2_UserIsoEigenSn(NLoad,NthBd)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
*===================================================
	Ex2_UserIsoEigenSn=0.d0
*===================================================
	return
	end
***************************************************
	Function Ex3_UserIsoEigenSn(NLoad,NthBd)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
*===================================================
	Ex3_UserIsoEigenSn=0.d0
!===================================================
	If(NthBd.eq.2.or.NthBd.eq.4) Ex3_UserIsoEigenSn=1.
*===================================================
	return
	end
***************************************************
	Function Ex4_UserIsoEigenSn(NLoad,NthBd)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
*===================================================
	Ex4_UserIsoEigenSn=0.d0
!===================================================
	If(NthBd.ge.2) Ex4_UserIsoEigenSn=-0.1
*===================================================
	return
	end