第7章流体流动问题的有限元法(13)

T T T T [N ] [N ] [N ] [N ] [N ] [N ] + ρv[N ] +µ +µ [K11 ] = ∫∫ 2 ρu[N ] x y x x y y ( e )

dxdy

T T [N ] { } U [N ] { } U T {F1 } = ∫ µ [N ] dx + ∫ µ [N ] [N ][N ] {P} dy y x (e) (e) S S

[K12 ] = ∫∫ ρu[N ] [N ] dxdy, y ( e )T

[C1 ] = ∫∫ [N ] [N ]T dxdy, y ( e )

[K 21 ] = ∫∫ ρv[N ] [N ] dxdy x ( e )T

T T [N ] { } U U [N ] { } T {F2 } = ∫ µ [N ] [N ][N ] {P} dx + ∫ µ [N ] dy y x S (e) S (e)

T T T T [N ] [N ] [N ] [N ] [N ] [N ] [K 22 ] = ∫∫ ρu[N ] +µ + 2 ρv[N ] +µ x x y y x y ( e ) [C 2 ] = ∫∫ [N ] [N ]T dxdy, y ( e )

dxdy