Research Papers: Design Automation

Thermal Performance Optimization in Electric Vehicle Power Trains by Locally Orthotropic Surface Layer Design

[+] Author and Article Information
Mario Petrovic

Department of Mechanical Engineering and
Science Graduate School of Engineering,
Kyoto University,
Kyoto 615-8540, Japan
e-mail: petrovic.mario.7e@kyoto-u.ac.jp

Tsuyoshi Nomura

Toyota Central R & D Labs,
Nagakute 480-1192, Japan
e-mail: nomu2@mosk.tytlabs.co.jp

Takayuki Yamada

Department of Mechanical
Engineering and Science,
Graduate School of Engineering,
Kyoto University,
Kyoto 615-8540, Japan
e-mail: takayuki@me.kyoto-u.ac.jp

Kazuhiro Izui

Department of Mechanical
Engineering and Science,
Graduate School of Engineering,
Kyoto University,
Kyoto 615-8540, Japan
e-mail: izui@me.kyoto-u.ac.jp

Shinji Nishiwaki

Department of Mechanical
Engineering and Science,
Graduate School of Engineering,
Kyoto University,
Kyoto 615-8540, Japan
e-mail: shinji@prec.kyoto-u.ac.jp

Contributed by the Design Automation Committee of ASME for publication in the JOURNAL OF MECHANICAL DESIGN. Manuscript received October 10, 2017; final manuscript received August 8, 2018; published online September 12, 2018. Assoc. Editor: James K. Guest.

J. Mech. Des 140(11), 111413 (Sep 12, 2018) (8 pages) Paper No: MD-17-1681; doi: 10.1115/1.4041220 History: Received October 10, 2017; Revised August 08, 2018

In this paper, the application of orthotropic material orientation optimization for controlling heat flow in electric car power trains is presented. The design process is applied to a case model, which conducts heat while storing heat-sensitive electronic components. The core of the case is designed using a low thermal conductivity material on order to focus the heat flow into the surface layer, which is designed using a high thermal conductivity material. Material orthotropy is achieved in the surface layer of the case by removing the material at points determined by the optimization analysis. For this purpose, an orthotropic material orientation optimization method was extended to calculate optimal material distribution. This is achieved by transforming the initially obtained optimal orientation vector field into a scalar field through the use of coupled time-dependent nonisotropic Helmholtz equations. Multiple parameters allow the control of the scalar field and therefore the control over material distribution in accordance to the optimal orientation. This allows the material distribution pattern to be scaled depending on the desired manufacturing method. The analysis method is applied to divert heat flow from a specific section of the model while focusing the heat flow to another section. The results are shown for a model with a 0.1 mm thick surface layer of copper and are compared to those results from several other materials and layer thicknesses. Finally, the manufactured design is presented.

Copyright © 2018 by ASME
Your Session has timed out. Please sign back in to continue.


Suzuki, K. , 1991, “ A Homogenization Method for Shape and Topology Optimization,” Comput. Methods Appl. Mech. Eng., 93(3), pp. 291–318. [CrossRef]
Gea, H. , and Luo, J. , 2004, “ On the Stress-Based and Strain-Based Methods for Predicting Optimal Orientation of Orthotropic Materials,” Struct. Multidiscip. Optim., 26(3–4), pp. 229–234. [CrossRef]
Hyer, M. W. , and Lee, H. H. , 1991, “ The Use of Curvilinear Fiber Format to Improve Buckling Resistance of Composite Plates With Central Circular Holes,” Compos. Struct., 18(3), pp. 239–261. [CrossRef]
IJsselmuiden, S. T. , Abdalla, M. M. , and Gürdal, Z. , 2010, “ Thermomechanical Design Optimization of Variable Stiffness Composite Panels for Buckling,” J. Therm. Stresses, 33(10), pp. 977–992. [CrossRef]
Peeters, D. M. J. , Hesse, S. , and Abdalla, M. M. , 2015, “ Stacking Sequence Optimisation of Variable Stiffness Laminates With Manufacturing Constraints,” Compos. Struct., 125, pp. 596–604. [CrossRef]
Li, Q. , Steven, G. P. , Querin, O. M. , and Xie, Y. M. , 1990, “ Shape and Topology Design for Heat Conduction by Evolutionary Structural Optimization,” Int. J. Heat Mass Transfer, 42(17), pp. 3361–3371. [CrossRef]
Gersborg-Hansen, A. , Bendsøe, M. P. , and Sigmund, O. , 2006, “ Topology Optimization of Heat Conduction Problems Using the Finite Volume Method,” Struct. Multidiscip. Optim., 31(4), pp. 251–259. [CrossRef]
Sigmund, O. , and Torquato, S. , 1999, “ Design of Smart Composite Materials Using Topology Optimization,” Smart Mater. Struct., 8(3), pp. 365–379. [CrossRef]
Kruijf, N. , Zhou, S. , Li, Q. , and Mai, Y.-W. , 2007, “ Topological Design of Structures and Composite Materials With Multiobjectives,” Int. J. Solids Struct., 44(22–23), pp. 7092 –7109. [CrossRef]
Dede, E. M. , 2010, “ Simulation and Optimization of Heat Flow Via Anisotropic Material Thermal Conductivity,” Comput. Mater. Sci., 50(2), pp. 510–515. [CrossRef]
Nomura, T. , Dede, E. M. , Lee, J. , Yamasaki, S. , Matsumori, T. , Kawamoto, A. , and Kikuchi, N. , 2015, “ General Topology Optimization Method With Continuous and Discrete Orientation Design Using Isoparametric Projection,” Int. J. Numer. Methods Eng., 101(8), pp. 571–605. [CrossRef]
Petrovic, M. , Nomura, T. , Yamada, T. , Izui, K. , and Nishiwaki, S. , 2017, “ Orthotropic Material Orientation Optimization Method in Composite Laminates,” Struct. Multidiscip. Optim., 57(2), pp. 815–828. [CrossRef]
Hiroshi, H. , and Minoru, T. , 1986, “ Equivalent Inclusion Method for Steady-State Heat Conduction in Composites,” Int. J. Eng. Sci., 24(7), pp. 1159 –1172. [CrossRef]
Turing, A. M. , 1952, “ The Chemical Basis of Morphogenesis,” Philos. Trans. R. Soc. London Ser. B, 237(641)., pp. 37–72. [CrossRef]
Svanberg, K. , 1987, “ The Method of Moving Asymptotes–a New Method for Structural Optimization,” Int. J. Numer. Methods Eng., 24(2), pp. 359–373. [CrossRef]


Grahic Jump Location
Fig. 2

Surface layer local material distribution

Grahic Jump Location
Fig. 3

Example for material distribution calculation: (a) results of Eq. (16) and (b) result of Eqs. (18) and (19)

Grahic Jump Location
Fig. 4

Case model: top side and bottom side views

Grahic Jump Location
Fig. 5

Temperature values: isotropic surface layer

Grahic Jump Location
Fig. 6

Temperature values: optimized orthotropic surface layer

Grahic Jump Location
Fig. 7

Local temperature values: (a) top side view and (b) bottom side view

Grahic Jump Location
Fig. 8

Comparison of temperature values for different surface layers

Grahic Jump Location
Fig. 9

Surface layer material distribution

Grahic Jump Location
Fig. 10

Local surface layer material distribution: (a) top side view and (b) bottom side view

Grahic Jump Location
Fig. 11

Manufactured prototype: (a) top side view and (b) bottom side view



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In