收藏 分享(赏)

外延薄膜生长过程的微观相场研究.doc

上传人:cjc2202537 文档编号:1529090 上传时间:2018-07-25 格式:DOC 页数:38 大小:71.73KB
下载 相关 举报
外延薄膜生长过程的微观相场研究.doc_第1页
第1页 / 共38页
外延薄膜生长过程的微观相场研究.doc_第2页
第2页 / 共38页
外延薄膜生长过程的微观相场研究.doc_第3页
第3页 / 共38页
外延薄膜生长过程的微观相场研究.doc_第4页
第4页 / 共38页
外延薄膜生长过程的微观相场研究.doc_第5页
第5页 / 共38页
点击查看更多>>
资源描述

1、固体力学专业毕业论文 精品论文 外延薄膜生长过程的微观相场研究关键词:相场模型 数值积分算法 半隐谱法 外延生长 台阶运动 外延薄膜 微观相场 材料生长摘要:材料生长过程的计算机模拟和仿真是当今国际公认的制造科学与材料科学的前沿领域之一,而相场法是一种用于描述非平衡态下复杂相界面演变的强有力工具。本文针对外延生长过程的介观尺度问题,基于与时间相关的Ginzburg-Landau 方程得到的纯物质相场模型,在 Neumann 边界条件下通过有限差分离散相场和温度场控制方程,构造该控制方程的显式数值积分方案,给出数值计算格式,编制相应模拟程序,通过适当的后处理再现了枝晶外延生长过程。 为克服显式积

2、分方案的不足,在考虑周期边界条件的情况下,利用半隐谱法结合快速傅里叶变换算法来求解相场和温度场的耦合控制方程,给出稳定收敛的有效积分方案。数值算例表明半隐谱法能在保证计算精度的条件下,降低对时间步长的依赖,同时扩大了相场法模拟的时间尺度,较有限差分法能有效提高计算效率。 通过系统的数值计算详细研究了材料凝固过程枝晶生长,再现了枝晶生长过程侧向分支的竞争生长机制,着重研究了初始晶核的尺寸、形状、Stefen 数等关键因素以及时间和空间步长选取等对外延生长行为的影响。结果表明,随着过冷度的增加,枝晶生长加快,同时界面变得不稳定;界面厚度取值在保证合理的情况下应尽可能小;参数的计算应与材料的热物性取

3、值相关等。 在对枝晶生长的相场模型详细研究的基础上,将相场模型应用到台阶式外延生长过程的数值研究中去,采用一个与时间相关的 Ginzburg-Landau 方程与活原子浓度扩散方程相耦合,从而对岛的形核和生长进行描述。利用该模型主要针对螺形和梯形这两种台阶模拟含有三层台阶的外延薄膜生长过程。数值结果表明,本文提出的方法完全可以再现台阶生长的演变趋势,并得出与试验结果一致的变化趋势。 最后将弹性应变能引入到相场方程中,并结合半隐谱法进行相关方面的理论推导工作,初步探讨了材料生长过程中弹性应力的作用。正文内容材料生长过程的计算机模拟和仿真是当今国际公认的制造科学与材料科学的前沿领域之一,而相场法是

4、一种用于描述非平衡态下复杂相界面演变的强有力工具。本文针对外延生长过程的介观尺度问题,基于与时间相关的Ginzburg-Landau 方程得到的纯物质相场模型,在 Neumann 边界条件下通过有限差分离散相场和温度场控制方程,构造该控制方程的显式数值积分方案,给出数值计算格式,编制相应模拟程序,通过适当的后处理再现了枝晶外延生长过程。 为克服显式积分方案的不足,在考虑周期边界条件的情况下,利用半隐谱法结合快速傅里叶变换算法来求解相场和温度场的耦合控制方程,给出稳定收敛的有效积分方案。数值算例表明半隐谱法能在保证计算精度的条件下,降低对时间步长的依赖,同时扩大了相场法模拟的时间尺度,较有限差分

