November 22, 2022 at 6:28 pmMani BhumaSubscriber
I am fairly new at Ansys and have been having some issues with accurately running simulations for a composite honeycomb material. Once I enter the orthotropic properties for the specific material I find the static structural results are fine; however, the modal results give ridiculous shapes as deformation results.
I have seen some of @peteroznewman posts on this topic, but am having some issues replicating the concepts.
I would appreciate help/feedback on this please!
November 23, 2022 at 1:49 ampeteroznewmanSubscriber
In Mechancial, right click on the Mesh branch and Clear Generated Data, then close Mechancial. In Workbench, save the Project. Then do File, Archive and save the .wbpz file. Upload that file (not the .wbpj file) to a sharing site such as Google Drive, OneDrive or Jumpshare. Copy the link to that file and paste it in your reply. Say what version of Ansys you are using in your reply. I will take a look at it.
November 23, 2022 at 8:42 amErik KostsonAnsys Employee
Perhaps this helps also.
The actual absolute values of the displacement results (mode shape) do not mean anything in modal analysis, as these are scaled/normalised eigenvectors/mode shapes.
Only the frequency and the overall mode-shape are of interest (so if the mode shape is bending, torsional, etc), but not the absolute/actual values of the displacements.
See the help manual (search for mass normalized), and our course for more info on this https://courses.ansys.com/index.php/courses/modal-analysis-in-ansys-mechanical/lessons/how-to-perform-modal-analysis-lesson-1-2/
Also this helps: https://www.eng-tips.com/viewthread.cfm?qid=83495
All the best
November 23, 2022 at 11:17 pmMani BhumaSubscriber
I followed the steps you mentioned. Here is the link:
I am using Ansys Workbench 2022 R1.
When you open the file you will see a 2 inch X 2 inch composite plate with an aluminum honeycomb center. My goal is to gather material information about the honeycomb core. In the future I will running a much bigger simulation with a larger composite plate. Due to the size of this new simulation it would be unreasonable to model the honeycomb core. So, using the material properties I find with this miniature simulation I will transfer the information to a model with 3 basic solid bodies. 2 rectangular solid bodies for the carbon fiber skins and 1 rectangular solid body for the honeycomb core.
Once I verify the static structural and modal results on this simulation is accurate would it be possible to extract information such as young's modulus in each axis from the honeycomb body for this purpose?
November 24, 2022 at 12:20 pmpeteroznewmanSubscriber
You mentioned orthotropic material properties in your original post, but when I look at your materials, they are all isotropic.
When I look at your Ansys project, you have modeled the hexagonal honeycomb cells using surface geometry and assigned a wall thickness and material properties for solid aluminum. By modeling it this way, you may not need an orthotropic material model. However, your model has not captured an important feature of honeycomb core, the double wall thickness created as a result of the manufacturing process. Foil sheets of thickness t are bonded together with strips of adhesive of width h and the strips of adhesive are staggered on the odd and even layers of foil. When that stack of foil sheets is pulled apart, a honeycomb cell is created, but notice in the figure below that some of the walls of the hexagonal cell have a wall thickness of t while others have a wall thickness of 2t.
Your model does not represent that.
A simpler way to model a honeycomb core sandwich is to use orthotropic material properties. Instead of modeling the hexagonal cells with surfaces, the entire core is represented by a solid body assigned an orthotropic material. The manufacturer provides the critical material properties of this core on their webpage. The properties are orthotropic relative to directions in the core called L and W. The L direction is horizontal in the figure above, and that is the length of the foil sheets. The W direction is vertical and that is the direction the sheets were pulled to open up the cells. Once face sheets are bonded to the top and bottom of the core, and the sandwich panel has some bending applied, you would expect a higher modulus in the L direction because that is shearing the walls with the 2t thickness, while the lower modulus in the W direction is shearing the angled walls with the thickness t. Below is a copy of a manufacturers data sheet. Notice that the Shear Modulus in the “L” direction is about twice as large as the “W” direction. The other critical property is the compressive modulus.
Ansys Orthotropic material models require more constants than those three, but those are the critical ones.
Composite facesheets can also be orthotropic, but if many layers are stacked at multiple angles, the sheet can end up with nearly isotropic properties, especially if the laminar are fabric and not uni-directional tapes.
The surfaces representing the composites would be modeled coincident with the top and bottom faces of the solid core. Use Shared Topology so that the mesh of shell elements shares the nodes on the solid elements that mesh the core. This avoids using any bonded contact.
November 28, 2022 at 6:26 pmMani BhumaSubscriber
I originally attempted a solid body with orthotropic materials as you mentioned. Similarly, to the datasheet you provided I have the exact same datasheet of the material I am using. However, as you mentioned, the datasheet only provides three constants of the many. I tried to use an educated guess to estimate the other values, but I kept getting ridiculous mode shapes in the modal analysis. My solution to this issue was to recreate the composite with the honeycomb completely modelled, as I know the isotropic material the honeycomb itself is made of. Then, use the results of this new model to adjust the material values on the solid body until I got similar results. Please let me know if there is a simpler solution than this.
Thank you for catching that important feature I was missing in the honeycomb. Is there a way to model that into my simulation with surfaces?
November 28, 2022 at 8:55 pmpeteroznewmanSubscriber
You opened this topic with a statement that you were getting ridiculous mode shapes from an orthotropic material. I was curious to see those shapes, but you took the orthotropic material out when you shared your model so I have not yet seen what you meant by that. Will I get a chance to see that?
You moved on to explicitly modeling the honeycomb cells. You can assign double the thickness to those faces that need it.
You say you plan to use the results of the honeycomb cell model to adjust the values of the orthotropic core material until the results match. You make it sound simple, but what you are saying is actually complicated.
What loads and boundary conditions you will use on the two models? What result will you pull from each model?
There are nine material constants in the orthotropic material model, how will you adjust those nine numbers based on a comparison of the results between the two models? Will you keep the honeycomb vendor values constant? There is a constraint equation involving the nine constants that must be satisfied, so you are not free to choose values that don't satisfy the constraint. How will you incorporate the constraint equation in your value update method?
November 29, 2022 at 1:52 amMani BhumaSubscriber
Here is the link to the simulation of solid bodies made of orthotropic materials:
In the honeycomb cells model, do I assign double the thickness of those faces by using the 'surface coating' feature?
Please correct me if I am mistaken in the following. Out of the nine constants, the vendor has provided three for me (Young's modulus Z direction, shear modulus YZ, and shear modulus XZ). Then, in each direction a constraint equation of E=2G(1+v) must be maintained. E is Young's modulus, G is shear modulus, and v is Poisson's ratio. Currently, in each direction there are two unknowns in the equation whether it is E and v or G and v. If I was able to extrapolate information on poisson's ratio in each direction I should be able to get the remaining material properties using the equation. Now, this is clearly the difficult part. I thought of running three separate static structurals where there is a fixed support and compressive force in each axis. Then, making estimated guesses on the poisson's ratios until I found relatively similar results in the honeycomb cell model and the solid body model. If I was to find a combination of poisson's ratios that satisify all three static structurals, I'd be quite confident with the material properties.
However, before I was able to test this theory, I received an issue with modeling the honeycomb cell model itself. As you may have seen in the file I sent, the modal results were unrealistic. I would appreciate any advice or suggestions on how I could model the solid body model accurately without the honeycomb cell model. From the link in the beginning of this message, are you able to see any missteps I took in this modelling approach?
November 29, 2022 at 12:06 pmpeteroznewmanSubscriber
The solid orthotropic materials give a perfectly good mode shape. Did you turn down the Deformation Scale Factor as shown in the image below? Displacements in a Modal analysis are arbitrary, so you might need to turn this down until you get normal looking displacements.
I suppressed the Fixed Support and requested 9 modes instead of 6 modes. With no support, the first six modes are rigid body modes with no deformation and a practically zero frequency. Modes 7 and up are bending modes.
Orthotropic materials are not required to satisfy the Isotropic relation you showed above.
Orthotropic materials are required to satisify this requirement:
The way I would get the 2t thickness is by doing a split body on the two planes that bound a column of 2t faces so that you have separate bodies for the t thickness walls from the 2t thickness walls. In SpaceClaim, on the Workbench tab, use the Share button to have the mesh connect the separate bodies together.
I don’t think it is necessary to build a honeycomb cell model to get good orthotropic material constants. What would be more useful is if you can do a 3 or 4 point bending test on a physical honeycomb sandwich sample, or do a shear modulus test on a physical sample of the core.
One way to make the model solve faster is to replace solid elements on the face sheets with shell elements, however, you will have to do some work to align the material properties to a coordinate system. The reason is solid elements automatically align the material direction to the global coordiantes by default. However, shell elements automatically align the material direction to the element coordinates by default, and that can result in the material direction randomly pointing in different directions depending on how the face was meshed. You are almost required to override the default direction and define the material direction you want.
November 29, 2022 at 10:29 pmMani BhumaSubscriber
Thank you for the detailed response. I have a lot more clarity now!
Is there a reason you chose a very specific deformation scale factor of 7.9288e-3?
I reran my simulation to adjust the deformation scale and take a closer look at the results, but I found my results were quite different from yours. Did you change any of the material properties? I noticed you also used a couple 'Sweep Method' features in your mesh. Does this just increase the speed of the mesh, or does this significantly change the formation of the mesh? I attached images of my results below.
As you can see from the images, it seems the face sheets are not experiencing any deformation and the bending modes' frequencies are quite different from your results.
How would I use the results from a 3 point bending test to estimate good material constants in modal analysis?
November 30, 2022 at 12:37 ampeteroznewmanSubscriber
Ansys Mechanical offers you an automatically calculated Display Scale Factor which is called Auto, Auto 0.5x and Auto 2x. I choose one of those that looked good. If the Scale Factor is too large, the display looks bad, which is what you were initially observing.
I couldn’t solve your mesh because it exceeded the Student License Limits, so I had to use larger elements. However, I still wanted 4 linear elements through the facesheet thickness. I solved using Linear element order. I don’t know why your deformed shapes are so different to mine. There must be a significant difference because my mode 7 frequency was 13.7 kHz while yours was 392 Hz so your and my materials have a dramatically different stiffness or density or both. Check what material was assigned to each layer.
If you have a physical test of a 3-point bending test, you will know the stiffness of the sample in lbf per inch of center deformation. You can build a model of the same test and tune the shear modulus until the model gives the same stiffness.
- You must be logged in to reply to this topic.
Boost Ansys Fluent Simulations with AWS
Computational Fluid Dynamics (CFD) helps engineers design products in which the flow of fluid components is a significant challenge. These different use cases often require large complex models to solve on a traditional workstation. Click here to join this event to learn how to leverage Ansys Fluids on the cloud, thanks to Ansys Gateway powered by AWS.
Earth Rescue – An Ansys Online Series
The climate crisis is here. But so is the human ingenuity to fight it. Earth Rescue reveals what visionary companies are doing today to engineer radical new ideas in the fight against climate change. Click here to watch the first episode.
Subscribe to the Ansys Blog to get great new content about the power of simulation delivered right to your email on a weekly basis. With content from Ansys experts, partners and customers you will learn about product development advances, thought leadership and trends and tips to better use Ansys tools. Sign up here.
- Solver Pivot Warning in Beam Element Model
- Saving & sharing of Working project files in .wbpz format
- Understanding Force Convergence Solution Output
- An Unknown error occurred during solution. Check the Solver Output…..
- What is the difference between bonded contact region and fixed joint
- User manual
- The solver engine was unable to converge on a solution for the nonlinear problem as constrained.
- whether have the difference between using contact and target bodies
- material damping and modal analysis
- Colors and Mesh Display
© 2023 Copyright ANSYS, Inc. All rights reserved.