基于COMSOL的简化对流传热结构拓扑优化方法

提出了一种基于COMSOL平台的简化对流传热结构拓扑优化实现方法。简化对流传热拓扑优化可用于获得性能更优的散热结构概念设计,但缺乏可行性强和适用范围广的实现方法,导致工程应用存在较高的技术门槛。为解决该问题,本研究采用在经典传热方程中引入对流换热项的方式,建立了可用于COMSOL平台简化对流换热分析的控制方程。进一步的,结合COMSOL平台提供的传热分析模块以及拓扑优化模块,搭建了简化对流传热结构拓扑优化计算框架。数值算例考虑了经典二维和三维的简化对流传热结构优化设计问题。结果显示了本文方法的有效性及可行性。     

多孔材料/结构尺度关联的一体化拓扑优化技术

针对多孔材料的尺度效应和微结构构型的可设计性,提出以宏观结构最大刚度为目标,材料表征体胞构型为变量的材料／结构尺度关联一体化设计新方法．采用有限元超单元技术,验证了材料表征体胞尺度、边长比、平移、对称周期分布方式对构型设计结果的影响,实现了材料宏观布局设计、材料表征体胞构型精细设计以及多尺度均匀化设计的统一。基于凸规划对偶优化求解技术与二次型周长控制约束,完成了快速设计与材料分布棋盘格效应的控制。计算结果表明,在给定材料用量的情况下,该方法能有效地实现蜂窝结构的拓扑优化设计,设计结果充分反映了蜂窝夹层结构的尺度效应,为轻质结构设计提供了新的设计方法。     

基于仿真和数据驱动的先进结构材料设计

先进结构材料近年来受到材料和结构设计领域的广泛关注,这些材料一般通过多个尺度的结构设计实现各种卓越的性能.在早期的材料设计中,有的基于设计者的丰富经验,从天然拓扑结构中抽象出合理的数学力学模型;有的基于生物系统的结构和功能特点提取出仿生力学模型.然而,仅依靠经验性的巧妙设计很难得到最优的设计方案,通过反复迭代设计和试验...     

生物结构自适应性的力学研究及其应用

以生物的骨组织、竹子以及植物根系、动物血管等生物分枝系统为例, 从微观组织结构和宏观形态, 以力学的观点, 对生物形态组织结构的自适应性进行了分析研究, 并概括地介绍了这种自适应性在工程材料及工程结构设计中的应用, 如基于骨组织重建机理的自修复材料、基于生物分枝系统形态生成机理的薄板加强筋以及热传递系统中散热通道的分布设计等. 模拟生物结构和系统的工程材料和结构及其设计方法, 因具有与生物结构和系统相似的优越性能, 可突破传统产品和设计方法的一些局限性. 基于仿生技术的生物形态、材料与机能、强度等关系的进一步研究, 以及相关的理想工程材料和结构的设计与制造方法的研究, 将是当前和未来十分重要的研究课题.      

Concurrent Analysis and Design of Structure and Its Periodic Material

The specific good properties of cellular materials and composite materials, such as low density and high permeability, make the optimal design of such materials necessary and attractive. However, the given materials for the structures may not be optimal or suitable, since the boundary condition and applied loads vary in practical applications; hence the macro-structure and its material micro-structure should be considered simultaneously. Although abundant studies have been reported on the structural and material optimization at each level, very few of them considered the mutual coordination on both scales. In this paper, two FE models are built for the macro-structure and the micro-structure, respectively; and the effective elastic properties of the periodic micro-structure are blended into the analysis of macro-structure by the homogenization theory. Here, a topological optimum is obtained by gradually re-distributing the constituents within the micro-structure and updating the topological shape at the macro-structure until converges are achieved on both scales. The mutual coordination between the roles of micro-scale and macro-scale is considered. Some numerical examples are presented, which illustrate that the proposed optimization algorithm is effective and highly efficient for the micro-structure design and macro-structure optimization. For the composite design, one can see reasonable effects of the stiffness of base materials on the resultant topologies.     

Planar solidification of a warm flowing-liquid under different boundary conditions

Planar solidification of a warm flowing liquid with the convective heat transfer to the growing solid layer, has been analysed for the boundary conditions of constant temperature, constant heat flux and convective heat flux at the surface respectively. The mathematical formulation of the problem resulted in a coupled set of two differential equations in temperature and solid thickness as function of position, time and the problem parameters. Analytical expressions for the temperature distribution within the growing solid layer, the rate of solidification and the solidification time are obtained. The perturbation techniques employed here is simple and straight forward in contrast with the earlier techniques. Good agreement between the experimental results and the present solutions is obtained for the convective heat flux boundary condition. The results of this analysis are useful in the design and analysis of experiments dealing with freezing/melting in one dimension. The role of the parameter Stefan number which is small for phase change materials, is discussed in context with the storage of thermal energy.     

