阅读说明：本技术 一种泡沫材料的多尺度结构反向设计方法 (Reverse design method for multi-scale structure of foam material ) 是由 李洋 林长华 陈红伟 苏文涛 马贵阳 于 2021-07-08 设计创作，主要内容包括：本发明公开了一种泡沫材料的多尺度结构反向设计方法,属于泡沫材料制造技术领域。本发明解决现有陶瓷或金属泡沫材料多尺度结构与其高温热应用的匹配、控制、优化等关键问题。本发明以泡沫材料高温应用中的辐射性能需求为牵引,基于光谱辐射特性信息的跨尺度传递规律,开展泡沫材料内辐射能量的跨尺度梯级调控机理与方法研究,依托3D打印等技术开展典型泡沫材料的多尺度结构反向设计理论与方法研究,突破相关高温热利用系统设计、优化和过程控制中的热辐射技术瓶颈。(The invention discloses a reverse design method for a multi-scale structure of a foam material, and belongs to the technical field of foam material manufacturing. The invention solves the key problems of matching, controlling, optimizing and the like of the multi-scale structure of the existing ceramic or metal foam material and the high-temperature heat application thereof. The invention takes the radiation performance requirement in the high-temperature application of the foam material as traction, develops the research of the cross-scale gradient regulation and control mechanism and method of radiation energy in the foam material based on the cross-scale transmission rule of spectral radiation characteristic information, develops the research of the multi-scale structure reverse design theory and method of typical foam materials by relying on the technologies such as 3D printing and the like, and breaks through the bottleneck of the heat radiation technology in the design, optimization and process control of related high-temperature heat utilization systems.)

1. A method for reversely designing a multi-scale structure of a foam material is characterized by comprising the following steps:

step 1, scanning by adopting an SEM (scanning Electron microscope) technology and a mu-CT (micro-computed tomography) technology to obtain pore morphology characteristics of the foam material on macro scale, pore scale and micro scale, and establishing a multi-scale structural parameter database of the foam material;

step 2, developing mathematical description of a multi-scale structure of the foam material, establishing a micro-scale-pore scale, pore scale-macro scale radiation characteristic transmission data/correlation, constructing a cross-scale radiation characteristic database of a typical foam material corresponding to the multi-scale structure, programming a user Application Program Interface (API) based on the mathematical description of the multi-scale structure, and controlling a modeling software SolidWorks to perform automatic simulation reconstruction of the multi-scale structure of the foam material;

step 3, performing three-level step classification on the structural parameters from the macro-scale, pore-scale and micro-scale levels of the foam material respectively;

step 4, determining the priority of radiation energy regulation and control based on the structural parameter step classification by combining the radiation performance requirements of the foam material in the actual engineering, and establishing a cross-scale step regulation and control method of the spectral radiation energy of the foam material;

step 5, reversely designing a multi-scale structure of the foam material with specific radiation performance aiming at different engineering practical application scenes of the foam material on the basis of the cross-scale radiation characteristic of the foam material and a cross-scale step regulation and control method;

step 6, combining the 3D printing technology and the thermal forming characteristics of the base material to adjust the multi-scale structural parameters of the foam material reversely designed in the step 5 until the actual application scene of the engineering is met;

and 7, programming an application program interface API of the user according to the structural parameters adjusted in the step 6, controlling the modeling software SolidWorks to carry out automatic simulation reconstruction on the multi-scale structure of the foam material, and exporting STL format data for 3D printing to obtain the foam material meeting the actual application scene of the project.

2. The method of claim 1, wherein the multi-scale structure parameter database of the foam material in step 1 comprises:

appearance data: length L, width W, and height D;

pore data: porosity phi, cell diameter d_cRib shape control parameters;

the rib shape control parameter comprises the average cell diameter d of the rib framework₀The rib longitudinal shape parameter t, the rib section shape parameter k and the rib hollowness parameter h;

③ microcosmic data: surface of ribRoughness R_aAnd R_yPore diameter d of micro-pores in the rib_p0Rib porosity phi₀。

3. The method for reverse design of multi-scale structure of foam material as claimed in claim 2, wherein the mathematical description of the foam rib shape control parameters is established by:

(I) porosity phi and cell diameter d_cDetermine the average cell diameter d of the rib skeleton₀：

In the formula, N_cIs the total number of cells; n is a radical of₀The total number of rib skeletons; n is a radical of_vIs the node total number of the non-repeated ribs; l is_0,jIs the length of the jth rib cage;

