An isogeometric numerical study of partially and fully implicit schemes for transient adjoint shape sensitivity analysis

doi:10.1007/s11465-019-0575-5

Frontiers of Mechanical Engineering

2020, Vol. 15

Issue (2): 279-293 https://doi.org/10.1007/s11465-019-0575-5

本期目录

An isogeometric numerical study of partially and fully implicit schemes for transient adjoint shape sensitivity analysis

Zhen-Pei WANG^1,², Zhifeng XIE³, Leong Hien POH¹(

)

¹. Department of Civil and Environmental Engineering, National University of Singapore, Singapore 117576, Singapore
². Institute of High Performance Computing (IHPC), Agency for Science, Technology and Research (A*STAR), Singapore 138632, Singapore
³. China Academy of Launch Vehicle Technology, Beijing Institute of Astronautical Systems Engineering, Beijing 100076, China

全文: PDF(1546 KB) HTML

Abstract：

In structural design optimization involving transient responses, time integration scheme plays a crucial role in sensitivity analysis because it affects the accuracy and stability of transient analysis. In this work, the influence of time integration scheme is studied numerically for the adjoint shape sensitivity analysis of two benchmark transient heat conduction problems within the framework of isogeometric analysis. It is found that (i) the explicit approach ( $β$ = 0) and semi-implicit approach with $β$ <0.5 impose a strict stability condition of the transient analysis; (ii) the implicit approach ( $β$ =1) and semi-implicit approach with $β$ > 0.5 are generally preferred for their unconditional stability; and (iii) Crank–Nicolson type approach ( $β$ =0.5) may induce a large error for large time-step sizes due to the oscillatory solutions. The numerical results also show that the time-step size does not have to be chosen to satisfy the critical conditions for all of the eigen-frequencies. It is recommended to use $β ≈ 0.75$ for unconditional stability, such that the oscillation condition is much less critical than the Crank–Nicolson scheme, and the accuracy is higher than a fully implicit approach.

Key words： isogeometric shape optimization design-dependent boundary condition transient heat conduction implicit time integration adjoint method

收稿日期: 2019-08-17 出版日期: 2020-05-25

Corresponding Author(s): Leong Hien POH

引用本文:

. [J]. Frontiers of Mechanical Engineering, 2020, 15(2): 279-293.
Zhen-Pei WANG, Zhifeng XIE, Leong Hien POH. An isogeometric numerical study of partially and fully implicit schemes for transient adjoint shape sensitivity analysis. Front. Mech. Eng., 2020, 15(2): 279-293.

链接本文:

https://academic.hep.com.cn/fme/CN/10.1007/s11465-019-0575-5
https://academic.hep.com.cn/fme/CN/Y2020/V15/I2/279

Fig.1

Fig.2

Fig.3

Fig.4

Fig.5

Fig.6

$I (i, j)$	Location		Weight	? $I (i, j)$	Location		Weight
$I (i, j)$	$x 1 I$	$x 2 I$	Weight	? $I (i, j)$	$x 1 I$	$x 2 I$	Weight
(1, 1)	0.0100	0.0000	1.00	?(4, 2)	0.0091	0.0121	0.85
(2, 1)	0.0100	0.0026	0.90	?(5, 2)	0.0039	0.0150	0.90
(3, 1)	0.0080	0.0061	0.85	?(6, 2)	0.0000	0.0150	1.00
(4, 1)	0.0061	0.0080	0.85	?(1, 3)	0.0200	0.0000	1.00
(5, 1)	0.0026	0.0100	0.90	?(2, 3)	0.0200	0.0100	1.00
(6, 1)	0.0000	0.0100	1.00	?(3, 3)	0.0238	0.0213	1.00
(1, 2)	0.0150	0.0000	1.00	?(4, 3)	0.0213	0.0238	1.00
(2, 2)	0.0150	0.0039	0.90	?(5, 3)	0.0100	0.0200	1.00
(3, 2)	0.0121	0.0091	0.85	?(6, 3)	0.0000	0.0200	1.00

Tab.1

