Thermo-elasticity in Time-dependent p... Log Out | Topics | Search
Moderators | Register | Edit Profile

FlexPDE User's Forum » User Postings » Thermo-elasticity in Time-dependent problem - Part 2 « Previous Next »

Author Message
Top of pagePrevious messageNext messageBottom of page Link to this message

fmpgomes
Member
Username: fmpgomes

Post Number: 7
Registered: 04-2010
Posted on Tuesday, August 31, 2010 - 08:39 am:   

Dear Mr. Nelson,
Thank you for your previous explanations. In my last question (07//07/2010) about the numeric problem when calculating the tension, it was sugested to me to add arbitrarily 1000 to the value of the elasticity modulus, which varies from zero to 30 thousand. I would like to ask if there is no other option, because the phenomenum I'm studying is exactly on the early ages.

What is interesting is that the numeric problem did not occur when I used the non-linear heat equation, which means, with a internal heat generation rate depending upon temperature. Somehow, the Flex PDE became more carefull (accurate) in the calculus.

I'd like to take this oportunity to say that I'm very close to developing something really interesting with the FlexPDE and the company in which I work is interested in acquiring two licenses of the software, once the one we have is not with me (I use the student version).

Thank you for your time,

Flávio mamede
application/octet-streamTérmico
SeparateTens-R0d.pde (11.7 k)
application/octet-streamTermoquímico
TermoQuimico-R0.pde (12.3 k)
application/vnd.ms-excelDiference between two models
Diferença Termico-Termoquímico.xls (183.8 k)
Top of pagePrevious messageNext messageBottom of page Link to this message

rgnelson
Moderator
Username: rgnelson

Post Number: 1402
Registered: 06-2003
Posted on Tuesday, August 31, 2010 - 02:04 pm:   

You'll have to be clearer about your problem.
Which script fails?
At what problem time does it fail?
What is the character of the failure?

I ran SeparateTens for 275 timesteps and it seemed to be happy.
Top of pagePrevious messageNext messageBottom of page Link to this message

fmpgomes
Member
Username: fmpgomes

Post Number: 8
Registered: 04-2010
Posted on Thursday, September 02, 2010 - 07:18 am:   

SeparateTens shows a really big oscilation in stress plots, how you can see in EXCEL plots. I don't know the character of the failure. Maybe an insuficient mesh because I use student version... But I will try with professional version of. prof. Gitirana.
Thank you.
Top of pagePrevious messageNext messageBottom of page Link to this message

jerrybrown11743
Member
Username: jerrybrown11743

Post Number: 64
Registered: 03-2004
Posted on Friday, September 03, 2010 - 06:51 pm:   

I found the same thing as Robert. If I prevented E from starting at zero by adding 1000 to the expression for E in region 2, I saw no problems.

When E is zero, the displacement PDEs collapse to zero = zero. I can't imagine why you would want to solve an elastic solid problem with zero elasticity. It's meaningless.
Top of pagePrevious messageNext messageBottom of page Link to this message

jerrybrown11743
Member
Username: jerrybrown11743

Post Number: 65
Registered: 03-2004
Posted on Friday, September 03, 2010 - 06:52 pm:   

By the way, I'm running version 5.1.
Top of pagePrevious messageNext messageBottom of page Link to this message

fmpgomes
Member
Username: fmpgomes

Post Number: 9
Registered: 04-2010
Posted on Friday, September 10, 2010 - 06:59 am:   

Hi, Mrs. Nelson & Jerry .
E = 0 to 30000 MPa because material is a cement composite. Water-cement mixture become an exothermic reaction. Inicially is a single fluid (in general, a Bingham's fluid), after few hours a viscoelastic solid (or viscoelastoplastic). After few days or months, most of cases is sufficient consider the concrete an elastic solid.
I posted part 2 fo anterior question because:
1) the problem was not correctly answered in part 1; the problem are not E=0, but variable 'temperature' error
2) the solution E'=E+1000 I didn't like

I try run this problem with professional version and run ok, both scritps. FlexPDE worked very well with E ramp funtion from 0 to 30000.
Thank you everyway.
Flávio Mamede

Add Your Message Here
Post:
Username: Posting Information:
This is a private posting area. Only registered users and moderators may post messages here.
Password:
Options: Enable HTML code in message
Automatically activate URLs in message
Action:

Topics | Last Day | Last Week | Tree View | Search | Help/Instructions | Program Credits Administration