(II) longitudinal rib shape parameter t: the size change of the rib along the length direction of the rib is characterized, namely the degree of the rib deviating from the equal-diameter cylinder, and the longitudinal shape parameter t is defined as follows:

in the formula (d)_minIs the diameter of the thinnest part of the middle section of the rib; d_maxIs the diameter of the thickest part at the two ends of the rib;

assuming that the outer edge of the longitudinal section conforms to the secondary distribution, the change of the rib diameter d along the length direction of the rib can be expressed as:

in the formula, L₀Is the length of the rib; taking the center point of the rib as the origin of coordinates, and taking the coordinate value of the point as-0.5L₀≤l≤0.5L₀；

At this time, the volume V of the single deformed rib₀' is represented as:

while the equivalent mean volume V of a single rib₀Expressed as:

simultaneous equations (4) and (5) can be obtained:

simultaneous equations (2) and (6) can be obtained:

in the formula (d)_maxIs the diameter of the thickest part at the two ends of the rib; average cell diameter d₀(ii) a A longitudinal shape parameter t;

by constructing a curvature equation of the outer edge of the longitudinal section of the rib and substituting the point coordinates determined by the formulas (6) and (7), the curvature radius r of the outer edge of the longitudinal section of the rib can be obtained as follows:

in the formula (d)_maxIs the diameter of the thickest part at the two ends of the rib; d_minIs the diameter of the thinnest part of the middle section of the rib; l is₀Is the length of the rib;

(III) rib section shape parameter k: representing the change of the section shape of the rib, namely the deviation degree of the section of the rib from the circular shape;

the cross-sectional shape parameter k is defined as:

in the formula, R is the radius of a circumscribed circle of an equilateral triangle of the section of the rib; r is the curvature radius of a circumscribed or inscribed arc on one side of the equilateral triangle of the section of the rib,indicating the concave cross section, k>0 represents a convex cross section;

as can be seen from the formula (9), the radius R of the circumscribed circle of the equilateral triangle of the rib section is kr, and the side length a of the equilateral triangle of the rib section is

The central angle alpha corresponding to the external or internal circular arc on one side of the equilateral triangle of the section of the rib is as follows:

sector area S corresponding to external or internal circular arc on one side of equilateral triangle of rib section_fanComprises the following steps:

triangle area S corresponding to external or internal circular arc on one side of rib section equilateral triangle_fan-ΔComprises the following steps:

the simultaneous formulas (11) and (12) can be used for corresponding to a circumscribed or inscribed circular arc on one side of an equilateral triangle with a rib sectionArcuate area S of_bowComprises the following steps:

section area S of rib₀Comprises the following steps:

at this time, the volume V of the single deformed rib₀Expressed as:

V₀＝S₀L₀ (15)

the curvature radius r of the outer edge of the section of the rib obtained by the simultaneous formulas (5) and (15) is as follows:

in the formula (d)₀Is the average cell diameter; k is a section shape parameter;

(IV) rib hollowness parameter h: characterizing the hollowness degree of the rib;

the rib hollowness parameter h is defined as:

in the formula, S_hIs a hollow cross-sectional area; s₀Is the total area of the cross section;

the hollow cross-sectional area S can be obtained from the formula (17)_hComprises the following steps:

order S_h＝S₀Substituting the formula (14) can obtain the outer edge curvature radius of the concave hollow sectionr is:

in the formula, k<0，d₀Is the average cell diameter; k is a section shape parameter; h is the degree of hollowness.

4. The method for reversely designing the multi-scale structure of the foam material according to claim 3, wherein BRDF, direction-hemisphere reflectance and specular reflectance data of a representative microstructure on the surface of a rib framework are obtained by simulation through an FDTD method, and radiation characteristic transfer data of a micro-scale-hole scale are obtained; if the foam material is semitransparent, acquiring equivalent attenuation coefficient, scattering albedo and scattering phase function data of microporosities/particle groups in the rib framework; the radiation characteristic data obtained from the micro-scale simulation is transmitted to the hole scale simulation for use, and is used for assigning the rib radiation characteristic, and the obtained micro-scale-hole scale radiation characteristic transmission data/correlation in the step 2 is as follows:

β₀＝-6.4339λ³+33.571λ²-60.097λ+40.893 (20)

in the formula: beta is a₀Is the foam material attenuation coefficient, and λ is the spectral wavelength;

ω₀＝0.0004λ³-0.0018λ²+0.0011λ+0.9997 (21)

