Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- TITLE 'New Problem' { the problem identification }
- VARIABLES { system variables }
- Temp Volt
- ux,uy
- DEFINITIONS
- v=0.3
- scale=1e3
- E
- CTE
- G=E/((1+v)*(1-2*v))
- C11=G*(1-v)
- C12=G*v
- C21=G*v
- C22=G*(1-v)
- C33=G*(1-2*v)/2
- b=g*cte*(1+v)
- lambda ro
- k=1000
- Vg=100 Va=0
- Ta=273
- DelT=Temp-Ta
- j=-ro*grad(Volt)
- q=-lambda*grad(Temp)
- heat=dot(j,j)/ro
- EQUATIONS
- Volt: div(j)=0
- ux:dx(C11*dx(ux)+C12*dy(uy)-b*DelT) + dy(C33*(dy(ux)+dx(uy)))= 0
- uy: dx(C33*(dy(ux)+dx(uy))) +dy(C21*dx(ux)+C22*dy(uy)-b*DelT) = 0
- Temp: div(q)=heat
- BOUNDARIES
- REGION 1
- E=10e9
- CTE=1000e-6
- lambda=100 ro=100
- START(0,0)
- natural(temp)=k *(temp-Ta)natural(volt)=0
- line to(5,0)
- natural(temp)=0 value(volt)=Va
- line to(5,1)
- natural(temp)=k *(temp-Ta)natural(volt)=0
- line to(0,1)
- natural(temp)=0 value(volt)=Vg
- LINE TO CLOSE
- region 2
- E=1e9
- CTE=1000e-6
- lambda=100 ro=100
- start(-4,-4)
- line to(8,-4)
- line to(8,0)
- line to(-4,0)
- line to close
- region 3
- E=100e9
- cte=1e-5
- lambda=100 ro=100
- start (0,1)
- line to(0,7)
- line to (1,7)
- line to (1,2.5)
- line to (2,2.5)
- line to (2,7)
- line to (3,7)
- line to (3,2.5)
- line to (4,2.5)
- line to (4,7)
- line to (5,7)
- line to (5,1)
- line to close
- PLOTS
- CONTOUR(Volt) vector(j)
- CONTOUR(Temp) vector(q)
- grid(x,y)
- grid(x+scale*ux,y+scale*uy)
- END
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement