为了使用它，您分配节点坐标 (NC) 并定义有限元模型中节点之间的连接。

图 9.1 节理岩体模型概念图示 节理岩体模型是一个各向异性的理想弹塑性模型，特别适用于模拟层理和节 理岩石层的行为。在这个模型中，假设原状岩体具有一个可能地层理方向和主要 节理方向。这种岩体被认为具有横观各向同性的弹性材料行为，由五个参数和一 个方向来量化。各向异性可能来自于层理现象或者其他现象。在主要节理方向上， 根据库伦准则，假设剪应力是有限的。当在这个方向达到 大剪应力时，塑性滑 移就会发生。 多可定义三个滑移方向（‘平面’），其中第一个平面假定为和弹 性各向异性方向重合。每一个平面都有不同的抗剪强度性质。除了塑性剪切，根 据预定义的拉力强度（拉力‘截断’），垂直于三个平面的拉应力都是有限的。 当呈现出节理簇和节理集时，节理岩石模型就能被证明是合理的。这些节理 集必须平行，中间不能填有断层泥，相对于结构的特征尺度来说，它们的间距必 须较小。 节理岩体模型的一些基本特征是： * 原状岩体的各向异性弹性行为 参数： E1，E2，v1，v2 ，G2 * 三个方向上根据库伦准则的剪切破坏 参数： c i ，ϕ i， ψi * 三个方向上的有限拉伸强度 参数：σt, j 9.1 各向异性弹性材料的刚度矩阵 节理岩体模型中的弹性材料行为是由弹性材料的刚度矩阵 D* 来描述的。与 Hooke 定律不同，节理岩体模型中的矩阵 *D 是横观各向同性的。在垂直于或者 平行于一个预先定义的方向（‘平面 1’）上的刚度可能是不同的。这个方向可以 对应着层理方向或者任何其他有显著不同的弹性刚度性质的方向。 比如考虑水平层理的情况，水平方向的刚度 E1 就不同于竖直方向的刚度 85

当受拉刚度退化系数取为零时：

基于weber能量法求解齿轮时变啮合刚度的程序，能够跑出刚度图，通过求解轮齿部分变形、基体变形及局部接触变形这三部分的变形，进而求得综合弹性变形，最终求出时变啮合刚度

使用局部刚度矩阵和连接矩阵组装全局刚度矩阵。

用apdl语言提取ansys有限元模型的整体、单元刚度和质量矩阵，感兴趣的可以学一下

使用石川公式计算齿轮时变啮合刚度

实验三 考虑结构刚度时的直流电动机-负载建模及仿真试验，包括word成品报告，matlab编辑的.m文件，Simulink搭建的模型.slx文件，MathType编辑的公式文件。 实验要求： 1） 根据给定的动力学方程及方块图，求出从电机电枢电压到转角的传递函数模型，并求其频率特性； 2）求出从 到 的传递函数模型，并求其频率特性和根轨迹； 3）分别取 和 ，编制MATLAB或Simulink程序，比较刚度系数不同时，电机-负载模型的频率特性（从 到转角 ）。 ......

金属材料刚度,强度及钢板知识.doc

此函数计算任意维度的常规 8 节点六面体有限元的 24x24 刚度矩阵。 数值积分是用 8 个高斯点进行的。 因此权重因子是一，它们不包含在代码中。