in the formula: omega₀Is the scattering reflectance of the foam material and λ is the spectral wavelength.

5. The method for reverse design of a multi-scale structure of a foam material according to claim 3, wherein based on a radiation transmission free path and a scattering distribution statistical model, MCRT method or DO method is adopted to simulate and obtain equivalent attenuation coefficient, scattering albedo and scattering phase function data of a foam pore simulation structure and a mu-CT scanning structure, and radiation characteristic transfer data of pore scale-macro scale in step 2 is obtained; the radiation characteristic data obtained from the pore scale simulation is transmitted to a macro scale simulation for use, and is used for medium radiation characteristic assignment of a macro equivalent medium, and the obtained pore scale-macro scale radiation characteristic transmission data/correlation in the step 2 is as follows:

in the formula: beta is the attenuation coefficient of the foam material,is porosity, d_cThe average cell diameter is shown, t is a rib longitudinal shape parameter, and k is a rib section shape parameter;

in the formula: omega is the scattering albedo, tau, of the foam₀The optical thickness of the rib is the optical thickness of the rib,β₀is the attenuation coefficient of the substrate; d₀Is the average rib diameter; h is the rib hollowness; omega₀Is the scattering albedo, n, of the substrate₀The refractive index of the base material, t, and k are the longitudinal shape parameters of the rib and the section shape parameters of the rib;

in the formula: g is the asymmetry factor of the foam material, t is the longitudinal shape parameter of the rib, k is the section shape parameter of the rib, tau₀The optical thickness of the rib is the optical thickness of the rib,ω₀is the scattering albedo, n, of the substrate₀Is the refractive index of the substrate.

6. The method of claim 2, wherein the three-level steps of the structural parameters in step 3 are classified as: the first level is appearance data; the second level is pore data; the third level is microscopic data.

7. The method as claimed in claim 5, wherein the step 4 of controlling the spectral radiant energy of the foam material in a cross-scale step is to sequentially change the first-stage, second-stage and third-stage structural parameters of the foam material to control the radiant properties of the foam material.

8. The method of claim 1, wherein the foam substrate is metal or ceramic.

Technical Field

The invention relates to a reverse design method for a multi-scale structure of a foam material, and belongs to the technical field of foam material manufacturing.

Background

The porous foam materials represented by ceramic foam, metal foam and the like have typical three-dimensional multi-scale pore structure characteristics and show excellent characteristics of designability, high temperature resistance, corrosion resistance, thermal shock resistance, large specific surface area and the like. The complex thermal radiation transport mechanisms and laws in foam materials involved in related high temperature applications continue to be the leading topic of the academic field. The previous work mainly focuses on the research of the positive problem of light-heat-radiation transmission of the existing foam materials, and the functional foam materials with specific radiation performance are rarely designed reversely from the radiation characteristic requirements. The complexity of the structure of the foam material, the cross-scale effect of spectral radiation transmission, the lack of the design theory of the multi-scale structure of the foam material and the like cause great difficulty in the regulation and optimization of radiation energy, so that the radiation energy step distribution mechanism related to the regulation and control process is lack of deep understanding, and the reverse design theory and method of the multi-scale structure of the foam material are not mastered, thereby restricting the high-temperature application of the foam material and the development of related new technologies. Therefore, it is necessary to provide a multi-scale structure reverse design method for foam materials to break through the thermal radiation technology bottleneck in the design, optimization and process control of related high-temperature heat utilization systems.

Disclosure of Invention

The invention provides a reverse design method of a multi-scale structure of a foam material, aiming at solving the key problems of matching, control, optimization and the like of the multi-scale structure of the existing ceramic or metal foam material and high-temperature heat application of the multi-scale structure.

The technical scheme of the invention is as follows:

a method for reverse design of a multi-scale structure of a foam material, the method comprising the steps of:

step 1, scanning by adopting an SEM (scanning Electron microscope) technology and a mu-CT (micro-computed tomography) technology to obtain pore morphology characteristics of the foam material on macro scale, pore scale and micro scale, and establishing a multi-scale structural parameter database of the foam material;

step 2, developing mathematical description of a multi-scale structure of the foam material, establishing a micro-scale-pore scale, pore scale-macro scale radiation characteristic transmission data/correlation, constructing a cross-scale radiation characteristic database of a typical foam material corresponding to the multi-scale structure, programming a user Application Program Interface (API) based on the mathematical description of the multi-scale structure, and controlling a modeling software SolidWorks to perform automatic simulation reconstruction of the multi-scale structure of the foam material;

