ANSYS几何非线性和线性对比分析的一个工程实例

1、ANSYS非线性和线性对比分析的一个工程实例二力杆几何非线性分析土木工程中，钢模板由于制作不精细或搬运模板过程受到碰撞或者挤压等外力作用常常会造成模板某处凸起，在活、恒载作用下或搬运过程中，该处常常会突然从凸起变成凹进状态。这一现象被称为油罐效应，通常采用桁架的失稳模型进行几何非线性简化分析因为也称为桁架的经典跳越问题。分析模型如图1所示，采用LINK1单元构成二力杆，两端完全约束，中间节点加集中力。图1荷载与顶点位移理论关系为：P=EA（）x（x-sin）（x2sin）式中，E为弹性模量，Ao为杆件初始截面积，x=V/Lo,V为顶点的竖向位移，Lo为杆长，6为杆件倾角。分析中取A0=10mm

2、2,E=200Gpa,Lo=100mm,臼=6。分析所用的非线性命令流如下。*二力杆几何非线性分析*L, FINI/CLEAR/FILENAME,NONLINEAR_ER-LI-GAN/TITLE,TheAnalysisofNONLINEAR_ER-LI-GAN/PREP7/PNUM,LINE,1/PNUM,KP,1LO=100CTA=6*AFUN,DEGL1=2*LO*COS(CTA)H1=LO*SIN(CTA)AA=10EM=2E5ET,1,LINK1MP,EX,1,EMR,1,AAK,1K,2,0.5*L1,H1K,3,L11,2$L,2,3LESIZE,ALL,1LMESH,ALLFI

3、NI/SOLUDK,1,ALLDK,3,ALLFK,2,FY,-1200ANTYPE,0NLGEOM,1NSUBST,100OUTRES,ALL,ALLARCLEN,ONSOLVEFINI/POST1/ESHAPE,1EPLOTSET,LASTPLDISP1!绘制变形图PRRSOLFINI/POST26NSOL,2,2,U,Y,DISPLACEMENTABS,3,2RFORCE,4,1,F,Y,F/AXLAB,X,DISPLACEMENT/AXLAB,Y,F/GRID,1XVAR,3PLVAR,4!绘制节点2位移和节点1的竖向支反力关系FINI通用后处理中可得到最后一个荷载步时节点1和节点3受

4、到的水平方向支座反力分别为-4759.0N和4759.0N,双杆均受拉。同时可得到结构变形图，如图2所示，并且由图可知此时节点2的位移最大为23.001mm,基本与实际情况一致。AIMJUN232-01021sIS;Z7._iaeJDrSPLACEMENTSTEP=1TTME=J.DMX=23.TheAnalysisofNONLZNEAR_ER-LZ-GAN图2由时间一历程后处理分析可知随着时间的增长，荷载的增加，节点2位移变大，节点1的竖向反力也不断增加，二者之间关系见图3。由图3可知位移一荷载曲线与理论公式分析基本符合，但随着杆件倾角B的增大，曲线两个顶点附近的误差会越来越大，这与理论公式

5、的假定有关，所以，该类模型只能模拟倾角比较的构件。PO3T26JUN23201021:19:IDTheAnalysisofNONLIWEAR_ER-LI-GAW以上问题是一个纯粹的几何非线性问题，分析可知在整个结构的变形过程中，尽管结构的位移很大，但是杆件的应力和应变之间一直保持着线性关系，即正比例关系，这也同实际相同，故该分析所得结果比较合理。二力杆线性分析分析所用的非线性命令流如下*以下为线性分析命令流*FINI/CLEAR/FILENAME,ER-LI-GAN/TITLE,TheAnalysisofer-li-gan/PREP7/PNUM,LINE,1/PNUM,KP,1LO=100CT

6、A=6*AFUN,DEGL1=2*LO*COS(CTA)H1=LO*SIN(CTA)AA=10EM=2E5ET,1,LINK1MP,EX,1,EMR,1,AAK,1K,2,0.5*L1,H1K,3,L11,2L,2,3LESIZE,ALL,1LMESH,ALLDK,1,ALLDK,3,ALLFK,2,FY,-1200FINI/SOLUDK,1,ALLDK,3,ALLFK,2,FY,-1200ANTYPE,0SOLVEFINI/POST1/ESHAPE,1EPLOTPLDISP.1!绘制变形图PRRSOL!列出节点力FINI通用后处理中可得到节点1和节点3受到的水平方向支座反力分别为5708.6N和-5708.6N。同时可得到结构变形图，如图4所示，并且由图可知此时节点2的位移最大为2.75mm,节点2还停留在X轴上方，远远小于非线性分析时的位移，同时杆件仍然受压，这与实际情况不符，所以对于这种经典的跳越式二力杆问