Bi-functional optimization of actively cooled, pressurized hollow sandwich cylinders with prismatic cores

All metallic, hollow sandwich cylinders having ultralight two-dimensional (2D) prismatic cores are optimally designed for maximum thermo-mechanical performance at minimum mass. The heated cylinder is subjected to uniform internal pressure and actively cooled by forced air convection. The use of two different core topologies is exploited: square- and triangular-celled cores. The minimum mass design model is so defined that three failure modes are prevented: facesheet yielding, core member yielding, and core member buckling. The intersection-of-asymptotes method, in conjunction with the fin analogy model, is employed to build the optimization model for maximum heat transfer rate. A non-dimensional parameter is introduced to couple the two objectives—structural and thermal—in a single cost function. It is found that the geometry corresponding to maximum heat transfer rate is not unique, and square-celled core sandwich cylinders outperform those having triangular cells. The eight-layered sandwich cylinders with square cells have the best overall performance in comparison with other core topologies. Whilst a sandwich cylinder with shorter length is preferred for enhanced thermo-mechanical performance, the influence of the outer radius of the cylinder is rather weak.     

Size effects of basic cell in static analysis of sandwich beams

In this paper, multilayered sandwich beam structures are considered. Within the scope of static analyses and stiffness design of such type of lightweight and functional structures, size effects of the basic cell are studied both theoretically and numerically in a systematic way for the first time. The direct FE discretization method, the homogenization method and the classical beam theory are examined systematically to reveal, on one hand, the existence of the size effect, and on the other hand, the ability of each method in capturing the size effect upon the static stress distribution and structural deflection. Particularly, limitations of the homogenization method are clarified although the latter is widely applied today in the equivalent modeling and topology design of cellular materials of sandwich structures. By means of the above methods, bending problems of multilayered beams and cellular core sandwiches are solved to illustrate variations of the deflection, stress as well as the computing accuracies in terms of the size of the basic cell. It is shown that the size effect is important when the basic cell has a considerable dimension relative to the structural size and that this effect decreases rapidly with the size reduction of the basic cell. Theoretically, the homogenized result corresponds to the limit solution when the size of the basic cell tends to be infinitely small.     

Topological design of structures and composite materials with multiobjectives

This paper studies the influence of heat conduction in both structural and material designs in two dimensions. The former attempts to find the optimal structures with the maximum stiffness and minimum resistance to heat dissipation and the latter to tailor composite materials with effective thermal conductivity and bulk modulus attaining their upper limits like Hashin–Shtrikman and Lurie–Cherkaev bounds. In the part of structural topology optimization of this paper solid material and void are considered respectively. While in the part of material design, two-phase ill-ordered base materials (i.e. one has a higher Young’s modulus, but lower thermal conductivity while another has a lower Young’s modulus but higher conductivity) are assumed in order to observe competition in the phase distribution defined by stiffness and conduction. The effective properties are derived from the homogenization method with periodic boundary conditions within a representative element (base cell). All the issues are transformed to the minimization problems subject to volume and symmetry constraints mathematically and solved by the method of moving asymptote (MMA), which is guided by the sensitivities with respect to the design variables. To regularize the problem the SIMP model is explored with the nonlinear diffusion techniques to create edge-preserving and checkerboard-free results. The illustrative examples show how to generate Pareto fronts by means of linear weighting functions, which provide an in-depth understanding how these objectives compete in the topologies.     

Free Convection Boundary Layer Flow from a Heated Upward Facing Horizontal Flat Plate Embedded in a Porous Medium Filled by a Nanofluid with Convective Boundary Condition

The steady laminar incompressible free convective flow of a nanofluid over a permeable upward facing horizontal plate located in porous medium taking into account the thermal convective boundary condition is studied numerically. The nanofluid model used involves the effect of Brownian motion and the thermophoresis. Using similarity transformations the continuity, the momentum, the energy, and the nanoparticle volume fraction equations are transformed into a set of coupled similarity equations, before being solved numerically, by an implicit finite difference numerical method. Our analysis reveals that for a true similarity solution, the convective heat transfer coefficient related with the hot fluid and the mass transfer velocity must be proportional to x ^−2/3, where x is the horizontal distance along the plate from the origin. Effects of the various parameters on the dimensionless longitudinal velocity, the temperature, the nanoparticle volume fraction, as well as on the rate of heat transfer and the rate of nanoparticle volume fraction have been presented graphically and discussed. It is found that Lewis number, the Brownian motion, and the convective heat transfer parameters increase the heat transfer rate whilst the thermophoresis decreases the heat transfer rate. It is also found that Lewis number and the convective heat transfer parameter enhance the nanoparticle volume fraction rate whilst the thermophoresis parameter decreases nanoparticle volume fraction rate. A very good agreement is found between numerical results of the present article for special case and published results. This close agreement supports the validity of our analysis and the accuracy of the numerical computations.     