step 3, performing three-level step classification on the structural parameters from the macro-scale, pore-scale and micro-scale levels of the foam material respectively;

step 4, determining the priority of radiation energy regulation and control based on the structural parameter step classification by combining the radiation performance requirements of the foam material in the actual engineering, and establishing a cross-scale step regulation and control method of the spectral radiation energy of the foam material;

step 5, reversely designing a multi-scale structure of the foam material with specific radiation performance aiming at different engineering practical application scenes of the foam material on the basis of the cross-scale radiation characteristic of the foam material and a cross-scale step regulation and control method;

step 6, combining the 3D printing technology and the thermal forming characteristics of the base material to adjust the multi-scale structural parameters of the foam material reversely designed in the step 5 until the actual application scene of the engineering is met;

and 7, programming an application program interface API of the user according to the structural parameters adjusted in the step 6, controlling the modeling software SolidWorks to carry out automatic simulation reconstruction on the multi-scale structure of the foam material, and exporting STL format data for 3D printing to obtain the foam material meeting the actual application scene of the project.

Further defined, the multi-scale structural parameter database of the foam material in step 1 comprises:

appearance data: length L, width W, and height D;

pore data: porosity phi, cell diameter d_cRib shape control parameters;

the rib shape control parameter comprises the average cell diameter d of the rib framework₀The rib longitudinal shape parameter t, the rib section shape parameter k and the rib hollowness parameter h;

③ microcosmic data: surface roughness R of rib_aAnd R_yPore diameter d of micro-pores in the rib_p0Porosity of rib

Further defined, the process of establishing the mathematical description of the foam rib shape control parameters is as follows:

(I) porosity phi and cell diameter d_cDetermine the average cell diameter d of the rib skeleton₀：

In the formula, N_cIs the total number of cells; n is a radical of₀The total number of rib skeletons; n is a radical of_vIs the node total number of the non-repeated ribs; l is_0,jIs the length of the jth rib cage;

(II) longitudinal rib shape parameter t: the size change of the rib along the length direction of the rib is characterized, namely the degree of the rib deviating from the equal-diameter cylinder, and the longitudinal shape parameter t is defined as follows:

in the formula (d)_minIs the diameter of the thinnest part of the middle section of the rib; d_maxIs the diameter of the thickest part at the two ends of the rib;

assuming that the outer edge of the longitudinal section conforms to the secondary distribution, the change of the rib diameter d along the length direction of the rib can be expressed as:

in the formula, L₀Is the length of the rib; taking the center point of the rib as the origin of coordinates, and taking the coordinate value of the point as-0.5L₀≤l≤0.5L₀

At this time, volume V 'of the single deformed rib'₀Expressed as:

while the equivalent mean volume V of a single rib₀Expressed as:

simultaneous equations (4) and (5) can be obtained:

simultaneous equations (2) and (6) can be obtained:

in the formula (d)_maxIs the diameter of the thickest part at the two ends of the rib; average cell diameter d₀(ii) a A longitudinal shape parameter t;

by constructing a curvature equation of the outer edge of the longitudinal section of the rib and substituting the point coordinates determined by the formulas (6) and (7), the curvature radius r of the outer edge of the longitudinal section of the rib can be obtained as follows:

in the formula (d)_maxIs the diameter of the thickest part at the two ends of the rib; d_minIs the diameter of the thinnest part of the middle section of the rib; l is₀Is the length of the rib;

(III) rib section shape parameter k: representing the change of the section shape of the rib, namely the deviation degree of the section of the rib from the circular shape;

the cross-sectional shape parameter k is defined as:

in the formula, R is the radius of a circumscribed circle of an equilateral triangle of the section of the rib; r is rib sectionThe radius of curvature of the circumscribed or inscribed arc on one side of the equilateral triangle,k<0 denotes the concave cross section, k>0 represents a convex cross section;

as can be seen from the formula (9), the radius R of the circumscribed circle of the equilateral triangle of the rib section is kr, and the side length a of the equilateral triangle of the rib section is

The central angle alpha corresponding to the external or internal circular arc on one side of the equilateral triangle of the section of the rib is as follows:

sector area S corresponding to external or internal circular arc on one side of equilateral triangle of rib section_fanComprises the following steps:

triangle area S corresponding to external or internal circular arc on one side of rib section equilateral triangle_fan-ΔComprises the following steps:

