泡沫铝的工程应用日益广泛,其表面相互作用问题的精细计算需求也日渐突出。泡沫铝表面相互作用问题在宏观上可视为经典力学问题,但在微观上则是涉及到原子间相互作用的物理问题。为此,本文研究泡沫铝表面相互作用问题的多尺度模拟方法。对于泡沫铝的宏、细观组织结构,首先建立与材料制备工艺相应的数值模型,作为多尺度模拟的前提。通过具体设计与发泡法和熔模铸造法相应的蒙卡方法,模拟得到了相应的数值试样。通过计算模拟和统计分析大量数值试样的力学性能,得到了不同制备工艺制得的泡沫铝的弹性模量的分散性特点。借鉴代表性体积单元方法(RVE)和完全多格方法(FMG)的思想,采用基于“巨原子”模型的多尺度关联方法,并以简化的物理势场描述巨原子间的关系,设计了原子—巨原子—有限元多尺度算法。以文中已建立的泡沫铝宏、细观组织结构模拟方法生成数值试样。采用IFEPG生成弹性的FEM计算程序,以GID进行相关计算的前、后处理。将FEM计算程序与原子—巨原子模拟程序整合,完成了泡沫铝与理想刚性平面纯弹性接触问题的一个多尺度数值模拟算例。通过算例探讨了与理想刚性平面接触时,泡沫铝接触表面的应力集中情况。计算结果显示泡沫铝在外部接触摩擦下可能出现表面局部应力集中,但是这一情况的出现与材料在原子和巨原子尺度的变形没有显著联系,而主要与泡沫铝内部所具有的宏观可见的孔洞有关。采用IFEPG生成弹塑性的FEM计算程序,对同一个接触算例进行原子—巨原子—有限元多尺度计算,这一模拟在宏观尺度上考虑了材料的弹塑性变形。其结果明显表现出这样的趋势:因塑性变形区的出现,接触面上的切向应力最大值有所减小,切向应力随空间坐标的变化也变得更为平缓。但由于引入了宏观层次的塑性计算,使得计算效率变得很低。在这一点上,泡沫铝弹塑性接触的原子—巨原子—有限元多尺度模拟方法还需大力改进。