useri_initidmg.for

两维弹性边界元程序

!For User
! to prescribe nodal INITIAL DAMAGE PARAMETER,
! which only applies to elements of crack and interface.
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!NN,NE: Nth node of Nth element
!NthGrp: Nth group
!XNod(1/2): nodal coordinates
!UserInitWDmg: damage parameter
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!BE SURE to be consistent with UserH_EKCohes.for
!For perfect bonding, set UserInitWDmg to be a tiny number (1.d-9 or so)
*************************************************************
	Function UserInitWDmg(NN,NE,NthGrp,XNod)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
	Dimension XNod(2)
	Common/Example/NthExample
*================================================================
	If(NthExample.eq.1) then
	 UserInitWDmg=Ex1_UserInitWDmg(NN,NE,NthGrp,XNod)
	 return
	elseif(NthExample.eq.2) then
	 UserInitWDmg=Ex2_UserInitWDmg(NN,NE,NthGrp,XNod)
	 return
	elseif(NthExample.eq.3) then
	 UserInitWDmg=Ex3_UserInitWDmg(NN,NE,NthGrp,XNod)
	 return
	elseif(NthExample.eq.4) then
	 UserInitWDmg=Ex4_UserInitWDmg(NN,NE,NthGrp,XNod)
	 return
	endif
*==================================================================
	If(IABS(NthGrp).le.100) then
	 UserInitWDmg=1.d-9
	elseif(IABS(NthGrp).gt.100) then
	 UserInitWDmg=1.d6
	endif
*==================================================================
	return
	end
*************************************************************
	Function Ex1_UserInitWDmg(NN,NE,NthGrp,XNod)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
	Dimension XNod(2)
*==================================================================
	Ex1_UserInitWDmg=1.d9
*==================================================================
	return
	end
*************************************************************
	Function Ex2_UserInitWDmg(NN,NE,NthGrp,XNod)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
	Dimension XNod(2)
*==================================================================
	Ex2_UserInitWDmg=1.d9
*==================================================================
	return
	end
*************************************************************
	Function Ex3_UserInitWDmg(NN,NE,NthGrp,XNod)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
	Dimension XNod(2)
*==================================================================
	Ex3_UserInitWDmg=1.d-7
*==================================================================
	return
	end
*************************************************************
	Function Ex4_UserInitWDmg(NN,NE,NthGrp,XNod)
	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
	Dimension XNod(2)
*==================================================================
	If(IABS(NthGrp).le.100) then
	 UserInitWDmg=1.d-7
	elseif(IABS(NthGrp).gt.100) then
	 UserInitWDmg=1.d6
	endif
*==================================================================
	Ex4_UserInitWDmg=UserInitWDmg
*==================================================================
	return
	end