the simultaneous formulas (11) and (12) can obtain the arch area S corresponding to the external or internal circular arc on one side of the equilateral triangle of the rib section_bowComprises the following steps:

section area S of rib₀Comprises the following steps:

at this time, the volume V of the single deformed rib₀Expressed as:

V₀＝S₀L₀ (15)

the curvature radius r of the outer edge of the section of the rib obtained by the simultaneous formulas (5) and (15) is as follows:

in the formula (d)₀Is the average cell diameter; k is a section shape parameter;

(IV) rib hollowness parameter h: characterizing the hollowness degree of the rib;

the rib hollowness parameter h is defined as:

in the formula, S_hIs a hollow cross-sectional area; s₀Is the total area of the cross section;

the hollow cross-sectional area S can be obtained from the formula (17)_hComprises the following steps:

order S_h＝S₀Substituting the formula (14) can obtain the outer edge curvature radius r of the concave hollow section as follows:

in the formula, k<0，d₀Is the average cell diameter; k is a section shape parameter; h is the degree of hollowness.

Further limiting, adopting an FDTD method to simulate and obtain BRDF, direction-hemisphere reflection ratio and specular reflection ratio data of a representative microstructure on the surface of the rib framework, and obtaining radiation characteristic transfer data of microscale-hole scale; if the foam material is semitransparent, acquiring equivalent attenuation coefficient, scattering albedo and scattering phase function data of microporosities/particle groups in the rib framework;

the radiation characteristic transfer data/correlation formula of the microscale-hole scale in the step 2 is as follows:

β₀＝-6.4339λ³+33.571λ²-60.097λ+40.893 (20)

in the formula: beta is a₀Is the foam material attenuation coefficient, and λ is the spectral wavelength;

ω₀＝0.0004λ³-0.0018λ²+0.0011λ+0.9997 (21)

in the formula: omega₀Is the scattering reflectance of the foam material and λ is the spectral wavelength.

Further limiting, based on a radiation transmission free path and a scattering distribution statistical model, adopting an MCRT method or a DO method to simulate and obtain equivalent attenuation coefficient, scattering albedo and scattering phase function data of a foam pore simulation structure and a mu-CT scanning structure, and obtaining the radiation characteristic transfer data of the pore scale-macro scale in the step 2;

the radiation characteristic transfer data/correlation of the meso-macro scale in the step 2 is as follows:

in the formula: beta is the attenuation coefficient of the foam material,is porosity, d_cThe average cell diameter is shown, t is a rib longitudinal shape parameter, and k is a rib section shape parameter;

in the formula: omega is a foam materialScattering albedo, τ of₀The optical thickness of the rib is the optical thickness of the rib,β₀is the attenuation coefficient of the substrate; d₀Is the average rib diameter; h is the rib hollowness; omega₀Is the scattering albedo, n, of the substrate₀The refractive index of the base material, t, and k are the longitudinal shape parameters of the rib and the section shape parameters of the rib;

in the formula: g is the asymmetry factor of the foam material, t is the longitudinal shape parameter of the rib, k is the section shape parameter of the rib, tau₀The optical thickness of the rib is the optical thickness of the rib,ω₀is the scattering albedo, n, of the substrate₀Is the refractive index of the substrate.

Further limiting, the three-level step classification of the structural parameters in step 3 is: the first level is appearance data; the second level is pore data; the third level is microscopic data.

Further defined, the three-level step classification of the structural parameters in step 3 is: the first level is appearance data; the second level is pore data; the third level is microscopic data.

Further defined, the foam substrate is a metal or a ceramic.

The invention has the following beneficial effects: the invention takes the radiation performance requirement in the high-temperature application of the foam material as traction, develops the research of the cross-scale gradient regulation and control mechanism and method of radiation energy in the foam material based on the cross-scale transmission rule of spectral radiation characteristic information, develops the research of the multi-scale structure reverse design theory and method of typical foam materials by relying on the technologies such as 3D printing and the like, and breaks through the bottleneck of the heat radiation technology in the design, optimization and process control of related high-temperature heat utilization systems. The reverse design method provided by the invention can meet the functional requirements of high-temperature application of the foam material, provides theoretical support for thermal design, optimization and process control of a high-temperature application system of the foam material, promotes rapid development of the field of novel material manufacturing represented by 3D printing and micro-nano scale processing technology, and can be applied to the technical fields of solar energy photo-thermal utilization, spacecraft thermal insulation, spacecraft thermonuclear propulsion, efficient combustion, electronic device cooling, enhanced heat exchange, high-temperature heat storage, high-temperature catalysis/thermochemical reaction and the like.

