ansys屈曲分析详细过程.doc

题目：跨径 L=89m，矢跨比 f/L=1/5 的圆弧拱，梁高 h/L=1/30,梁宽 b/L=1/15 求：1.弹性屈曲荷载； 2.非线性极限承载能力。 1、 线性屈曲荷载理论计算 在理论计算时，先根据圆弧拱的矢跨比查出稳定系数 K ： 2 表 1 圆弧拱理论计算的稳定系数 f/l 0.1 0.2 0.3 0.4 0.5 K 58.9 90.4 93.4 86.7 64.0 2 根据表 1 查得： K2 ? 90.4 故其理论弹性屈曲荷载为： 1 3.25?104 ? ?5933.33? 2966.673 EI N 12 q K l m ? x ? 90.4? ? 5.38?107 cr 2 3 3 89000 2、拱的弹性屈曲与非线性屈曲 对于一般的特征值屈曲分析，主要是在平衡状态，考虑到轴向力或者中面内 力对弯曲变形的影响，由最小势能原理，结构弹性屈曲分析归结为求解特征值问 题： ?K ?? ??K ?? 0 D G 通过特征值分析求得的解有特征值和特征向量，特征值就是临界荷载系数， 特征向量是临界荷载系数对应的屈曲模态。 特征值屈曲分析的流程图如下： 1 图 1 弹性屈曲分析流程图 非线性屈曲分析是考虑结构平衡受扰动（初始缺陷、荷载扰动）的非线性静 力分析，该分析是一直加载到结构极限承载状态的全过程分析，分析中可以综合 考虑材料塑性和几何非线性。结构非线性屈曲分析归结为求解矩阵方程： ??K ???K ???????F? D G 非线性屈曲分析的流程图如下： 图 2 非线性屈曲分析流程图 2 3、非线性方程组求解方法 （1）增量法 增量法的实质是用分段线性的折线去代替非线性曲线。增量法求解时将荷载 分成许多级荷载增量，每次施加一个荷载增量。在一个荷载增量中假定刚度矩阵 保持不变，在不同的荷载增量中，刚度矩阵可以有不同的数值，并与应力应变关 系相对应。 （2）迭代法 迭代法是通过调整直线斜率对非线性曲线的逐渐逼近。迭代法求解时每次迭 代都将总荷载全部施加到结构上，取结构变形前的刚度矩阵，求得结构位移并对 结构的几何形态进行修正，再用此时的刚度矩阵及位移增量求得内力增量，并进 一步得到总的内力。 （3）混合法 混合法是增量法和迭代法的混合使用。在一般的非线性分析中常采用增量迭 代混合法，将荷载分成若干级增量，在每一荷载增量上进行多次迭代。混合法综 合了增量法、迭代法的优点，并且与单纯的迭代法相比，混合法并不增加太多的 迭代次数。 4、曲梁加载问题 曲梁径向和切分布荷载可在圆柱坐标系下直接施加，而非径向和切向的分布 荷载可以将荷载等效到沿曲梁轴向分布，然后分解为径向和切向分布两部分施 加，其分解后为： Q Q / (1 y ) 和Q Q y / (1 y ) ? ? ?2 ? ? ? ?2 SV H SP H ? ? ? ? ? 其中，y =df (x) / dx,当单元足够小时，可采用y y / y 图 3 曲梁均布荷载等效与分解 3 5、问题的描述 矩形截面圆弧拱桥，在梁上受均布荷载，荷载及截面尺寸如图 1 所示，要求 用 ANSYS 进行拱圈的非线性屈曲分析，求极限荷载。本作业进行了圆弧拱的弹 性屈曲分析和非线性屈曲分析，分线性分析时引入初始几何缺陷。 图 4 荷载及拱的几何尺寸示意 几何尺寸： f/L=1/5， H/L=1/30，B/L=1/15 L=89m，f=17.8m，H=2.97m，B=5.93m，R=64.5m 材料性能：拱圈采用 C40 混凝土，弹性模量为 32500MPa，泊松比为 0.3。 在 ANSYS 中进行建模分析，初始荷载为 q=100000kN/m 其过程如下： 基本步骤： ? 建模； ? 获得静力解； ? 线性屈曲分析； ? 打开大变形开关,进行非线性分析； ? 查看结果。 Ansys 命令流如下： finish /clear /prep7 /units,si l=89 !定义圆弧拱尺寸 r=64.525 b=l/15 4 h=l/30 et,a,beam189 mp,ex,1,3.25e10 !材料 C40 mp,prxy,1,0.3 sectype,1,beam,rect secdata,b,h !定义梁截面几何数据 k,1 k,2,89 k,3,44.5,-1 larc,1,2,3,r !建立圆弧拱 dk,1,all !约束拱角位移 dk,2,all latt,1,,1,,3,,1 lesize,all,,,60 lmesh,all !划分网格 p=1e8 !初始荷载 *get,ne,elem,,count *do,i,1,ne ni=nelem(i,1) nj=nelem(i,2) dy=ny(nj)-ny(ni) dx=nx(nj)-nx(ni) dyx=dy/dx dyx2=s