降维变分方法.pdfVIP

第25卷增刊 2008年12月 计算力学学报 Chinese Journal ofCompt啦fional Mechanics V01．25，Sup． December 2008 文章编号：1007—4708(2008)Sup-0001—05 降维变分方法 唐小行 (中国矿业大学机电学院，徐州221008) 摘要：依据能量原理建立降维变分法，该方法将物体离散为条，在一定约束下，一条独立求解能量极值，可 有效简化数值计算，本文研究分条离散方法，极值约束，边界条件处理并给出证明．更重要的，该方法揭示了 新的物体形变运动规律．在此也做了分析． 关键词：降维变分法；徼条；组合变形；组合位移：塑性滑移 中图分类号：029：0344 文献标识码：A 1 引 言 有限元方法广泛地用于求解大型科学与工程 问题，取得了巨大的成就。现代科学技术面临求 解高维、非线性、奇异解、复杂结构等困难问题， 要求计算规模越来越大，有限元方法因需在整个 区域上进行离散化，网格剖分复杂、繁琐，未知 量存储多，计算量大。如流体运动数值模拟，目 前只能就低雷诺数条件下进行湍流直接数值模 拟，因为网格最小尺度难以达到湍流最小尺度Il J。 在求解金属冲压变形，高速撞击，裂纹动态扩展 等涉及大变形时，有限元网格可能会产生严重扭 曲，不仅需要网格重构，而且严重影响解的精度【2】。 降维变分计算将物体离散为条，每条在一定约束 下独立求解能量极值，可有效简化数值计算。 2位移降维变分法 2．1必要性证明 2．1．1 内部一条的独立计算 图l中立方体材料均匀，坐标原点与立方体 形心重合，而，jc2，毛方向边长分别为2“272， 2厶。边界位移条件给定，将其分为平行毛轴微条 (见图)，沿毛向排列为肌，而向为行的微条记为 仰，”)，取位移变分万虬，依据最小势能原理，在 满足位移边界条件的各组位移中，实际存在的位 移应使物体形变能取最小值。设”，满足毛=一“ 收稿日期：2008-08-08． 作者简介：唐小行(1957一)，男．副教授 (E--lail：dataagl i@126．COg) 屯=-f2，毛=±‘边界位移条件 ％=瓦 瓦给定，形变满足几何方程 毛=jI(坼．J+“川) (2) 圈I将立方体分为平行毛轴微条 Fig．1 Cube divided into a sd of infinitesimal蛳parallel to 与axis 为分条独立计算，玑对分条截面的法向导数 变分万％．1，艿％应有约束：使扰，所含变量组合，在 本条形成 8uj=0 (3) 如此，对分出的内部任一微条(历，聆=3，4，5?) 8A=0 以下周向约束据边界条件而定，例如，假定 外力作用与变形关于毛=0，x2=0，毛=0三面 对称，有约束 fUj．。咖=cf。 I，％．：咖=cj： (4，5) 对万％：不受(3)式约束，(4)、(5)式改为： Ldi_uflv=0 (6) 式中：A为微条形变势能，矿为微条所占区域。 2 计算 力 学 学 报 第25卷 2．1．2起始端微条共取极值分析 以五=-，l端为例，因取应变变分，左边界位 移条件应在(1，刀)，(2，刀)(n=l，2，3，?)每两 条形变能极值的分部积分中对两条积分 。L(q』6uj)，1 dv中各项的约束，因此要求端部每条 与右邻条共取极值。毛=‘端同理，xl=‘， 毛=‘相交处头4条(1，1)-(1，2)，(2，1)，(2， 2)也共求极值。 2．1．3周向约束 同时可将物体分为平行五轴微条的类似计算 (参见以下计算公式)，同理而=‘，毛=．，3相 交处头4条为最大的共求极值组，此4条左右端 边界位移同时取值，可确定一个线性组合。假定 外力作用与变形关于xl=0，x2=0，暑=0三面 对称，可将头4条左右端边界条件写成一线性组 合(4条积分)：L珥，。av=c,i，头4条以上共取极 值的不向每两条应有相同组合，这对平行墨轴微 条计算的头4条，(研，1)，(m，2)(历=3，4，5，?) (以下未加说明t?／、刀、所’、刀’均依次取3，4， 5，?)每两条分别有 I。，玛1dv=ql (7) 同理，对(1，疗)，(2，丹)每两条及内部各 条有约束(7)。 同时可将物体分为平行矗轴微条计算，与上 同理，岛=坫边界位移条件组合，对(m，1)，(m， 2)每两条，(1。刀)，(2，刀)每两条，头4条， 其余内部各条分别要求 lvu,-2dv=c『2 (8) 依据对(1，刀)，(2，玎)的约束(7)，而=-fl 边界位移条件，左端第二排(历--2)微条的而=故 边界位移条件，对平行矗轴的(2，1)，(2，2)两 条底端两排4微元，再往上每排两微元应满足 L坼．2av=C,2 (9) 对(所，1)，(肌，2)每两条她．2，占％的约 束：据约束(7)，