{"id":160172,"date":"2022-09-26T10:00:48","date_gmt":"2022-09-26T10:00:48","guid":{"rendered":"\/knowledge\/forums\/topic\/discovery-aim-tutorial-analysis-of-cantilever-beam-with-i-cross-section\/"},"modified":"2023-08-16T06:33:38","modified_gmt":"2023-08-16T06:33:38","slug":"discovery-aim-tutorial-analysis-of-cantilever-beam-with-i-cross-section","status":"publish","type":"topic","link":"https:\/\/innovationspace.ansys.com\/knowledge\/forums\/topic\/discovery-aim-tutorial-analysis-of-cantilever-beam-with-i-cross-section\/","title":{"rendered":"Discovery AIM tutorial &#8211; Analysis of cantilever beam with I cross-section"},"content":{"rendered":"<p><strong>This example is taken from <u><a href=\"https:\/\/confluence.cornell.edu\/display\/SIMULATION\/ANSYS+AIM+-+I+Beam\" target=\"_blank\" rel=\"nofollow noopener noreferrer\">Cornell University&#8217;s ANSYS AIM Learning Modules<\/a><\/u><\/strong><\/p>\n<hr \/>\n<nav class=\"toc\">Contents<\/p>\n<ol class=\"toc__section -lev0\">\n<li class=\"toc__item -lev0\">Learning Goals<\/li>\n<li class=\"toc__item -lev0\">Problem Specification<\/li>\n<li class=\"toc__item -lev0\">Pre-Analysis<\/li>\n<li class=\"toc__item -lev0\">Geometry<\/li>\n<li class=\"toc__item -lev0\">Mesh<\/li>\n<li class=\"toc__item -lev0\">Physics Setup<\/li>\n<li class=\"toc__item -lev0\">Results Evaluation<\/li>\n<li class=\"toc__item -lev0\">Verification and Validation<\/li>\n<\/ol>\n<\/nav>\n<hr \/>\n<h4  id=\"LEARNING-GOALS\"><strong>Learning Goals<\/strong><\/h4>\n<p>The purpose of this tutorial is to showcase, in a relatively simple situation, where simple beam theory is no longer as valid as it is in the limit of a long and slender beam geometry. \u00a0In some commercial codes, simple one-dimensional cubic beam elements for bending deflection, do not capture shear deflection when the beam is no longer slender. Alternatively in ANSYS, if shear deflection is accounted for in the 1D element formulation, results for the beam\u2019s tip deflection will not agree with tip deflections predicted\u00a0by simple Euler-Bernoulli beam theory. Again, attempts to capture this effect by using more elements will ultimately fail. Either the necessary physics is not contained in the element formulation or it is and the results are compared to simpler theory. Either way, using more elements captures \u201cno more\u201d of the solution than does a coarser discretization.<\/p>\n<p>This tutorial is meant to highlight where it is relatively straightforward to apply 3D FEA and resolve a correct solution, which contradicts analytical treatment with simple formulae such as bending tip deflection =<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-157375\" src=\"\/knowledge\/wp-content\/uploads\/sites\/4\/2022\/08\/hg301.png\" alt=\" width=\"54\" height=\"46\" \/><\/p>\n<hr \/>\n<h4 content_id=\"problem-specification\" class=\"toc__permalink\" content_id=\"problem-specification\" class=\"toc__permalink\"  id=\"PROBLEM-SPECIFICATION\"><strong>Problem Specification<\/strong><\/h4>\n<p>Consider a fixed-end, aluminum cantilever I-beam point-loaded at its tip as shown in the figure below. \u00a0We will be solving for directional deformations and normal stresses in this tutorial.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-157376\" src=\"\/knowledge\/wp-content\/uploads\/sites\/4\/2022\/08\/hg302.png\" alt=\" width=\"459\" height=\"208\" \/><\/p>\n<p>P = 1000 lbf, L = 24 in, c = 4 in<\/p>\n<p>The cross-sectional area and second moment of inertia are A = 9.45 in^2 \u00a0and I = 112.3 in^4 respectively. \u00a0These values correspond to the I-beam cross-section shown in the next figure along with a fully three-dimensional solid model of the beam for purposes of visualization. \u00a0Use the following material properties for the aluminum beam:<\/p>\n<p>Young&#8217;s Modulus = 1e7 psi, \u00a0Poisson&#8217;s Ratio = 0.33<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-157377\" src=\"\/knowledge\/wp-content\/uploads\/sites\/4\/2022\/08\/hg303.png\" alt=\" width=\"599\" height=\"244\" \/><\/p>\n<hr \/>\n<h4 content_id=\"pre-analysis\" class=\"toc__permalink\" content_id=\"pre-analysis\" class=\"toc__permalink\"  id=\"PRE-ANALYSIS\"><strong>Pre-Analysis<\/strong><\/h4>\n<p>You are told that experiments have been performed that confirm that the tip deflection under the load is approximately<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-157378\" src=\"\/knowledge\/wp-content\/uploads\/sites\/4\/2022\/08\/hg304.png\" alt=\" width=\"410\" height=\"43\" \/><\/p>\n<p>Generally, in long, slender beams, the transverse displacement due to pure bending dominates when compared with the accompanying shear deflection. In relatively short beams, however, displacement contributions from shear can be large enough as to be non-negligible compared with those from pure bending. In the present case, the shear contribution is no longer a &#8220;small percentage&#8221; of that from bending:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-157381\" src=\"\/knowledge\/wp-content\/uploads\/sites\/4\/2022\/08\/hg305.png\" alt=\" width=\"296\" height=\"294\" \/><\/p>\n<p>where the relevant shear area is approximately the area of the beam web only. This is as prescribed in most Mechanics of Materials textbooks because any load applied to the flanges tend to deform the flanges locally and do not contribute to the deflection of the neutral axis. You can see this if you attempt to perform this analysis and distribute the end load over the entire cross sectional area including the flanges.<\/p>\n<p>Assuming shear modulus (G) = 3770 ksi:<br \/>\nBending deformation= 0.0041033 in<br \/>\nShear deformation = 0.0030557 in<\/p>\n<p>Here the bending deformation alone varies enough from the total deformation to question whether the difference might be a discretization error when , in fact, it is due to a difference in physical modeling assumptions. To perform a proper validation of any numerical results post-mortem, one must know what theory is embedded within the finite element model.<\/p>\n<p>When the beam is long and slender, transverse displacement due to bending is the dominant contribution to the tip deflection and the pure bending stress component is the dominant normal stress component acting at the wall fixture:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-157389\" src=\"\/knowledge\/wp-content\/uploads\/sites\/4\/2022\/08\/hg306.png\" alt=\" width=\"284\" height=\"63\" \/><\/p>\n<p><strong>What to Expect<\/strong><\/p>\n<p>Performing the three-dimensional analysis results in tip deflections in agreement with theory including the shear deformation but does not render normal stresses at the wall given by simple beam theory. \u00a0Below are the hand calculations for comparison in the Verification and Validation \u00a0section of the tutorial:<\/p>\n<p>&nbsp;<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-157393\" src=\"\/knowledge\/wp-content\/uploads\/sites\/4\/2022\/08\/hg307.png\" alt=\" width=\"603\" height=\"72\" \/><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-157395\" src=\"\/knowledge\/wp-content\/uploads\/sites\/4\/2022\/08\/hg308.png\" alt=\" width=\"410\" height=\"43\" \/><\/p>\n<hr \/>\n<h4 content_id=\"geometry\" class=\"toc__permalink\" content_id=\"geometry\" class=\"toc__permalink\"  id=\"GEOMETRY\"><strong>Geometry\u00a0<\/strong><\/h4>\n<p><a style=\"color:#1E6DDC;font-weight:bold;text-decoration:none;\"  href=\"https:\/\/ansys13.ansys.com\/KnowledgeArticles\/Discovery\/Analysis-of-cantilever-beam-with-I-cross-section.zip\">Analysis-of-cantilever-beam-with-I-cross-section.zip<\/a><\/p>\n<p>The video below shows the steps necessary to import the provided geometry into Discovery AIM.<\/p>\n<p><iframe loading=\"lazy\" class=\"vidyard_iframe\" src=\"https:\/\/play.vidyard.com\/WoTXBRpMe85ESt22bDUB74.html?v=3.1.1&amp;\" width=\"700\" height=\"400\" frameborder=\"0\" scrolling=\"no\" allowfullscreen=\"allowfullscreen\" data-mce-fragment=\"1\"><\/iframe><\/p>\n<hr \/>\n<h4 content_id=\"mesh\" class=\"toc__permalink\" content_id=\"mesh\" class=\"toc__permalink\"  id=\"MESH\"><strong>M<\/strong><strong>e<\/strong><strong>sh<\/strong><\/h4>\n<p>In this tutorial, we will be using\u00a0<strong>Physics-Aware Meshing<\/strong>. Physics-aware meshing helps automate and simplify your problem setup. With physics-aware meshing, the computational mesh is generated automatically based on the solution fidelity setting and the physics inputs.<\/p>\n<hr \/>\n<h4 content_id=\"physics-setup\" class=\"toc__permalink\" content_id=\"physics-setup\" class=\"toc__permalink\"  id=\"PHYSICS-SETUP\"><strong>Physics Setup<\/strong><\/h4>\n<p>In this video, you will learn how to:<\/p>\n<ul>\n<li>Assign material to the I-beam<\/li>\n<li>Assign support and force to the I-beam<\/li>\n<\/ul>\n<p>Also, in this video you will notice that the load is applied to the central web.\u00a0The reason behind this is to focus on the deflection that occurs across the entire span of the I-Beam.\u00a0 This tutorial isn&#8217;t concerned with the deflections that occur to the top and bottom flanges. The problem focuses on the increasing impact which shear deflection has on shorter, non-slender beams.\u00a0 This shear causes a difference between the analytical Euler-Bernoulli solution and the 3D finite element model.<\/p>\n<p><iframe loading=\"lazy\" class=\"vidyard_iframe\" src=\"https:\/\/play.vidyard.com\/e3UHMqwLtJ5565N3VJEjwP.html?v=3.1.1&amp;\" width=\"700\" height=\"400\" frameborder=\"0\" scrolling=\"no\" allowfullscreen=\"allowfullscreen\" data-mce-fragment=\"1\"><\/iframe><\/p>\n<hr \/>\n<h4 content_id=\"results-evaluation\" class=\"toc__permalink\" content_id=\"results-evaluation\" class=\"toc__permalink\"  id=\"RESULTS-EVALUATION\"><strong>Results Evaluation<\/strong><\/h4>\n<p>In this video, you will learn how to view normal stress and displacement.<\/p>\n<p><iframe loading=\"lazy\" class=\"vidyard_iframe\" src=\"https:\/\/play.vidyard.com\/cKKDfiM27uLg4rT5kLZQwk.html?v=3.1.1&amp;\" width=\"700\" height=\"400\" frameborder=\"0\" scrolling=\"no\" allowfullscreen=\"allowfullscreen\" data-mce-fragment=\"1\"><\/iframe><\/p>\n<hr \/>\n<h4 content_id=\"verification-and-validation\" class=\"toc__permalink\" content_id=\"verification-and-validation\" class=\"toc__permalink\"  id=\"VERIFICATION-AND-VALIDATION\"><strong>Verification and Validation<\/strong><\/h4>\n<p>&#8220;Verification and validation&#8221; can be thought of as a formal process for checking results. Validation consists of assuring oneself that the solution is, in fact, correct. This consists of making sure the discretization error is minimized by performing a convergence study to assure oneself that the results are insensitive to the chosen mesh. This can be done by adding Size Controls on the Mesh Task.<\/p>\n<p>Once this is accomplished, one needs to compare the converged results with experimental data or directly applicable theory. In this case, we have experimental results for the deflection and simple Euler-Bernoulli theory for comparison with the stress results. The deflection results are obtained within the uncertainty of the measured tip deflection:<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-157398\" src=\"\/knowledge\/wp-content\/uploads\/sites\/4\/2022\/08\/hg309.jpg\" alt=\" width=\"403\" height=\"31\" \/><\/p>\n<p><img decoding=\"async\" src=\"https:\/\/s3-us-west-2.amazonaws.com\/media.forumbee.com\/i\/71f13239-e768-4071-94d1-49d70e27db4a\/547.jpg\" width=\"464\" \/><\/p>\n<p>Performing the three-dimensional analysis results in tip deflections in agreement with theory including the shear deformation but\u00a0<strong>does not<\/strong>\u00a0render normal stresses at the wall given by simple beam theory. \u00a0Below are the hand calculations from (one-dimensional) Euler-Bernoulli beam theory and the predictions for the normal stress at the fixed wall obtained from a fully three-dimensional analysis that includes Poisson effects from the dimensions in the cross section. Recall, from calculation in the Pre-Analysis section, the normal stress at the wall due to bending deformation alone is 854.9 psi whereas that predicted by the FEA analysis is more than double at 1816 psi.<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-157399\" src=\"\/knowledge\/wp-content\/uploads\/sites\/4\/2022\/08\/hg311.png\" alt=\" width=\"933\" height=\"578\" \/><\/p>\n<p>The Euler-Bernoulli theory predicts normal stresses in the absence of out=of=plane stresses (which are presumed negligibly small. Here the out-of-plane ZZ component of stress is predicted to be approximately 816 psi, i.e on the order of normal stresses due to bending alone. These substantial out-of-plane stresses result in Poisson effects on the normal wall stress, increasing its value substantially.<\/p>\n","protected":false},"template":"","class_list":["post-160172","topic","type-topic","status-publish","hentry","topic-tag-aim-tutorial","topic-tag-discovery-aim","topic-tag-structures"],"aioseo_notices":[],"acf":[],"custom_fields":[{"0":{"_wp_page_template":["default"],"_bbp_last_active_time":["09-13-2022  20:20:20"],"_bbp_forum_id":["159552"],"_btv_view_count":["8757"],"family":[""],"application_name":[""],"product_version":[""],"_bbp_likes_count":["1"]},"test":"solution"}],"_links":{"self":[{"href":"https:\/\/innovationspace.ansys.com\/knowledge\/wp-json\/wp\/v2\/topics\/160172","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/innovationspace.ansys.com\/knowledge\/wp-json\/wp\/v2\/topics"}],"about":[{"href":"https:\/\/innovationspace.ansys.com\/knowledge\/wp-json\/wp\/v2\/types\/topic"}],"version-history":[{"count":0,"href":"https:\/\/innovationspace.ansys.com\/knowledge\/wp-json\/wp\/v2\/topics\/160172\/revisions"}],"wp:attachment":[{"href":"https:\/\/innovationspace.ansys.com\/knowledge\/wp-json\/wp\/v2\/media?parent=160172"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}