Design Sensitivity Analysis of Thermomechanical Problems under Thermal and Mechanical loads with Complex Variables Method

Document Type : Original Research

Authors
M. Sheikhi Azqandi 1
M. Hasan zadeh 2
1 University of Birjand
2 Department of Mechanical Engineering, Gorgan Branch, Islamic Azad University, Gorgan, Iran
Abstract
On the mechanics of solids, in many cases, the body is affected by mechanical force and heat load simultaneously. One of the most widely used and important parameters in the optimal design of engineering problems is the utilization of the sensitivity numbers with respect to the values of its design variables. The present paper analyzes sensitivity in thermomechanical problems using the complex variables method. The most important drawback of past methods for the analysis of sensitivity values such as finite difference and semi-analytical method is the dependence of the calculated values of sensitivity numbers on the step size, which leads to unreliable answers. Complex variables method (CVM) is a new robust approach that uses a Taylor series extension in a complex space and which is not sensitive to step size, disturbance and given that during the computational process, subtraction of two numbers is almost equally avoided, so does not influence from error caused by removing meaningful figures that commonly occur in finite difference method due to subtracting of two near the number. As a result, it is possible to achieve accurate answers by choosing a small step size. Also, it has some potential advantages over other methods. At first order sensitivities using the complex variables method, the implementation is straightforward, only requiring a perturbation of the finite element mesh along the imaginary axis. Implementation of sensitivity using CVM requires complex variable finite element code such that complex nodal coordinates can be used to implement a perturbation in the shape of interest in the complex domain. All resulting finite element outputs such as temperature, displacements, strains, stresses, etc. become complex and accurate derivatives of all finite element outputs to the shape parameter of interest are available. So, unlike most methods, using the proposed method does not require calculations to select the appropriate change in step size. However, the value of the change in step size is not known in advance and changes from one problem to another Therefore, it is necessary to repeat the simulation to ensure the results. The study found that the source of the rounding error was due to the rotation of the element, this also makes it difficult to use two methods, semi-analytical and finite-difference. In this case, the use of CVM is not limited and can be achieved by effective answers without worrying about the small step size. In the present study, to examine how efficient and valid the proposed method is, how it works to solve several thermomechanical problems is described. The methodologies are demonstrated using two-dimensional finite element models. Obtained sensitivity derivatives are compared to the semi-analytical solution and also finite-difference solution and it is shown that the proposed method is effective and can predict stable and accurate sensitivity results. The complex variable method presented in this study can be used to solve a wide range of engineering problems.

Keywords

Subjects

[1] Sheikhi Azqandi, M., Ghoddosian, A. Classical optimization methods for engineering design, Bozorgmehr university of Qaenat press, 2020 (in Persian).
[2] Ghoddosian, A., Sheikhi Azqandi, M., Meta-heuristic optimization Methods in engineering, Semnan university press, 2013(in Persian).
[3] Haftka, Raphael T., Gürdal, Z., Elements of Structural Optimization, Springer, 1992.
[4] Sadd, Martin, Elasticity Theory, Applications, and Numerics, Third Edition) Pages 283-321, chapter 10 (Complex Variable Methods), Academic Press, 2014.
[5] Anderson, WK, Newman, J.C, Whitfield D.L, Nielsen E.J. Sensitivity analysis for Navier–Stokes equations on unstructured meshes using complex variables. AIAA J, 39(1):56–63. 2001.
[6] Newman J.C, Whitfield, D.L, Anderson, W.K., Step-size independent approach for multidisciplinary sensitivity analysis, J Aircraft, 40(3):566–73, 2003.
[7] Burg C.O.E, Newman J.C, Efficient Computationally. Numerically exact design space derivatives via the complex Taylor’s series expansion method. Comput Fluids, 32(3): 373–83, 2003.
[8] Martins, J., Kroo, I. M., Alonso, J. J., An automated method for sensitivity analysis using complex variables, AIAA paper, 689: 2000-2012, 2001.
[9] Martins, J., P. Sturdza, J., Alonso. J., The complex-step derivative approximation, ACM Transactions on Mathematical Software (TOMS), 29(3): 245-262, 2003.
[10] Wang B.P, Apte A.P. Complex variable method for eigen solution sensitivity analysis, AIAA J., 44(12): 2958–61, 2006.
[11] Hassanzadeh, M., Kashani, S., Computation of first and second-order sensitivities for steady state incompressible laminar flow using extended complex variables method, Modares Mechanical Engineering, 19(1): 237-246, 2019, (in Persian).
[12] Gao, X.W, He, M.C., A new inverse analysis for multi-region heat conduction BEM using complex-variable-differentiation method, Eng Anal Boundary Elem, 29(8): 788–95, 2005.
[13] Araújo, A.L., Soares, C. M., Herskovits, J., Pedersen, P., Development of a finite element model for the identification of mechanical and piezoelectric properties through gradient optimisation and experimental vibration data, Composite Structures, 58(3): 307-318, 2002.
[14] Haukaas, T., Scott, M. H., Shape sensitivities in the reliability analysis of nonlinear frame structures, Computers structures, 84(15): 964-977, 2006.
[15] Voorhees, A., Millwater H., Bagley, R., Golden, P., Fatigue sensitivity analysis using complex variable methods, International Journal of Fatigue, 40: 61–73, 2012.
[16] Kim, J., Bates D.G., Postlethwaite, I., Nonlinear Robust performance analysis using complex-step gradient approximation. Automatica, 42(1):177–82, 2006.
[17] Tamayo, D. R., Montoya, A., Millwater, H, A virtual crack extension method for thermoelastic fracture using a complex-variable finite element method, Engineering Fracture Mechanics, 192: 328-342, 2018.
[18] Randal Fielder, R., Millwater, H., Montoya, A., Golden, P., Efficient estimate of residual stress variance using complex variable finite element methods, International Journal of Pressure Vessels and Piping, 173:101-113, 2019.
[19] Voorhees, A., Millwater, H., Bagley, R., Complex variable methods for shape sensitivity of finite element models, Finite Elements in Analysis and Design, 47:1146–1156, 2011.
[20] Hassanzadeh, M., Safarpour, P., Employing a new method for computation of design sensitivities for nonlinear FEM, Modares Mechanical Engineering, 18(4): 645-654, 2018 (in Persian).
[21] Yang, R., Shape design sensitivity analysis of thermo elasticity problems, Computer methods in applied mechanics and engineering,102(1): 41-60, 1993.
[22] Lund, E., Finite element based design sensitivity analysis and optimization, Institute of Mechanical Engineering, Aalborg University, Denmark, 1994.
[23] Lyness, J.N., Moler, C.B., Numerical differentiation of analytic functions, SIAM Journal on Numerical Analysis, 4(2): 202-210, 1967.
[24] Sheikhi Azqandi, M., Hasanzade, M., Arjmand, M., Sensitivity analysis based on complex variables in FEM for linear structures, Advances in Computational Design, 4(1): 15-32, 2019.
[25] Koutromanos, I., Fundamentals of Finite Element Analysis: Linear Finite Element Analysis, Wiley, 2018.
Modares Civil Engineering journal
Volume 21, Issue 4
September and October 2021
Pages 7-17