几何非线性微分刚度有限元分析.docVIP

几何非线性微分刚度有限元分析 　　摘 要：本文对大位移几何非线性问题作有限元分析，将简单情形分析的基本结论推广到一般情形，给出一般情形几何非线性问题微分刚度的计算方法，并对微元体、杆单元及等截面梁单元的微分刚度进行计算，这种微分刚度计算对于大位移几何非线性问题的有限元分析具有非常重要的理论价值及应用价值。 　　关键词：有限元方法；大位移；广义力；微分刚度 　　中图分类号：TU32 文献标志码：A文章编号：1671-7953(2009)01-0103-03 　　 　　Finite Element Analysis of Differential Rigidity?? 　　About Geometry Nonlinearity?? 　　HU Ke?? 　　（Guangdong Institute Of Science And Technology Zhuhai Guangdong 519090，China） 　　Abstract: The article has made a finite element analysis about large displacement geometry nonlinearity question，spreading simple analyzing result to an ordinary way，calculating common method has been givens out about geometry nonlinearity differential rigidity，as an example it has been set up differential rigidity of small element，bar element and constant cross section beam element，in theory and use it is of great value to that calculating differential rigidity about large displacement geometry nonlinearity question. 　　Key words: Finite element；Large displacement；Generailzed forces；Differential rigidity. 　　 　　对几何非线性问题的有限元分析，计算大位移引起的附加刚度效应，比较简单的方法是微分刚度计算法，如图1所示，刚性杆AB长为l，A端铰链连接并有一个扭转弹簧（扭转刚度为K），B端有一水平荷载P。?? 　　图1 　　当杆由AB位置转过θ角到达AB′位置时，P力做功为W=Pu??B=Pl（cosθ-1），弹性势能为V=1/2Kθ 2，以转角θ为广义坐标，相应于弹性势能的广义力为?? 　　Q??θ=-??V?功?=-Kθ（1） 　　而相应于P力做功的广义力为?? 　　Q??θ=??W?功?=-Plsinθ（2） 　　若考虑小转动，即θ为小量（sinθ≈θ），则有Q??θ=-Plθ，由（1）式与（2）式可见，Pl与刚度K相当，这里Pl就是P力在小转动过程中附加的刚度，称为微分刚度Kd，通常线刚度K由（1）式得?? 　　K=-??Q??θ?功? 　　?ね?样，微分刚度由（2）式得?? 　　Kd=-??Q??θ?功?=-??2W?功?2=Pl 　　将上述结论推广到多自由度弹性系统，设想将一般弹性结构离散成具有n个自由度弹性系统，其上作用m个力P??i（i=1，2，…，m），沿作用力方向位移为u??i，则m个力总功为?? 　　W=??mi=iP??iu??i 　　?び氲?r个广义坐标q??r相应的广义力为?? 　　Q??r=??W??q??r=??mi=iP??i??u??i??q??r；（r=1，2，…，n） 　　式中假定P??i不随位移而变化，但位移u??i可以是广义坐标的函数，第r个自由度对第s个自由度的微分刚度为?? 　　Kdrs??=-??Q??r??q??s=-??2W??q??r??q??s；（r，s=1，2，…，n）（3） 　　这就是一般弹性体出现大位移时，在线性刚度阵上附加的微分刚度。显然它是一个对称矩阵，即Kdrs??=Kdsr??，式中功W本来是外荷载的功，但把弹性结构离散成n个自由度系统后，可用与外荷载等效的内力的功来代替。 　　1 微元体的微分刚度 　　微元体单元如图2所示，它有三个移动自由度和三个转动自由度（ω??x、ω??y、ω??z），在小应变条件下，应力在移动位移中所做功之和与转动位移相比可以略去不计。?? 　　图2 　　 　　几