几何非线性与屈曲分析.ppt

1,第五章非线性问题与屈曲分析,一、几何非线性问题基本理论二、几何非线性有限元方程的建立三、非线性结构分析简介四、非线性分析的基本过程五、几何非线性六、屈曲分析七、材料非线性,第7讲,2,一、几何非线性问题基本理论,3,1有限应变与应力,4,2变形率和本构关系,5,2变形率和本构关系,6,二、几何非线性有限元方程的建立,7,二、几何非线性有限元方程的建立,8,9,2结构稳定性和屈曲问题,10,2结构稳定性和屈曲问题,11,1线性稳定分析,12,2.非线性稳定分析,SeeP128-130,13,三、非线性结构分析简介,1非线性结构,非线性,14,2非线性分析的特殊性,（1）逐步递增载荷和平衡迭代,非线性问题需要一系列带校正的线性近似来求解,15,3非线性求解的组织级别,非线性求解级别：,载荷步子步平衡迭代,16,4收敛容差,当确定收敛准则时，可以选择将收敛检查建立在力、力矩、位移、转动或这些项目的任意组合上，且每一个项目可以有不同的收敛容差值。对多自由度问题，同样也有收敛准则的选择问题。,17,5保守行为与非保守行为过程依赖性,18,6子步7自动时间步长8载荷和位移方向9非线性瞬态过程的分析,（1）子步数或时间步长（2）自动时间步长,19,四、非线性分析的基本过程,非线性分析的步骤,1建模2加载并求解3检查结果,20,非线性瞬态分析实例,21,/prep7et,1,visco106keyopt,1,3,1mp,ex,1,117e9mp,dens,1,8930mp,nuxy,1,0.35tb,biso,1,1,0tbmodif,2,1,4e8tbmodif,3,1,le8tbplot,biso,1rectng,0,0.0032,0,0.0324lesize,2,20,1lesize,1,4,1mshape,0,2dmshkey,1amesh,1fini/soluantype,4trnopt,fulllumpm,onnlgeom,1ic,all,uy,-227d,all,0,uznsel,s,loc,x,0d,all,0,uxnsel,s,loc,y,0d,all,0,uy,allseltime,8e-5autots,1deltim,4.4e-7kbc,1outres,all,4solvefini/postlset,lasttop_node=26*get,deform,node,top_node,u,y/dscale,1,1pldisp,2plnsol,epto,eqvfini/post26nsol,2,top_node,u,y,displacederiv,3,2,1,velocity,1/axlab,y,velocityplvar,3fini,建模,网格,22,等效应变,合位移m,23,X方向位移,等效应力Pa,24,25,不同时刻等效应力及变形,26,不同时刻等效应力及变形,27,28,五、几何非线性,(一)、简介,如果结构经受大变形，变化的几何形状可能会引起结构的非线性响应。,29,(二)、几何非线性分析,30,1大应变效应,一个结构的总刚度,当一个单元的节点经历位移后，那个单元对总体结构刚度的贡献可以以两种方式改变。,31,2关于大应变的特殊建模讨论,（1）应力应变,在大应变求解中，所有应力应变输入和结果将采用真实应力和真实（或对数）应变。,32,(2)单元的形状,在大变形分析的任何迭代中粗劣的单元形状（即大的横纵比，过度的顶角及具有负面积的已扭曲单元）是不好的，所以必须象注意单元的原始形状一样注意单元已扭曲后的形状。,33,3.具有小应变的大挠度（转角）,大挠度效应,在一个大挠度分析中，单元的转动可以任意地大，但是应变假定是小的。,用NLGEOM，ON命令来激活那些支持这一特性的单元中的大挠度效应。,GUI：MainMenuSolutionAnalysisOptions,34,4.应力刚化,它在薄的、高应力的结构中，如缆索或薄膜中，是最明显的。一个鼓面，当它绷紧时会产生垂向刚度，这是应力刚化结构的一个普通的例子。,35,5旋转软化,在小位移分析中这种调整近似于由于大的环形运动而导致几何形状改变的效应,36,六、屈曲分析,1.概念,ANSYS提供两种分析结构屈曲载荷和屈曲模态的技术,特征值（或线性）屈曲分析,37,2非线性屈曲分析的过程,（1）施加载荷增量（2）自动时间步长功能（3）重要注意事项（4）施加初始缺陷（扰动）（5）弧长法,38,3特征值屈曲分析的过程,39,如图5一15所示，一根长为L，两端铰支的细长杆受到轴向加载。此杆截面的高度为h，面积为A。由于对称性，只给杆的上端建模，则上半部分的边界条件变为自由一固支。为了描述屈曲模态，在x方向取10个主自由度。杆的惯性矩为I=Ah2120.0052083in4详细参数如下：材料性质：模型尺寸：载荷。EX=30E6psiL=200inF=1lbfA=0.25in2h0.5in,特征值屈曲分析实例,40,/PREP7/TITLE,BUCKLING.ET,1,BEAM3R,1,0.25,52083E-7,0.5MP,EX,1,30E6N,1N,11,100FILLE,1,2EGEN,10,1,1FINISHSave,buckle1,dbResume,buckle1,db/SOLUANTYPE,STATICPSTRES,OND,1,ALL,11,41,F,11,FY,-1SOLVEFINISHSAVE,BUCKLE1,DBRESUME,BUCKLE1,DB/SOLUANTYPE,BUCKLEBUCOPT,REDUC,1MXPAND,1M,2,UX,11,1SOLVE*get,FCR,MODE,1,FREQ!FcR-曲屈分析中第1阶频率为失稳载荷NFINISH/POST1SET,FIRSTPLDISP,1FINISH,42,43,fini/clear/prep7antype,transnlgeom,onet,1,beam3,1r,1,161.29,2167.86,12.7mp,ex,1,2.07e5mp,prxy,0.3!n,1n,11,2540fille,1,2egen,10,1,1finish/ESHAPE,10/solusolcontrol,0neqit,250!outpr,basic,lastoutres,all,alld,1,allfcr=-171.5pi=3.1415926!time,1DELT,0.05,0.005,0.5AUTO,ONkbc,0f,11,fy,fcr*1.884f,11,fx,1solve!,超过失稳载荷的动力计算,44,七、材料非线性,(一)、弹塑性分析,1有关塑性的概念,(1)路径相关性(2)率相关性(3)工程应力、应变与真实的应力、应变,ANSYS的材料非线性分析,45,2塑性选项,1).经典双线性随动强化BKIN2).双线性等向强化BISO3).多线性随动强化MKIN4).多线性等向强化MISO,注意，使用MP命令来定义弹性模量；弹性模量也可以是与温度相关的；切向斜率Et不可以是负数，也不能大于弹性模量。,46,3.塑性分析实例,一个周边简支的圆盘，在其中心受到一个冲杆的周期作用。由于冲杆被假定是刚性的，因此在建模时不考虑冲杆，而将圆盘上和冲杆接触的结点的Y方向上的位移耦合起来。由于模型和载荷都是轴对称的，因此用轴对称模型来进行计算。求解通过四个载荷步实现。问题详细说明：材料性质：EX70000MPa(杨氏模量)NUXY：0.325(泊松比)塑性时的应力一应变关系如下：应变应力MPa0.0007857550.005751120.029251720.1241加载历史：时间s载荷N001-600027503-6000,47,