^{1}

^{*}

^{2}

^{1}

This article proposes a synthesis and contribution at three levels: generation of dynamic equations of shell structures interacting with fluids, reduction of implicit resolution, and cross-applications to aerospace tanks and living systems. The synthesis of the equations is proposed around the four principles of thermodynamics at the level of discrete, structural and digitized systems. The implicit approach envisages an innovative analysis in terms of condensation and digitization, with in particular a perspective towards singular and integral methods. Some illustrations are proposed, in the field of performed research models and also in the fields of educational applications in biodynamics. The proposed bridge links, on one hand, the analytical Lagrange-Feynman’s approach, and on the other hand experimental results obtained in laboratory and numerical experiments obtained with multiphysics software. Finally, the realized models concern conservative and dissipative models for the active and passive control of complex systems, in a unified approach.

Living cells are complex, controlled and dynamic systems. They are governed by the laws of the universe, i.e. physics, and therefore fundamentally by the four principles of thermodynamics. The development, life and death of cells are comparable to their counterparts in structures, especially thin shells interacting with fluids.

The purpose of this paper is to present in parallel, in a synoptic way, the reversible and irreversible phenomena observed at the cellular, structural and organic scales (see [

In terms of observability and controllability, complex systems give rise to a variety of models which—in the end—are all based on the four principles of thermodynamics. This approach has made it possible to model several biological systems, including vestibular elements, the basilar membrane, the organ of Corti, the cerebrospinal system in its traumatological approach (see [

The simplest dynamic modeling consists of declining the first two principles—based on the two state functions temperature and internal energy—according to exchanges between kinetic and potential energies. Basically, this amounts to expressing Newton’s equivalence principle digitally in matrix form, via the Lagrange-Hamilton equations:

m O P → ° ° − f → ≡ O → ⇒ [ M i j ] ( q j ) ° ° + [ K i j ] ( q j ) ≡ ( O ) (Feynman) (1)

where

m: punctual mass (particle);

O: affine Galilean origin (euclidian);

P: point associated with m (trajectory);

°: particle derivative with respect to time;

f → force associated with m;

( q j ) : column vector of n Lagrange displacement parameters ( j = 1 , ⋯ , n ∈ N );

[ M i j ] : n × n matrix of system inertias, associated with the q j parameters ( i , j = 1 , ⋯ , n );

[ K i j ] : n × n matrix of system stiffnesses, associated with the q j parameters.

To be more precise, the central arrow of Equation (1) contains all the dynamic principles, whose synoptic presentation could be summarized as shown hereunder; it can be observed that the structure of all these equations is similar, as shown hereunder.

· First step: Equivalent enthalpy

Equivalence (closed system)

v → ° − g → ≡ 0 →

(Newton) 1st principle

Dynamics (digital system)

〈 L j 〉 ( E l a g ) − Q j d y n ≡ Q j s t a t

(Lagrange) 4th principle,

with 〈 L j 〉 ≡ d d t ( ∂ ∂ q j ° ) − ∂ ∂ q j

· Second step: Kinetic exchanges & exergy

Continuous 1D U ° ° – c N 2 U ″ ≡ 0 with c N 2 ≡ E ρ (Mach-Hooke)

Continuous 1.5 - 2.5D ρ S V ° ° + E I V ′ ′ ′ ′ ≡ 0 ; ρ h W ° ° + D Δ Δ W ≡ 0

(Mach-Hooke)

Discrete 3D h 2 Δ ψ + 8 π 2 m E h a m ψ ≡ 0 (Schrödinger)

Feynman’s equivalence formalism is valid in the linear domain with defined [K] and [M] matrices, but remains applicable in the non-linear domain with living matrices. Thus, the fundamental approach works step by step.

To be complete, we can observe that Feynman’s equivalence is based on the four principles, by truncature at resonance. In fact, the complete equation, shown here in the general formalism of Lagrange and in the linear digital frame, is in correspondance with the principles as proposed hereunder:

〈 L j 〉 ( E k i n ) − 〈 L j 〉 ( E o p t ) ≡ Q j f e y n + Q j s t a t (2)

[ M ] ( q ) ° ° + [ K ] ( q ) ≡ ( Q ) f e y n + ( Q ) s t a t i c l o a d s { P r 2 ; P r 1 } { P r 3 ; P r 4 }

where Q j f e y n ≡ Q j d y n p o t + Q j (dynamic loads & dissipation).

It should be noted at this level that Maxwell-Betti’s reciprocity theorem implies the symmetry of these matrices (constant or not) and the principle of inertia the fact that they are real positive semi-definite matrices. They will even be defined positive ( m i i > 0 , k i i > 0 ) in the case of structures which admit a static reference configuration.

However, eigenmodes can be explained algebraically in the linear domain because the whole system can be resolved on a harmonic basis. Indeed, the diagonalization allows to make an equivalence of the whole dynamics with systems represented by a single Degree of Freedom (DoF).

Thus, it is possible to define an implicit analytical impedance matrix, in which the time has explicitly disappeared:

[ Z i j ] ≡ − ω 2 [ M i j ] + [ K i j ] ; i , j = 1 , ⋯ , n ⇒ [ Z i j ] ( q j ) ≡ ( 0 ) ⇒ − ω i 2 μ i + χ i ≡ 0 (3)

where [ Z i j ] structural harmonic impedance matrix of the system, associated with the q_{j} parameters.

ω ∈ ℝ + ∗ : system pulsation (rd/s with classical Lagrange parameters);

ω i ∈ ℝ + ∗ : eigenpulsation of the system (natural or circular pulsation);

( μ i , χ i ) ∈ ℝ + ∗ 2 : effective modal inertia & stiffness of the system.

This approach is the basis for all classical dynamic models, including biodynamic systems. In particular, it allows tracing a first genome of the systems with linear tools (see [

Of particular note are the important possibilities for thin structures (beams and shells) in vibroacoustic interaction with their fluid environment, which lend themselves very well to simplectic and numerical analytical models, thus providing an efficient modal representation and a powerful tool for first approach theoretical observability and controllability.

It is the fundamental tool for topological design support, initial dynamic diagnosis and Generalized Predictive Control (see [

The advantage of the implicit approach is to include time in the harmonic response; it is also its weak point.

The first two principles of thermodynamics respond perfectly one to another, and it is therefore an excellent classical model. It covers perfectly the classical motion of an electron around a proton, the vibrations of a Lagrange beam, or the classical Maxwell equations in the Minkowski-Einstein Euclidean universe. Constitutively, however, it is euclidian and does not assume the history of the system.

Even if it can be solved explicitly (Newmark’s scheme) this system of representation is thus intrinsically adapted to an implicit resolution, the Kramer-Gauss modal one. However, the latter has been most powerful in the representation by truncated modal bases in cyclic dynamics and by shock spectra in transient dynamics. It can thus be very useful to represent the severity of a dynamic environment in fatigue and structural seismology (see [

It does not mean that the causality does not exist, but that it is only implicit. But the implicit time is not the time of the history of the structure, it is simply the number of cycles, or the implied frequencies, but fundamentally the system is built around 2 base units of length and mass (m and kg in ISO system, see

If we wish to go further, we must introduce the other two principles of thermodynamics, which are inseparable from one another (as are the first two). The following example will be used in the linear framework, but the principle extends to the general case.

We therefore introduce an explicit function of time which is neither inertia nor generalized potential, and which is sometimes called excitation, but which we will call external force and which will be represented canonically—according to Fourier and Kronecker—by an elementary harmonic function.

This induces an explicit response (called particular solution of the complete equation) emphasizing the dynamic amplification factor & resonance:

μ p e ° ° + χ p e ≡ cos ( ω e t ) ⇒ p e ≡ β e cos ( ω e t ) ⇒ α e ≡ χ β e = 1 1 − γ e 2 (4)

where ( μ , χ ) ∈ ℝ + ∗ 2 modal inertia and stiffness;

p e ( t ) ∈ ℝ main parameter in the canonical response to canonical cos ( ω e t ) external dynamic load;

β e ∈ ℝ amplitude of the canonical response;

α e ∈ ℝ dynamic amplification ratio (DAR);

γ e ≡ ω e ω ∈ ℝ + ∗ dimensionless external pulsation (frequency or angular rate)

( ω e ∈ ℝ + ∗ ).

The α_{e}DAR factor introduces causality into the model, but in a particular solution. It raises the concept of resonance, which will require an additional (fourth) term of dissipation. This makes it possible to envisage explicit models incorporating structural damage.

In linear models, the dissipation is the viscosity (solid and fluid) but as this is only very approximate in a solid (and in a non-Newtonian fluid) we will rather introduce an equivalent dimensionless viscous damping (integrating all the dissipations—linear or not—around the resonance) which truncates the amplitudes (or operates a frequency shift) in the vicinity of the resonances.

The vicinity of the resonances is moreover precisely defined by the bandwidth:

Frequency amplitude-resonance bandwidth

{ γ e ∈ [ 1 − ζ ; 1 + ζ ] ⊂ ℝ α e max = ± 1 2 ζ (5)

where ζ dimensionless dissipation; for structures ζ ≅ 0.01 ;

α e max dynamic quality factor or resonance surtension ratio.

The ζ factor from which α e max is derived integrates the nonconformities of the structure; it should be noted that it depends—in non-linear models—on absolute frequency values, and must therefore be associated with each resonance. In an n-mode model, this set is therefore represented by then ζ i values ( i = 1 , ⋯ , n ) which define a truncation or frequency slip envelope that contains non-conformity or system damage information.

It should be noted that this entropic representation had been proposed by Feynman in the linear framework, indicating that at resonance the terms viscous and excitation cancel each other out, leaving the two Newtonian terms free. In other words, at resonance, the external force—by definition resonance—compensates the dissipation, and maintains the perpetual motion of the perfect model. This observation—extremely pertinent—is generalized to all dynamic systems thanks to dissipation envelopes.

It applies as much in observability as in modal controllability (see [

Then we have a complete mechanical model, which can be declined quite simply to thin solid elements (straight beams, curved beams, plates, shells) in their Newtonian interaction (incompressible perfect fluid, compressible perfect fluid, Newtonian fluid). They have given rise to advanced structural modelling involving solid-fluid metamodes (see

The models have been performed with COMSOL Multiphysics® software with the shell elements, and with Newtonian and non-newtonian fluids. This numerical approach has been performed both for the cryotechnic tank and the organic systems; the aim here is precisely to show that this kind of model, with similar equations, can be applied to both kinds of systems.

Indeed, the general matrix model applies as well to Lagrangian displacements in the case of a solid as to Eulerian pressure fluctuations in the case of a fluid, which makes it possible to offer hybrid models combining, for example, very

simply shells and perfect fluids (see [

These models are transposable to living systems such as a cell (membrane with fluid interaction) or an organic structure: shell in interaction with soft materials, and in a specific environment. For instance, in biodynamics, the tonopic coupling in the cochlea, with the basila membrane, is modelised by the plate-fluid theory (

Last but not least, in the dissipation itself, a linear term of viscous damping can be distinguished from a contribution related to other frictions, both in solids and fluids. In solids, the non-linear part represents precisely the damage or defects (delaminations, cracks, discontinuities) and in the fluid, the realization of tests with equivalent granular fluids in similarity to Froude makes it possible to observe this type of irreversibility (see

On the basis of the resonance zone—delimited by dissipation in the probabilistic sense—both observability and controllability of the models can be fairly easily shown. Indeed, this zone of unknowingness (in the Heisenberg sense) is the intersection between an observable and non-controllable past and a non-observable & controllable future. And what guarantees this possibility is precisely that it is based on principles three and four, and that there are in fact conservations that address both observable kinetic quantities (Routh conservations) and controllable energy quantities (Hamilton-Painlevé conservation).

This synthesis intrinsically covers the non-linear domain (see [

Remaining in the field of solid-fluid mechanics, the set of the 4 terms of dynamic equations thus constitutes an important potential source that can feed both cellular simulators in their interactions and organic simulators (see [

On the basis of the resonance zone—delimited by dissipation in the probabilistic sense—both observability and controllability of the models can be fairly easily shown. Indeed, this zone of unknowingness (in the Heisenberg sense) is the intersection between an observable and non-controllable past and a non-observable and controllable future. And what guarantees this possibility is precisely that it is based on principles three and four, and that there are in fact conservations that address both observable kinetic quantities (Routh conservations) and controllable energy quantities (Hamilton-Painlevé conservation).

The lagrangian bended plate moves with the 2.5D classical equation, given by principles 2 and 1, inducing W displacement along Z (X, Y, t field):

D Δ Δ W + ρ h W ° ° ≡ p Z (6)

D ≡ E h 3 12 ( 1 − ν 2 ) ∈ ℝ + ∗ lineic bending stiffness;

( ρ , h ) ∈ ℝ + ∗ 2 platevolumic mass & thickness;

p Z ( X , Y , t ) ∈ ℝ external normal surface static load applied on the plate.

It is pertinent to observe that this equation is similar to the prestressed membrane law. This is also the case for the 1D rod equation, similar to the rope law.

In parallel, we could consider a perfect compressible fluid in a tank, precisely in a rigid tank (galilean solid reference) closed by a flexible plate (see

This demonstrator is the typical flexible structure for the basic hybrid dynamic representation:

1) W ( X , Y , t ) lagrangian relative dynamic Z-displacement 2D field in the plate/membrane;

2) p ( X , Y , Z , t ) eulerian pressure fluctuation 3D field in the perfect fluid.

In the fluid, the dynamic behaviour is given by the Euler’s equation (Navier-Stokes 2 without viscous stress):

( ρ V → ) ° + d i v → ( ρ V → ⊗ V → + p I ¯ ¯ ) ≡ 0 → (7)

⊗ dyadic product

completed with the continuity condition (N-S 1):

ρ ° + d i v ( ρ V → ) ≡ 0 (8)

Thus, Equations (5) to (7) drive the motion of the system (plate or fluid).

We observe that implicitly in the linear modal frame, the Euler-Helmholtz equations in the (perfect) fluid can be written as:

{ a 2 Δ p − ω i 2 p = 0 in Ω ∂ Z p − ρ f ω i 2 W = 0 on ∂ Ω ∭ Ω ( fluid ) p d X d Y d Z + ∬ ∂ Ω ( plate ) ρ f a 2 W d X d Y = 0 ( Stokesresultant ) (9)

( a , ρ , p ) ∈ ℝ + ∗ 3 fluid sound celerity, volumic mass, eulerian pressure fluctuation;

Ω ⊂ ℝ 3 volumic fluid domain;

∂ Ω ⊂ ℝ 2 surfacic plate interface (plate surface).

It is pertinent to observe that this equation is similar to the prestressed membrane law. This is also the case for the 1D rod equation, similar to the rope law. Thus, Equations (5) to (7) drive the motion of the system (plate or fluid).

Then, the 2D plate will be discretized, as well as the 3D fluid volume. At each node of the system will be attached: a Lagrangian W displacement at each node of the shell (plate) and an Eulerian p pressure fluctuation at each node of the (perfect) fluid. In total, we have nodes: n_{s} structural (solid) nodes in the plate and n_{f} nodes in the fluid.

n = n s + n f (10)

The dynamic digital modes of the system will be valid for these n hybrid Degrees of Freedom (DoF).

At this level—and in order to reinforce the full mechanical pertinence of these hybrid modes—it is useful to observe that the Mach-Hooke equation previously mentioned in a rod has exactly the same expression in displacement and in stress:

U ° ° − c N 2 U ″ = 0 ⇔ σ X X ° ° − c N 2 σ ″ X X = 0 (11)

The reason is simply the rheological relation between σ X X tension/compression stress and U deflection in a rod:

σ X X = E ε X X = E U ′ (12)

This allows to represent the global dynamic system (plate + fluid) with n hybrid DoFs: n_{s} W Lagrangian normal deflections and n_{f} pressure Eulerian fluctuations.

Thus, n_{s} is the number of nodes in the plate, and n_{f} the number of nodes in the fluid. The hybrid impedance matrix becomes natural to model the motion of the system:

[ Z i j ( ω k 2 ) ] ( r j ) ≡ ( 0 j ) , i , j , k = 1 , ⋯ , n { [ Z α β ] s [ C β δ ] c [ C β δ ] t c [ Z γ δ ] f } ( ( W β ) s ( p δ ) f ) ≡ ( ( 0 β ) ( 0 δ ) ) | α , β ≡ 1 , ⋯ , n s γ , δ ≡ 1 , ⋯ , n f (13)

with [ Z α β ] s solid n_{s}× n_{s} structural impedance matrix of the solid plate 2D domain;

[ Z γ δ ] f fluid n_{f}_{ }× n_{f} acoustic (Helmholtz) impedance matrix of the fluid 3Ddomain;

[ C β δ ] c coupling impedance n_{s}× n_{f} matrix structure-fluid 2.5D domain;

( W β ) s plate structural n_{s} vector of W lagrangian deflections (plate DoFs);

( p δ ) f fluid n_{f} vector of eulerian pressure fluctuations (Helmholtz DoFs).

By the principles 1 and 2 of thermodynamics (Galilean reversibility & Maxwell’s reciprocity) all [Z] matrices are symmetrical, real positive semi-definite in classical physics.

In addition, as the fluid is here considered perfect, the Helmholtz matrix is diagonal in any base (for any paving/meshing of the fluid).

Thus, for the considered mesh/pavage, the information contained in the [Z] matrices are summarized as follows:

1) the [Z]_{S} solid/plate matrix is defined by its upper n_{s}_{ }× n_{s} triangle:

[ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ s y m ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ∗ ] ;

2) it is thus defined by n s + ( n s − 1 ) + ( n s − 2 ) + ⋯ + 1 = n s ( n s + 1 ) / 2 independent scalar numbers;

3) the [Z]_{f} fluid matrix is n_{f} diagonal [ ∗ 0 ∗ ∗ 0 0 0 0 0 0 s y m ∗ 0 0 ∗ ∗ 0 ∗ ∗ ∗ ] representable by the n_{f} vector: ( H 1 H 1 ⋮ H n f ) ;

4) it is thus defined by n_{f}independent scalar numbers;

5) then the [C]_{sf} solid-fluid coupling matrix is a full n_{s }× n_{f}one;

6) it is thus defined by n_{s }× n_{f} independent scalar numbers.

Consequently, the global {Z} matrix (n× n) which models the complete plate-fluid dynamics, is represented by ( n s + 1 ) n s / 2 + n f + n s n f = ( n s + 1 ) ( 2 n s + n f ) / 2 scalar numbers. That is the implicit dynamic representation of the global system, valid for any method of digitization (FEM, DEM, BEM—F, D, B).

This preliminary approach can be extended to large & complex systems, for which n_{s}& n f ≫ 1 . For this purpose, let’s consider that the plate (which represents a generic thin structure) is linked with a cubic volume. With these assumptions, if the mesh is coherent between the fluid and the solid: n f ~ n s 3 / 2 . The previous expression of representative scalar numbers is then depending only on the structural representation ( n s + 1 ) ( 2 n s + n s 3 / 2 ) / 2 ~ 0.5 n s 5 / 2 .

To conclude this approach of the digitization of a shell—perfect fluid n_{s }× n_{f}structure, the expression of the N D o F number of scalars required (necessary and sufficient) for a complete implicit digital representation—i.e. the number of digital hybrid Degrees of Freedom—is given by:

2 N D o F = ( n s + 1 ) ( 2 n s + n f ) ; 2 N D o F ¯ ~ n s 5 / 2 (14)

That is the implicit dynamic representation of the global system. The application field is here the domain of thin shells with gaseous or ultra-light fluid (LH2).

If we consider an optimized/local mesh, focusing consequently on a FEM-type model, the matrices are not only symmetrical, but also hollowed. With standard quadrangular/hexaedral elements, the number of representative scalar numbers is reduced to 8 n s + 24 n f + 12 n s = 8 ( 3 n s + 4 n f ) .

In explicit solution, the corresponding [M] and [K] matrices are well adapted to Cholesky’s direct inversion.

Thus, the standard FEM cost in terms of scalar number (memory size for Z/M/K matrices) is given by:

S F E M = 8 ( 3 n s + 4 n f ) ; S F E M ¯ ~ 32 n s 5 / 2 (15)

This figure is valid to size the memory and power of the digital process, and is thus relevant in terms of internal (for computation) representation of the global system, whereas the N D o F number represents the physical external approach.

The symmetrical impedance matrix, including the structural and fluid dynamics, gives implicitly—by diagonalization—the eigenvalues & eigenshapes, Lagrange-Euler metamodes (solid + fluid). That is the modal dynamic representation of the global system. Inventoriating the information in the complete external (meta-)modal base, we obtain:

1) n ω i eigenvalues (pulsations);

2) n ( p j ) i eigenvectors (of W or p amplitudes for each mode).

Each eigenvector is defined with a multiplicative constant, and so represented by n − 1 independent scalars. Thus, the whole modal base requires n + n ( n − 1 ) = n 2 independent scalars. That is the complete dynamic external signature (or “DNA”) of the solid-fluid system.

2 N m o d = n 2 ≡ ( n s + n f ) 2 ; 2 N m o d ¯ ~ n s 3 (16)

In dynamic engineering, it is often possible to truncate the base. A performing strategy consists of considering the three first modes (the three lowest frequencies) or the three first ones plus a fourth representative of a specific phenomenon (optional local mode or sensitive mode, considering the external loads).

For a well-conditioned system, this approach allows a precision better than 1%. Consequently, the truncated dynamic model of the global system will be represented by 4n independent scalars, which is a very simple representation, the basic dynamic DNA of the system.

2 N t r u n c 4 = 8 n ≡ 8 ( n s + n f ) ; 2 N t r u n c 4 ¯ ~ 8 n s 3 / 2 (17)

Through the boundary-integral methods, it is allowed to represent a plate (2D) with a mesh on its 1D frontier, and the volumic fluid by its frontier with the plate. Thus, with this ultimate condensation, the [M] and [K] global matrices (plate + fluid) are reduced to a n_{s }× n_{s}one. These matrices are full and non-symmetrical but solvable explicitly by iterative methods. Provided that M and K switch, then they admit a common basis for diagonalization and the implicit approach can be applied to this integral model.

Considering perfect fluids, the previous fluid equations are unchanged in the general case of a volumic tank, both for internal and external fluids, the X, Y, Z deflection on the shell being referenced as U, V, W:

δ P ¯ → ≡ ( U V W ) X , Y , Z

If n_{n} is the number of nodes in the shell (structural nodes) then the number of shell structural DoFs is n s = 3 n n . The fluid equations are then connected with the Reissner-Naghdi shell equations where the lagrangian deflection of mid-surface is expressed in the local curvature axes (see [_{XZ}, Q_{YZ} hybrid variables, but U, V, V are actually the structural ones for the shell dynamics.

The previous system shown on the plate example remains valid, the solid mesh being on curvilinear coordinates on the shell. The unique difference is, for the external fluid only, the Stoke’s resultant condition being replaced by Sommerfeld limit condition:

p → 0 when ‖ O P → ‖ → ∞ (Sommerfeld) (18)

Connected with the Reissner-Naghdi U, V, V equations of shells, the previous Euler-Stokes-Sommerfeld system generates the global dynamics of the shell in a perfect fluid.

A viscous newtonian fluid requires the third Navier-Stokes equation, including the shear stress. Consequently, the standard shell digital dynamic implicit matrix equation of a Reissner-Naghdi shell in a newtonian fluid keeps the shape defined in (13); the shapes of the implied block-matrices are as follows:

[ Z α β ] s = [ Z 11 Z 12 ∗ 0 0 ∗ Z 22 ∗ ∗ 0 ∗ ∗ ∗ s y m ∗ ∗ Z n s n s ] ;

[ Z γ δ ] f = [ H 11 S 12 S 13 ∗ 0 H 22 S 23 ∗ 0 ∗ ∗ ∗ s y m ∗ S ( n f − 1 ) n f H n f n f ] ;

( W β ) s = ( ( U 1 V 1 W 1 ) ( U 2 V 2 W 2 ) ∗ ∗ ∗ ( U n s V n s W n s ) )

The field application is here thick shells in viscous fluids, such as hydraulic systems or submarine skins, but these models can also be combined with the previous ones.

To determine the level of the model (viscous/non-viscous, shell/membrane, static/dynamic) it is necessary to examine several adimensional numbers, considering the dimensions and characteristics of the tank (

For the fluid model, it is possible to consider Stoke’s equation valid for Newtonian fluids:

ρ ∂ V → ∂ t + ρ ( V → ⋅ g r a d → ) V → + g r a d → p − μ Δ V → − ( μ ′ + μ 3 ) g r a d → ( d i v V → ) − ρ f → v o l ≡ 0 → (19)

Σ ¯ ¯ f ≡ 2 μ E ¯ ¯ f ° + μ ′ t r ( E ¯ ¯ f ° ) I ¯ ¯ (20)

where V → Eulerian velocity in the fluid;

μ , μ ′ primary & secondary Newtonian viscosities of the fluid;

Σ ¯ ¯ f , E ¯ ¯ f ° stress & strain-rate tensors in the fluid;

I ¯ ¯ unitary Kronecker’s tensor.

It is thus possible to observe different adimensional numbers:

1) Acoustic Reynolds Number R e a c ≡ ω c L c 2 ν c | ν c kinematic viscosity ( ν c ≡ μ ρ )

High: model with perfect fluid. Low: model with viscous fluid.

2) Oscillatory Froude Number F r o s c ≡ L c ω c 2 g |g equivalent gravity, for static

loads.

High: dynamics independent from statics. Low: coupled model.

3) Flow Detachment Number F d s ≡ V t ω c R e x t |tangential velocity on the shell;

R e x t shell curvature radius, on external normal.

High: turbulent model. Low: laminar Coanda model.

4) Weber’s Number W e f ≡ ρ ω c L c 3 σ |σ surface tension in the fluid.

High: surface tension neglected. Low: with surface tension.

Important parameter for multiphasic systems.

The topology, geometry and singularity of the structure and its control-model can be used as an indicator of its entropy. That approach is also an occasion to mention the modern methods (in progress) whose potentiality open wide perspectives in the domain of complex structure thermodynamics and, consequently, in the field of cellular and organic modeling.

In comparison with FEM and DEM, the Boundary Elements Method (BEM) is the most promizing one for such combined problems. In the case of an uncompressible fluid, the Equation (6) is reduced to anholonomic condition between the boundary (plate) DoFs. In the general case of the shell with a newtonian fluid, the condensation leads to a ns dimension problem. In a certain way, it is equivalent to condensate the behavior of the cellular tank to its boundary represented by the shell. The global BEM can be used for multiphasic 3D structures, with 2D external and internal interfaces (see [

In the context of BEM, particular boundary elements can be developed, in similarity with Robinson’s finite element for a crack (see [

But we can offer to go further. Indeed, the shell itself can be meshed according to its interfaces, which are reduced to the singular discrete lines between panels (see [

Thus, a cubic shell will be meshed by its arches, which is an important advance in terms of digital condensation. In this perspective a complex assembly of

shells will be represented by the intersection lines, or even the points (nodes) interfacing these lines. We could thus obtain a very concentrated thermodynamic model of the structural cell. Its implicit resolution gives metamodes, and its explicit integration allows to explore transient evolutions. And the ratio between the explicit and implicit results could even constitute an indicator of the entropy of the controlled system.

The methods presented in this article bring together achieved and present work on current aerospace structures and research on innovative methods in two directions. Vibratory dynamics thus open tracks that integrate both the modeling of structures and the modal representation of information.

The first explored direction concerns active and passive control: Generalized Predictive Control (see [

The second direction concerns the optimal hybrid modeling of thin structures coupled with Newtonian fluids. The proposed synthesis, in the presented structuration of the equations, integrates fluid-solid meta-modes in Gauss’ sense, as well as proposals for condensation by integral methods with perspectives on the two-levels condensed representation: by boundary conditions and modal extraction, leading to vibratory optimal compression of models.

Jean-Emmanuel Chambe, Miguel Charlotte, Jean-Emmanuel Alibert, Matthieu Hospital, Yulia Akisheva; ISAE-SUPAERO. Vincent Lapoujade, Charlotte Michel, Tess Legaud, Thomas Pierrot, Antoine Guilpin; DynaS+. ValiaFascio, Mathieu Delorme; Ateca.

The authors declare no conflicts of interest regarding the publication of this paper.

Gourinat, Y., Christon, L. and Lachaud, F. (2021) Structural Dynamics in Biology: A Bridge Given by Implicit Vibratory Crossed Models. Engineering, 13, 237-256. https://doi.org/10.4236/eng.2021.135018

Below is proposed a sequential and synoptic presentation of fourteen dynamic equations with comparable 2 + 2 term structures.

Thus is highlighted the explicit relation with the four principles of thermodynamics, and thus open the possibility of both implicit and explicit resolutions.

Mechanics

V → ° − g → ≡ 0 → ; m O P → ° ° − f → ≡ 0 → ⇒ [ M ] ( q ) ° ° + [ K ] ( q ) = ( 0 )

(Newton-Feynman,pr. 1-2) (1.1)

Kinetic constant ∂ E l a g ∂ q k ≡ 0 ⇒ d d t ( ∂ E l a g ∂ q k ° ) ≡ 0

(Lagrange-Routh,pr. 2) (1.2)

Control energy ∂ E l a g ∂ t ≡ 0 ⇒ d d t ( E h a m − ∑ j ∂ E p o t ∂ q j ° q j ° ) ≡ 0

(Hamilton-Painlevé,pr. 3) (1.3)

Dynamic digital system 〈 L 〉 j ( E l a g ) − Q j f e y n ≡ Q j s t a t

(Lagrange-Feynman,pr. 1-4) (1.4)

Structures

Continuous 1D solid U ° ° − c N 2 U ″ ≡ 0 , c N 2 ≡ E ρ

(Hooke-Mach,pr. 1-2) (2.1)

Continuous 1.5D solid ρ S V ° ° + E I V ′ ′ ′ ′ ≡ λ Y s t a t

(Bresse-Lagrange,pr. 1-2) (2.2)

Continuous 2.5D euclidian ρ h W ° ° + D Δ Δ W ≡ p Z s t a t , D ≡ E h 3 12 ( 1 − ν 2 )

(Kirchoff-Lagrange,pr. 1-2) (2.3)

Cont. 2.5D non-euclidian D X , Y , Z − F X , Y , Z = p X , Y , Z s t a t ; T Y , X + Q Y Z , X Z = 0

(Reissner-Mindlin,pr. 1-3) (2.4)

Solid 3D ρ v → ° − d i v → Σ ¯ ¯ − f → v o l d y n = f → v o l s t a t

(Navier-Lamé,pr. 1-4) (2.5)

Physics

Newtonian fluid { ρ ° + d i v → ( ρ v → ) = 0 ( ρ v → ) ° + g r a d ¯ ¯ ( ρ v → ) ⋅ v → + ρ v → d i v v → − d i v → Σ ¯ ¯ − f → v o l = O → ( ρ e t o t ) ° + d i v ( ρ e t o t v → − Σ ¯ ¯ v → − λ T g r a d → T ) − f → v o l ⋅ v → = q v o l

(Navier-Stokes pr. 1-3) (3.1)

Newtonian gas { e g ≡ c v T h ≡ e g + p ρ ≡ c p T T d s g ≡ d e g + p d ( 1 ρ )

(Joule-Thomson pr. 1-3) (3.2)

Mechanical discrete 3D 8 π 2 m E h a m ψ + h 2 Δ ψ ≡ 0

(Schrödinger pr. 2, 1) (3.3)

Discrete 3D half-spin 2 π c 0 ( m c 0 [ α 0 ] + α → ⋅ p → ) ψ ≡ i h ψ °

(Dirac-Schrödinger pr. 3, 2, 1) (3.4)

General continuous 4.5D 8 π G 0 T ¯ ¯ − c 0 4 R ¯ ¯ ≡ c 0 4 ( Λ c − t r R ¯ ¯ 2 ) μ ¯ ¯

(Einstein-Hilbert pr. 2, 1, 4, 3) (3.5)