ansys分析钢结构稳定问题

ANSYS软件分析轴压和压弯构件的稳定性问题2摘要轴心受压杆件和压弯杆件广泛应用于工程中，本文通过ANSYS软件对该两种杆件进行分析，对于轴心受压杆件，运用BEAM189、SOLID95、SHELL65单元，进行弹性稳定分析和非线性分析，得到其屈曲荷载和变形情况；对于压弯杆件，在集中荷载和分布荷载的条件下，运用BEAM3单元进行非线性分析，得到其最大弯矩值，通过和理论值相比较，验证其正确性。关键词ANSYS；轴心受压杆件；压弯杆件；非线性分析ABSTRACTAXIALSTRUTPIECESANDBENDINGRODSAREWIDELYUSEDINENGINEERINGTHISPAPER,USINGANSYSSOFTWARE,ANALYZESTHETWORODSFORCENTRALLYCOMPRESSEDMEMBERS,THISPAPERUSINGBEAM189,SOLID95,SHELL65UNIT,CARRIESOUTELASTICSTABILITYANALYSISANDNONLINEARANALYSIS,GETTINGTHEBUCKLINGLOADANDDEFORMATIONFORTHEBENDINGRODUNDERCONDITIONSOFCONCENTRATEDLOADSANDDISTRIBUTEDLOADS,NONLINEARANALYSISWASCONDUCTEDUSINGBEAM3UNIT,GETTINGITSGREATESTMOMENT,ANDWASCOMPAREDTOTHEORETICALVALUETOVERIFYITSCORRECTNESSKEYWORDSANSYSCENTRALLYCOMPRESSEDMEMBERSTHEBENDINGRODMEMBERNONLINEARANALYSIS钢材具有高强度、质轻、力学性能良好的优点，是制造结构物的一种极好的建筑材料，所以广泛运用于工程实例中，它和钢筋混凝土结构相比，对于充任相同受力功能的构件，具有截面轮廓尺寸小、构件细长和构件柔薄的特点。对于因受压、受弯和受剪等存在受压受压区的构件或板件，如果技术上处理不当，可能使钢结构出现整体失稳或局部失稳。失稳前结构物的变形可能很微小，突然失稳使结构物的几何形状急剧改变而导致结构物完全丧失抵抗能力，以致整体塌落。钢结构的稳定性能是决定其承载力的一个特别重要的因素1。对于钢结构稳定性的研究也就极其重要。而轴压杆件和压弯杆件是钢结构的基础，对此杆件进行稳定性分析也就是不可避免的和尤为重要的。所以，非常有必要利用大型通用ANSYS软件对这两类杆件进行分析，得到一系列的研究成果。一、基本理论结构在荷载作用下由于材料的弹性性能而发生变形，若变形后结构上的荷载保持平衡，这种状态称为弹性平衡。如果结构在平衡状态时，受到扰动而偏离平衡位置，当扰动消除后仍能恢复到原来平衡状态的，这种平衡状态称为稳定平衡状态。根据失稳的性质，结构的稳定问题可以分为平衡分岔失稳，极值点失稳和跃越失稳三种情况。结构的弹性稳定分析属于平衡分岔失稳，在ANSYS中对应的分析类型是特征值屈曲分析（BUCKLINGANALYSIS）2。3关于特征值屈曲分析有以下说明（1）该分析对结构临界失稳力的预测往往要高于结构实际的临界失稳力，因此在实际的工程结构分析时一般不用特征值屈曲分析；（2）特征值屈曲分析能够预测临界失稳力的大致所在，因此在做非线性屈曲分析时所加力的大小便有了依据；（3）特征值屈曲分析所预测的结果我们只取最小的第一阶；（4）特载值分析得到的是第一类稳定问题的解，只能得到屈曲荷载和相应的失稳模态，它的优点就是分析简单，计算速度快。基本步骤（1）创建模型；（2）获得静力解；（3）获得特征值屈曲解；（4）查看结果。注意事项为（1）在建模时,仅考虑线性行为；定义材料的弹性模量或某种形式的刚度；另外,单元网格对屈曲荷载系数影响很大；（2）在获得静力解时,必须激活预应力效应由屈曲分析得到的特征值是屈曲荷载系数,在此直接施加单位荷载,得到的屈曲荷载系数即屈曲荷载；（3）若想用命令流获取第N阶模态的特征值（屈曲荷载系数）直接采用以下命令流GET,FREQN,MODE,N,FREQ，其中FREQN为用户定义的变量，存放第N阶模态的屈曲荷载系数。结构的非线性问题可以分为几何非线性、材料非线性、状态非线性三种情况。在此题中我们主要考虑几何非线性。其基本步骤如下（1）创建模型；（2）设置求解控制参数，包括设置分析类型和分析选项，设置时间和时间步，设置输出控制，设置求解器选项，设置重启动控制等；（3）加载求解，注意变形前后荷载的方向；（4）查看结果。因几何变形引起结构刚度改变的问题属于几何非线性问题。通常分为大应变，大位移和应力刚化三类。导致结构刚度3变化的原因如下（1）单元形状改变（如面积、厚度等），导致单刚变化；（2）单元形状改变（如大转动），导致单刚向总体系坐标系下转换时发生变化；（3）单元较大的应变使得单元在某个面内具有较大的应力状态，从而显著的影响面外的刚度。在分析时，应该注意单元选择，单元形状，网格密度，荷载和边界条件等问题2。二、轴心受压杆件分析采用两端简支的受压柱，设截面尺寸BH003M005M,柱长L3M,弹性模量E210GPA根据欧拉临界公式4，其临界荷载为2217234052903EIPCRKNL4采用BEAM189单元时，需要约束绕单元轴的转动自由度，否则虽可进行静力分析，但是会出现异常屈曲模态（模态分析时会出现零值）。BEAM189是3D二次有限应变梁。BEAM188/189是不支持跨间集中荷载和跨间部分分布荷载，仅支持在整个单元长度上分布的荷载。采用SOLID95和SHELL63模拟此模型时，仅在下端截面中心约束Y方向的平动自由度，而不能约束整个截面，否则与简支约束条件不符。在SOLID95单元上，施加的为面荷载，在SHELL63上施加的为线荷载。其杆件模型如图一所示图一，轴心受压杆件21采用BEAM189单元其弹性稳定分析的命令流如下FINISH/CLEAR/PREP7B003H005L3E21E11ET,1,BEAM189MP,EX,1,EMP,PRXY,1,03SECTYPE,1,BEAM,RECTSECDATA,B,HK,1K,2,LK,10,0,L/2,L/2L,1,2DK,1,UX,UY,UZ,ROTYDK,2,UX,UZ,ROTYLATT,1,1,10,1LESIZE,ALL,20LMESH,ALLFINISH/SOLUFK,2,FY,1PSTRES,ONSOLVEFINISH/SOLU5ANTYPE,BUCKLEBUCOPT,LANB,5MXPAND,5OUTRES,ALL,ALLSOLVEFINISH/POST1SET,LIST得出以下图形图二,杆件一阶变形图和位移矢量图图三,杆件二阶、三阶变形图图四，杆件四阶、五阶变形图非线性分析的命令流FINISH/CLEAR/PREP76B003H005L3E21E11ET,1,BEAM189MP,EX,1,EMP,PRXY,1,03SECTYPE,1,BEAM,RECTSECDATA,B,HK,1K,2,LK,10,0,L/2,L/2L,1,2LATT,1,1,10,1LESIZE,ALL,20LMESH,ALLFINISH/SOLUDK,1,UX,UY,UZ,ROTYDK,2,UX,UZ,ROTYFK,2,FY,27000F,10,FX,50PSTRES,ONSOLVEFINISH/SOLUANTYPE,STATICNLGEOM,ONOUTRES,ALL,ALLNSUBST,50AUTOTS,ONLNSRCH,ONSOLVEFINISH/POST26/GROPT,DIVY,10/COLOR,AXES,8/COLOR,CURVE,2/AXLAB,X,DEFLECTION/AXLAB,Y,FORCERFORCE,3,1,F,YNSOL,4,10,U,XXVAR,4PLVAR,3/AXLAB,X,FORCE/AXLAB,Y,DISPLACEMENT3RFORCE,5,1,F,YNSOL,6,2,U,YXVAR,5PLVAR,6图五,杆件中点荷载挠度曲线从上图可以得知，当荷载达到欧拉临界荷载时，挠度突然增加。图六，杆件顶点荷载位移曲线从上图可以看出，当荷载达到临界荷载时，该点位移突然增加，向下移动。把27000改为56000，得到屈曲前后的图像，注意，改后的值不能过大。FINISH/CLEAR/PREP7B003H005L3E21E11ET,1,BEAM1897MP,EX,1,EMP,PRXY,1,03SECTYPE,1,BEAM,RECTSECDATA,B,HK,1K,2,LK,10,0,L/2,L/2L,1,2LATT,1,1,10,1LESIZE,ALL,20LMESH,ALLFINISH/SOLUDK,1,UX,UY,UZ,ROTYDK,2,UX,UZ,ROTYFK,2,FY,56000F,10,FX,50PSTRES,ONSOLVEFINISH/SOLUANTYPE,STATICNLGEOM,ONOUTRES,ALL,ALLNSUBST,50AUTOTS,ONLNSRCH,ONSOLVEFINISH/POST26/GROPT,IVVY,10/COLOR,AXES,8/COLOR,CURVE,2/AXLAB,X,DEFLECTION/AXLAB,Y,FORCERFORCE,3,1,F,YNSOL,4,10,U,XXVAR,4PLVAR,3/AXLAB,X,FORCE/AXLAB,Y,DISPLACEMENTRFORCE,5,1,F,YNSOL,6,2,U,YXVAR,5PLVAR,6图七，杆件中点屈曲前后挠度荷载曲线图八，杆件顶点位移荷载曲线22采用SOLID95单元SOLID95弹性稳定分析的命令流FINISH/CLEAR/PREP7B003H005L3E21E11ET,1,SOLID958MP,EX,1,EMP,PRXY,1,03BLC4,B,L,HWPOFF,B/2,H/2VSBW,ALLWPROTA,90VSBW,ALLWPCSYS,1ESIZE,3/20VMESH,ALLDK,KPB/2,0,H/2,UYASEL,S,LOC,Y,0ASEL,A,LOC,Y,LDA,ALL,UXDA,ALL,UZASEL,S,LOC,Y,LSFA,ALL,1,PRES,1/B/HALLSEL,ALL/SOLUPSTRES,ONSOLVEFINISH/SOLUANTYPE,BUCKLEBUCOPT,LANB,5MXPAND,5OUTRES,ALL,ALLSOLVEFINISH/POST1SET,LISTSET,1,1PLDISP,1PLVECT,UPLNSOL,U,X,1图九，一阶变形图、位移矢量图图图十，二阶、三阶变形图23采用SHELL65单元SHELL65弹性稳定分析的命令流FINISH/CLEAR/PREP7B003H005L3E21E11ET,1,SHELL63MP,EX,1,EMP,PRXY,1,03R,1,BWPROTA,90BLC4,H,LWPCSYS,1WPOFF,H/29ASBW,ALLESIZE,3/20AMESH,ALLLSEL,S,LOC,Y,0LSEL,A,LOC,Y,LDL,ALL,UXDL,ALL,UZDK,KP0,0,H/2,UYLSEL,S,LOC,Y,LSFL,ALL,PRES,1/HALLSEL,ALL/SOLUPSTRES,ONSOLVEFINISH/SOLUANTYPE,BUCKLEBUCOPT,LANB,5MXPAND,5OUTRES,ALL,ALLSOLVEFINISH/POST1SET,LIST图十一，一阶变形图、位移矢量图三、压弯杆件分析31两端铰接横向荷载下压弯构件分析图十二，杆件简图、力学分析为计算分析方便，在此选用之前所建模型的数据，B003M,H005M,L3M所以，该构件的欧拉临界力为259077KN。由教材P78P801可知，跨中最大弯矩，其中MAX00178/12/EEPPMM04QLM10在此取P10000KN，Q500KN,计算得15167。375KN。所以，在构件MA0M中点处，有最大弯矩1516737556877KN。接着利用ANSYS软件建0MAM立此模型，并进行非线性分析，得出构件的最大弯矩，以此和计算所得数据相比较，验证软件分析的正确性。由于采用BEAM189模拟比较复杂，在此采用BEAM3模型，该单元为2D梁单元，分析比较简单，也能够较好的模拟此状态。命令流FINISH/CLEAR/PREP7CSYS,0B003H005L3E21E11A0BHI1BH3/12I2HB3/12MP,PRXY,03MP,EX,1,EET,1,BEAM3R,1,A0,I2,HK,1K,2,LK,3,L/2L,1,3L,2,3DK,1,UX,UYDK,2,UXLATT,1,1,1LESIZE,ALL,20LMESH,ALLFINISH/SOLUFK,2,FY,10000FK,3,FX,500ANTYPE,0NLGEOM,ONSUBST,50AUTOTS,ONLNSRCH,ONSOLVEFINISH/POST1ETABLE,MI,SMISC,6ETABLE,MJ,SMISC,12PLLS,MI,MJ图十三，弯矩图通过上图可知，构件最大弯矩位于中点处，其值为566385NM,通过公式计算为56877NM。二者相差不大。11图十四，变形图、位移矢量图32横向均布荷载下的压弯构件图十五，杆件简图、力学分析图构件尺寸保持不变，横向荷载500N/M，可知，其中通过计算可得MAX00128/1/EEPMMP208QLM5625NM,1646，926NM00/AMAX命令流FINISH/CLEAR/PREP7CSYS,0B003H005L3E21E11A0BHI1BH3/12I2HB3/12MP,PRXY,03MP,EX,1,EET,1,BEAM3R,1,A0,I2,HK,1K,2,LL,1,2DK,1,UX,UYDK,2,UXLATT,1,1,1LESIZE,ALL,20LMESH,ALLFINISH/SOLUFK,2,FY,10000SFBEAM,ALL,1,PRES,500ANTYPE,0NLGEOM,ONSUBST,50AUTOTS,ONLNSRCH,ONSOLVEFINISH/POST1ETABLE,MI,SMISC,6ETABLE,MJ,SMISC,12PLLS,MI,MJ12图十六，弯矩图通过ANSYS软件分析可知，在构件中点处获得弯矩最大，其值为922837NM,和公式计算所得相差不大。图十七，杆件变形图、位移矢量图33集中力作用下的压弯构件图十八，杆件简图构件尺寸保持不变Q300N，MAX0105/13/EEPQLMMP经计算得，498206NM。MAX命令流如下FINISH/CLEAR/PREP7CSYS,0B003H005L3E21E11A0BHI1BH3/12I2HB3/1213MP,PRXY,03MP,EX,1,EET,1,BEAM3R,1,A0,I2,HK,1K,2,LK,3,L/3K,4,2L/3L