小弟最近在做一个铝制穹顶结构，跨度125.3m，矢高44.5m，截面为工字型梁，有限元计算模型见图1，荷载模式为1.0倍恒载+1.0倍活载
（1）采用beam188单元，计算时建立表面效应单元surf154覆于结构之上，以便与施加荷载；
（2）进行特征值屈曲分析，第一阶特征值屈曲模态见图2；
（3）几何非线性分析：以第一阶屈曲模态为基础对结构引入L/300的初始缺陷，进行几何非线性分析，荷载因子随分析步的变化曲线见图3；
（4）双非线性：施加初始缺陷，选取铝的材料特性模型为双线性等向强化模型，设置其屈服强度200GPa，切线模量假设为0，荷载因子随分析步的变化曲线见图4

问题1：为什么特征值屈曲因子小于非线性分析中得到的荷载因子？是结构对所施加的缺陷不敏感吗？
2：非线性分析中得到的荷载因子相对于所有荷载吗？就是说极限承载力是荷载因子*（恒载+活载）还是荷载因子*活荷载？
3：进行非线性分析时，我没有用弧长法（总是不收敛），用的是N-R法，但是没有得到下降段，请问采用N-R法是否合理呢？
4：对于铝材料的本构模型的选取是否合理，其本构关系曲线为标准弹塑性本构关系，而铝的实际本构关系曲线并不是这样的


希望前辈们不吝赐教，在此先谢过。

1.材料非线性会使结构的极限承载力系数降低。
2.相对于所有的荷载。
3.荷载-位移曲线，曲线处于明显的上升段即突然中止，这往往是计算过程中由于各种原因导致的不收敛使计算无法继续，将该点对应的荷载作为稳定承载力并不合适；尽管由于采用Newton-Raphson法无法跟踪曲线的下降段，但合理的曲线应是逐渐接近水平，意味着刚度矩阵接近奇异，结构临近失稳。建议作者对上述情况重新分析计算，也可考虑改用弧长法，确保分析结果的正确。
4.影响不大。

楼主的surf154单元是怎么加上去的？