5、法能有效提高计算效率。 通过系统的数值计算详细研究了材料凝固过程枝晶生长,再现了枝晶生长过程侧向分支的竞争生长机制,着重研究了初始晶核的尺寸、形状、Stefen 数等关键因素以及时间和空间步长选取等对外延生长行为的影响。结果表明,随着过冷度的增加,枝晶生长加快,同时界面变得不稳定;界面厚度取值在保证合理的情况下应尽可能小;参数的计算应与材料的热物性取值相关等。 在对枝晶生长的相场模型详细研究的基础上,将相场模型应用到台阶式外延生长过程的数值研究中去,采用一个与时间相关的 Ginzburg-Landau 方程与活原子浓度扩散方程相耦合,从而对岛的形核和生长进行描述。利用该模型主要针对螺形和梯形这

6、两种台阶模拟含有三层台阶的外延薄膜生长过程。数值结果表明,本文提出的方法完全可以再现台阶生长的演变趋势,并得出与试验结果一致的变化趋势。 最后将弹性应变能引入到相场方程中,并结合半隐谱法进行相关方面的理论推导工作,初步探讨了材料生长过程中弹性应力的作用。材料生长过程的计算机模拟和仿真是当今国际公认的制造科学与材料科学的前沿领域之一,而相场法是一种用于描述非平衡态下复杂相界面演变的强有力工具。本文针对外延生长过程的介观尺度问题,基于与时间相关的 Ginzburg-Landau 方程得到的纯物质相场模型,在 Neumann 边界条件下通过有限差分离散相场和温度场控制方程,构造该控制方程的显式数值积

7、分方案,给出数值计算格式,编制相应模拟程序,通过适当的后处理再现了枝晶外延生长过程。 为克服显式积分方案的不足,在考虑周期边界条件的情况下,利用半隐谱法结合快速傅里叶变换算法来求解相场和温度场的耦合控制方程,给出稳定收敛的有效积分方案。数值算例表明半隐谱法能在保证计算精度的条件下,降低对时间步长的依赖,同时扩大了相场法模拟的时间尺度,较有限差分法能有效提高计算效率。 通过系统的数值计算详细研究了材料凝固过程枝晶生长,再现了枝晶生长过程侧向分支的竞争生长机制,着重研究了初始晶核的尺寸、形状、Stefen 数等关键因素以及时间和空间步长选取等对外延生长行为的影响。结果表明,随着过冷度的增加,枝晶生

8、长加快,同时界面变得不稳定;界面厚度取值在保证合理的情况下应尽可能小;参数的计算应与材料的热物性取值相关等。在对枝晶生长的相场模型详细研究的基础上,将相场模型应用到台阶式外延生长过程的数值研究中去,采用一个与时间相关的 Ginzburg-Landau 方程与活原子浓度扩散方程相耦合,从而对岛的形核和生长进行描述。利用该模型主要针对螺形和梯形这两种台阶模拟含有三层台阶的外延薄膜生长过程。数值结果表明,本文提出的方法完全可以再现台阶生长的演变趋势,并得出与试验结果一致的变化趋势。 最后将弹性应变能引入到相场方程中,并结合半隐谱法进行相关方面的理论推导工作,初步探讨了材料生长过程中弹性应力的作用。材

9、料生长过程的计算机模拟和仿真是当今国际公认的制造科学与材料科学的前沿领域之一,而相场法是一种用于描述非平衡态下复杂相界面演变的强有力工具。本文针对外延生长过程的介观尺度问题,基于与时间相关的 Ginzburg-Landau 方程得到的纯物质相场模型,在 Neumann 边界条件下通过有限差分离散相场和温度场控制方程,构造该控制方程的显式数值积分方案,给出数值计算格式,编制相应模拟程序,通过适当的后处理再现了枝晶外延生长过程。 为克服显式积分方案的不足,在考虑周期边界条件的情况下,利用半隐谱法结合快速傅里叶变换算法来求解相场和温度场的耦合控制方程,给出稳定收敛的有效积分方案。数值算例表明半隐谱法

10、能在保证计算精度的条件下,降低对时间步长的依赖,同时扩大了相场法模拟的时间尺度,较有限差分法能有效提高计算效率。 通过系统的数值计算详细研究了材料凝固过程枝晶生长,再现了枝晶生长过程侧向分支的竞争生长机制,着重研究了初始晶核的尺寸、形状、Stefen 数等关键因素以及时间和空间步长选取等对外延生长行为的影响。结果表明,随着过冷度的增加,枝晶生长加快,同时界面变得不稳定;界面厚度取值在保证合理的情况下应尽可能小;参数的计算应与材料的热物性取值相关等。在对枝晶生长的相场模型详细研究的基础上,将相场模型应用到台阶式外延生长过程的数值研究中去,采用一个与时间相关的 Ginzburg-Landau 方程

