Abstract:<p>We present a Bernstein–Bézier finite element method with local mesh refinement to simulate the heat transfer problem in the welding process. A novel constraint method is formulated to handle the incompatible degrees of freedom associated with the transitional elements during the local quad-tree/oct-tree refinement. A dynamic refinement strategy is proposed with a local <em>L</em><sub>2</sub> projection approach to transfer the control variables between different meshes, which enables the reduction of the model size, on the premise that good accuracy of the solution around the heat source is obtained. The cladding process during welding is simulated by the element addition technique. Exactness of the geometry is preserved by use of the rational form of the basis functions, which avoids frequent communication between the CAD and simulation when refining the mesh locally. The proposed method is investigated through examples, including convergence studies and industrial applications. The numerical results obtained using the proposed method and serendipity finite elements are compared in terms of accuracy and efficiency, revealing the potential of the proposed method for welding analysis.</p>
computer science, interdisciplinary applications,engineering, multidisciplinary, software engineering
What problem does this paper attempt to address?
This paper aims to solve several key problems in heat transfer analysis during the welding process:
1. **High - temperature gradients caused by local heating**: During the welding process, due to high - intensity local heating, the material will exhibit complex nonlinear behavior near the heating area. This requires that the mesh division near the heat source must be dense enough to ensure the convergence of the problem and the rationality of the solution. However, such dense mesh division will significantly increase the computational burden, especially when dealing with complex geometric structures.
2. **Adaptive mesh refinement**: In order to effectively deal with the above problems, the paper proposes an adaptive mesh refinement technique based on the Bernstein - Bézier finite - element method. This method can dynamically refine the mesh near the heat source while maintaining the accuracy of the geometric shape, avoiding frequent communication between CAD and simulation in traditional methods.
3. **Treatment of transition elements**: During the local mesh refinement process, "hanging nodes" will be generated, and these nodes will cause mesh incompatibility. The paper proposes a new constraint method to deal with these incompatible degrees of freedom, ensuring the quality of the mesh and the accuracy of the solution.
4. **Simulation of the surfacing process**: During the welding process, material addition, that is, surfacing, usually occurs. The paper simulates this process by directly adding elements, effectively dealing with the material addition problem in the surfacing process by using the \(C^0\) continuity characteristics of Bernstein - Bézier elements.
### Formula summary
- **Bernstein basis function**:
\[
B_{i,p}(\xi)=\frac{p!}{i!(p - i)!}\xi^i(1 - \xi)^{p - i},\quad i = 0,1,\ldots,p,\quad\xi\in[0,1]
\]
- **Bézier surface**:
\[
S(\xi,\eta)=\sum_{i = 0}^{p}\sum_{j = 0}^{q}B_{i,p}(\xi)B_{j,q}(\eta)P_{i,j}
\]
- **Rational Bézier surface**:
\[
R(\xi,\eta)=\frac{\sum_{i = 0}^{p}\sum_{j = 0}^{q}w_{i,j}B_{i,p}(\xi)B_{j,q}(\eta)P_{i,j}}{\sum_{i = 0}^{p}\sum_{j = 0}^{q}w_{i,j}B_{i,p}(\xi)B_{j,q}(\eta)}
\]
- **Heat conduction equation**:
\[
\rho c\frac{\partial u}{\partial t}-\nabla\cdot(k\nabla u)=Q(x,t)
\]
- **Boundary conditions**:
\[
\begin{aligned}
&u(x,t)=\bar{u}(x,t),\quad x\in\Gamma_d\\
&-k\nabla u\cdot\mathbf{n}=q(x,t),\quad x\in\Gamma_n\\
&-k\nabla u\cdot\mathbf{n}=h_c(u - u_a),\quad x\in\Gamma_c
\end{aligned}
\]
- **Discrete weak form**:
\[
\mathbf{K}\mathbf{u}+\mathbf{C}\frac{\partial\mathbf{u}}{\partial t}=\mathbf{Q}
\]
- **Time integration**:
\[
\mathbf{u}^{t+\Delta t}=\mathbf{u}^t+\Delta t\left((1-\alpha\right)