Drawings

FIG. 1 is a technical roadmap for the reverse engineering process of the present invention;

FIG. 2 is a schematic cross-scale transfer of radiation characteristic information within a foam material;

FIG. 3 is a foam structure of varying porosity;

FIG. 4 is a foam structure of different cell diameters;

FIG. 5 is a longitudinal shape configuration of the rib;

FIG. 6 shows rib-shaped structures with different longitudinal shape parameters t;

FIG. 7 is a cross-sectional configuration of a rib;

FIG. 8 is a cross-sectional structure of a rib with different cross-sectional shape parameters k;

FIG. 9 is a cross-sectional hollow structure of a rib;

FIG. 10 shows the cross-sectional structure of ribs with different hollowness parameters h;

Detailed Description

In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.

The experimental procedures used in the following examples are conventional unless otherwise specified. The materials, reagents, methods and apparatus used, unless otherwise specified, are conventional in the art and are commercially available to those skilled in the art.

Example 1:

characterizing and reconstructing a multi-scale (macro-scale, pore-scale and micro-scale) structure of a typical foam material (nickel foam, alumina foam, silicon carbide foam and carbon foam) by adopting an SEM (scanning electron microscope) technology and a mu-CT (micro-computed tomography) technology, and establishing a multi-scale structure database of the typical foam material;

the multi-scale structural parameter database of foam materials comprises:

appearance data: length L, width W, and height D;

pore data: porosity phi, cell diameter d_cRib shape control parameters;

the rib shape control parameter comprises the average cell diameter d of the rib framework₀The rib longitudinal shape parameter t, the rib section shape parameter k and the rib hollowness parameter h;

③ microcosmic data: surface roughness R of rib_aAnd R_yPore diameter d of micro-pores in the rib_p0Porosity of rib

Secondly, developing a mathematical description method of a multi-scale structure of the foam material, wherein the macro-scale of the foam plate is described by taking engineering application as guidance, the pore scale mainly quantitatively describes the characteristics of the porosity, the pore size, the pore cellular configuration (Lord Kelvin type, Weiire-Phelan type, Voronoi mosaic type and the like), the rib size, the longitudinal shape (spindle shape, straight column shape and the like) of the rib, the section shape (concave shape, convex shape and the like) of the rib, the hollow shape and the like of the rib and the like, and the micro-scale mainly considers the surface roughness of the rib, the particle diameter/distribution/packing density and the like in the rib;

the structure of the macro-scale and pore-scale is mathematically described as follows:

(I) as shown in FIGS. 3 and 4, the porosity φ and the cell diameter d_cDetermine the average cell diameter d of the rib skeleton₀：

In the formula, N_cIs the total number of cells; n is a radical of₀The total number of rib skeletons;N_vis the node total number of the non-repeated ribs; l is_0,jIs the length of the jth rib cage.

Theoretically, t ∈ (0, 1), and t ═ 1 indicates ribs of equal diameter, but t is not usually too small, and the ribs t ∈ [0.5,1] of common metal foams and ceramic foams are found statistically.

(II) longitudinal rib shape parameter t: characterizing the dimensional change of the rib along its length direction, i.e. the degree of deviation of the rib from the cylinder of constant diameter, as shown in fig. 5, the longitudinal shape parameter t is defined as:

in the formula (d)_minIs the diameter of the thinnest part of the middle section of the rib; d_maxIs the diameter of the thickest part at the two ends of the rib;

assuming that the outer edge of the longitudinal section conforms to the secondary distribution, the change of the rib diameter d along the length direction of the rib can be expressed as:

in the formula, L₀Is the length of the rib; taking the center point of the rib as the origin of coordinates, and taking the coordinate value of the point as-0.5L₀≤l≤0.5L₀；

At this time, volume V 'of the single deformed rib'₀Expressed as:

while the equivalent mean volume V of a single rib₀Expressed as:

simultaneous equations (4) and (5) can be obtained:

simultaneous equations (2) and (6) can be obtained:

in the formula (d)_maxIs the diameter of the thickest part at the two ends of the rib; average cell diameter d₀(ii) a A longitudinal shape parameter t;

by constructing the curvature equation of the outer edge of the longitudinal section of the rib (as shown by the black line in fig. 5) and substituting the point coordinates determined by the formulas (6) and (7), the curvature radius r of the outer edge of the longitudinal section of the rib can be obtained as follows:

