Update TSHIKWAND_JABER_FULLY_COUPLED_STRAIN_BASED_SMA_MODEL.f (7040f17c) · Commits · FunctionalMaterials / Coupled Finite Element Simulation of Shape Memory Bending Microactuator

TSHIKWAND_JABER_FULLY_COUPLED_STRAIN_BASED_SMA_MODEL.f

+0 −64

Original line number	Diff line number	Diff line
		@@ -355,13 +355,6 @@ C*************************************************************************

		E1_MA=EPSIL1EQ+C_MA*(T2-T1)/YOUNG
		E2_MA=EPSIL2EQ
		C****************************************************************************
		C VERIFY START CONDITIONS OF AUSTENITE TO TWINNED MARTENSITE TRANSFORMATION
		C AND UPDATE THE VALUES OF VOLUME FACTIONS OF MARTENSITE
		C****************************************************************************
		CALL FORWARD_TRANS_ATM(T1,T2,ES_AM,XCI,XCI_T,
		+ XCI_D,H,DXCI,EPSIL2EQ)

		C*************************************************************************
		C VERIFY START CONDITIONS OF THE TRANSFORMATION OF AUSTENITE OR/AND
		C TWINNED MARTENSITE TO DETWINNED MARTENSITE
		@@ -419,63 +412,6 @@ C*************************************************************************
		END SUBROUTINE

		C##############################################################################
		C##############################################################################
		SUBROUTINE FORWARD_TRANS_ATM(T1,T2,ES_AM,XCI,XCI_T,
		+ XCI_D,H,DXCI,EPSIL2EQ)

		IMPLICIT NONE
		REAL*8 T1,T2,H,ES_AM,EPSIL2EQ,DELTA
		REAL*8 XCI,XCI_T,XCI_D,DXCI,DXCI_T,DXCI_D
		C****************************************************************************
		C MATERIAL PARAMETERS
		C****************************************************************************
		REAL*8 EPS_L,EA,EM,NUA,NUM,C_MA,C_AM,MS,MF,AS,AF
		REAL*8 SS_CR,SF_CR,ALPHAM,ALPHAA
		COMMON / MATERIAU1/ EPS_L,EA,EM,NUA,NUM
		COMMON / MATERIAU2/C_MA,C_AM,MS,MF,AS,AF
		COMMON /MATERIAU3/ SS_CR,SF_CR,ALPHAM,ALPHAA
		C##############################################################################

		IF(T2.LT.T1)THEN
		IF(XCI.LT.1.D0)THEN
		IF(T2.LT.MS)THEN
		IF(EPSIL2EQ.LT.ES_AM)THEN
		H=1.D0
		IF(T2.LT.MF)THEN
		IF(T1.LT.MS)THEN
		DELTA=(T2-MS)/(MF-MS)
		ELSE
		DELTA=(T2-T1)/(MF-T1)
		ENDIF
		ELSE
		DELTA=1.D0
		ENDIF
		IF(DELTA.GE.1.D0)THEN
		DXCI_T=1.D0-XCI
		DXCI_D=0.D0
		DXCI=1.D0-XCI
		ELSE
		DXCI_T=DELTA*(1.D0-XCI)
		DXCI_D=0.D0
		DXCI=DELTA*(1.D0-XCI)
		ENDIF
		XCI_T=XCI_T+DXCI_T
		XCI_D=XCI_D+DXCI_D
		XCI=XCI+DXCI
		ENDIF
		ENDIF
		ENDIF
		ENDIF

		RETURN
		END SUBROUTINE
		C##############################################################################


		C#########################################################################
		C#########################################################################
		C#########################################################################
		C#########################################################################
		C#########################################################################

		SUBROUTINE FORWARD_TRANS(T3,E1_AM,E2_AM,ES_AM,