WarningYour internet explorer is in compatibility mode and may not be displaying the website correctly.
You can fix this by pressing 'F12' on your keyboard, Selecting 'Document Mode' and choosing 'standards' (or the latest version
listed if standards is not an option).

Structural Analysis with Thin Elastic Layers

Some structural applications involve thin or high aspect ratio structures sandwiched between other relatively low aspect ratio structures. For example, if a piezoelectric transducer is glued on the surface of a mechanical system, the thickness of the adhesive layer is very small in comparison to the two structures it glues together. Numerical modeling of such a thin layer in two or three dimensions requires resolving it with an appropriate finite element mesh. This can result in a large concentration of finite elements near the adhesive layer, leading to high computational cost and time. To avoid this issue, one common assumption many numerical simulations make is to assume perfect bonding between the two structures. By making this assumption, the numerical model ignores the effect of the flexibility of the adhesive layer. This will lead to inaccurate results because the adhesive layer is not infinitely stiff. In this blog post, we will show how to model such thin layers using COMSOL Multiphysics — without the need to draw the layer’s thickness explicitly, while still accounting for the effect of the thin elastic layer. This can lead to a more efficient structural analysis with significant reduction in computational cost and time, without sacrificing the accuracy of the simulation.

About the Thin Elastic Layer Boundary Condition

The Thin Elastic Layer boundary condition is suitable for modeling systems involving very thin elastic layers. This boundary condition has both elastic and damping properties and acts between two parts to model a thin elastic layer with specified stiffness and damping properties. On an interior boundary, the thin elastic layer decouples the displacements between two sides of the boundary. The two boundaries are then connected by elastic and viscous forces with equal size but opposite directions, proportional to the relative displacements and velocities. This means that you can replace any thin elastic layer in your geometry with just an interior surface with zero thickness, then use the Thin Elastic Layer boundary condition on the interior surface and assign it the properties of the elastic layer. You can include stiffness and damping values of the elastic layer. The Thin Elastic Layer boundary condition can be used in any type of analysis, including stationary, time dependent, frequency domain, and eigenfrequency.

Thin Elastic Layer boundary condition settings.

Structural Analysis of a Composite Piezoelectric Transducer with an Adhesive Layer

Next, let’s analyze an application that can benefit from the use of the Thin Elastic Layer boundary condition. Consider a composite piezoelectric transducer, for example. The composite piezoelectric transducer consists of piezoceramic, aluminum, and adhesive layers. The adhesive layer binds the piezoceramic and aluminum layers. In the simulation considered here, the length of the longest side of the adhesive layer is 27.5 mm and it can be 0.02 mm thick. Thus, the aspect ratio of the adhesive layer is quite high — it’s equal to 1.38 X 103. The piezoceramic and aluminum layers have the same length for the longest side, but are far thicker, with values equal to 10 mm and 5 mm, respectively. An AC potential is applied on the electrode surfaces of both sides of the piezoceramic layer.

The objective of the following simulation is to find the first six eigenfrequencies of this composite piezoelectric transducer without explicitly drawing and meshing the thin adhesive layer. Instead we will consider its effect by using the Thin Elastic Layer boundary condition available in the MEMS, Structural Mechanics, and Acoustics Modules. This transducer model is taken from the Model Library of the MEMS Module where the adhesive is modeled explicitly for a thicker adhesive layer. The figure below shows that the geometry containing a very thin adhesive layer is replaced by a geometry with just an interior boundary, which is then assigned a Thin Elastic Layer boundary condition.

The Piezoelectric Device interface is used to simulate the composite structure. The boundary conditions consist of an applied potential and ground boundary conditions on both sides of the piezoceramic layer. The inner boundary representing the adhesive layer is assigned the Thin Elastic Layer boundary condition. Given the Young’s modulus (E), Poisson’s ratio (v), and the thickness (t) of the thin elastic layer, the spring constant per unit area in directions normal (kn) and tangential (kt) to the boundary can be estimated as: and , where G is the Shear modulus. Note the two asymptotes for kn: if v is low or closer to 0, , and if v is high or closer to 0.5, , where K is the bulk modulus. The values of normal and tangential stiffness per unit area are used for spring constant per unit area in the Thin Elastic Layer boundary condition.

An eigenfrequency calculation is performed for values of the thickness of the adhesive layer ranging between 0.01 mm and 0.38 mm.

Results of our Structural Analysis

The following results show a comparison of six eigenfrequencies calculated using the Thin Elastic Layer boundary condition against the eigenfrequencies calculated by modeling the adhesive layer explicitly. The eigenfrequencies calculated using the Thin Elastic Layer boundary condition are indicated by an asterisk marker on solid lines, and the eigenfrequencies calculated by modeling the adhesive layer explicitly are indicated by dashed lines without any markers.

The above results show that for lower adhesive layer thickness values the solution obtained by using the Thin Elastic Layer boundary condition closely matches the solution obtained by modeling the adhesive explicitly. For large adhesive layer thickness values, the differences between the solution obtained from two calculations will increase with increased adhesive layer thickness. A sufficiently thick adhesive or elastic layer can be modeled explicitly. However, for very low values of the adhesive layer thickness or for simulating any thin elastic layer, the Thin Elastic Layer boundary condition is a more efficient way of modeling the effect of the thin layer without modeling the thin layer explicitly, while still accounting for its thickness and elastic properties.

There are several other important uses of the Thin Elastic Layer boundary condition. One such application is for large CAD models represented by a COMSOL Assembly, where the Thin Elastic Layer boundary condition can be used to approximate mechanical contact conditions for frequency-domain studies. Other uses, typically in combination with using a nonlinear stiffness, include modeling of fracture zones in geomechanics, simplified mechanical contact, and delamination of adhesive layers.

Comments

Hi,
I would like to know if I can use the thin elastic layer property to model compliant, ie. elastic electrodes for a capacitor. I should be able to specify the conductivity also.

Thank you

Chandan Kumar April 11, 2014 at 4:27 pm

Dear George,

Thin elastic layer can be used for structural analysis of thin layers. Similar boundary conditions exist for other types of physics like, Distributed impedance in Electric currents physics or Distributed capacitance in Electrostatics physics. If you have further questions about the use of COMSOL in the analysis you are doing, please contact us by emailing your question to us at support@comsol.com.

Best regards,
Chandan

Andrew Melro September 3, 2014 at 5:35 pm

Hi,

I was just wondering where the formulations for the spring constant per area came from? I would like to get some more understanding about how they are formulated

Thanks,
Andrew

Peter York January 2, 2015 at 3:52 pm

Great post.

Joseph O’Day April 11, 2015 at 1:59 pm

I’m also interested in the spring constant formulation. It seems as though it would be possible to simulate a pseudo-plasticity effect if they could be reformulated to be functions of extension.

Chandan Kumar April 13, 2015 at 8:42 am

Dear Andrew and Joseph,

There is no rigorous theoretical derivation for the spring constants that can be referenced here. The equations used are based on the asymptotic values explained in the article and its behavior for other Poisson’s ratio values.

Best regards,
Chandan

Joseph O’Day April 24, 2015 at 1:42 pm

It seems as though it’s just Hooke’s Law in 3D. Then the strain terms are effectively replaced with the spring deformation divided by the thickness layer.