CN112149310A - 一种空间质子辐照环境下材料表面brdf建模方法 - Google Patents

一种空间质子辐照环境下材料表面brdf建模方法 Download PDF

Info

Publication number
CN112149310A
CN112149310A CN202011074446.0A CN202011074446A CN112149310A CN 112149310 A CN112149310 A CN 112149310A CN 202011074446 A CN202011074446 A CN 202011074446A CN 112149310 A CN112149310 A CN 112149310A
Authority
CN
China
Prior art keywords
model
brdf
proton irradiation
space
material surface
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011074446.0A
Other languages
English (en)
Other versions
CN112149310B (zh
Inventor
侯晴宇
谭凡教
苏金宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN202011074446.0A priority Critical patent/CN112149310B/zh
Publication of CN112149310A publication Critical patent/CN112149310A/zh
Application granted granted Critical
Publication of CN112149310B publication Critical patent/CN112149310B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/26Composites
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Abstract

本发明公开了一种空间质子辐照环境下材料表面BRDF建模方法。步骤1:选择一种经过验证且没有考虑质子辐照的空间目标材料表面BRDF模型作为基础模型;步骤2:根据步骤1所述基础模型建立空间质子辐照环境下材料表面BRDF模型;步骤3:利用质子辐照量及对应的BRDF测量数据对步骤2的模型进行拟合,确定模型中的待定参数集合Ω中的所有待定参数,即完成了空间质子辐照环境下材料表面BRDF的建模。本发明能有效分析不同在轨时间内不同量的质子辐照对BRDF特性的影响规律,支撑空间辐照环境下空间目标反射特性仿真、分析及相关应用。

Description

一种空间质子辐照环境下材料表面BRDF建模方法
技术领域
本发明属于材料表面散射特性领域;具体涉及一种空间质子辐照环境下材料表面BRDF建模方法。
背景技术
空间目标表面的包覆材料一般为聚酰亚胺镀铝薄膜,也称Kapton/Al薄膜二次表面镜,是一种复合表面,常用作空间目标表面包覆的热控材料。
材料表面BRDF是影响空间目标反射特性的主要因素,是空间目标探测、识别的数据基础及源头,是该领域的核心研究内容。近些年来,在BRDF建模方面开展了诸多研究工作,但主要还是在地面光照条件下的BRDF测量与建模。
然后,空间目标处于太空辐照环境下,质子辐照对空间目标表面BRDF特性的影响不可忽略。因此,为了体现质子辐照对材料表面BRDF的影响,需要将质子辐照注量作为BRDF模型中的变量,建立质子辐照对BRDF的影响模型,进而对空间目标光度信号的仿真、分析和应用提供有效支撑。
发明内容
本发明提供了一种空间质子辐照环境下材料表面BRDF建模方法,有效分析不同在轨时间内不同量的质子辐照对BRDF特性的影响规律,支撑空间辐照环境下空间目标反射特性仿真、分析及相关应用。
本发明通过以下技术方案实现:
一种空间质子辐照环境下材料表面BRDF建模方法,所述建模方法包括以下步骤:
步骤1:选择一种经过验证且没有考虑质子辐照的空间目标材料表面BRDF模型作为基础模型;
步骤2:根据步骤1所述基础模型f1irir1)建立空间质子辐照环境下材料表面BRDF模型,
f(θirir,Ω)=f1irir1)f2(IPr2)
其中,θi
Figure BDA0002716252570000011
分别为入射天顶角、入射方位角;θr
Figure BDA0002716252570000012
分别为反射天顶角、反射方位角;Ω1为基础模型的待定参数集合;
Figure BDA0002716252570000013
表示质子辐照对BRDF的影响模型;IP表示质子辐照注量;Ω2为影响模型的待定参数集合,Ω2={A,B,C,D,E,F},Ω为整体模型的待定参数集合;
步骤3:利用质子辐照量及对应的BRDF测量数据对步骤2的模型进行拟合,确定模型中的待定参数集合Ω中的所有待定参数,即完成了空间质子辐照环境下材料表面BRDF的建模。
进一步的,所述步骤1的基础模型具体可为以下模型中的一种,Torrance-Sparrow五参数模型、Phong模型、Cook/Torrance模型、Ward模型或经过验证的基础模型。
进一步的,所述步骤1的基础模型选择为
Figure BDA0002716252570000021
其中,θi
Figure BDA0002716252570000022
分别为入射天顶角、入射方位角;θr
Figure BDA0002716252570000023
分别为反射天顶角、反射方位角;Ω1为基础模型的待定参数集合;D为微面元分布函数,D=exp[c(1-cosγ)d],
Figure BDA0002716252570000024
为入射方位角与反射方位角的夹角,即
Figure BDA0002716252570000025
G为几何衰减因子,
Figure BDA0002716252570000026
F(θi,λ)为菲涅尔反射系数,F(θi)=exp[a(b-cosθi)2];ks为镜面反射系数;kd为漫反射系数;可调参数集合Ω1={a,b,c,d,ks,kd}。
进一步的,所述步骤2中的空间质子辐照环境下材料表面BRDF模型具体为,
Figure BDA0002716252570000027
其中,Ω为整体模型的待定参数集合,Ω={a,b,c,d,ks,kd,A,B,C,D,E,F}。
本发明的有益效果是:
(1)本发明考虑空间质子辐照,建立了含有质子辐照注量的空间目标材料表面的BRDF模型,可描述质子辐照条件下材料表面BRDF的演变规律。
(2)本发明的模型包含两部分,一部分为基础模型,一部分为辐照影响模型,二者融合即可建立质子辐照条件下的材料表面BRDF模型。
(3)本发明的基础模型可以采用任何一种能够表征非质子辐照条件下的材料表面BRDF的典型模型形式,因此融合模型是一种通用BRDF表征模型。
(4)本发明中引入的质子辐照对BRDF的影响模型,采用多待定参数形式,可以有效拟合测量数据。
附图说明
附图1本发明的建模流程。
附图2本发明入射角30度时BRDF测量结果与建模结果的比对示意图,图2-(a)入射角30度时BRDF测量结果与建模结果未辐伏示意图,图2-(b)入射角30度时BRDF测量结果与建模结果辐照注量为5×1014cm-2示意图,图2-(c)入射角30度时BRDF测量结果与建模结果辐照注量5×1015cm-2示意图,图2-(d)入射角30度时BRDF测量结果与建模结果辐照注量为1×1015cm-2示意图,图2-(e)入射角30度时BRDF测量结果与建模结果辐照注量1×1016cm-2示意图。
具体实施方式
下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
一种空间质子辐照环境下材料表面BRDF建模方法,所述建模方法包括以下步骤:
步骤1:选择一种经过验证且没有考虑质子辐照的空间目标材料表面BRDF模型作为基础模型;
步骤2:根据步骤1所述基础模型f1irir1)建立空间质子辐照环境下材料表面BRDF模型,
f(θirir,Ω)=f1irir1)f2(IPr2)
其中,θi
Figure BDA0002716252570000031
分别为入射天顶角、入射方位角;θr
Figure BDA0002716252570000032
分别为反射天顶角、反射方位角;Ω1为基础模型的待定参数集合;
Figure BDA0002716252570000033
表示质子辐照对BRDF的影响模型;IP表示质子辐照注量;Ω2为影响模型的待定参数集合,Ω2={A,B,C,D,E,F},Ω为整体模型的待定参数集合;
步骤3:利用质子辐照量及对应的BRDF测量数据对步骤2的模型进行拟合,确定模型中的待定参数集合Ω中的所有待定参数,即完成了空间质子辐照环境下材料表面BRDF的建模。
进一步的,所述步骤1的基础模型具体可为以下模型中的一种,如Torrance-Sparrow五参数模型、Phong模型、Cook/Torrance模型、Ward模型或经过验证的基础模型。
进一步的,所述步骤1的基础模型选择为
Figure BDA0002716252570000041
其中,θi
Figure BDA0002716252570000042
分别为入射天顶角、入射方位角;θr
Figure BDA0002716252570000043
分别为反射天顶角、反射方位角;Ω1为基础模型的待定参数集合;D为微面元分布函数,D=exp[c(1-cosγ)d],
Figure BDA0002716252570000044
为入射方位角与反射方位角的夹角,即
Figure BDA0002716252570000045
G为几何衰减因子,
Figure BDA0002716252570000046
F(θi,λ)为菲涅尔反射系数,F(θi)=exp[a(b-cosθi)2];ks为镜面反射系数;kd为漫反射系数;可调参数集合Ω1={a,b,c,d,ks,kd}。
进一步的,所述步骤2中的空间质子辐照环境下材料表面BRDF模型具体为,
Figure BDA0002716252570000047
其中,Ω为整体模型的待定参数集合,Ω={a,b,c,d,ks,kd,A,B,C,D,E,F}。
实施例2
选160keV质子进行辐照试验,选择试验注量为5×1014个/cm2、1×1015个/cm2、5×1015个/cm2及1×1016个/cm2。测量得到材料表面的光谱BRDF,针对532nm波段的测量结果进行建模,采用神经网络拟合方法求解待定参数集合中的所有待定参数,求解结果见表1所示,入射角30度时BRDF测量结果与建模结果的比对示意图如图2所示,可见建模方法具有较高的精度。
可见,辐照注量增大时,材料表面BRDF下降;相对无质子辐照条件,辐照注量为1×1016cm-2时,BRDF的峰值下降超过50%,说明了质子辐照对材料表面的BRDF有较严重的影响。同时,说明了本发明建立的质子辐照条件下的BRDF参数模型可较好的表征材料表面BRDF特性。
表1空间质子辐照条件下的材料表面BRDF建模结果
Figure BDA0002716252570000048

