


{"id":31453,"date":"2018-11-12T10:44:01","date_gmt":"2018-11-12T10:44:01","guid":{"rendered":"\/forum\/forums\/topic\/progressive-damage-analysis-using-continuum-damage-mechanics-cdm\/"},"modified":"2018-11-12T10:44:01","modified_gmt":"2018-11-12T10:44:01","slug":"progressive-damage-analysis-using-continuum-damage-mechanics-cdm","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/progressive-damage-analysis-using-continuum-damage-mechanics-cdm\/","title":{"rendered":"Progressive Damage Analysis using Continuum Damage Mechanics (CDM)"},"content":{"rendered":"<p>Hi everyone,<\/p>\n<p>Currently, in my PhD thesis, I&#8217;m dealing with prediction of burst pressure of composite overwrapped pressure vessels (COPVs). For that I&#8217;ve implemented, Hashin or Puck Failure Criteria for damage initiation law. And for damage&nbsp;evolution &#8220;Material Property Degradation&#8221; is selected. Results are OK for GF overwrapped steel liners, but not satisfactory for CF overwrapped Al liners.&nbsp;<\/p>\n<p>So, I decided that I should go for Continuum Damage Mechanics (CDM) for damage evolution. I&#8217;m aware that is restricted with Hashin failure criteria but Hashin is OK for me now.<\/p>\n<p>I&#8217;ve selected simpler case for experimenting with CDM. [+-45]ns laminate with applied tensile load (specifically ASTM D3518 test). I&#8217;m also interested in Open Hole Tests in similar manner.<\/p>\n<p>I&#8217;ve modelled the [+-45]ns case for non-linear shear behavior of the laminate. I&#8217;ve found the constants needed for CDM. Gft = around 80-130 N\/mm (or kJ\/m2), Gfc = 80-110 N\/mm, Gmt= 0.2-1 N\/mm, Gmc=0.5 &#8211; 1 N\/mm and viscous damping coefficients are&nbsp;from 1e-4 to 0.002 for GFRP or CFRP materials.<\/p>\n<p>I&#8217;ve found several articles about the case and material properties are provided.&nbsp;<\/p>\n<p>But in the end, my analysis doesn&#8217;t converge at all, at least for expected loads or displacements.&nbsp;<\/p>\n<p>I tried several things, lets recap those&#8230;<\/p>\n<p>1. Element types: hex20, hex8, quad4, quad8. Best convergence acquired with hex20 but there is still an element violation at about %10-15 of the desired final displacement\/force value.<\/p>\n<p>2. Mesh density: Mesh was mapped homogeneously both on D3518 specimen and open hole specimen. Started with 3 mm element size and gone down to 0.5 mm max size. It doesn&#8217;t help at all.<\/p>\n<p>3. Boundary conditions. I&#8217;ve used both face or edge BCs using&nbsp;ANSYS mechanical interface and also tried &#8220;direct FE&#8221; option with nodal BCs with symmetry and no symmetry. Result is always the same. I can give more detailed information about my BCs if you want.<\/p>\n<p>4. I&#8217;ve tried some Non-linear diagnostics with newton-rhapson residuals and HDST (high distortion) elements. Generally element violation occurs at near BCs because of Poisson&#8217;s ratio (I guess). But it happens too quickly!<\/p>\n<p>5. Large Deflections ON. I also turned on Line Search and when I apply force rather than displacement, I also use Newton-Raphson Unsymmetric.&nbsp; Stabilization option is also tried with constant but didn&#8217;t go well. But I didn&#8217;t try too much values on that option.<\/p>\n<p>I always get an element violation error. Sometimes only 1, sometimes several elements. The solid\/shell body seems to be distorted around that element. Sometimes that is too much that you cant see the body itself.&nbsp;<\/p>\n<p>From trial-and-errors that I&#8217;ve conducted, I understand that CDM (at least ANSYS Workbench implemented version) may not be suitable for non-linear shear behavior of composite materials. For that, custom material models might be modelled through USERMAT subroutines. I&#8217;ll try that way but it seems very complicated for me, at least for now.<\/p>\n<p>Do you have any&nbsp;suggestions&nbsp;about this problem? I can provide my project, BCs, mat properties, schematics, journal articles etc.<\/p>\n<p>Regards&#8230;<\/p>\n","protected":false},"template":"","class_list":["post-31453","topic","type-topic","status-closed","hentry"],"aioseo_notices":[],"acf":[],"custom_fields":[{"0":{"_bbp_old_topic_id":["3625"],"_bbp_old_topic_author_name_id":["Anonymous"],"_bbp_old_is_topic_anonymous_id":["false"],"_bbp_old_closed_status_id":["publish"],"_bbp_author_ip":[null],"_bbp_old_sticky_status_id":["normal"],"_bbp_likes_count":["0","0","0","0","0"],"_btv_view_count":["1536"],"_bbp_subscription":["237810"],"_bbp_topic_status":["unanswered"],"_bbp_status":["publish"],"_bbp_topic_id":["31453"],"_bbp_forum_id":["27791"],"_bbp_engagement":["240","160292"],"_bbp_voice_count":["2"],"_bbp_reply_count":["2"],"_bbp_last_reply_id":["70167"],"_bbp_last_active_id":["70167"],"_bbp_last_active_time":["2018-11-22 07:47:27"]},"test":"obiforeva"}],"_links":{"self":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/31453","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics"}],"about":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/types\/topic"}],"version-history":[{"count":0,"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/31453\/revisions"}],"wp:attachment":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/media?parent=31453"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}