线性和非线性有限元分析 4（）.pptVIP

常识 常用烤箱的温度是 180-350 C 在250C 左右紅花油开始汽化 牛肉80C就可以煮熟 有限元（ANSYS)分析结果 用ANSYS和简化手算结果在2mm处温度和时间曲线比较 局部边界载荷向量(Neumann) 局部边界载荷向量(Robin ) 以上双线性热学矩形元的局部边界矩阵类似于线性三角元的边界局部矩阵的详细计算公式请参考我的课本APENDIX B。 由于线性三角元和以上双向线性矩阵元在公共边上都是线性的，故它们可以搭配使用： 三角元和矩形元搭配使用例子 (见附件） 有限元热瞬态分析简介 这部分我们使用课本：的第九章和 附录 C 温度和应力的单向耦合焊接案例 One-D wielding problem Wilding for 1s, and then let it cooling in the water for 60s. Number of elements The Corresponding Shape Functions 2 quadratic element and 2 linear element Shape functions Finite element function Quadratic element: Linear element: 自编以上问题的MATLAB 程序 one_D_welding(2, 2, 1000, 14, 1, 0.01, 7840, 498, [1,1500], 25, 175, 1, 2/3,60) %Where QNE: 2 quadratic elements; %LNE: 2 linear elements; %h: convection coefficient=1000 ; k: conductivity= 14; L: total length=1m; Dens: density=7848 %d: is the diameter=0.01; c: specific heat=498 %T_BC: boundary temperature which is a 1 by 2 matrix contains the time (s) and %temperature. [1s, 1500℃] %T_AMB: ambient temp=25℃; T0: initial temp=175℃. det_t: is the time for one step=1s %theta: 2/3; res_t: is the time specimen when cooling=60s; The power of wielding gun at the wielding step Q = 2.0045e+008 W Out put We can not see the temperature distribution near the tip clearly by just 4 element, so I calculate it again just concerned about the area near the tip, just like zoom in to the tip. My setting is divided the 0.01m tip area into 10 linear element, and the time step is 0.1s. The result shows below: In put one_D_welding(0, 10, 1000, 14, 0.01, 0.01, 7840, 498, [1,1500], 25, 175, 0.1, 2/3,10) Out put Power of the wielding gun at each wielding step Q = 1.0e+008 *(0.1632 0.8036 0.9596 1.0333 1.0742 1.1000 1.1176 1.1303 1.1399 1.1474) W Axisymmetric 2D FEA Model – by Abaqus 10mm 1000mm Circular Bar Model Diameter: 20 mm Length : 1000mm Element type: DCAX4 Axisymmetric element O O Material Information Thermal Analysis: Density: 7890 kg/m3 Specific heat: 0.498 kJ/(kg oC) Conductivity: 14 W/(m oC