in the formula (d)_maxIs the diameter of the thickest part at the two ends of the rib; d_minIs the diameter of the thinnest part of the middle section of the rib; l is₀Is the length of the rib.

To this end, the dimensional change of the rib along its length is characterized fully parametrically. d_min、d_maxR are basic input quantities required in parametric modeling, and can be calculated based on the longitudinal shape parameter t. Fig. 6 shows the rib shapes corresponding to different longitudinal shape parameters t, and it can be seen that the parameters t change the longitudinal shape of the rib significantly.

(III) rib section shape parameter k: representing the change of the section shape of the rib, namely the deviation degree of the section of the rib from the circular shape; the section of the actual foam rib is more non-circular, if the section of the polyurethane rib is more concave triangular, the section of the aluminum rib is more convex triangular, and the section of the aluminum oxide rib is close to circular, as shown in fig. 7, the definition of the section shape parameter k is as follows:

in the formula, R is the radius of a circumscribed circle of an equilateral triangle of the section of the rib; r is the curvature radius of a circumscribed or inscribed arc on one side of the equilateral triangle of the section of the rib.

In theory, it is possible to use,k<0 denotes the concave cross section, k>0 represents the convex cross section of the outer shell,the concave triangular section representing the maximum curvature, the constant-straight triangular section represented by k-0 and the circular section represented by k-1 hardly appear in the actual foam ribThe statistics show that the k E of the common foam ribs is [ -0.3,1 [ ]]。

As can be seen from the formula (9), the radius R of the circumscribed circle of the equilateral triangle of the rib section is kr, and the side length a of the equilateral triangle of the rib section is

The central angle alpha corresponding to the external or internal circular arc on one side of the equilateral triangle of the section of the rib is as follows:

sector area S corresponding to external or internal circular arc on one side of equilateral triangle of rib section_fanComprises the following steps:

triangle area S corresponding to external or internal circular arc on one side of rib section equilateral triangle_fan-ΔComprises the following steps:

the simultaneous formulas (11) and (12) can obtain the arch area S corresponding to the external or internal circular arc on one side of the equilateral triangle of the rib section_bowComprises the following steps:

section area S of rib₀Comprises the following steps:

at this time, the volume V of the single deformed rib₀Expressed as:

V₀＝S₀L₀ (15)

the curvature radius r of the outer edge of the section of the rib obtained by the simultaneous formulas (5) and (15) is as follows:

in the formula (d)₀Is the average cell diameter; k is a cross-sectional shape parameter.

To this end, the shape variation of the rib section is characterized by a complete parameterization. a. r are the basic input quantities required in parametric modeling, and they can all be calculated based on the section shape parameter k. Fig. 8 shows the rib cross-sectional shapes corresponding to different parameters k, and it can be seen that the parameter k significantly changes the rib cross-sectional shape.

(IV) rib hollowness parameter h: characterizing the hollowness degree of the rib; in practical foam ribs there may be cavities which have a negligible effect on the radiation transmission of opaque metal foam but a non-negligible effect on the radiation transmission of translucent ceramic foam, so that a parametric characterization of the foam rib void is necessary. As shown in fig. 9, the rib hollowness parameter h is defined as:

in the formula, S_hIs a hollow cross-sectional area; s₀Is the total area of the cross section.

Theoretically, h ∈ [0,1), h ═ 0 means no hollowness (solid); statistics shows that the hollowness of the real ribs does not exceed 0.3 generally, so the application takes h epsilon [0,0.3 ].

The hollow section of the ceramic rib mostly presents a concave triangular shape, so the shape parameter k defined by the formula (9) is adopted to describe the shape of the hollow section. The hollow cross-sectional area S can be obtained from the formula (17)_hComprises the following steps:

order S_h＝S₀Substituting the formula (14) can obtain the outer edge curvature radius r of the concave hollow section as follows:

in the formula, k<0，d₀Is the average cell diameter; k is a section shape parameter; h is the degree of hollowness.

To this end, the hollowness of the ribs is characterized parametrically. a. r are the basic input quantities required for parametric modeling, which can be derived from the hollow cross-sectional shape parameter k and the degree of hollowness h.

The statistics show that the hollow section shape parameter k of different ribs does not change greatly, and is about-0.25, and the parameters are obtained by substituting the following formula (19):