Claims (4)

1.一种空间质子辐照环境下材料表面BRDF建模方法,其特征在于,所述建模方法包括以下步骤:
步骤1:选择一种经过验证且没有考虑质子辐照的空间目标材料表面BRDF模型作为基础模型;
步骤2:根据步骤1所述基础模型f1irir1)建立空间质子辐照环境下材料表面BRDF模型,
f(θirir,Ω)=f1irir1)f2(IPr2)
其中,θi
Figure FDA0002716252560000011
分别为入射天顶角、入射方位角;θr
Figure FDA0002716252560000012
分别为反射天顶角、反射方位角;Ω1为基础模型的待定参数集合;
Figure FDA0002716252560000013
表示质子辐照对BRDF的影响模型;IP表示质子辐照注量;Ω2为影响模型的待定参数集合,Ω2={A,B,C,D,E,F},Ω为整体模型的待定参数集合;
步骤3:利用质子辐照量及对应的BRDF测量数据对步骤2的模型进行拟合,确定模型中的待定参数集合Ω中的所有待定参数,即完成了空间质子辐照环境下材料表面BRDF的建模。
2.根据权利要求1所述一种空间质子辐照环境下材料表面BRDF建模方法,其特征在于,所述步骤1的基础模型具体可为以下模型中的一种,Torrance-Sparrow五参数模型、Phong模型、Cook/Torrance模型、Ward模型或经过验证的基础模型。
3.根据权利要求2所述一种空间质子辐照环境下材料表面BRDF建模方法,其特征在于,所述步骤1的基础模型选择为
Figure FDA0002716252560000014
其中,θi
Figure FDA0002716252560000015
分别为入射天顶角、入射方位角;θr
Figure FDA0002716252560000016
分别为反射天顶角、反射方位角;Ω1为基础模型的待定参数集合;D为微面元分布函数,D=exp[c(1-cosγ)d],
Figure FDA0002716252560000017
Figure FDA0002716252560000018
为入射方位角与反射方位角的夹角,即
Figure FDA0002716252560000019
G为几何衰减因子,
Figure FDA00027162525600000110
F(θi,λ)为菲涅尔反射系数,F(θi)=exp[a(b-cosθi)2];ks为镜面反射系数;kd为漫反射系数;可调参数集合Ω1={a,b,c,d,ks,kd}。
4.根据权利要求1所述一种空间质子辐照环境下材料表面BRDF建模方法,其特征在于,所述步骤2中的空间质子辐照环境下材料表面BRDF模型具体为,
Figure FDA0002716252560000021
其中,Ω为整体模型的待定参数集合,Ω={a,b,c,d,ks,kd,A,B,C,D,E,F}。
CN202011074446.0A 2020-10-09 2020-10-09 一种空间质子辐照环境下材料表面brdf建模方法 Active CN112149310B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011074446.0A CN112149310B (zh) 2020-10-09 2020-10-09 一种空间质子辐照环境下材料表面brdf建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011074446.0A CN112149310B (zh) 2020-10-09 2020-10-09 一种空间质子辐照环境下材料表面brdf建模方法

