WOLFRAM

Wolfram言語 ＆ システムドキュメントセンター
詳細
    • Details and Options Details and Options
    例題  
    例  
    関連項目
    NDSolve`FEM`

    ElementMesh::femimq

    詳細

    • このメッセージは，有効ではないメッシュ要素が指定された場合に生成される．
    • Off[message]はメッセージを非表示にし，On[message]はメッセージを表示させる．例えば，Off[ElementMesh::femimq]のように使う．

    例題

    例  (2)

    このメッセージは，指定されたメッシュ要素が自己交差する場合，誤った順序を持つ場合，あるいは領域を持たない場合に生成される：

    Wolfram Language code: Needs["NDSolve`FEM`"]
    Wolfram Language code: ToElementMesh["Coordinates" -> {{0, 0}, {1, 0}, {1, 1}, {0, 1}}, "MeshElements" -> {QuadElement[{{3, 2, 1, 4}}]}]

    これを修正するには，正しい順序を指定する必要がある：

    Wolfram Language code: ToElementMesh["Coordinates" -> {{0, 0}, {1, 0}, {1, 1}, {0, 1}}, "MeshElements" -> {QuadElement[{{1, 2, 3, 4}}]}]

    このメッセージは，補間の際に生成されることがある：

    Wolfram Language code: data = {{2.62713, 145.509, 1.}, {2.62713, 155.609, 2.}, {4.50713, 154.346, 2.}, {6.38713, 153.084, 2.}};
    Wolfram Language code: Interpolation[data, InterpolationOrder -> 1];

    この問題は，数値的に非厳密であるために，メッシュの作成中に，領域を持たない三角形が生成されたために起こるものである：

    Wolfram Language code: m = ToElementMesh[data[[All, 1 ;; 2]]]

    3つの三角形が生成されたが，その1つの質が悪い：

    Wolfram Language code: m["Quality"]

    生成されたメッシュを可視化する：

    Wolfram Language code: m["Wireframe"]

    質の閾値を下回る要素を求めて可視化する：

    Wolfram Language code: pos = Position[m["Quality"], _ ? (# <= 0.1&)]; Show[ m["Wireframe"], m["Wireframe"[pos, "MeshElement" -> "MeshElements", "ElementMeshDirective" -> Directive[EdgeForm[Red], FaceForm[]]]]]

    入力データに少し摂動を与えることにより，この問題を修正することができる：

    Wolfram Language code: data2 = data; data2[[All, 1 ;; 2]] += RandomReal[10 ^ -12 * {-1, 1}, Dimensions[data2[[All, 1 ;; 2]]]];
    Wolfram Language code: if = Interpolation[data2, InterpolationOrder -> 1];

    質をチェックする：

    Wolfram Language code: if["ElementMesh"]["Quality"]