Main Research Contributions

Mathematical Modeling

   Beam Models:

1.      Linear elastic beam  model (with D. Russell, 1996)   

It is known that the classical Euler–Bernoulli beam theory is valid only for long span (equivalently, thin) beams. In 1921, S. P. Timoshenko proposed a dynamic beam theory with two generalized displacements; i.e., the deflection w(x) and the transverse shear deformation  v(x).

In 1996, Gao and Russell proposed an extended beam model which allows the shear deformation vary in y-direction, v(x,y).

 

Static model

Dynamical model

 

Finite element Simulations extended beam model

 

Numerical simulations and comparison of Euler- Bernoulli beam, Timoshenko beam with Gao-Russell beam model

 

Clamped beam

Clearly, these numerical simulations show that by allowing the shear deformation, the Gao-Russell beam is much weaker than both Euler-Bernoulli beam and Timoshenko beam models.  This extended beam model can be used for studying contact problems with friction.

 

 

 

2. Elasto-perfectly plastic beam model:

The dual problem of this model is equivalent to a quadratic minimization problem subjected to linear and nonlinear inequality constraints given below

 

 3. Large deformation nonlinear elastic beam models (GAO, 1996)

  wtt  + k wxxxx + (ν p – a wx2 ) wxx  = f(x,t)

where k,  a and ν are given constants, p is a given axial load, f(x,t) is a given distributed input.

 This equation has at most three solutions w(x,t) at each given (x, t), represents the two possible buckling states (upper and down positions), and one unbuckled (unstable) position. These three solutions could lead to chaos in dynamical vibration for certain given data k, a, ν, the axial load p and the distributed load f(x,t)

This equation is also equivalent to one-dimensional Landau-Ginzburg equation in phase transitions.

In literature, this beam model is known as the Gao beam

By using this beam model, the following interesting phenomena in chaotic vibration have been discovered in 2003.

·         Tri-chaos: for a given periodic load f(x, t) and a linear axial load p = po  + t k, where  po  and k are positive constants, the beam may experiences three chaotic vibrations before final collapse (see the following figures).

 

·         Meta chaos: there exists a very short transition period before the chaos occurred, see the following figures.

 

chaos2.jpg

chaos3.jpg

                                Meta Chaos                                                                                      Chaos vase:   A closed vision

              

 

4. Some more general beam models are proposed in 2000

·         Fourth-order large deformation dynamical beam model with shear effects:

where a = 3 h (1 - n2 ), and the shear deformation v(x,y) must be an old function of y.   Ignoring the shear effects  (v = 0 ), this beam model is reduced to the beam model given above.

·         Second-order large deformation dynamical beam model with shear effects:

 

where the shear deformation x could be any function of (x,y), so it can be used for studying frictional contact problems

 

·         Second-order static beam model with shear effects:

The total potential energy is a nonconvex functional

Its canonical dual is given by

Which is concave if the axial stress  is positive for all x.

 

 

Multi-scale modeling for super-conductivity materials (click here for details)

 

 

 

 

 

 

 

 

 

 

Eigenvalue Problem on Optimal Surface (Gao and S.-T Yau, 1991, published in Gao’s book)

Let Ω be a domain in R2 with boundary. A regular parametrical surface S  in R3  is a C2 mapping

 c : W           S

The eigenvalue problem on the optimal surface is a minimization problem

Where the contrain sets are

The criticality condition leads to the following coupled nonlinear partial differential equations

Where  is the Laplace-Beltrami operator, I3  = det C , and H(c) is the mean curvature of the surface

And  is the second fundamental form of the surface with N = { Ni } being the unit vector normal to the surface.

This problem plays an important role in optimal design of many engineering structures. To solve this coupled nonlinear partial differential equations is fundamentally difficult. By using numerical simulation, some very interesting singularities on both eigenfunction and optimal surface are discovered (see below)

Singularities on eigenvalue problem over the optimal surface.

Some movies about finite element simulations for dynamical post-buckling analysis of a rubber diaphragm created by my Ph.D. student Axinte Ionita:
 

Movie 1(5.2MB)

Movie 2(2.32MB)

Movie 3(8.9MB)

Movie 4(3.26MB)


 

 


Optimal Shape Design for extended beam model.
 

Back to the home of David Yang Gao

E-mail to me at gao at vt.edu