本论文以Breit-Pauli哈密顿的球张量形式为基础,借助角动量耦合理论和不可约张量理论,导出了铍原子非相对论能量、相对论修正（含质量修正、达尔文修正、自旋-自旋接触修正、轨道-轨道相互作用修正）以及精细结构（含自旋-轨道相互作用、自旋-其它轨道相互作用、自旋-自旋相互作用）能级的实用计算式,其中完成了所有的角向积分和自旋求和计算,将相关能量表示为一系列径向积分之和。在此基础上,利用Mathematica数学软件编写了计算铍原子能级的计算程序,包含Racah表象与Slater表象之间的变换、铍原子的非相对论能量及其相对论修正的计算、铍原子精细结构能级的计算等子程序。具体计算了类铍体系电子组态1s22s2和1s22snp（n=2-6,Z=4-8）的非相对论能量、相对论修正、以及1s22snp3P谱项的精细结构能量,并将计算结果与实验数据进行了比较。探索了提高铍原子能量计算精度的两种方案。其一是重新构造单电子径向波函数的方案,其要点是以类氢型径向波函数为原形,将其中的广义拉盖尔多项式中的每一项前面添加一个变分系数。利用此方案具体计算了铍原子1s22sns（n=3-6）组态的能级结构,计算结果与实验值的相对误差很小。其二是采用WBEPM理论计算铍原子高里德堡态能级的方案,利用此方案具体计算了铍原子1s22sns1S0（n=3-50）、1s22sns1S0（n=3-50）和1s22snd 1D2 （n =3-50）里德堡系列的能级和量子亏损。计算结果与已知的实验数据符合得非常好,同时预言了很多实验上尚未观察到的能级数据。