Abstract
在有限單元分析中, 數值計算之準確度及可靠度對單元形狀的扭曲相當敏感. 因此, 扭曲參數的定義以及單元敏感度的分析, 實為一重要研究課題.由微分幾何中的測地線理論, 已成功地導出平面4節點, 平面8節點和8節點六面體單元等參映射的反函數關係, 在這種反函數關係裡假設了測地線座標無限級數, 此座標可證明為由等參映射時在原點的Jacobian值所決定而得. 若將測地線座標以等參座標的多項式表示, 則多項式中的係數可定義成元素的扭曲參數, 而這些扭曲參數能完全地描述等參映射的反函數關係及其Jacobian值.本文基於上述之推廣, 對於6節點三角形和10節點四面體選擇其形心位置為面積座標和體積座標之原點. 亦可導出其等參映射的反函數關係由此, 則6節點三角形可得6個扭曲參數; 10節點四面體則有18個扭曲參數.此扭曲參數可做為單元形狀扭曲之測度.本文定義扭曲參數的方法具有數學理論上之一致性與一般性, 並可應用於其它有限單元之相關研究.In analyzing the finite element, The accurance and realibity ofnumerical calculation are rather sensitive to the distortionof element shape. so, defining distortion parameter andanalyzing sensitivity of element are an important researchcourse. By using the theory of geodesics in differentialgeometry, inverse relations of the mapping for 4-node elements,8-node quadrilaterals with curved boundaries and 8-nodehexahedra can be successfully derived and expressed in termsof the element coordinates defined at the local origin. Suchinverse relations assume the form of infinite power series inthe element geodesic coordinates, which are shown to be theskew Cartesion coordinates determined by the Jacobian of themapping evaluated at the origin By expressing the geodesiccoordinates in turn in terms of the isoparametric coordinates,the coefficients in the resulted polynomials are suggested tobe the distortion parameters of the element. These distortionparameters can be used to completely describe the inverserelations and the determinant of the of the mapping. In thisstudy, for quadratic triangular element in two dimensions andquadratic tetrahedral element in three dimensions, choose theircentroids as the origins of the area and volume coordinates.Therefore, we can also derived the inverse relations ofmapping. For quadratic triangular element, there are 6parameters; For quadratic tetrahedral element there are 18parameters. From these distortion parameters, they can be themeasures of the distortion of element shape. These methods ofdefining the distortion measures and deriving the inverserelations of the mapping are completely general and can beapplied to any other two- or three- isoparametric elements.