11、与活原子浓度扩散方程相耦合,从而对岛的形核和生长进行描述。利用该模型主要针对螺形和梯形这两种台阶模拟含有三层台阶的外延薄膜生长过程。数值结果表明,本文提出的方法完全可以再现台阶生长的演变趋势,并得出与试验结果一致的变化趋势。 最后将弹性应变能引入到相场方程中,并结合半隐谱法进行相关方面的理论推导工作,初步探讨了材料生长过程中弹性应力的作用。材料生长过程的计算机模拟和仿真是当今国际公认的制造科学与材料科学的前沿领域之一,而相场法是一种用于描述非平衡态下复杂相界面演变的强有力工具。本文针对外延生长过程的介观尺度问题,基于与时间相关的 Ginzburg-Landau 方程得到的纯物质相场模型,在 N

12、eumann 边界条件下通过有限差分离散相场和温度场控制方程,构造该控制方程的显式数值积分方案,给出数值计算格式,编制相应模拟程序,通过适当的后处理再现了枝晶外延生长过程。 为克服显式积分方案的不足,在考虑周期边界条件的情况下,利用半隐谱法结合快速傅里叶变换算法来求解相场和温度场的耦合控制方程,给出稳定收敛的有效积分方案。数值算例表明半隐谱法能在保证计算精度的条件下,降低对时间步长的依赖,同时扩大了相场法模拟的时间尺度,较有限差分法能有效提高计算效率。 通过系统的数值计算详细研究了材料凝固过程枝晶生长,再现了枝晶生长过程侧向分支的竞争生长机制,着重研究了初始晶核的尺寸、形状、Stefen 数等

13、关键因素以及时间和空间步长选取等对外延生长行为的影响。结果表明,随着过冷度的增加,枝晶生长加快,同时界面变得不稳定;界面厚度取值在保证合理的情况下应尽可能小;参数的计算应与材料的热物性取值相关等。在对枝晶生长的相场模型详细研究的基础上,将相场模型应用到台阶式外延生长过程的数值研究中去,采用一个与时间相关的 Ginzburg-Landau 方程与活原子浓度扩散方程相耦合,从而对岛的形核和生长进行描述。利用该模型主要针对螺形和梯形这两种台阶模拟含有三层台阶的外延薄膜生长过程。数值结果表明,本文提出的方法完全可以再现台阶生长的演变趋势,并得出与试验结果一致的变化趋势。 最后将弹性应变能引入到相场方程

14、中,并结合半隐谱法进行相关方面的理论推导工作,初步探讨了材料生长过程中弹性应力的作用。材料生长过程的计算机模拟和仿真是当今国际公认的制造科学与材料科学的前沿领域之一,而相场法是一种用于描述非平衡态下复杂相界面演变的强有力工具。本文针对外延生长过程的介观尺度问题,基于与时间相关的 Ginzburg-Landau 方程得到的纯物质相场模型,在 Neumann 边界条件下通过有限差分离散相场和温度场控制方程,构造该控制方程的显式数值积分方案,给出数值计算格式,编制相应模拟程序,通过适当的后处理再现了枝晶外延生长过程。 为克服显式积分方案的不足,在考虑周期边界条件的情况下,利用半隐谱法结合快速傅里叶变

15、换算法来求解相场和温度场的耦合控制方程,给出稳定收敛的有效积分方案。数值算例表明半隐谱法能在保证计算精度的条件下,降低对时间步长的依赖,同时扩大了相场法模拟的时间尺度,较有限差分法能有效提高计算效率。 通过系统的数值计算详细研究了材料凝固过程枝晶生长,再现了枝晶生长过程侧向分支的竞争生长机制,着重研究了初始晶核的尺寸、形状、Stefen 数等关键因素以及时间和空间步长选取等对外延生长行为的影响。结果表明,随着过冷度的增加,枝晶生长加快,同时界面变得不稳定;界面厚度取值在保证合理的情况下应尽可能小;参数的计算应与材料的热物性取值相关等。在对枝晶生长的相场模型详细研究的基础上,将相场模型应用到台阶