$C I$	Component	FD			Referential FD
$C I$	Component	$β = 0.5$	$β = 0.75$	$β = 1$	Referential FD
$C 1$	1	135874.5067	135876.4566	135887.5961	135879.5198
$C 2$	1	243335.5767	243333.5103	243350.4778	243339.8549
$C 2$	2	–37170.7865	–37162.2082	–37166.9958	–37166.6635
$C 3$	1	897323.9019	897326.7759	897308.0330	897319.5703
$C 3$	2	306189.0384	306194.4444	306182.1844	306188.5557
$C 4$	1	306189.1548	306194.2916	306176.7857	306186.7440
$C 4$	2	897323.3635	897312.5441	897306.3377	897314.0818
$C 5$	1	–37165.4896	–37167.1631	–37163.8816	–37165.5115
$C 5$	2	243336.8791	243328.4244	243347.0800	243337.4611
$C 6$	2	135881.9864	135871.1306	135886.7594	135879.9588

Tab.2

Fig.7

Tab.3

Fig.8

Fig.9

$I (i, j)$	Location		Weight	? $I (i, j)$	Location		Weight
$I (i, j)$	$x 1 I$	$x 2 I$	Weight	? $I (i, j)$	$x 1 I$	$x 2 I$	Weight
(1, 1)	0.00	100.00	1.00	?(5, 2)	90.00	20.00	1.00
(2, 1)	0.00	80.00	1.00	?(6, 2)	145.00	20.00	1.00
(3, 1)	0.00	30.00	1.00	?(7, 2)	200.00	20.00	1.00
(4, 1)	0.00	0.00	0.71	?(1, 3)	30.00	100.00	1.00
(5, 1)	30.00	0.00	1.00	?(2, 3)	30.00	80.00	1.00
(6, 1)	140.00	0.00	1.00	?(3, 3)	30.00	65.00	1.00
(7, 1)	200.00	0.00	1.00	?(4, 3)	30.00	45.00	1.00
(1, 2)	15.00	100.00	1.00	?(5, 3)	70.00	45.00	1.00
(2, 2)	15.00	80.00	1.00	?(6, 3)	120.00	45.00	1.00
(3, 2)	15.00	65.00	1.00	?(7, 3)	200.00	45.00	1.00
(4, 2)	15.00	20.00	1.00

Tab.4

Fig.10

Fig.11

$C I$	Component	FD			Referential FD
$C I$	Component	$β = 0.5$	$β = 0.75$	$β = 1$	Referential FD
$C 1$	1	–6.2906×10⁵	–6.2891×10⁵	–6.2885×10⁵	–6.2894×10⁵
	2	?2.0000×10	?3.0000×10	1.0000×10	2.0000×10
$C 2$	1	–3.2439×10⁵	–3.2434×10⁵	–3.2422×10⁵	–3.2432×10⁵
	2	???3.4358×10⁵	???3.4351×10⁵	???3.4350×10⁵	???3.4353×10⁵
$C 3$	1	???2.0760×10⁶	???2.0759×10⁶	???2.0758×10⁶	???2.0759×10⁶
	2	???1.1210×10⁶	???1.1210×10⁶	???1.1210×10⁶	???1.1210×10⁶
$C 4$	1	???8.6030×10⁴	???8.6030×10⁴	???8.6010×10⁴	???8.6020×10⁴
	2	–6.2714×10⁵	–6.2713×10⁵	–6.2714×10⁵	–6.2714×10⁵
$C 5$	1	?2.0000×10	?1.0000×10	1.0000×10	?2.0000×10
	2	???1.5853×10⁶	???1.5852×10⁶	???1.5849×10⁶	???1.5851×10⁶

Tab.5

Fig.12

Tab.6

Fig.13

