<< Click to Display Table of Contents >> transfer_import |
{ TRANSFER_IMPORT.PDE
This problem demonstrates the use of the TRANSFER facility to import
temperatures from TRANSFER_EXPORT.PDE as the source of thermal expansion
driving a stress calculation.
TRANSFER_EXPORT.PDE must be run before running this problem.
}
title 'Testing the TRANSFER input function'
select
painted { paint all contour plots }
variables
U
V
definitions
nu = 0.3 { define Poisson's Ratio }
E = 21 { Young's Modulus x 10^-11 }
G = E/(1-nu^2)
C11 = G
C12 = G*nu
C22 = G
C33 = G*(1-nu)/2
alpha = 1e-3
b = G*alpha*(1+nu)
{ HERE IS THE TRANSFER INPUT FUNCTION: }
transfer('transfer_export_output/transfer.xfr',Temp,Kxfer)
Sxx = C11*dx(U) + C12*dy(V) - b*Temp
Syy = C12*dx(U) + C22*dy(V) - b*temp
Sxy = C33*(dy(U) + dx(V))
initial values
U = 0
V = 0
equations
U: dx(Sxx) + dy(Sxy) = 0
V: dy(Syy) + dx(Sxy) = 0
constraints
integral(u) = 0
integral(v) = 0
integral(dx(v)-dy(u)) = 0
boundaries
Region 1
start "OUTER" (0,0)
natural(U)=0 value(V)=0 line to(1,0)
natural(U)=0 natural(V)=0
arc (center=0,0) to (0,1) { free outer boundary }
value(U)=0 natural(V)=0 line to close
{ free inner boundary }
start "INNER" (0.4,0.2)
natural(U)=0 natural(V)=0
arc (center=0.4,0.4)
to (0.6,0.4)
to (0.4,0.6)
to (0.2,0.4)
to close
monitors
grid(x+100*U,y+100*V)
plots
contour(Temp) report(Kxfer)
grid(x+100*U,y+100*V)
vector(U,V) as "Displacement"
contour(U) as "X-Displacement"
contour(V) as "Y-Displacement"
contour(Sxx) as "X-Stress"
contour(Syy) as "Y-Stress"
surface(Sxx) as "X-Stress"
surface(Syy) as "Y-Stress"
end