16、式外延生长过程的数值研究中去,采用一个与时间相关的 Ginzburg-Landau 方程与活原子浓度扩散方程相耦合,从而对岛的形核和生长进行描述。利用该模型主要针对螺形和梯形这两种台阶模拟含有三层台阶的外延薄膜生长过程。数值结果表明,本文提出的方法完全可以再现台阶生长的演变趋势,并得出与试验结果一致的变化趋势。 最后将弹性应变能引入到相场方程中,并结合半隐谱法进行相关方面的理论推导工作,初步探讨了材料生长过程中弹性应力的作用。材料生长过程的计算机模拟和仿真是当今国际公认的制造科学与材料科学的前沿领域之一,而相场法是一种用于描述非平衡态下复杂相界面演变的强有力工具。本文针对外延生长过程的介观尺度

17、问题,基于与时间相关的 Ginzburg-Landau 方程得到的纯物质相场模型,在 Neumann 边界条件下通过有限差分离散相场和温度场控制方程,构造该控制方程的显式数值积分方案,给出数值计算格式,编制相应模拟程序,通过适当的后处理再现了枝晶外延生长过程。 为克服显式积分方案的不足,在考虑周期边界条件的情况下,利用半隐谱法结合快速傅里叶变换算法来求解相场和温度场的耦合控制方程,给出稳定收敛的有效积分方案。数值算例表明半隐谱法能在保证计算精度的条件下,降低对时间步长的依赖,同时扩大了相场法模拟的时间尺度,较有限差分法能有效提高计算效率。 通过系统的数值计算详细研究了材料凝固过程枝晶生长,再现

18、了枝晶生长过程侧向分支的竞争生长机制,着重研究了初始晶核的尺寸、形状、Stefen 数等关键因素以及时间和空间步长选取等对外延生长行为的影响。结果表明,随着过冷度的增加,枝晶生长加快,同时界面变得不稳定;界面厚度取值在保证合理的情况下应尽可能小;参数的计算应与材料的热物性取值相关等。在对枝晶生长的相场模型详细研究的基础上,将相场模型应用到台阶式外延生长过程的数值研究中去,采用一个与时间相关的 Ginzburg-Landau 方程与活原子浓度扩散方程相耦合,从而对岛的形核和生长进行描述。利用该模型主要针对螺形和梯形这两种台阶模拟含有三层台阶的外延薄膜生长过程。数值结果表明,本文提出的方法完全可以

19、再现台阶生长的演变趋势,并得出与试验结果一致的变化趋势。 最后将弹性应变能引入到相场方程中,并结合半隐谱法进行相关方面的理论推导工作,初步探讨了材料生长过程中弹性应力的作用。材料生长过程的计算机模拟和仿真是当今国际公认的制造科学与材料科学的前沿领域之一,而相场法是一种用于描述非平衡态下复杂相界面演变的强有力工具。本文针对外延生长过程的介观尺度问题,基于与时间相关的 Ginzburg-Landau 方程得到的纯物质相场模型,在 Neumann 边界条件下通过有限差分离散相场和温度场控制方程,构造该控制方程的显式数值积分方案,给出数值计算格式,编制相应模拟程序,通过适当的后处理再现了枝晶外延生长过

20、程。 为克服显式积分方案的不足,在考虑周期边界条件的情况下,利用半隐谱法结合快速傅里叶变换算法来求解相场和温度场的耦合控制方程,给出稳定收敛的有效积分方案。数值算例表明半隐谱法能在保证计算精度的条件下,降低对时间步长的依赖,同时扩大了相场法模拟的时间尺度,较有限差分法能有效提高计算效率。 通过系统的数值计算详细研究了材料凝固过程枝晶生长,再现了枝晶生长过程侧向分支的竞争生长机制,着重研究了初始晶核的尺寸、形状、Stefen 数等关键因素以及时间和空间步长选取等对外延生长行为的影响。结果表明,随着过冷度的增加,枝晶生长加快,同时界面变得不稳定;界面厚度取值在保证合理的情况下应尽可能小;参数的计算

