有限元分析在快速设计中的应用.docxVIP

有限元分析在快速设计中的应用 1 弹性力学有限单元法 在计算鞋底和鞋顶的设计中，通常需要处理一些商标的艺术变形，并对一些图案进行特殊的变形。一般的做法是通过控制文字串或图形的矩形包围框的变形, 然后根据合理的假定, 推算出包围框内文字串或图形的变形。虽然CorelDRAW等优秀的设计绘图软件有提供文字串的某些艺术变形功能, 但其算法未见报道。我们在开发鞋底和鞋帮款式快速设计系统过程中, 尝试把弹性力学有限单元法的形函数概念, 引入到图形和文字的特殊变形处理中来, 取得良好的效果, 该功能已经加入到我们开发的制鞋CAD系统中。通过计算机实现和对一些特殊变形的限制处理, 表明这种方法具有对各种变形算法统一、操作方便、令二维图形有立体感等特点。 2 有限元算法介绍 有限元法广泛应用于解决弹性力学、固体力学、流体力学、电磁学和传热学等方面的问题。然而, 在图形图像处理方面的应用还少见报道。根据单元选取的不同可得到定义不同、公式不同、适用范围不同的有限元算法。本文为解决文字和图形的变形, 采用平面应力问题中12结点的矩形单元。有限元的单元分析中有一个较为重要的步骤:即由结点位移求内部任一点的位移。要假设单元的位移分布形式, 使其满足结点和单元边界的位移连续条件, 而形函数就是这种位移变化的形状函数。比如, 点i的形函数表示当单元结点i发生单元位移, 单元其它结点位移分量为零时, 在单元内部产生相应的位移的分布形状。利用这种形函数的定义, 我们就可以通过改变包围框上结点的位置, 来改变框内文字或图形上结点的位置, 并保持框内的结点始终位于框内, 得到有效变形的结果。 3 变形算法用于计算函数 3.1 边界框架的处理 3.1.1 y右上上,y右下x,y 当选择一些图形或文字对象后, 生成一个矩形窗口, 落在窗口内的图形和文字的结点为欲参加变形的结点, 该窗口称为这些结点的包围框, 如图1所示。其左上角的结点坐标为 (x左上, y左上) , 右下角的结点坐标为 (x右下, y右下) 。 由于算法的要求, 采用无因次的自然坐标系, 必须把原坐标 (x, y) 根据公式 (1) 进行坐标规格化变换。公式如下: ξ=(x-x0)/a,η=(y-y0)/b(1)其中∶x0=(x右下+x左上)/2,y0=(y右下+y左上)/2a=(x右下-x左上)/2,b=(y右下-y左上)/2(2)ξ=(x?x0)/a,η=(y?y0)/b(1)其中∶x0=(x右下+x左上)/2,y0=(y右下+y左上)/2a=(x右下?x左上)/2,b=(y右下?y左上)/2(2) 3.1.2 becier曲线拟合 将矩形包围框的四条边分别看成特殊的三次Bezier曲线, 如图2所示。 Bezier曲线是由一组折线集 (也称之为特征多边形) 来定义的, 其中第一点和最后一点在曲线上, 第一条和最后一条折线表示曲线在起点、终点处的切线方向, 曲线形态趋于折线集的形态。Bezier曲线广泛应用于艺术设计中, 而三次参数曲线是得到曲线之间光滑连接的最低次曲线。用三次Bezier曲线来表示包围框, 可以较好地拟合边界的光滑程度, 其参数方程为: γ(t)=(1-t)3γ0+3t(1-t)2γ1+3t2(1-t)γ2+t3γ3(3)γ(t)=(1?t)3γ0+3t(1?t)2γ1+3t2(1?t)γ2+t3γ3(3) 包围框的四边在变形前为直线段, 当把它们分别看成特殊的三次Bezier曲线时, 取线段的三等分点作为曲线的特征点, 有: {γ1=(2γ0+γ3)/3γ2=(γ0+2γ3)/3(4){γ1=(2γ0+γ3)/3γ2=(γ0+2γ3)/3(4) 这样一来, 四条边共有12个控制点, 在图1中用圆点来表示控制点, 处于四个角上的控制点是相邻两边共用的。移动这些控制点 (总共有24个自由度) , 可使包围框变成所需要的形状。现在问题归结为已知这12个控制点变形后的位置, 如何求出包围框内结点 (x, y) 的新位置 (x′, y′) 。 3.2 求框内结点的新位置 首先, 在矩形包围框的边界上设定12个结点, 其位置与变形前控制点的位置重合, 如图3所示, 我们给它们注上编号。 每条边内部的两个结点, 恰好三等分该边。设控制点移动后的新位置分别为P1′, P2′, …, P12′, 令矩形框上12个结点的新位置分别为Q1′, Q2′, …, Q12′。于是, 由 (3) 式可知: Q1′=P1′,Q4′=P4′,Q7′=P7′,Q10′=P10′ (5) 而 类似地有: Q5′=(1-13)3Ρ4′+3×(13)(1-13)2Ρ5′+3×(13)2(1-13)Ρ6′+(13)3Ρ7′Q6′=(1-23)3Ρ4′+3×(23)(1-23)2Ρ5′+3×(23)2(1-23