Publications (2)

Publication Number Publication Date
CN112149310A true CN112149310A (zh) 2020-12-29
CN112149310B CN112149310B (zh) 2022-08-02

Family

ID=73952639

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011074446.0A Active CN112149310B (zh) 2020-10-09 2020-10-09 一种空间质子辐照环境下材料表面brdf建模方法

Country Status (1)

Country Link
CN (1) CN112149310B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070018996A1 (en) * 2005-07-25 2007-01-25 Microsoft Corporation Real-time rendering of partially translucent objects
CN105243289A (zh) * 2015-11-17 2016-01-13 上海无线电设备研究所 复杂目标红外辐射特性一体化建模方法
CN109932341A (zh) * 2019-03-11 2019-06-25 北京环境特性研究所 野外环境下典型目标的双向反射分布函数测量方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070018996A1 (en) * 2005-07-25 2007-01-25 Microsoft Corporation Real-time rendering of partially translucent objects
CN105243289A (zh) * 2015-11-17 2016-01-13 上海无线电设备研究所 复杂目标红外辐射特性一体化建模方法
CN109932341A (zh) * 2019-03-11 2019-06-25 北京环境特性研究所 野外环境下典型目标的双向反射分布函数测量方法

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
ABDELAZIZ KALLEL: ""FluLCVRT: Reflectance and fluorescence of leaf and canopy modeling based on Monte Carlo vector radiative transfer simulation"", 《JOURNAL OF QUANTITATIVE SPECTROSCOPY AND RADIATIVE TRANSFER》 *
HOU QING-YU等: ""Inversion and reconstruction of the macroscopic photometric characterization model for on-orbit space object"", 《ACTA PHYSICA SINICA》 *
HOU QINGYU等: ""Measurement of Equivalent BRDF on the Surface of Solar Panel with Periodic Structure"", 《COATINGS》 *
Z.M.ZHANG等: ""Optical and Thermal Radiative Properties of Semiconductors Related to Micro/Nanotechnology"", 《ADVANCES IN HEAT TRANSFER》 *
杨本永等: "光学传感器星上定标漫射板的特性测量", 《光学精密工程》 *
石良臣等: "空间飞行器红外光学信号的模拟", 《宇航学报》 *
谭凡教 等: ""基于时谱信号分析的在轨空间目标姿态感知"", 《物理学报》 *
赵艳华等: "漫反射板全光路全视场全口径在轨辐射定标技术", 《航天返回与遥感》 *
韩意等: "空间目标光学散射特性研究进展", 《红外与激光工程》 *

