Thermoelasticity Log Out | Topics | Search
Moderators | Register | Edit Profile

FlexPDE User's Forum » User Postings » Thermoelasticity « Previous Next »

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

Nir Pony (pony)
New member
Username: pony

Post Number: 1
Registered: 04-2009
Posted on Sunday, April 12, 2009 - 10:45 pm:   

Hi all,
I am trying to solve the set of coupled thermoplastic equations in 2D for a beam using version 5.1.0:
{equation of motion}
Dis: Mu* div(grad (Dis))+(Lambda+mu)*grad(div(Dis))+X=gama*grad(Temp)+rho*dt(dt(Dis)
{Heat conduction eqation}
Temp: div(grad (Temp))-1/aleph*dt(Temp)-etha*div(dt(Dis)=0

Where T is a scalar, and Dis is a 2D vector (x,y)

IC:

Temp=0
Dis=(0,y+0.1*x^2)
dt(dis)=(0,0)

The first approach was to break the equations along x and y. when I run it I got an error “memory protection fault”.
The solution did not converge (the number of unknown).

The second approached was to define Dis as a vector variable “dis = VECTOR (disx, disy)
I got an error “illegal item encountered in input”.
I get the same error when I run the example 'Vector Variables'(http://www.flexpde.com/help/index.html?vector_variables2.html)


What will be the optimal approach to solve this problem?
Why I get the error when using the vector example?
How can I define initial condition (displacement) for a vector?

Pleas guide me how to address this problem.

Thank you,

pony
Top of pagePrevious messageNext messageBottom of page Link to this message

Robert G. Nelson (rgnelson)
Moderator
Username: rgnelson

Post Number: 1243
Registered: 06-2003
Posted on Monday, April 13, 2009 - 02:07 pm:   

1.
VECTOR variables are a capability of FlexPDE version 6, and are not supported in version 5. You will need to declare each component as a separate variable and provide a scalar equation for each.

The HELP system we show on our website is for version 6 (the current shipping version). For version 5 Help, see the Help system installed on your computer by the version 5 installation.

2.
FlexPDE will not accept second derivatives in time. You will need to declare an additional variable for the first time derivatives. See "Samples | Time_Dependent | Stress | Vibrate.pde" for a similar use of second-order time derivatives.

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