×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Are you an
Engineering professional?
Join Eng-Tips Forums!
• Talk With Other Members
• Be Notified Of Responses
• Keyword Search
Favorite Forums
• Automated Signatures
• Best Of All, It's Free!

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

#### Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

# CURRENT TEMPERATURE IN ABAQUS

## CURRENT TEMPERATURE IN ABAQUS

(OP)
Dear all,

I am doing a phase transformation analysis in Abaqus . I tried to use UEXPAN subroutine to include the volumetric strain due to the coefficient of thermal expansion. But when I am using the Koistinen-Marburger equation for the martensite formation I want to call the current temperature in Abaqus. Abaqus defies the current temperature as TEMP(1). But when I gave the TEMP(1) it is not picking the right temperature. Can anyone help?

### RE: CURRENT TEMPERATURE IN ABAQUS

Indeed TEMP(1) is current temperature (at the end of the increment). There's also TEMP(2) which is temperature increment. But if your subroutine doesn't work properly then there might be some error in the code. It's really hard to find exemplary Abaqus subroutines, especially for phase transformation, but there is one code attached to the paper titled "Two-dimensional simulation of thermal cutting of low-alloyed steels" by A. Laitinen from Tampere University of Applied Sciences.

### RE: CURRENT TEMPERATURE IN ABAQUS

(OP)
Dear FEA way,
In the code, I tried with the thermal stress alone and it gives me correct output. Also, the code runs fine with austenite transformation. For checking the code I just created a simple model with 1500 predefined temperature and tried to cool down to room temperature. The expansion coefficient was given as subroutine (ie the equation for martensite expansion)In the subroutine if I change the temp(1) to a numerical value say for example 200/300 degree Celcius I will get exact analytical stress at this temperatures.

SUBROUTINE UEXPAN(EXPAN,DEXPANDT,TEMP,TIME,DTIME,PREDEF,
1 DPRED,STATEV,CMNAME,NSTATV,NOEL)
C
INCLUDE 'ABA_PARAM.INC'
C
CHARACTER*80 CMNAME
C
DIMENSION EXPAN(*),DEXPANDT(*),TEMP(2),TIME(2),PREDEF(*),
1 DPRED(*),STATEV(NSTATV)
REAL alpha,P,G,L,K,E,MS

MS=375

P=TEMP(1)-MS
K=0.011*P
E= EXP(K)
L=(-0.011)*E

G=(3.75)*(10.**(-3))*1.17*L*TEMP(2)

EXPAN(1)= G

RETURN
END

#### Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

#### Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

#### Resources

White Paper - The Evolving Landscape of Commercial Battery-Powered Trucks
Whatâ€™s driving the evolving landscape of truck electrification? What are the barriers, motivators and strategies for accelerating the electric transition? What insights and resources are available for todayâ€™s design engineers working to achieve industry disruption and evolution? For answers to these and other pertinent questions, read this white paper. Download Now
A lot has changed since the 90s. You don't surf the Web using dial-up anymore, so why are you still using a legacy PLM solution that's blocking your ability to innovate? To develop and launch products today, you need a flexible, cloud-based PLM, not a solution that's stuck in the past. Download Now

Close Box

# Join Eng-Tips® Today!

Join your peers on the Internet's largest technical engineering professional community.
It's easy to join and it's free.

Here's Why Members Love Eng-Tips Forums:

• Talk To Other Members
• Notification Of Responses To Questions
• Favorite Forums One Click Access
• Keyword Search Of All Posts, And More...

Register now while it's still free!