Advertisement
Guest User

Untitled

a guest
Apr 25th, 2019
99
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.32 KB | None | 0 0
  1. TITLE 'New Problem' { the problem identification }
  2. VARIABLES { system variables }
  3. Temp Volt
  4. ux,uy
  5. DEFINITIONS
  6. v=0.3
  7. scale=1e3
  8. E
  9. CTE
  10. G=E/((1+v)*(1-2*v))
  11. C11=G*(1-v)
  12. C12=G*v
  13. C21=G*v
  14. C22=G*(1-v)
  15. C33=G*(1-2*v)/2
  16. b=g*cte*(1+v)
  17. lambda ro
  18. k=1000
  19. Vg=100 Va=0
  20. Ta=273
  21. DelT=Temp-Ta
  22. j=-ro*grad(Volt)
  23. q=-lambda*grad(Temp)
  24. heat=dot(j,j)/ro
  25. EQUATIONS
  26. Volt: div(j)=0
  27. ux:dx(C11*dx(ux)+C12*dy(uy)-b*DelT) + dy(C33*(dy(ux)+dx(uy)))= 0
  28. uy: dx(C33*(dy(ux)+dx(uy))) +dy(C21*dx(ux)+C22*dy(uy)-b*DelT) = 0
  29. Temp: div(q)=heat
  30. BOUNDARIES
  31. REGION 1
  32. E=10e9
  33. CTE=1000e-6
  34. lambda=100 ro=100
  35. START(0,0)
  36. natural(temp)=k *(temp-Ta)natural(volt)=0
  37. line to(5,0)
  38. natural(temp)=0 value(volt)=Va
  39. line to(5,1)
  40. natural(temp)=k *(temp-Ta)natural(volt)=0
  41. line to(0,1)
  42. natural(temp)=0 value(volt)=Vg
  43. LINE TO CLOSE
  44. region 2
  45. E=1e9
  46. CTE=1000e-6
  47. lambda=100 ro=100
  48. start(-4,-4)
  49. line to(8,-4)
  50. line to(8,0)
  51. line to(-4,0)
  52. line to close
  53.  
  54. region 3
  55. E=100e9
  56. cte=1e-5
  57.  
  58. lambda=100 ro=100
  59. start (0,1)
  60. line to(0,7)
  61. line to (1,7)
  62. line to (1,2.5)
  63. line to (2,2.5)
  64. line to (2,7)
  65. line to (3,7)
  66. line to (3,2.5)
  67. line to (4,2.5)
  68. line to (4,7)
  69. line to (5,7)
  70. line to (5,1)
  71. line to close
  72.  
  73. PLOTS
  74. CONTOUR(Volt) vector(j)
  75. CONTOUR(Temp) vector(q)
  76. grid(x,y)
  77. grid(x+scale*ux,y+scale*uy)
  78. END
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement