初始应力刚度矩阵
初始应力刚度矩阵的理论基础
初始应力刚度矩阵
老师，初始应力刚度矩阵和几何刚度矩阵是同一个东西吗？
是同一概念的不同名称。$[K_\sigma]$ 也被称为“初始应力刚度”、“几何刚度”或“应力刚度”。用于屈曲分析和振动分析。
构成
$[K_\sigma]$ 由当前应力状态 $\{\sigma\}$ 构成。对于梁单元是轴力 $N$，对于壳/实体单元则是膜应力 $N_{xx}, N_{yy}, N_{xy}$。
$[G]$ 是位移梯度矩阵，$[S]$ 是应力矩阵。
总结
欧拉屈曲与初始应力刚度矩阵的发现
将压缩载荷下的屈曲作为线性代数问题来公式化的是Leonhard Euler在1744年的工作。他所推导的“临界载荷”概念实际上对应着现代表述中的“初始应力刚度矩阵失去正定性（特征值变为零）”。1960年代，Turner等人将初始应力刚度矩阵引入FEM，使得任意形状的屈曲分析首次成为可能。
数值解法与实现
$[K_\sigma]$的计算
FEM中的计算：
1. 通过静力分析（或非线性分析）求得应力分布
2. 根据各单元应力计算$[K_\sigma]_e$
3. 组装成整体$[K_\sigma]$
所有求解器均可自动计算（屈曲分析、预应力模态）。NLGEOM非线性分析内部也会使用。
总结
初始应力刚度矩阵的组装方法
初始应力刚度矩阵Kg通过单元内应力σij与形函数导数矩阵的积分∫[G]ᵀ[σ][G]dV计算。当压缩应力占主导时，Kg具有负定值成分，(K+Kg)的行列式为0时的载荷即为线性屈曲载荷Pcr。FEM的线性屈曲分析归结为求解(K+λKg)φ=0的广义特征值问题，最小特征值λmin表示安全系数（λmin>1为安全）。
初始应力刚度矩阵初始应力刚度矩阵实践指南
实务检查清单
钢框架结构的几何非线性设计
AISC 360-10推荐在高层钢结构建筑设计中采用考虑初始应力的“直接分析法（DAM）”。DAM能通过初始应力刚度矩阵自动考虑柱承受压缩载荷状态下的水平刚度降低（P-δ效应），相比对虚拟长细比0.2进行换算的“有效长度法”，虽然保守约3~5%但精度更高。在超过400米的超高层建筑中，DAM已成为标准分析方法。
初始应力刚度矩阵软件与求解器比较
工具
所有求解器都支持。概念共通但名称不同。
SAP2000直接分析法的设置
CSI SAP2000的Direct Analysis Method（DAM）设置中，可以自动应用初始应力刚度矩阵和柱刚度折减系数（0.8×EI），执行考虑P-δ·P-Δ效应的非线性分析。作为AISC 360认证软件，在美国结构设计事务所的采用率很高，特别是在超高层建筑的风·地震响应分析中，DAM的标准设置已经确立。
尖端技术
尖端研究
初始应力刚度：预应力混凝土的设计
初始应力刚度的概念由1940年代的预应力混凝土技术工程师弗莱西奈在桥梁设计中实用化。PC钢绞线的张拉预应力（通常1,200〜1,400 MPa）赋予混凝土初始压缩应力刚度，使承载能力倍增的原理，在现代CAE分析中也作为`PRESTRESS`选项被各求解器继承。
初始应力刚度矩阵常见问题与调试
故障
线性屈曲特征值小于1的不合理结果
线性屈曲分析中即使出现特征值λ<1（在设计载荷下已屈曲）的结果，实际上也可能不发生屈曲。因为线性屈曲假设为无初始缺陷的“完美形状”，对于像膜压力容器这样几何刚度化效应较大的结构，会变得非保守。对于壳体，建议通过添加初始缺陷（板厚0.1〜1倍的几何误差）进行非线性屈曲分析来确认。