Also Published As

Publication number Publication date
CN112149310B (zh) 2022-08-02

Similar Documents

Publication Publication Date Title
McNeil et al. A validation of a ray-tracing tool used to generate bi-directional scattering distribution functions for complex fenestration systems
Andersen et al. Bi-directional transmission properties of Venetian blinds: experimental assessment compared to ray-tracing calculations
AU2014346447B2 (en) Texture analysis of a coated surface using kepler's planetary motion laws
Andersen et al. Comparison between ray-tracing simulations and bi-directional transmission measurements on prismatic glazing
CN109932341B (zh) 野外环境下典型目标的双向反射分布函数测量方法
CN108955608A (zh) 绝缘子表面rtv涂层涂覆效果评估方法、装置及系统
CN103646175A (zh) 一种目标光谱辐射亮度的计算方法
Zinzi et al. Directional and angular response of construction materials solar properties: Characterisation and assessment
CN114216559B (zh) 一种星上定标机构的部分孔径因子测量方法和装置
CN112149310B (zh) 一种空间质子辐照环境下材料表面brdf建模方法
Pfreundt et al. Rapid calculation of the backsheet coupling gain using ray groups
Riedel‐Lyngskær et al. Interlaboratory comparison of angular‐dependent photovoltaic device measurements: Results and impact on energy rating
CN104807616B (zh) 一种Spectralon漫反射板校正方法
Boppana et al. Impact of uncertainty in IAM measurement on energy predictions
Balling A comparative study of the bidirectional reflectance distribution function of several surfaces as a mid-wave infrared diffuse reflectance standard
Freeman et al. Using measured plane-of-array data directly in photovoltaic modeling: Methodology and validation
Aldaz et al. The influence of optical characterization at different angles of incidence on optical efficiency calculation of a novel small-size parabolic trough collector for process heat applications
Jonsson et al. Light-scattering properties of a woven shade-screen material used for daylighting and solar heat-gain control
Pal et al. Understanding and Benchmarking Ground Reflectors for Bifacial Photovoltaic Yield Enhancement
CN107247038B (zh) 一种获取河流冰凌红外波段散射特性的方法
CN104008268A (zh) 空间目标可见光散射特性分析模型校验方法
Schregle et al. Physical validation of global illumination methods: measurement and error analysis
Qin et al. Application of small angle approximation in circumsolar irradiance modelling
Riedel et al. Incident Angle Modifier (IAM) Round Robin Updates
CN111207912A (zh) 光学元件散射光束的空间分布检测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant