基于面向对象的程序设计方法，以C++语言构造了自适应有限元P加密程序，它涉及到形函数的构造、高阶单元的形成、自升阶的算法以及误差估计一系列类及成员函数．解线性方程组时，根据P加密的特点，提出了一种节约存储和运算空间的数据结构．所给算例的计算结果表明，理论推导和程序实现是正确的．
采用Biot提出的流固两相介质耦合的波动方程来研究半空间饱和土在垂直集中力作用下的瞬态问题．利用Laplace-Hankel积分变换得到变换域内的解析表达式，并利用Laplace-Ha nkel数值逆变换得到半空间饱和土在时域内的数值解．退化到线弹性中的解与Pekeri s的闭合解进行比较，验证了本文结果的正确性和数值积分的可靠性．
