自然单元法（Natural Element Method,简称NEM）是一种近20年才发展起来的数值算法。它采用基于Voronoi图的Sibson或non-Sibonian插值方法构造全域近似函数和试函数。自然单元法具有无网格不依赖网格的优势，同时形函数构造简单，在边界上满足Delta插值性质，是一种发展前景广阔的求解力学和其他工程技术中偏微分方程的数值方法针对自然单元法无法直接求解节点应力且应力解精度不高的缺陷以及采用高阶连续形函数恢复应力计算效率不高的问题，通过引入应力杂交的思想，采用对应力和位移分别插值，提出了杂交自然单元法。利用杂交自然单元法对弹性力学问题、弹塑性力学问题、粘弹性力学问题以及弹性大变形问题进行了分析。具体工作如下：将复变量移动最小二乘法引入到自然单元法应力恢复算法中，提出了改进的自然单元法应力恢复算法。与传统应力恢复算法相比，复变量移动最小二乘法待定系数少，可以提高应力恢复算法的计算效率。采用自然邻点插值对应力和位移分别进行逼近，结合弹性力学的多变量Hellinger-Reissner变分原理，提出了弹性力学的杂交自然单元法。与传统自然单元法相比，杂交自然单元法有着更高的应力计算精度并能够直接求解节点应力。同时，在相同精度条件下，杂交自然单元法比应力恢复算法有着更高计算效率。在杂交自然单元法的基础上，通过对应力和位移增量分别进行插值，提出了适用于非线性问题的增量形式杂交自然单元法。结合弹塑性问题的Hellinger-Reissner变分原理，提出了弹塑性力学问题的杂交自然单元法，并推导了计算公式。与自然单元法相比，杂交自然单元法有着更高的计算精度，比应力恢复算法更加节省计算时间。在详细分析了线性粘弹性材料的蠕变特性的基础上，采用各个时间节点上的变形来描述蠕变过程。通过杂交自然单元法与在各个时间节点上建立的Hellinger-Reissner变分原理，推导了适用于线性粘弹性材料蠕变分析的杂交自然单元法。本文方法能够准确的描述材料的蠕变性能；并能够直接求解得到准确的节点应力。结合大变形问题的特性，推导了增量型杂交自然单元法分析弹性大变形问题的Total Lagrangian格式下表达方式，数值算例中采用了Newton-Raphson增量迭代法，提出了适用于弹性大变形问题分析的杂交自然单元法。算例表明杂交自然单元法有着更高的计算精度和更快的求解效率。利用增量形式的Hellinger-Reissner变分原理结合岩土工程的特点，推导了适用于岩土工程开挖问题的模拟的杂交自然单元法；同时为了模拟隧道的支护过程提出了适用于本文方法的添加节点的方法。为了证明本文提出的杂交自然单元法的有效性，采用Matlab编制了数值算例。通过数值算例分析可以看出本文方法的有效性和优越性。