有限元_4-薄板弯曲问题

《有限元》讲义第4章弹性薄板弯曲问 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 的有限元法板和壳是指厚度比其他尺寸要小得多的平面或曲面构件。由于它的这种几何特点，三维单元并不适合用来 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 它们的变形。因为三维单元在三个方向的尺寸应尽量接近，否则求解精度由于“剪切自锁”(shearlocking)或系统矩阵病态而大大降低，甚至得到错误的结果。所以必须采用很细密的网格来适应板和壳的几何特征，但是这将导致有限元模型的自由度疯狂地增长。仿照根据梁理论建立梁单元的思路，自然想到根据板理论建立板单元。这里讨论两种板理论，一是薄板理论，也被称为Kirchhoff板理论，它忽略了板的横向剪切变形；另一种是Mindlin板理论，它考虑了板的横向剪切变形的影响，适合于板的厚跨比较大的情形。后者也常被称为Reissner板理论[8]或中厚板理论。根据这两种理论可以建立不同的板单元。薄板弯曲问题在理论上和应用上都具有重要意义,并有专门著作加以论述（如杨耀乾《平板理论》）。象其它弹性力学问题一样，用微分方程、差分法等经典方法所能求解的薄板问题很有限，一般只能解决等厚、小孔口、支承情况较简单的单跨板。故工程设计中以往多采用简化、近似、图 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 等方法来解决板的设计问题。在板的分析中，常取板的中面为xoy平面(如图)。平板结构按其厚度t与短边a的比值大小而分为:厚板（Thickplate）和薄板（Thinplate)两种。当ta<<1时称为薄板平板上所承受的荷载通常有两种:1.面内拉压荷载。由面内拉压刚度承担,属平面应力问题。2.垂直于板的法向荷载,弯扭变形为主,具有梁的受力特征,即常说的弯曲问题。平板在垂直于板面的荷载作用下产生挠度W。当最大挠度w远小于t时,称为小挠度问题(or刚性板)(stiffnessplate)当最大挠度w与t相差不大时,称为大挠度问题(or柔性板)(flexureplate)(工程定义:15£wt为刚性板；15£wt£5为柔性板；>5为绝对柔性板。)4.1wt基本理论一、基本假定１、略去垂直于中面的法向应力。（sz=0），即以中面上沿Z方向的挠度W代表板的挠度)２、变形前垂直中面的任意直线，变形后仍保持为垂直中面的直线。（─法向假定tzx=0,tzy=0）３、板弯曲时，中面不产生应力。(─中面中性层假定)上述假定常称为薄板小挠度问题假定（or柯克霍夫假定）。符合上述假定的平板即为刚性板。二、基本方法1¾将¾U.¾代¾®¶x2ïìexüï{e}＝íeyý=í-z2ýïeïïïîxyþï¶wï¶x2¶y2¶x¶y《有限元》讲义以上述假定为基础，板分析中常用挠度w作为基本未知量,下面介绍以w为基本未知量所导出的有关方程。１、几何方程（应变─挠度关系）①弹性曲面沿x,y方向的倾角从中面取出一微小矩形ABCD，如图所示，设其边长为dx,dy，变形后弯曲成曲面A'B'C'D'设A点挠度w,则沿x方向倾角(绕y轴)¶w¶xqy=¶w¶x(B’点绕度w+dx)沿y方向倾角(绕x轴)¶w¶yqx=¶w¶y(D’点绕度w+dy)②沿x,y方向位移作平行于xoz平面,设中面上点A到A1的距离为Z，变形后,A点有挠度W,同时发生弯曲，曲面沿x方向的倾角为¶w¶x,根据法线假定,则A1点沿x方向的位移:¶xu=-z¶w(负号为方向与x相反)同理取yoz平面得：¶yv=-z¶w(4-1-1)③Z平面的应变分量和曲、扭率基本假定，由于sz=tzx=tzy=0,故板内任意点的应变与平面问题相同:¶v¶x+¶u¶x¶v¶y¶u¶yex=ey=exy=V入îþì¶2wüï-zïïïï¶2wï¶y2ï-2z¶x¶yï(4-1-2)¶2w¶2w¶2w,,为曲面在X,Y方向的曲、此为Z平面的应变─挠度度几何方程。上式中的扭率，记为：2¶2wìcxüï-{c}=ïícyïý=ïí-¶w¶y2ïýïcïïîxyþï-2¶2wïê1êm1-mú《有限元》讲义ìü2¶x2ïïî¶x¶yþ(4-1-3)所以,{e}=z{c}î２、物理方程（应力─挠度关系）由于忽略σz对变形的影响,因此z平面的应力─应变关系具有与平面问题相同的形式:ìEEïEï2(1-m)将(4-1-2)代入得:2ç2ìsxüï22îxyþïïïïîþ或简写为：{s}=z[D0]{x}(4-1-4)式中弹性矩阵:ù0ú0úú2ûm10éE1-m2êê0ëD0=３、内力方程（内力─挠度关系）从板内取微元体tdxdy,由其上正应力sx,sy和剪应力txy,可在截面上合成合力矩:Mx(y0z面上由sx产生的绕Y轴弯矩)3{F}=t[D0]{c}[s]=123z{F}w(x,y)=ååAmnsinq(x,y)=ååqmnsinçç2+2÷÷í绕x轴转角qx(上节为沿y方向倾角《有限元》讲义My(x0z面上由sy产生的绕X轴弯矩)扭矩:假定Mxy（由剪应力产生，如图）Mx,My,Mxy分别表示单位宽度上的内力矩。如是，内力矩阵：üìMxüï¶x2ïtt3212¶y22ïïî简写成312(4-1-5)比较(4-1-4)和(4-1-5)可得用内力矩表示的平板应力:t由此可见，平板上、下表面处的应力最大：{s}z=±t=±6{F}22以上是薄板弯曲问题中的基本公式,从中可见其挠度W是弯曲问题中的基本未知函数。且由于忽略了z方向的变化,因此它只是x，y的函数:w=w(x,y)。若w已知，则位移，内力、应力均可按上述相应公式求出。在经典解析法中，W(x,y)常设为三角级数形式。例如，四边简支矩形板的W(x,y)设为:(纳维尔解)sin¥¥m=1n=1npybmpxa式中Amn为待定系数。假定荷载sin¥¥m=1n=1npybmpxa则可得位移函数:qmnnpybmpxasinsin1Dpw(x,y)=2æm2n2öèabø4åå4.2有限元分析方法一、矩形单元的典型形式将图示矩形薄板沿x,y方向划分成若干小矩形(常取等分)从中取出一小矩形(单元)，共有四个结点，每个结点有三个位移分量：挠度w,绕x、y轴转角))îì挠度wïï绕y轴转角qy(上节为沿x方向倾角即结点i的位移4ìwiüïìüïïï¶w{di}=íqxiý=íïqïïîyiþï-¶w{d}e=[d1Ld4]=[wqx1qy1Lwqx4qy4]LF4]=[f1Mx1My1Lf4Mx4My4]在上式中,前10项取到了三次项的全部,最后两项则是从五个四次项(x第二项qy=qx=¶wqy=-a2ïý表示刚体转动ïþ《有限元》讲义þw¶yýïî¶xïiïï(i=1,K4)同理,相应的结点力ìfiü竖向力ïïîyiþxT符号重新定义是为了有限元表示的方便，由此得单元结位移向量T节点力TT{F}e=[F1二、位移模式（函数）１、位移模式的选取插值多项式取为：w(x,y)=a1+a2x+a3y+a4x2+a5xy+a6y2+a7x3+a8x2y+a9xy2+a10y3+a11x3y+a12xy3(4-2-1)4x3yx2y2xy3y4)中选用了两个。没2选xy2是因为它没有多一项与其配对，没选x4,y4它们在边界上结出的挠度函数是四次的，比x3y和xy3要高一次,较之更难满足边界的协调和条件。２、位移模式的检验（三个基本要求：刚体位移，常应变，尽可能的边界协调）①前三项含单元的刚体位移状态：第一项a1与坐标x,y无关,表示z方向的挠度是─常量,刚体移动第三项üïqx=a3ï-¶w¶x¶y②二次项代表均匀变形状态:曲率¶2w¶x2=-2a4，¶2w¶y2=-2a6，=-2a5¶2w¶x¶y③能保证相邻单元在公共边界上挠度的连续性。5{w}=[N]{d}[N]=ê¶N1/¶y¶Ny4/¶yúNxi=-bhi(1+xix)(1+hih)1-h2Nyi=axi(1+xix)(1+hih)(1-x2)6《有限元》讲义④不能保证相邻单元在公共边界上法线转角的连续性。y=-b=常量，代入位移模式4-2-1，可知边界上的挠度W是x的以单元１～２边界为例，在此边界上三次函数，合并整理后可得：w1-2=c1+c2x+c3x2+c4x3两个端点共有4个边界条件,(结点1,2的挠度W1,W2,和转角qy1,qy2。利用他们可唯一确定四个常数C1～C4。因为相邻单元在结点1,2的W,θy对应相同，则两个单元依据四个条件得到的C1～C4亦相同，即两单元在边界具有同一挠度函数W。⑤法线转角仍以1-2边界为例，将y=-b代入后，此时qx=d1+d2x+d3x2+d4x3但对θx来讲，1,2结点只能提供2个已知条件，不能完全确定上式，故边界的法线转角不能保证连续性。因此，这种单元是非协调元，但可以验证这种非协调远是能通过分片试验的。（即当单元划分不断缩小时，计算结果仍能收敛于精确解。）三、形函数和形函数矩阵。¶w¶y分别将单元结点1,2,3,4的坐标值代入(4-2-1)，并事先求出qx=¶w¶x，qy=-，便可得到各结点的位移值。一共可得12个关于ai的方程组，联立求解可得：e(4-2-2)形函数矩阵:ëNy4ùú¶Ny4/¶xúûéN1êê¶N1/¶xNx4LLNx1LLNy1LN4LLLLLL式中形函数:18(1+xix)(1+hih)(2+xix+hih-x2-h2)Ni=)1818(4-2-3)(i=1234)在上面的推导中,我们仍然选用了局部坐标(无因次坐标)。局部坐标与整体坐标的关系为:1a16(y-y0)h=é¶2Nê¶x2[B]=-êê¶N1¶y2ê2¶N1ë¶2Ny4¶Ny4¶2Ny4¶x¶x¶yû7《有限元》讲义V=z四、单元的几何矩阵[B]和内力矩阵[S]1．几何矩阵[B]由前可知{e}=z{c},将（4－2－3）代入(4-2-4)得到几何矩阵：¶x¶y¶xêê222ù2úúú¶y2úúú2¶x21¶2Ny1LL¶2Nx1LLLLL2（4-2-5）或以子块形式表示：[B]=[B1B2B3B4]。式中：2.内力矩阵[S]由基本方程（4-2-5）可得到：{F}=[D][B]{d}e=[S]{d}e（4-2-6）òòt[B][D][B]dxdyT《有限元》讲义[S]称为内力矩阵，把单元的四个结点坐标分别代入4-2-4，求得[B]后，即可获得[S]，各节点内力矩阵[S]的显式：五、单元刚度矩阵ab-a-b由一般公式得:[K]=。将几何矩阵[B]和弹性矩阵[D]的表达式代入，积分可得薄板弯曲问题矩形单元的单元刚度矩阵的显示：8{}()òò--=abTdxdyyxqyxNR,)],([ïïïïaNïïbïï3òòí=qR}{ýí=ýqabdxdy1ïïN3ïïbïïxN3ïïïïN3ïï-þyN4ïïa3《有限元》讲义9六、荷载等效变换由荷载等效变换的一般公式可得ab1．法向均布荷载q代入上述公式得：ï-ïïïî-a-bî3þï-ïïïï3ïïNïï3ïì1üïbïïNy3ïïaïï4ïï1ïïNx4ïïbïabïNy2ïïx2ïïaïïïï1ïïN2ïìN1üï3ïïx1ïï-ïïNy1ïï3ïîïïïï《有限元》讲义102．单元中心点受法向集中力P代入上述公式可得：七、位移边界条件对称、固定边和简支边上支点的已知位移条件如下：对称轴:法线转角=0固定边:挠度=0(或已知值)边线转角=0(或已知值)法线转角=0(或已知值)简支边:挠度=0(或已知值)边线转角=0(或已知值)自由边上节点的挠度、边线和法线转角均为特定参数，同内部节点一样。与板铰接的固定立柱，其节点挠度Ｗ=0，也可以是已知值。八、计算例题例题1:计算图示四边固定方板方板的边长为l，厚度为t，弹性模型量为E，波松比μ=0.3，全板承受均布法向荷载ｑ，求薄板中的挠度和内力。单元划分：为了说明解题方法，采用最简单的网络2×2，即把方板分成四个矩形单元。由于对称性，只需计算一个单元，例如，计算图中有阴影的单元，单元的节点编号为１，２，３，４。此时，单元的a,b是l4a=b=计算节点荷载:由前面的均布荷载计算公式得：ql219212-l-l]T[12l-l12ll12-ll{R}=边界条件：边界23和34为固定边，因此节点2,3,4的挠度、边线和法线转角均为零。边界12和14为对称轴，因此θx1=0、θy1=0。于是，在4个节点和12个位移分量中,只有一个待求的未知量w1。结构的代数方程组：这是一个单元的计算题目，单元刚度矩阵在此处即为总刚度矩阵。引入支承条件后，在总刚度矩阵中只取第一行、列元素，在方程组右端项中只保留第一个元素。于是结构的代数方程为：ql2168D015l28D015l2k1w1=(81-6m)w1={P}=2üb-aïï2ïbïaïý2ï-bïïa2ïï-bï-aïþìïïïïïïïPïí8ïïïïïïïï同此解出w1=0.00148ql。其中《有限元》讲义4D0=0.09158Et3D0=Et312(1-m2)内力:利用式(4-2-6)可求得方板中点力矩为：由表看出，网格越密，计算结果越接近于精确 答案 八年级地理上册填图题岩土工程勘察试题省略号的作用及举例应急救援安全知识车间5s试题及答案 。还可看出，位移的精度一般比内力的精度高，这是因为在位移法中，位移是由基本方程直接求出的，而内力则是根据位移间接求出的。11《有限元》讲义4.3薄板有限元程序设计一、总框图根据弯曲板有限元分析方法的解题过程，可写出其总框图如下：┌───────┐│输入原始数据││orCAI│└───┳───┘↓┌──────┐┌──┤算等效结点力│┌───┻───┐││形成荷载列阵├←┘└──────┘┌─────┐│├←───┤形成单元│└───┳───┘┌─┤定位向量├─┐└─────┘↓┌───┻───┐│││││形成总刚├←─┘┌─────┐││├←───┤单刚││└───┳───┘└─────┘│↓┌───┻────┐│解方程输出位移│└───┳────┘↓┌───┻───┐││││┌──────┐││几何矩阵[B]│││├←───┤弹性矩阵[D]│││计算单元内力等│└──────┘││├←───────────┘└───┳───┘↓┌──┻──┐│结束│└─────┘下面结合程序对框图中的 内容 财务内部控制制度的内容财务内部控制制度的内容人员招聘与配置的内容项目成本控制的内容消防安全演练内容 加以说明。二、子框图１、单元坐标结点编号及单刚形式。为了取挠度向下为正，又能与前述坐标系统统一，特将前述坐标前翻180°(如图)为了能适用板的弹型性分析,程序采用了应力元和弯曲元的组合形式，即每个结点考虑5个位移分量:U,V,W,θx,θy,前2个为平面应力问题的未知量,后3个为弯曲板的结点未知量。当只作弹性分析时，平面应力元和弯曲元是非藕连的，即单刚的两个副块垣为0,单刚的形式为：12《有限元》讲义w1θx1θy4…w4θx4θy4u1v1…u4v4┌┐│平面应力元│0│[K]e=│（8×8）││├──────┼──────────────┤││弯曲元││0│（12×12）│└┘程序中单刚数组为DK(20,20),子程序:SubroutineDG(A,B,E,T,U)为其形成单刚的子程序。２、自动形成单元编号信息（单元信息数组：[IB]）。３、结点定位向量。４、形成荷载列向量。（a.结点力；b.非结点力(只考虑均布力)）５、总刚，SubroutineZG(M,N,LD,A,B,E,T,U)６、解方程。FJZG(),HUD()７、算单元力。SubroutineDYL()８、算等效结点力。９、弹性矩阵[D]。１０、几何矩阵[B]。三、输入数据说明１、总信息。共11个(见程序)２、结点约束信息数组[JB]JB(I,1)──i结点的结点号JB(I,2)──i结点的约束分量号(1～5)结点约束信息应根据支承条件或对称条件决定，如算例中所给出的四边简支方板，承受满布均布力，此时可只取板的1/4作为分析对象，如下图只取右上角1/4板，采用6×6网络，则每个单元的边长为1米(A=0.5,B=0.5)。设结点编号如图示:在y=0的边界上(1-7结点):挠度w=0(第3个分量)绕y轴转角θy=0(第5个分量)同理,在平形于y轴的x=6m边界上:w=0,θx=0(3,4分量)在对称轴x=0边界上u=0,θy=0在对称轴y=6边界上v=0,θx=0中点(43点)除W外,其余均=0同时，在简支边上,也可设u和v均为0这样便共有77个约束。另外，也可通过改变约束信息来改变计算简图，如同样网格数的1/2、1/4板，固支边界板等。３、结点荷载信息pp(I,1)荷载值pp(I,2)结点号·位移分量号(如竖向力的位移分量为3)４、非结点荷载信息13《有限元》讲义14（程序中均为考虑了满布均布力）PF(I,1)荷载值PF(I,2)荷载作用单元号四、输出信息1.结点位移；2.单元力（王勖成P144,5章2节应力计算结果的处理与改善）。