练习13：支架屈曲分析
本案例是一个支架面内受载屈曲问题，介绍引入初始扰动至后屈曲分析中的新方法。
问题描述
支架一端固定，受集中荷载作用；尺寸如图，厚0.004m。
材料信息
Aluminum；杨氏模量70GPa；泊松比0.3.
工作目录
几何模组
上图的分割是为了布局网格，对于涉及弧形或者曲边几何，我们考虑其角度不超过90°（即圆形考虑四等分），否则加以分割。按照上图分割，有利于布置弯曲段的网格，保证网格均匀不歪曲。
分割有许多原因，但这里是出于网格考虑。软件自动将发了个以后的区域绑定起来，所以无需担心断开。
属性模组
装配模组
分析步模组
切换至Step模组，建立屈曲分析程序，请求输出5个模态。源文件已经帮大家定义好了，观察下图：
几何非线性是关闭的，且不允许我们改变，因为我们执行的是特征值屈曲分析，不考虑非线性因素。如果你给非线性因素，如非线性本构或者接触，软件仍然以线性求解。
定义载荷
定义边界
交互模组
网格模组
使用的单元类型为S8R5。
分析任务
可视化模组
切换至Visualization模组
整体位移
模态跳转
一阶由于方向与施加的方向相反才可能产生对应的模态，所以这里忽略一阶，而二阶刚好与一阶同阶（且阶次最低），所以选择二阶作为最低感兴趣的阶次。
载荷位移分析
引入初始扰动
在*Output关键字后（这个模块后，不是这个keyword紧挨着后面），使用Add After放置光标（你需要点击一下此项）。
输入如下字段：
*NODE FILEU注意，我们是修改特征值屈曲分析的模型。*NODE FILE用于选择节点变量，将其写入到结果文件。
修改载荷位移分析模型文件，目的是从结果文件读取节点位移。同样的进入方法，在边界*Boundary模块修改。*Step之前使用Add After放置光标位置，输入如下字段：
*IMPERFECTION, FILE=Eigen, STEP=12,0.000004
*IMPERFECTION用于引入几何扰动到后屈曲分析的模型中，Eigen为FILE的参数（特征值屈曲分析的结果文件名），STEP=1，即告诉软件具体的分析步。2，0.000004为屈曲模态的阶次和缩放系数（取第二节模态，缩放0.000004）。
屈曲模态缩放系数是怎么定的？设定缩放系数为0.1%的支架厚度（0.004X0.1%）。
以上就设置好了需要的内容，接下来就开始分析了。我们修改了特征值屈曲分析的模型，所以需要重新提交。点击对应的任务名，然后提交计算即可。
基于载荷位移分析模型，建立新的分析任务
可视化后处理
位移（U3）
面外位移很小，结构没有发生屈曲。可能的原因有：引入的扰动太小或者增量太大没有捕捉到屈曲（如果一步就跨过去呢），首先尝试缩小增量步。
设定初始增量为0.01，最大增量为0.05，重新计算。（注意，我们修改的是载荷位移分析，不是特征值屈曲分析）
位移大了很多，结构发生屈曲。有人说，为什么位移小或者没有位移就说结构没有屈曲，有相当的量就说结构发生屈曲。别忘了结构屈曲的定义，必须具有相当量的位移（这里相当量的量，是指与结构本身的特征尺寸相比较，不是你目测分析结果的大小）。
Force VS Displacement
为了进一步探究几何扰动的影响，这里将曲线重新命名为Imp_01P。
观察上图，载荷达到1.5E3之前，刚度极大。载荷位移分析预测的与特征值屈曲分析基本上一直（1444.4N），约误差3.7%，可以认为是一致的。
扰动影响
*IMPERFECTION, FILE=Eigen, STEP=12, 0.00004*IMPERFECTION, FILE=Eigen, STEP=12, 0.0004
模型文件下载：
引入初始扰动至后屈曲分析中的新方法，忒好奇的是初始屈曲模态和缩放因子是怎么确定呢，凭心而动？