降维变分法.pdfVIP

降维变分法 唐晓行 【中国矿业大学机电工程学院) ，摘要依据最小势能．余能原理，建立了静态降维芟分法，以此为基础， 由哈密尔顿原理推得动态降维变分计算方法．该方法将物体离散为条，在任一时 刻．每条在一定约束下独立求解一维能量极值，从而大大简化了工程力学数值计 算．同时谊方法揭示了新的物体形变、运动规律，本文对此作了分析． 。，0 关键词降维变分法，微条，合变形，舍位移 在固体力学数值计算等领域有限元法以分域多项式插值对结构完全离散化，使之能应 付各种复杂边界问题，有着广泛的适用性及简易性，但该方法又不可避免地带来自由度 多，计算量大的一面．随着有限元方法深入到非线性、高速冲击、流体、断裂、细观力学 等领域，深感其所需单元，内存与计算量十分大。造成工程应用上的困难．降维变分法将 物体分条单独求解能量极值，有效地简化了数值计算．本文研究物体的分条离散方法，对 极值的约束及边界条件处理并给出证明． ． 1静态位移降维变分法 I．I计算方法简介 图1中立方体整体位移边界条件给定，将其分为平行z轴微条沿x方向排别为第妙 向为』的微条记为({『)，取位移变分，对任一微条要求 万爿=0 一组周向约束(依据边界条件而定，见对极值的周向约束一节)，如： C3 (2) 展出《偿d：=c2协OWdz f考出=。』考出=G．f茜a：=G(3) 6个应变 矛H a“ a” 图1 5，2瓦 ～2万+瓦 f4) 上下位移边界条件 “=U V=V∞=曲 其中爿为微条形变势能，c为待定常量，=o、妒=0边界位移条件给定，计算沿枷，方向逐条进 行． 1．2证明方程(1)为最小势能原理的必要条件 在整体边界位移条件给定，满足相容条件下，最小势能原理要求物体形变能取最小 值。除一条外其余位移变分取0，得(1)、(4)、(5)式。 a起始端微条计算 以a-=O端为例，据(4)，’第一排每条右边界位移条件给定，使得本条一些应变变分为其 左截面位移变分，引起左邻条应变变化，因此要求端部每条与左邻条在(4)及上、下、右边 ．-■■■■r 界条件约束下共求极值．．，卸端同理， “头四条”11，12，21，22也有此要求． b对极值的周向约束 t一…同时将物体分为平行J轴微条沿y、：方向逐条计算(参见充分性证明一节)，同理。 轴处头四条是最大的一共求极值组，因此四条左右端部边界位移为最大的边界位移共取值 组．以“为例，据特性：两点确定一线性组合，三点定一面，四条左右位移边界条件可写 c取值其分析相同)，以下前后每两微条有相应组合．若￡脚对)鄙端平行z轴前后共取 极值的两条(头四条除外)要求 JJ昙蛐％同理JJ塞妙a：％jJ罢挑％ 头四条有 Ⅲ罢蛐d=铂。同理『』J詈出州：咱，『『『豢出挑铂： (6) 同理x=O端部左右两两微条(头四条除外) 』『塞她弘， 『『塞蛐铂。．『『芸她硝。 物体分为平行y轴微条计算：以“为例，z,--o端上下共取极值的每两条前后端边界位 移条件写成组合可考oI+P