[所有分类]网壳结构稳定基本理论—052.ppt

有实际工程意义的点： 上临界点、下临界点 上临界点：数值远高于 后面临界点，整条曲 线呈退化趋势。荷载 一旦超过上临界点， 结构一跨到底。从变形 角度，上临界点对应 的位移已经很大 下临界点：反应缺陷 敏感程度 Oran 沈院士 两个矩阵区别在于是否考虑两个主轴方向弯曲的耦合作用。 9.3.2 节点角位移的修正 对于大转角问题，刚性节点的最终转动位置取决于节点绕坐标轴的转动次序，即转动位移不满足交换律—— 增量计算中，角位移不能像线位移那样简单叠加，必须给予修正 在三维空间中，变形后的节点方向可以由与节点刚性相连的三条直线相对于坐标轴的方向余弦表示，即可唯一确定任意时刻的节点方向。 设定这三条直线平行于三个坐标轴，其方向余弦为3╳3的正交矩阵，随节点转动而变化——“节点方向矩阵” 假定一个节点绕三个坐标轴的转角 9.3.3 单元的转换矩阵 难点： 线性分析：局部坐标被定义为杆端截面主轴方向 非线性分析：每步增量计算中两个杆端截面相对位置变化 ——引入“单元方向矩阵”、“杆端方向矩阵”、“节点方向矩阵” “杆端方向矩阵”——确定杆端截面方向 [P(i)]、[P(j)]。规定矩阵列元素代表杆端截面的垂直轴和两个主轴对整体坐标轴的方向余弦 单元方向矩阵： 9.4 平衡路径的跟踪 9.4.1 全过程路径的计算方法及迭代策略 关键因素： 1、计算的累积误差问题——自由度数多，累积误差大；计算步数多，累积误差大。导致病态矩阵 2、计算的CPU问题——计算时间多长导致无法实施——自动变步长 实施： 1、第一步荷载增量法，得到ΔtU 2、第二部开始可以采用柱面弧长法，将ΔtU带入式（2.59）， 得到弧长增量Δl，再在每级荷载的计算开始首先计算增量 弧长： 方程（2.59）是二次方程，有两种情况： 9.4.2 临界点的判别准则 K的正负或0标志结构平衡状态的稳定性 如果k级荷载D矩阵主元皆为0，而k+1级荷载D矩阵主元 出现负值——k+1级荷载越过了临界点 进一步判别临界点类型： 1、PkPk+1——极值点 2、PkPk+1—— Pk+1Pk+2——极值点 Pk+1Pk+2——分枝点 9.4.3 分枝路径的跟踪——扰动荷载法 极值点屈曲：前、后位移性态一致——路径跟踪不需要采取措施 分枝屈曲：扰动荷载法（少自由度） 计算的累积误差（多自由度） 扰动荷载法：在分枝点采用非比例加载，即计算中引入与该临界点屈曲模态相一致的扰动荷载。 具体做法：第k+2级荷载判断为临界荷载——退回k级 荷载，子空间迭代法求出此时失稳模态 {Y}，再引入扰动荷载。 式（2.75）的各参数： 过渡点的个数N影响因素： （1）计算步长。当Fn偏向分枝路径时取消扰动荷载，否则跟踪会原路返回。 （2）屈曲前的Pk与Pc的接近程度（越接近N越小） 计算经验：N=2~5 β——扰动系数： 太小：计算返回 太大：计算发散 主要影响因素： Pk与Pc的接近程度 当Pk与Pc非常接近时，无论β取值怎样小，计算都难以收敛——由Pk再退回Pk-1，重新计算。 β取值根据经验。 自由度数超过1000的结构：计算累积误差可视为干扰。靠改变迭代收敛值的大小控制计算路径。要进入分枝路径，收敛值适当放宽。 例1： 采用全过程分析，无需设定 失稳类型，各种失稳形态可 自然表现。 例1：24单元穹顶 例2：342单元K6球壳 9.1.1 基本概念 1、失稳现象 几何失稳——几何可变 约束失稳——刚体位移 变形失稳——大变位或位移跳跃 2、屈曲类型 极值点屈曲 分枝点屈曲 第9章 网壳结构的稳定性分析 3、初始缺陷敏感性 经典试验：均匀轴压钢圆柱壳体 试验屈曲荷载为经典解的30%，离散性大。 根本原因：初始几何缺陷（会使分枝屈曲转化为极值屈曲） 9.1.2 网壳结构的失稳模态 失稳模态：失稳后因大变形而形成的新的几何形状 影响因素：网壳类型、几何形状、荷载条件、边界条件、 节点刚度 杆件失稳：单根杆件发生屈曲而结构其余部分不受任何 影响（两端铰接压杆模型） 点失稳：节点出现很大的几何变位、偏离平衡位置（某个节点受力偏大、某节点具有初始缺陷） 条状失稳：网壳结构的某个方向出现一条失稳带（柱壳沿一条母线、球壳的一个环带） 整体失稳：结构大部分发生很大的几何变位、偏离平衡位置。（薄膜应力到弯曲应力；呈波状，波长大于单根杆件长度；整体失稳往往从局部失稳开始） 耦合失稳：不同失稳模态间的相互作用严重——点失稳或整体失稳的临界荷载远大于单根杆件