userf_farss.for

两维弹性边界元程序

!For User
! to prescribe FAR-FIELD LOADING for THE infinite body, if any.
!It may be a function of loading step.
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!NthBd: Nth Boby
!NLoad: Nth loading step
!UserFarSs: value of far-field stress component (I,J)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!BE SURE that NthBdInfinite is given the correct number of the infinite body
! if there is one.
!ONLY ONE infinite body is allowed.
***************************************************
	Function UserFarSs(NLoad,I,J,NthBd)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
	Common/Example/NthExample
*================================================================
	If(NthExample.eq.1) then
	 UserFarSs=Ex1_UserFarSs(NLoad,I,J,NthBd)
	 return
	elseif(NthExample.eq.2) then
	 UserFarSs=Ex2_UserFarSs(NLoad,I,J,NthBd)
	 return
	elseif(NthExample.eq.3) then
	 UserFarSs=Ex3_UserFarSs(NLoad,I,J,NthBd)
	 return
	elseif(NthExample.eq.4) then
	 UserFarSs=Ex4_UserFarSs(NLoad,I,J,NthBd)
	 return
	endif
*===================================================
	NthBdInfinite=1
!===================================================
	If(NthBd.eq.NthBdInfinite) then
	 If(I.eq.1.and.J.eq.1) UserFarSs=0.d0
	 If(I.eq.1.and.J.eq.2) UserFarSs=0.d0
	 If(I.eq.2.and.J.eq.1) UserFarSs=0.d0
	 If(I.eq.2.and.J.eq.2) UserFarSs=0.d0
	else
	 UserFarSs=0.d0
	endif
*===================================================
	return
	end
***************************************************
	Function Ex1_UserFarSs(NLoad,I,J,NthBd)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
*================================================================
	Ex1_UserFarSs=0.d0
*===================================================
	return
	end
***************************************************
	Function Ex2_UserFarSs(NLoad,I,J,NthBd)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
*================================================================
	NthBdInfinite=1
!===================================================
	If(NthBd.eq.NthBdInfinite) then
	 If(I.eq.1.and.J.eq.1) Ex2_UserFarSs=0.d0
	 If(I.eq.1.and.J.eq.2) Ex2_UserFarSs=0.d0
	 If(I.eq.2.and.J.eq.1) Ex2_UserFarSs=0.d0
	 If(I.eq.2.and.J.eq.2) Ex2_UserFarSs=-1.d-1
	else
	 Ex2_UserFarSs=0.d0
	endif
*===================================================
	return
	end
***************************************************
	Function Ex3_UserFarSs(NLoad,I,J,NthBd)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
*================================================================
	Ex3_UserFarSs=0.d0
*===================================================
	return
	end
***************************************************
	Function Ex4_UserFarSs(NLoad,I,J,NthBd)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
*================================================================
	NthBdInfinite=1
!===================================================
	If(NthBd.eq.NthBdInfinite) then
	 If(I.eq.1.and.J.eq.1) Ex4_UserFarSs=0.d0
	 If(I.eq.1.and.J.eq.2) Ex4_UserFarSs=0.d0
	 If(I.eq.2.and.J.eq.1) Ex4_UserFarSs=0.d0
	 If(I.eq.2.and.J.eq.2) Ex4_UserFarSs=1.d-1
	else
	 Ex4_UserFarSs=0.d0
	endif
*===================================================
	return
	end