利用微分方程求解单元刚度矩阵.pdfVIP

利用微分方程求解单元刚度矩阵 1,2 1 1 甘荣飞 ，曹文龙 ，孙靖立 1 2 （ 吉林建筑大学城建学院 吉林长春 130111， 内蒙古工业大学 内蒙古呼和浩特 010051） Using DifferentialEquationMethodtoSolvethe Using DifferentialEquationMethodtoSolvethe UUssiinnggDDiiffffeerreennttiiaallEEqquuaattiioonnMMeetthhooddttooSSoollvveetthhee ElementStiffnessMatrix ElementStiffnessMatrix EElleemmeennttSSttiiffffnneessssMMaattrriixx a,b a a GanRong-fei CaoWen-long SunJing-li a ( TheCity College ofJinlin JIANZHU University ,Changchun130111,China; bInnerMongoliaUniversityofTechnology,Hohhot010051,China ) 摘要： 摘要： 摘摘要要：：求解单元刚度矩阵是有限元求解问题的必要的过程，通常采用虚位移原 理或是最小势能原理进行求解。本文将以一维的拉（压）杆为例采用求解微分方 程的方法推导其单元刚度矩阵。根据其原理也可推广到其他形式单元中求解单元 刚度矩阵。 Abstract: Abstract: AAbbssttrraacctt::It’s necessaryprocessforthefiniteelement solutiontosolvetheelement stiffnessmatrix,whichissolved bytheprincipleofvirtualdisplacement andthe principleofminimumpotentialenergy.Thepapertakes theone-dimensionalpullbar (pressure)tosolveitselement stiffnessmatrixbythedifferential equationmethod. Themethodcan bealso appliedtosolveotherformsoftheelementstiffnessmatrix. 关键字： 关键字： 关关键键字字：：单元刚度矩阵；微分方程；有限元 中图分类号：O122 文献标识码： A Keyword Keyword KKeeyywwoorrdd：TheElementStiffnessMatrix;Differential Equation;theFiniteElement 1 引言 1 引言 11 引引言言 利用有限元方法对结构进行分析，是将结构进行离散化，形成一系列的单元 集合，这些单元在节点与相邻单元连接在一起，要确定整个结构的特性，必须先 确定个别单元的的特性。通常根据虚位移原理或是最小势能原理来推导得到单元 刚度矩阵，单元刚度矩阵是将单元节点位移和单元节点力、单元等效节点荷载联 系起来的联系矩阵，是利用有限元求解问题的必要过程。本文将介绍除了根据虚 位移原理和最小势能原理推导单元刚度矩阵之外，也可以由求解微分方程来推导 单元的刚度矩阵。 2由微分方程推导单元刚度矩阵 2由微分方程推导单元刚度矩阵 22由由微微分分方方程程推推导导单单元元刚刚度度矩矩阵阵 按单元的应力 位移的微分关系，建立关于位移函数的微分方程，在一定的- 边界条件下求解，即可得到单元节点力-位移的关系式。 1 2