userg_fric.for

两维弹性边界元程序

!For User
! to prescribe FRICTION COEFFICIENT for each group of elements
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!NthGrp: Nth group
!FC: coefficient of Coulomb friction
!WTaoD: penalty parameter, normally taken 1/10000 of displacement magnitude
***************************************************
	Subroutine UserFrictParam(NthGrp,FC,WTaoD)
	Real*8 FC,WTaoD
	Integer NthGrp
	Common/Example/NthExample
*================================================================
	If(NthExample.eq.1) then
	 Call Ex1_UserFrictParam(NthGrp,FC,WTaoD)
	 return
	elseif(NthExample.eq.2) then
	 Call Ex2_UserFrictParam(NthGrp,FC,WTaoD)
	 return
	elseif(NthExample.eq.3) then
	 Call Ex3_UserFrictParam(NthGrp,FC,WTaoD)
	 return
	elseif(NthExample.eq.4) then
	 Call Ex4_UserFrictParam(NthGrp,FC,WTaoD)
	 return
	endif
*===================================================
	FC=0.0
	WTaoD=1.d-5
*===================================================
	return
	end
***************************************************
	Subroutine Ex1_UserFrictParam(NthGrp,FC,WTaoD)
	Real*8 FC,WTaoD
	Integer NthGrp
*===================================================
	FC=0.0
	WTaoD=1.d-5
*===================================================
	return
	end
***************************************************
	Subroutine Ex2_UserFrictParam(NthGrp,FC,WTaoD)
	Real*8 FC,WTaoD
	Integer NthGrp
*===================================================
	FC=0.3
	WTaoD=1.d-5
*===================================================
	return
	end
***************************************************
	Subroutine Ex3_UserFrictParam(NthGrp,FC,WTaoD)
	Real*8 FC,WTaoD
	Integer NthGrp
*===================================================
	FC=0.
	WTaoD=1.d-5
*===================================================
	return
	end
***************************************************
	Subroutine Ex4_UserFrictParam(NthGrp,FC,WTaoD)
	Real*8 FC,WTaoD
	Integer NthGrp
*===================================================
	FC=0.
	WTaoD=1.d-5
*===================================================
	return
	end