单自由度体系ansys与matlab阻尼对比

matlab中阻尼的计算：

单自由度结构体系采用钢筋混凝土材料，假定阻尼比ξ=0.05；

ξ=c2mω

推出c=2mωξ;并用该式子计算结构的整体的常数阻尼。经过验证，常数阻尼和变阻尼对结果的影响不大。

ansys中阻尼的计算：

依据上表得出瞬态分析（full）法中只能考虑α，β阻尼，材料阻尼和单元阻尼；其他定义的阻尼在瞬态分析中均被忽略，瞬态分析（full）法中阻尼矩阵的计算公式如下（若为单自由度体系，阻尼C为粘滞系数）：

式中α 为α阻尼，用ALPHAD命令定义。 β为β阻尼，用命令BETAD定义。

[M]为质量矩阵，[K]为刚度矩阵；

Nm为具有材料阻尼的材料数目，用命令MP,DAMP定义；

βjm为材料的刚度阻尼系数，用命令MP,DAMP定义

[Kj]与材料j相关的部分的结构刚度矩阵；

Nc为具有单元阻尼的单元数目，

CK为单元阻尼矩阵；

公式的前半部分为Rayleigh阻尼，可以用结构整体的阻尼。

Rayleigh阻尼是最常用的粘滞阻尼模型，也称比例阻尼；即：

式子中：

α：质量矩阵系数，用ALPHAD命令定义；

β：刚度矩阵系数，用BETAD命令定义；

[M]：质量矩阵；

[K]：刚度矩阵；

通常β与α未知，一般采用模态阻尼比来计算，用该方法来计算等效的结构整体阻尼比；一般假定各阶模态的阻尼比相同；

式子中：

ωiωj分别为第i阶，第j阶的自振频率；

为使得ansys软件的计算条件和matlab软件计算条件一致；单自由度体系在ansys软件的阻尼定义有两种类型：

α=0时，β=2ξωi；或者β=0时；α=2ωξ.。

Rayleigh阻尼即阻尼矩阵c的前半部分，为全局性的整个结构阻尼，适合所有的动力分析。当模型由许多不同材料部件或构件组成时,各部分的材料阻尼相差甚远,需要考虑阻尼的局部差异；在此种情况下，采用材料阻尼考虑结构的阻尼的一种附加手段。

验证结果下次再附上