Contents     Functions         Previous Next     PDF Index

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 [43]

  s 2 [
M0
CTKp 
] {
q
p
}+ [
K(s)C
0F 
] {
q
p
} = {
Fext
0
}     (8.2)

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

 


ΩS
σij(uij(δ udx⇒δ qT K q
 


ΩS
ρS u.δ u dx⇒δ qT M q
1
ρF
 


ΩF
∇ p ∇ δ p dx⇒δ pT F p
1
ρF c2
 


ΩF
p δ p dx⇒δ pT Kp p
 


Σ
p δ u.n dx⇒δ qT C p
    (8.3)

To assemble fluid/structure coupling matrix you should declare a set of surface elements (any topology) with property p_solid('dbval 1 fsc'). The C matrix (solid forces induced by pressure field) is assembled with the stiffness (matrix type 1), while the CT matrix (fluid pressure due to normal velocity of solid) is assembled with the mass (matrix type 2).

Some formulations, consider a surface impedance proportional to the pressure. This matrix can be computed by defining a group of surface elements with an acoustic material (see m_elastic 2) and a standard surface integration rule (p_solid('dbval 1 d2 -3')). This results in a mass given by

δ pT Kp p=
1
ρF c2
 


ΩF
δ p p dx     (8.4)

Follower force

One uses the identity

n   dS = 
∂ x
∂ r
 ∧ 
∂ x
∂ s
   drds,     (8.5)

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
 Π    
∂ x
∂ r
 ∧ 
∂ x
∂ s

·δv drds.     (8.6)

On thus must add the non-linear stiffness term:

dδ Wp = 
 


r,s
 Π   
∂ du
∂ r
 ∧ 
∂ x
∂ s
∂ x
∂ r
 ∧ 
∂ du
∂ s

·δv drds.     (8.7)

Using ∂ x/∂ r={x1,r   x2,r   x3,r}T (idem for s), and also

[Axr]= 


0x,r3x,r2
x,r30x,r1
x,r2x,r10



,   [Axs]= 


0x,s3x,s2
x,s30x,s1
x,s2x,s10



,

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.8)

Tests : fsc3 testsimple and fsc3 test.

In the RivlinCube test , the pressure on each free face is given by

Π1 = −
1+λ1
(1+λ2)(1+λ3)
 Σ11
onface(x1=l1)
Π2 = −
1+λ2
(1+λ1)(1+λ3)
 Σ22
onface(x2=l2)
Π3 = −
1+λ3
(1+λ1)(1+λ2)
 Σ33
onface(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-2012 by SDTools
Previous Up Next