21、应与材料的热物性取值相关等。在对枝晶生长的相场模型详细研究的基础上,将相场模型应用到台阶式外延生长过程的数值研究中去,采用一个与时间相关的 Ginzburg-Landau 方程与活原子浓度扩散方程相耦合,从而对岛的形核和生长进行描述。利用该模型主要针对螺形和梯形这两种台阶模拟含有三层台阶的外延薄膜生长过程。数值结果表明,本文提出的方法完全可以再现台阶生长的演变趋势,并得出与试验结果一致的变化趋势。 最后将弹性应变能引入到相场方程中,并结合半隐谱法进行相关方面的理论推导工作,初步探讨了材料生长过程中弹性应力的作用。材料生长过程的计算机模拟和仿真是当今国际公认的制造科学与材料科学的前沿领域之一,而

22、相场法是一种用于描述非平衡态下复杂相界面演变的强有力工具。本文针对外延生长过程的介观尺度问题,基于与时间相关的 Ginzburg-Landau 方程得到的纯物质相场模型,在 Neumann 边界条件下通过有限差分离散相场和温度场控制方程,构造该控制方程的显式数值积分方案,给出数值计算格式,编制相应模拟程序,通过适当的后处理再现了枝晶外延生长过程。 为克服显式积分方案的不足,在考虑周期边界条件的情况下,利用半隐谱法结合快速傅里叶变换算法来求解相场和温度场的耦合控制方程,给出稳定收敛的有效积分方案。数值算例表明半隐谱法能在保证计算精度的条件下,降低对时间步长的依赖,同时扩大了相场法模拟的时间尺度,

23、较有限差分法能有效提高计算效率。 通过系统的数值计算详细研究了材料凝固过程枝晶生长,再现了枝晶生长过程侧向分支的竞争生长机制,着重研究了初始晶核的尺寸、形状、Stefen 数等关键因素以及时间和空间步长选取等对外延生长行为的影响。结果表明,随着过冷度的增加,枝晶生长加快,同时界面变得不稳定;界面厚度取值在保证合理的情况下应尽可能小;参数的计算应与材料的热物性取值相关等。在对枝晶生长的相场模型详细研究的基础上,将相场模型应用到台阶式外延生长过程的数值研究中去,采用一个与时间相关的 Ginzburg-Landau 方程与活原子浓度扩散方程相耦合,从而对岛的形核和生长进行描述。利用该模型主要针对螺形

24、和梯形这两种台阶模拟含有三层台阶的外延薄膜生长过程。数值结果表明,本文提出的方法完全可以再现台阶生长的演变趋势,并得出与试验结果一致的变化趋势。 最后将弹性应变能引入到相场方程中,并结合半隐谱法进行相关方面的理论推导工作,初步探讨了材料生长过程中弹性应力的作用。材料生长过程的计算机模拟和仿真是当今国际公认的制造科学与材料科学的前沿领域之一,而相场法是一种用于描述非平衡态下复杂相界面演变的强有力工具。本文针对外延生长过程的介观尺度问题,基于与时间相关的 Ginzburg-Landau 方程得到的纯物质相场模型,在 Neumann 边界条件下通过有限差分离散相场和温度场控制方程,构造该控制方程的显

25、式数值积分方案,给出数值计算格式,编制相应模拟程序,通过适当的后处理再现了枝晶外延生长过程。 为克服显式积分方案的不足,在考虑周期边界条件的情况下,利用半隐谱法结合快速傅里叶变换算法来求解相场和温度场的耦合控制方程,给出稳定收敛的有效积分方案。数值算例表明半隐谱法能在保证计算精度的条件下,降低对时间步长的依赖,同时扩大了相场法模拟的时间尺度,较有限差分法能有效提高计算效率。 通过系统的数值计算详细研究了材料凝固过程枝晶生长,再现了枝晶生长过程侧向分支的竞争生长机制,着重研究了初始晶核的尺寸、形状、Stefen 数等关键因素以及时间和空间步长选取等对外延生长行为的影响。结果表明,随着过冷度的增加