Numerical study of natural convection flow in a vertical slot

The stability of convective motion, generated by a lateral temperature difference across a vertical slot, is studied numerically over a range ofGr=5000 to 1.5 × 10⁵,Pr=0.01 to 10, andA=8,16 and 20. Various cellular flow structures and temperature patterns are illustrated. Several branches of solutions characterized by different numbers of the cells in the flow patterns as well as by both steady and unsteady multicellular patterns are found for low-Prandtl-number fluid in the vertical slot. Meanwhile, the behaviors of the temperature variation and heat transfer are also discussed. The project supported by the National Natural Science Foundation of China (59776011) and by the Returnee from Abroad Funding of Academia Sinica.     

Microstructure and Tensile Mechanical Properties of Anisotropic Rigid Polyurethane Foam

An understanding of the mechanical properties of solid foams facilitates effective use of such materials, which are often deployed in load-bearing applications such as impact absorbers, cushioning and sandwich structures. This study is an experimental investigation that focuses on the deformation response of rigid polyurethane foam to tension. Microstructural features such as the size and geometry of constituent cells and the struts that define the cell edges, as well as their stiffness and tensile strength, are examined. The nature of cell deformation and fracture are identified through microscopy and the associated micromechanics analyzed. Results show that the cells are elongated and thus the foam exhibits anisotropic properties. Flexure of the struts that define the cell edges is the primary mechanism governing deformation and failure. Information on the mechanical, microstructural and deformation characteristics elicited through this investigation will facilitate formulation of idealized cell-based models to account for the mechanical response of rigid polymeric foams.     

????????????????????????

超轻多孔金属材料因其具有优良的结构效率和广泛的多功能应用前景,近年来已成为国内外 众多研究所关注的热点. 本文就学术界关注的多孔金属材料的宏观等效模型建模与材料微 结构的多功能优化设计进行了综述.     

机械密封动环外周表面织构换热机理及结构优化

本文中以接触式机械密封为研究对象,考虑端面摩擦热,建立了包含密封环及密封腔的三维轴对称传热模型,并采用Fluent软件进行了传热分析,得到了密封环及密封腔的流场、温度场、速度场和Nusselt Number(Nu数)分布情况. 通过分析织构及其周边的流场、流态和温度场,揭示了织构的对流换热机理,主要研究了等边三角形织构的深径比、排数、排间距、排列方式和方向角等几何因素对端面温度的影响,并对等边三角形的结构参数进行了优选. 结果表明:动环外周表面织构具有增强对流换热效果和降低密封环整体温度的能力;织构的型式、深径比、个数以及排列方式均会对其换热效果产生影响,其中织构个数和旋转角影响较大,且三角形织构具有较强的换热能力,为今后高性能机械密封的优化设计及应用提供了理论依据.      

圆柱夹层多孔材料的主动散热性能研究与优化设计

通过推导不同边界条件的圆柱状夹层多孔材料散热指标,研究了一种特殊结构的圆柱状夹层多孔材料主动散热问题。这种圆柱夹层材料的每一层胞体个数相同,胞体尺寸随着外径的增大而增大,从而保持每一层胞体的相对密度相同。通过计算两种不同换热边界条件下圆柱夹层多孔材料的散热性能,比较并分析了与夹层材料层数相对应的最大散热效率和最优相对密度等指标,并最终得到这些工况下的最优质量。通过分析得到,无论哪种边界条件,正六边形胞体的夹层结构散热性能优于其他构型。同时当为达到某一限定散热效率值时,正六边形胞体结构的质量最小,正六边形构型的多孔材料具有明显的综合性能优势。     

Thermoelastic Stress Analysis of Structures Under Natural Vibrations

The paper focuses on stress analyses of structures subjected to excitation forces operating at resonant frequencies. The structures are analysed experimentally using the Thermoelastic Stress Analysis (TSA) technique. Experiments are carried out for fixed-free beams of different dimensions and materials, and also for a steel rectangular plate with clamped edges. These structures are excited by a shaker via a stinger. For materials with low thermal conductivity, the agreement between the theory, numerical results and experimental results is excellent. As the thermal conductivity of the material is increased, the correspondence is not as close. This is because of non-adiabatic behaviour. The implications of these results are discussed in detail in the paper and a means of deriving the severity of heat transfer is provided. Other factors that influence the TSA results from structures under natural loading are also discussed.     

Experimental investigation of the coherent structures in a spirally corrugated pipe

Previous numerical and theoretical results (Chen et al., 2019; Liu et al., 2018; Zhao et al., 2019) based on the optimization theory of convective heat transfer reveal that the optimized flow structures in a straight circular pipe enhancing convective heat transfer are multiple longitudinal vortices. This conclusion encourages us to find out whether such flow structures really exist in some enhanced heat transfer pipes by means of advanced experimental techniques. Therefore, a typical enhanced heat transfer pipe was selected, namely a spirally corrugated pipe, and stereoscopic particle image velocimetry (SPIV) was employed to measure its internal instantaneous flow field. Moreover, the proper orthogonal decomposition (POD) method was used to extract the large-scale coherent structures from the measured instantaneous velocity fields. Besides the spirally corrugated pipe, the fully developed turbulent flow in a straight pipe was also analyzed as benchmark of the enhanced heat transfer pipes. The results reveal that longitudinal whirling flow with multi-vortices is formed in both the fully developed turbulent flow field of the straight pipe and the spirally corrugated one. It is thus easy to explain the heat transfer enhancement mechanism of the above flow structures from the perspective of momentum transfer. The flow structures of the fully developed turbulent flow in a straight pipe are quite similar to the optimal flow pattern from the optimization theory. More specifically, multiple longitudinal vortices are spontaneously generated due to turbulence without external heat transfer enhancement techniques. Furthermore, the flow structures similar to multiple longitudinal vortices also exist in the spirally corrugated pipe, although these flow structures deviate from symmetric multiple vortices. Moreover, the flow structures in the spirally corrugated pipe are much more energetic than those in the fully developed turbulent flow in a straight pipe. This is probably the reason why a spirally corrugated pipe can enhance heat transfer compared with a straight circular pipe.     

A Two-Equation Analysis of Convection Heat Transfer in Porous Media

This paper presents a two-equation analysis on the convection heat transfer in porous media based on the modeling developed by Carbonell and Whitaker (1984). The porous system under consideration is bounded by two parallel walls and heated uniformly from one side surface. The Darcy flow is imposed and the fully developed heat transfer is assumed. General solutions, which take into account the additional convective and conductive terms, are obtained for the temperature fields and the Nusselt number. The detailed studies are presented for the porous systems characterized by consolidated and unconsolidated circular unit cells. The results show that, for the consolidated unit cell case, a prediction without the additional convective term overestimates the heat transfer, while for the unconsolidated unit cell case, this effect is negligible. The additional conductive terms are also examined and found to act conventionally as part of the conductive terms.     

Laminar convective heat transfer of non-Newtonian nanofluids with constant wall temperature

Nanofluids are obtained by dispersing homogeneously nanoparticles into a base fluid. Nanofluids often exhibit higher heat transfer rate in comparison with the base fluid. In the present study, forced convection heat transfer under laminar flow conditions was investigated experimentally for three types of non-Newtonian nanofluids in a circular tube with constant wall temperature. CMC solution was used as the base fluid and γ-Al₂O₃, TiO₂ and CuO nanoparticles were homogeneously dispersed to create nanodispersions of different concentrations. Nanofluids as well as the base fluid show shear thinning (pseudoplastic) rheological behavior. Results show that the presence of nanoparticles increases the convective heat transfer of the nanodispersions in comparison with the base fluid. The convective heat transfer enhancement is more significant when both the Peclet number and the nanoparticle concentration are increased. The increase in convective heat transfer is higher than the increase caused by the augmentation of the effective thermal conductivity.     

A variational asymptotic micromechanics model for predicting thermoelastic properties of heterogeneous materials

A variational asymptotic micromechanics model has been developed for predicting effective thermoelastic properties of composite materials, and recover the local fields within the unit cell. This theory adopts essential assumptions within the concept of micromechanics, achieves an excellent accuracy, and provides a unified treatment for 1D, 2D, and 3D unit cells. This theory is implemented using the finite element method into the computer program, VAMUCH, a general-purpose micromechanics analysis code. Several examples are used to validate the theory and the code. The results are compared with those available in the literature and those produced by a commercial finite element package.