fig. 10 shows the rib hollow sections corresponding to different parameters h, and it can be seen that the parameter h changes the rib hollow significantly.

Therefore, the structural form of the foam rib is comprehensively represented by the longitudinal shape parameter t, the section shape parameter k and the hollowness parameter h.

Thirdly, simulating and acquiring BRDF (bidirectional reflectance-hemispherical reflectance), direction-hemispherical reflectance, specular reflectance and other data of the representative microstructure on the surface of the rib framework by adopting a relatively mature FDTD (fully drawn differential TD) method; if the foam rib substrate has translucency, characteristic data such as equivalent attenuation coefficient, scattering albedo and scattering phase function of internal micropores/particle groups of the foam rib substrate are required to be obtained. The radiation characteristic data obtained from the micro-scale simulation is transmitted to the hole scale simulation for use, and is used for assigning the rib radiation characteristic, and the obtained micro-scale-hole scale radiation characteristic transmission data/correlation formula is as follows:

β₀＝-6.4339λ³+33.571λ²-60.097λ+40.893 (20)

in the formula: beta is a₀Is the foam material attenuation coefficient, and λ is the spectral wavelength;

ω₀＝0.0004λ³-0.0018λ²+0.0011λ+0.9997 (21)

in the formula: omega₀Is the scattering reflectance of the foam material and λ is the spectral wavelength.

And fourthly, simulating and acquiring equivalent attenuation coefficient, scattering albedo and scattering phase function data of a foam pore simulation structure and a mu-CT scanning structure by adopting an MCRT method or a DO method with strong adaptability based on a radiation transmission free path and a scattering distribution statistical model, and combining the MCRT method and a space subdivision algorithm by considering the complexity and the randomness of the pore structure to accelerate the solving of a radiation transfer process and improve the calculation efficiency. The radiation characteristic data obtained from the hole scale simulation is transmitted to the macro scale simulation for use, and is used for medium radiation characteristic assignment of a macro equivalent medium, and the obtained hole scale-macro scale radiation characteristic transmission data/correlation formula is as follows:

in the formula: beta is the attenuation coefficient of the foam material,is porosity, d_cThe average cell diameter is shown, t is a rib longitudinal shape parameter, and k is a rib section shape parameter;

in the formula: omega is the scattering albedo, tau, of the foam₀The optical thickness of the rib is the optical thickness of the rib,β₀is the attenuation coefficient of the substrate; d₀Is the average rib diameter; h is the rib hollowness; omega₀Is the scattering albedo, n, of the substrate₀The refractive index of the base material, t, and k are the longitudinal shape parameters of the rib and the section shape parameters of the rib;

in the formula: g is the asymmetry factor of the foam material, t is the longitudinal shape parameter of the rib, k is the section shape parameter of the rib, tau₀The optical thickness of the rib is the optical thickness of the rib,ω₀is the scattering albedo, n, of the substrate₀Is the refractive index of the substrate.

And fifthly, classifying the radiation characteristic requirements of the foam material in high-temperature application into various situations and combination situations of strong absorption, strong scattering, gradual change of absorption/scattering, local strong absorption/strong scattering/low emission and the like, and using a cross-scale radiation characteristic database to iteratively optimize the radiation characteristic combination meeting the actual application scene of the engineering in turn from the macro-scale, pore-scale and micro-scale levels according to the priority order.

And sixthly, combining the limitations of the foam material on the conditions such as size, strength, pore space, weight and the like in practical application, reversely matching the cross-scale radiation characteristic parameters with the multi-scale structure data, and determining the corresponding foam material multi-scale structure parameters from the macro-scale, pore-scale and micro-scale levels in sequence.

And seventhly, fine adjustment is carried out on the structural parameters of the designed foam material by combining the 3D printing technology and the thermal forming characteristics of the base material (metal or ceramic), the radiation characteristic optimization evaluation function is called again, and the rest structural parameters are continuously adjusted and optimized to meet the actual application scene of the engineering.

Eighthly, according to the designed multi-scale structural scheme of the foam material, the foam material with the corresponding specification is manufactured by using a 3D printing technology, the multi-scale structural morphology data of the printed foam material is obtained by using SEM and mu-CT technologies, and the similarity degree of the obtained material structure and the design scheme is inspected.

Ninth, measuring the apparent spectral radiation characteristics of the foam material, mainly directional-directional transmittance, directional-directional reflectance and directional emission ratio data, and directly verifying whether the designed foam material has specific radiation performance.

The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.