26、,枝晶生长加快,同时界面变得不稳定;界面厚度取值在保证合理的情况下应尽可能小;参数的计算应与材料的热物性取值相关等。在对枝晶生长的相场模型详细研究的基础上,将相场模型应用到台阶式外延生长过程的数值研究中去,采用一个与时间相关的 Ginzburg-Landau 方程与活原子浓度扩散方程相耦合,从而对岛的形核和生长进行描述。利用该模型主要针对螺形和梯形这两种台阶模拟含有三层台阶的外延薄膜生长过程。数值结果表明,本文提出的方法完全可以再现台阶生长的演变趋势,并得出与试验结果一致的变化趋势。 最后将弹性应变能引入到相场方程中,并结合半隐谱法进行相关方面的理论推导工作,初步探讨了材料生长过程中弹性应力的

27、作用。材料生长过程的计算机模拟和仿真是当今国际公认的制造科学与材料科学的前沿领域之一,而相场法是一种用于描述非平衡态下复杂相界面演变的强有力工具。本文针对外延生长过程的介观尺度问题,基于与时间相关的 Ginzburg-Landau 方程得到的纯物质相场模型,在 Neumann 边界条件下通过有限差分离散相场和温度场控制方程,构造该控制方程的显式数值积分方案,给出数值计算格式,编制相应模拟程序,通过适当的后处理再现了枝晶外延生长过程。 为克服显式积分方案的不足,在考虑周期边界条件的情况下,利用半隐谱法结合快速傅里叶变换算法来求解相场和温度场的耦合控制方程,给出稳定收敛的有效积分方案。数值算例表明

28、半隐谱法能在保证计算精度的条件下,降低对时间步长的依赖,同时扩大了相场法模拟的时间尺度,较有限差分法能有效提高计算效率。 通过系统的数值计算详细研究了材料凝固过程枝晶生长,再现了枝晶生长过程侧向分支的竞争生长机制,着重研究了初始晶核的尺寸、形状、Stefen 数等关键因素以及时间和空间步长选取等对外延生长行为的影响。结果表明,随着过冷度的增加,枝晶生长加快,同时界面变得不稳定;界面厚度取值在保证合理的情况下应尽可能小;参数的计算应与材料的热物性取值相关等。在对枝晶生长的相场模型详细研究的基础上,将相场模型应用到台阶式外延生长过程的数值研究中去,采用一个与时间相关的 Ginzburg-Landa

29、u 方程与活原子浓度扩散方程相耦合,从而对岛的形核和生长进行描述。利用该模型主要针对螺形和梯形这两种台阶模拟含有三层台阶的外延薄膜生长过程。数值结果表明,本文提出的方法完全可以再现台阶生长的演变趋势,并得出与试验结果一致的变化趋势。 最后将弹性应变能引入到相场方程中,并结合半隐谱法进行相关方面的理论推导工作,初步探讨了材料生长过程中弹性应力的作用。特别提醒 :正文内容由 PDF 文件转码生成,如您电脑未有相应转换码,则无法显示正文内容,请您下载相应软件,下载地址为 http:/ 。如还不能显示,可以联系我 q q 1627550258 ,提供原格式文档。“垐垯櫃 换烫梯葺铑?endstream

30、endobj2x 滌?U 閩 AZ箾 FTP 鈦X 飼?狛P? 燚?琯嫼 b?袍*甒?颙嫯?4)=r 宵?i?j 彺帖 B3 锝檡骹笪 yLrQ#?0 鯖 l 壛枒l 壛枒 l 壛枒 l 壛枒 l 壛枒 l 壛枒 l 壛枒 l 壛枒 l 壛枒 l 壛枒 l 壛枒 l 壛渓?擗#?“?# 綫 G 刿#K 芿$?7. 耟?Wa 癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$Fb癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$Fb 皗 E|?pDb 癳$Fb 癳$Fb癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$Fb 癳$F?責鯻 0 橔 C,f 薍秾腵薍秾腵薍秾腵薍秾腵薍秾腵薍秾腵薍秾腵薍秾腵薍秾腵薍秾腵薍秾腵薍秾腵薍秾腵薍秾腵薍秾腵秾腵薍秾腵%?秾腵薍秾腵薍秾腵薍秾腵薍

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 学术论文 > 毕业论文

本站链接:文库   一言   我酷   合作


客服QQ:2549714901微博号:道客多多官方知乎号:道客多多

经营许可证编号: 粤ICP备2021046453号世界地图

道客多多©版权所有2020-2025营业执照举报