1	Q Li, G P Steven, O M Querin, et al. Shape and topology design for heat conduction by evolutionary structural optimization. International Journal of Heat and Mass Transfer, 1999, 42(17): 3361–3371 https://doi.org/10.1016/S0017-9310(99)00008-3
2	G Xie, Y Liu, B Sunden, et al. Computational study and optimization of laminar heat transfer and pressure loss of double-layer microchannels for chip liquid cooling. Journal of Thermal Science and Engineering Applications, 2013, 5(1): 011004 https://doi.org/10.1115/1.4007778
3	O Sigmund, S Torquato. Design of materials with extreme thermal expansion using a three-phase topology optimization method. Journal of the Mechanics and Physics of Solids, 1997, 45(6): 1037–1067 https://doi.org/10.1016/S0022-5096(96)00114-7
4	T Gao, W Zhang, J Zhu, et al. Topology optimization of heat conduction problem involving design-dependent heat load effect. Finite Elements in Analysis and Design, 2008, 44(14): 805–813 https://doi.org/10.1016/j.finel.2008.06.001
5	A Iga, S Nishiwaki, K Izui, et al. Topology optimization for thermal conductors considering design-dependent effects, including heat conduction and convection. International Journal of Heat and Mass Transfer, 2009, 52(11–12): 2721–2732 https://doi.org/10.1016/j.ijheatmasstransfer.2008.12.013
6	K Yaji, T Yamada, S Kubo, et al. A topology optimization method for a coupled thermal-fluid problem using level set boundary expressions. International Journal of Heat and Mass Transfer, 2015, 81: 878–888 https://doi.org/10.1016/j.ijheatmasstransfer.2014.11.005
7	Q Xia, L Xia, T Shi. Topology optimization of thermal actuator and its support using the level set based multiple-type boundary method and sensitivity analysis based on constrained variational principle. Structural and Multidisciplinary Optimization, 2018, 57(3): 1317–1327 https://doi.org/10.1007/s00158-017-1814-1
8	K K Choi, N H Kim. Structural Sensitivity Analysis and Optimization 1: Linear Systems. New York: Springer, 2005
9	K Dems, B Rousselet. Sensitivity analysis for transient heat conduction in a solid body-Part I: External boundary modification. Structural Optimization, 1999, 17(1): 36–45 https://doi.org/10.1007/BF01197711
10	K Dems, B Rousselet. Sensitivity analysis for transient heat conduction in a solid body-Part II: Interface modification. Structural Optimization, 1999, 17(1): 46–54 https://doi.org/10.1007/BF01197712
11	R T Haftka, C P Shore. Approximation Methods for Combined Thermal/Structural Design. NASA Technical Paper 1428. 1979
12	R T Haftka. Techniques for thermal sensitivity analysis. International Journal for Numerical Methods in Engineering, 1981, 17(1): 71–80 https://doi.org/10.1002/nme.1620170106
13	W H Greene, R T Haftka. Computational aspects of sensitivity calculations in transient structural analysis. Computers & Structures, 1989, 32(2): 433–443 https://doi.org/10.1016/0045-7949(89)90054-0
14	Z Y Gao, R V Grandhi. Sensitivity analysis and shape optimization for preform design in thermo-mechanical coupled analysis. International Journal for Numerical Methods in Engineering, 1999, 45(10): 1349–1373 https://doi.org/10.1002/(SICI)1097-0207(19990810)45:10<1349::AID-NME634>3.0.CO;2-1
15	R T Haftka, D S Malkus. Calculation of sensitivity derivatives in thermal problems by finite differences. International Journal for Numerical Methods in Engineering, 1981, 17(12): 1811–1821 https://doi.org/10.1002/nme.1620171206
16	Z P Wang, S Turteltaub, M M Abdalla. Shape optimization and optimal control for transient heat conduction problems using an isogeometric approach. Computers & Structures, 2017, 185: 59–74 https://doi.org/10.1016/j.compstruc.2017.02.004
17	P Michaleris, D A Tortorelli, C A Vidal. Tangent operators and design sensitivity formulations for transient non-linear coupled problems with applications to elastoplasticity. International Journal for Numerical Methods in Engineering, 1994, 37(14): 2471–2499 https://doi.org/10.1002/nme.1620371408
18	D A Tortorelli, R B Haber, S C Y Lu. Design sensitivity analysis for nonlinear thermal systems. Computer Methods in Applied Mechanics and Engineering, 1989, 77(1–2): 61–77 https://doi.org/10.1016/0045-7825(89)90128-X
19	D A Tortorelli, R B Haber. First-order design sensitivities for transient conduction problems by an adjoint method. International Journal for Numerical Methods in Engineering, 1989, 28(4): 733–752 https://doi.org/10.1002/nme.1620280402
20	Z P Wang, D Kumar. On the numerical implementation of continuous adjoint sensitivity for transient heat conduction problems using an isogeometric approach. Structural and Multidisciplinary Optimization, 2017, 56(2): 487–500 https://doi.org/10.1007/s00158-017-1669-5
21	J H Kane, B L Kumar, M Stabinsky. Transient thermoelasticity and other body force effects in boundary element shape sensitivity analysis. International Journal for Numerical Methods in Engineering, 1991, 31(6): 1203–1230 https://doi.org/10.1002/nme.1620310612
22	Y Jarny, M N Ozisik, J P Bardon. A general optimization method using adjoint equation for solving multidimensional inverse heat conduction. International Journal of Heat and Mass Transfer, 1991, 34(11): 2911–2919 https://doi.org/10.1016/0017-9310(91)90251-9
23	M Kleiber, A Służalec. Material derivative and control volume approaches to shape sensitivity analysis of nonlinear transient thermal problems. Structural Optimization, 1996, 11: 56–63 https://doi.org/10.1007/BF01279656
24	G A Dorai, D A Tortorelli. Transient inverse heat conduction problem solutions via Newton’s method. International Journal of Heat and Mass Transfer, 1997, 40(17): 4115–4127 https://doi.org/10.1016/S0017-9310(97)00044-6
25	R Korycki. Sensitivity analysis and shape optimization for transient heat conduction with radiation. International Journal of Heat and Mass Transfer, 2006, 49(13–14): 2033–2043 https://doi.org/10.1016/j.ijheatmasstransfer.2006.01.007
26	C H Huang, M T Chaing. A transient three-dimensional inverse geometry problem in estimating the space and time-dependent irregular boundary shapes. International Journal of Heat and Mass Transfer, 2008, 51(21–22): 5238–5246 https://doi.org/10.1016/j.ijheatmasstransfer.2008.03.019
27	A Służalec, M Kleiber. Shape optimization of thermo-diffusive systems. International Journal of Heat and Mass Transfer, 1992, 35(9): 2299–2304 https://doi.org/10.1016/0017-9310(92)90072-Z
28	Y X Gu, B S Chen, H W Zhang, et al. A sensitivity analysis method for linear and nonlinear transient heat conduction with precise time integration. Structural and Multidisciplinary Optimization, 2002, 24(1): 23–37 https://doi.org/10.1007/s00158-002-0211-5
29	B Chen, L Tong. Sensitivity analysis of heat conduction for functionally graded materials. Materials & Design, 2004, 25(8): 663–672 https://doi.org/10.1016/j.matdes.2004.03.007
30	R T Haftka, R V Grandhi. Structural shape optimization—A survey. Computer Methods in Applied Mechanics and Engineering, 1986, 57(1): 91–106 https://doi.org/10.1016/0045-7825(86)90072-1
31	F van Keulen, R T Haftka, N H Kim. Review of options for structural design sensitivity analysis, Part 1: Linear systems. Computer Methods in Applied Mechanics and Engineering, 2005, 194(30–33): 3213–3243 https://doi.org/10.1016/j.cma.2005.02.002
32	S Cho, S H Ha. Isogeometric shape design optimization: Exact geometry and enhanced sensitivity. Structural and Multidisciplinary Optimization, 2009, 38(1): 53–70 https://doi.org/10.1007/s00158-008-0266-z
33	X Qian. Full analytical sensitivities in nurbs based isogeometric shape optimization. Computer Methods in Applied Mechanics and Engineering, 2010, 199(29–32): 2059–2071 https://doi.org/10.1016/j.cma.2010.03.005
34	A P Nagy, M M Abdalla, Z Gürdal. Isogeometric sizing and shape optimisation of beam structures. Computer Methods in Applied Mechanics and Engineering, 2010, 199(17–20): 1216–1230 https://doi.org/10.1016/j.cma.2009.12.010
35	A P Nagy, M M Abdalla, Z Gürdal. Isogeometric design of elastic arches for maximum fundamental frequency. Structural and Multidisciplinary Optimization, 2011, 43(1): 135–149 https://doi.org/10.1007/s00158-010-0549-z
36	H Liu, D Yang, X Wang, et al. Smooth size design for the natural frequencies of curved Timoshenko beams using isogeometric analysis. Structural and Multidisciplinary Optimization, 2019, 59(4): 1143–1162 https://doi.org/10.1007/s00158-018-2119-8
37	O Weeger, B Narayanan, M L Dunn. Isogeometric shape optimization of nonlinear, curved 3D beams and beam structures. Computer Methods in Applied Mechanics and Engineering, 2019, 345: 26–51 https://doi.org/10.1016/j.cma.2018.10.038
38	A P Nagy, S T IJsselmuiden, M M Abdalla. Isogeometric design of anisotropic shells: Optimal form and material distribution. Computer Methods in Applied Mechanics and Engineering, 2013, 264: 145–162 https://doi.org/10.1016/j.cma.2013.05.019
39	T Hirschler, R Bouclier, A Duval, et al. Isogeometric sizing and shape optimization of thin structures with a solid-shell approach. Structural and Multidisciplinary Optimization, 2019, 59(3): 767–785 https://doi.org/10.1007/s00158-018-2100-6
40	H Lian, P Kerfriden, S Bordas. Implementation of regularized isogeometric boundary element methods for gradient-based shape optimization in two-dimensional linear elasticity. International Journal for Numerical Methods in Engineering, 2016, 106(12): 972–1017 https://doi.org/10.1002/nme.5149
41	H Lian, P Kerfriden, S Bordas. Shape optimization directly from CAD: An isogeometric boundary element approach using T-splines. Computer Methods in Applied Mechanics and Engineering, 2017, 317: 1–41 https://doi.org/10.1016/j.cma.2016.11.012
42	C Wang, S Xia, X Wang, et al. Isogeometric shape optimization on triangulations. Computer Methods in Applied Mechanics and Engineering, 2018, 331: 585–622 https://doi.org/10.1016/j.cma.2017.11.032
43	Z P Wang, L H Poh, J Dirrenberger, et al. Isogeometric shape optimization of smoothed petal auxetic structures via computational periodic homogenization. Computer Methods in Applied Mechanics and Engineering, 2017, 323: 250–271 https://doi.org/10.1016/j.cma.2017.05.013
44	Z P Wang, L H Poh. Optimal form and size characterization of planar isotropic petal-shaped auxetics with tunable effective properties using IGA. Composite Structures, 2018, 201: 486–502 https://doi.org/10.1016/j.compstruct.2018.06.042
45	D Kumar, Z P Wang, L H Poh, et al. Isogeometric shape optimization of smoothed petal auxetics with prescribed nonlinear deformation. Computer Methods in Applied Mechanics and Engineering, 2019, 356: 16–43 https://doi.org/10.1016/j.cma.2019.07.014
46	Y Wang, D J Benson. Geometrically constrained isogeometric parameterized level-set based topology optimization via trimmed elements. Frontiers of Mechanical Engineering, 2016, 11(4): 328–343 https://doi.org/10.1007/s11465-016-0403-0
47	Y Wang, H Xu, D Pasini. Multiscale isogeometric topology optimization for lattice materials. Computer Methods in Applied Mechanics and Engineering, 2017, 316: 568–585 https://doi.org/10.1016/j.cma.2016.08.015
48	X Xie, S Wang, M Xu, et al. A new isogeometric topology optimization using moving morphable components based on R-functions and collocation schemes. Computer Methods in Applied Mechanics and Engineering, 2018, 339: 61–90 https://doi.org/10.1016/j.cma.2018.04.048
49	Q X Lieu, J Lee. Multiresolution topology optimization using isogeometric analysis. International Journal for Numerical Methods in Engineering, 2017, 112(13): 2025–2047 https://doi.org/10.1002/nme.5593
50	W Hou, Y Gai, X Zhu, et al. Explicit isogeometric topology optimization using moving morphable components. Computer Methods in Applied Mechanics and Engineering, 2017, 326: 694–712 https://doi.org/10.1016/j.cma.2017.08.021
51	H Liu, D Yang, P Hao, et al. Isogeometric analysis based topology optimization design with global stress constraint. Computer Methods in Applied Mechanics and Engineering, 2018, 342: 625–652 https://doi.org/10.1016/j.cma.2018.08.013
52	P Hao, X Yuan, C Liu, et al. An integrated framework of exact modeling, isogeometric analysis and optimization for variable-stiffness composite panels. Computer Methods in Applied Mechanics and Engineering, 2018, 339: 205–238 https://doi.org/10.1016/j.cma.2018.04.046
53	Y Guo, M Ruess. Nitsche’s method for a coupling of isogeometric thin shells and blended shell structures. Computer Methods in Applied Mechanics and Engineering, 2015, 284: 881–905 https://doi.org/10.1016/j.cma.2014.11.014
54	S Y Cai, W H Zhang, J Zhu, et al. Stress constrained shape and topology optimization with fixed mesh: A B-spline finite cell method combined with level set function. Computer Methods in Applied Mechanics and Engineering, 2014, 278: 361–387 https://doi.org/10.1016/j.cma.2014.06.007
55	W Zhang, L Zhao, T Gao, et al. Topology optimization with closed B-splines and Boolean operations. Computer Methods in Applied Mechanics and Engineering, 2017, 315: 652–670 https://doi.org/10.1016/j.cma.2016.11.015
56	Y Wang, Z P Wang, Z Xia, et al. Structural design optimization using isogeometric analysis: A comprehensive review. Computer Modeling in Engineering & Sciences, 2018, 117(3): 455–507 https://doi.org/10.31614/cmes.2018.04603
57	L Xia, Q Xia, X Huang, et al. Bi-directional evolutionary structural optimization on advanced structures and materials: A comprehensive review. Archives of Computational Methods in Engineering, 2018, 25(2): 437–478 https://doi.org/10.1007/s11831-016-9203-2
58	L Meng, W Zhang, D Quan, et al. From topology optimization design to additive manufacturing: Today’s success and tomorrow’s roadmap. Archives of Computational Methods in Engineering, 2019 (in press) https://doi.org/10.1007/s11831-019-09331-1
59	W Kaminski. Hyperbolic heat conduction equation for materials with a nonhomogeneous inner structure. Journal of Heat Transfer, 1990, 112(3): 555–560 https://doi.org/10.1115/1.2910422
60	L Xia, P Breitkopf. Recent advances on topology optimization of multiscale nonlinear structures. Archives of Computational Methods in Engineering, 2017, 24(2): 227–249 https://doi.org/10.1007/s11831-016-9170-7
61	Z P Wang, S Turteltaub. Isogeometric shape optimization for quasi-static processes. International Journal for Numerical Methods in Engineering, 2015, 104(5): 347–371 https://doi.org/10.1002/nme.4940
62	Q Xia, T Shi, S Liu, et al. A level set solution to the stress-based structural shape and topology optimization. Computers & Structures, 2012, 90: 55–64 https://doi.org/10.1016/j.compstruc.2011.10.009
63	Q Xia, T Shi, L Xia. Stable hole nucleation in level set based topology optimization by using the material removal scheme of BESO. Computer Methods in Applied Mechanics and Engineering, 2019, 343: 438–452 https://doi.org/10.1016/j.cma.2018.09.002
64	J N Reddy, D K Gartling. The Finite Element Method in Heat Transfer and Fluid Dynamics. Boca Raton: CRC Press, 2001
65	J M Bergheau, R Fortunier. Finite Element Simulation of Heat Transfer. Hoboken: John Wiley & Sons, 2013
66	W C Carter. Lecture Notes on Mathematics for Materials Science and Engineers. MIT 3.016, 2012
67	D Ho Lee, B Man Kwak. Shape sensitivity and optimization for transient heat diffusion problems using the BEM. International Journal of Numerical Methods for Heat & Fluid Flow, 1995, 5(4): 313–326 https://doi.org/10.1108/EUM0000000004068

Viewed

Full text

Abstract

Cited

Shared

Discussed