Contents     Functions     Index     Previous Next     PDF

fsc



Purpose

Fluid structure/coupling with non-linear follower pressure support.



Description

Elasto-acoustic coupling is used to model structures containing a compressible, non-weighing fluid, with or without a free surface.

The FE formulation for this type of problem can be written as [41]
  s 2 [
M 0
CT Kp
] {
q
p
}+ [
K(s) -C
0 F
] {
q
p
} = {
Fext
0
}     (8.1)

with q the displacements of the structure, p the pressure variations in the fluid and Fext the external load applied to the structure, where
ó
õ
 


WS
ij(u)ij( u) dxÞ qT K q
ó
õ
 


WS
S u. u dxÞ qT M q
1
F
ó
õ
 


WF
Ñ p Ñ p dxÞ pT F p
1
F c2
ó
õ
 


WF
p p dxÞ pT Kp p
ó
õ
 


S
p u.n dxÞ qT C p
    (8.2)


Follower force

One uses the identity
n   dS =
x
r
Ù
x
s
  drds,     (8.3)
where (r,s) designate local coordinates of the face (assumed such that the normal is outgoing). Work of the pressure is thus:
Wp = - ó
õ
 


r,s
P   (
x
r
Ù
x
s
) ·v drds.     (8.4)
On thus must add the non-linear stiffness term:
-d Wp = ó
õ
 


r,s
P   (
du
r
Ù
x
s
+
x
r
Ù
du
s
) ·v drds.     (8.5)

Using x/ r={x1,r   x2,r   x3,r}T (idem for s), and also
[Axr]= æ
ç
ç
è
0 -x,r3 x,r2
x,r3 0 -x,r1
-x,r2 x,r1 0
ö
÷
÷
ø		 ,   [Axs]= æ
ç
ç
è
0 -x,s3 x,s2
x,s3 0 -x,s1
-x,s2 x,s1 0
ö
÷
÷
ø		 ,
this results in
     (
dx
r
Ù
x
s
+
x
r
Ù
dx
s
) ·v =
 
    { qik}T {Nk} ( Axrij {Nl,s}T - Axsij {Nl,r}T ) {d qj}.     (8.6)

Tests : fsc3 testsimple and fsc3 test.

In the RivlinCube test , the pressure on each free face is given by
P1 = -
1+1
(1+2)(1+3)
S11
on face (x1=l1)
P2 = -
1+2
(1+1)(1+3)
S22
on face (x2=l2)
P3 = -
1+3
(1+1)(1+2)
S33
on face (x3=l3).


Non-conform

SDT supports non conforming element for fluid/structure coupling Terms corresponding to the structure are computed using the classical elements of the SDT, and terms corresponding to the fluid are computed using the fluid elements (see flui4).

The coupling term C is computed using fluid/structure coupling elements (fsc elements).

Only one integration point on each element (the centre of gravity) is used to evaluate C.

When structural and fluid meshes do not match at boundaries, pairing of elements needs to be done. The pairing procedure can be described for each element. For each fluid element Fi, one takes the center of gravity Gf,i (see figure), and searches the solid element Si which is in front of the center of gravity, in the direction of the normal to the fluid element Fi. The projection of Gf,i on the solid element, Pi, belongs to Si, and one computes the reference coordinate r and s of Pi in Si (if Si is a quad4, -1<r<1 and -1<s<1). Thus one knows the weights that have to be associated to each node of Si. The coupling term will associate the DOFs of Fi to the DOFs of Si, with the corresponding weights.


See also

flui4, m_elastic©1991-2007 by SDTools
Previous Up Next