


{"id":356897,"date":"2024-03-13T17:59:29","date_gmt":"2024-03-13T17:59:29","guid":{"rendered":"\/forum\/?post_type=topic&#038;p=356897"},"modified":"2024-03-13T17:59:29","modified_gmt":"2024-03-13T17:59:29","slug":"about-the-energy-plot-in-ls-dyna","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/about-the-energy-plot-in-ls-dyna\/","title":{"rendered":"About the energy plot in LS-DYNA"},"content":{"rendered":"<p>Hello,<\/p>\n<p>I&#8217;m using the SPH method in LS-DYNA to simulate the high-speed impact problem (Using a metal ballistic projectile to impact a metal plate). I find that the internal energy is not consistent with different numbers of particles (different meshes) for the plate and the projectile. The results are listed below, which are the coarse mesh (70*70*19=93100 nodes), medium mesh (107*107*28=320572 nodes), and fine mesh (151*151*40=912040 nodes).<\/p>\n<p><img decoding=\"async\" src=\"\/forum\/wp-content\/uploads\/sites\/2\/2024\/03\/13-03-2024-1710350564-mceclip0.png\"><\/p>\n<p><img decoding=\"async\" src=\"\/forum\/wp-content\/uploads\/sites\/2\/2024\/03\/13-03-2024-1710350740-mceclip2.png\"><\/p>\n<p><img decoding=\"async\" src=\"\/forum\/wp-content\/uploads\/sites\/2\/2024\/03\/13-03-2024-1710350662-mceclip1.png\"><\/p>\n<p>All the other setups are kept the same for the three cases. As can be seen, the internal energy is not convergent for these cases. For example, the internal energy at 1e-5 s is 2.8e9 (coarse mesh), 3e9 (medium mesh), and 3.5e9 (fine mesh) respectively. These results are calculated by the R 13.0 solver. I also tried to use R 14.0 solver but the results still show a different trend. Even for some cases, the job will stuck at 33% or 46% and I don&#8217;t know why.<\/p>\n<p><img decoding=\"async\" src=\"https:\/\/attachments.office.net\/owa\/rong.jin%40uky.edu\/service.svc\/s\/GetAttachmentThumbnail?id=AAMkAGUyZGVkMDY5LWUwODMtNDdmMy04NTdmLTI0YTlhMTQ5YWFiNwBGAAAAAACAWHUrmz4KSZLe4xEuOJELBwA7jrtVU0w5SaQRNJN%2BH1DOAAAAAAEJAAA7jrtVU0w5SaQRNJN%2BH1DOAAHtoN1WAAABEgAQADeGVUm%2F9KdMpNMvcqATlyQ%3D&amp;thumbnailType=2&amp;token=eyJhbGciOiJSUzI1NiIsImtpZCI6IkU1RDJGMEY4REE5M0I2NzA5QzQzQTlFOEE2MTQzQzAzRDYyRjlBODAiLCJ0eXAiOiJKV1QiLCJ4NXQiOiI1ZEx3LU5xVHRuQ2NRNm5vcGhROEE5WXZtb0EifQ.eyJvcmlnaW4iOiJodHRwczovL291dGxvb2sub2ZmaWNlLmNvbSIsInVjIjoiNWQ0Y2I1ODY5MGMyNGUxNjg5YzgxMzgxOWY2YTY0MzMiLCJzaWduaW5fc3RhdGUiOiJbXCJpbmtub3dubnR3a1wiXSIsInZlciI6IkV4Y2hhbmdlLkNhbGxiYWNrLlYxIiwiYXBwY3R4c2VuZGVyIjoiT3dhRG93bmxvYWRAMmIzMDUzMGItNjliNi00NDU3LWI4MTgtNDgxY2I1M2Q0MmFlIiwiaXNzcmluZyI6IldXIiwiYXBwY3R4Ijoie1wibXNleGNocHJvdFwiOlwib3dhXCIsXCJwdWlkXCI6XCIxMTUzODAxMTIxODkyMTAyOTM3XCIsXCJzY29wZVwiOlwiT3dhRG93bmxvYWRcIixcIm9pZFwiOlwiODkyZDZhOWYtYmNiYy00MThkLWFhZTEtMTlhMmRkM2E0NmU2XCIsXCJwcmltYXJ5c2lkXCI6XCJTLTEtNS0yMS0zNTkwODQzODE4LTMwMDQ5NjAxNjMtNDIzMzM2NjMwOS0zMzY3MDEyNFwifSIsIm5iZiI6MTcxMDM1MjcxNCwiZXhwIjoxNzEwMzUzMDE0LCJpc3MiOiIwMDAwMDAwMi0wMDAwLTBmZjEtY2UwMC0wMDAwMDAwMDAwMDBAMmIzMDUzMGItNjliNi00NDU3LWI4MTgtNDgxY2I1M2Q0MmFlIiwiYXVkIjoiMDAwMDAwMDItMDAwMC0wZmYxLWNlMDAtMDAwMDAwMDAwMDAwL2F0dGFjaG1lbnRzLm9mZmljZS5uZXRAMmIzMDUzMGItNjliNi00NDU3LWI4MTgtNDgxY2I1M2Q0MmFlIiwiaGFwcCI6Im93YSJ9.by2w_OD3o8WmX2VQhXlep_gA5tVSziW6Xc_XjTNZ0DiasIVQwx_S_t7RhORik9BYLMHGhITt8FYMN6Uq_ye62B8W8i9XEwIUJqHxQKtzJxTBa634lSjkzoP_6SnN7hCkeECDrQiPBX9wB-3CAbql4MTqIZo1qnuipJ1K4NMnlThS0KwQ8Ac2g3DPKdZm5u_CgkQgAfOor1-vRUhJcZqrgtzgLg5LyR59sEgVAxBYvCIednRCchvv6W10f7z5LCGRdGCuukA0YZ3Z91HMueDiT-Ui4Omyvg6LZkaW8Pmq5c6QTM_HvdmsM9llpcpxoLPtL_tKNR66fcz2D7SBugPTAQ&amp;X-OWA-CANARY=X-OWA-CANARY_cookie_is_null_or_empty&amp;owa=outlook.office.com&amp;scriptVer=20240229049.16&amp;clientId=D9F310B6B8724C7FA227DA946AD8500E&amp;animation=true&amp;persistenceId=726052c4-4327-40ac-9f9e-582b7740ed30\"><\/p>\n<p>I would like to know how LS-DYNA calculates the internal energy and how I can verify it. We know that the internal energy can be divided into two parts. One is the strain energy, which can be calculated by particle strain and stress. The other is heat, which can be calculated by the particle temperature change. For strain energy, I plan to output all the stress and strain of all the particles in each timestep, then I can calculate the increment of strain &delta;&epsilon;, then time &sigma;, and the total volume and sum, which is the internal energy. For the heat, I plan to output all the temperature change, then times the total mass and specific heat capacity, which is the total heat generated. However, in some cases, I find that the sum of total strain energy calculated and total heat calculated does not equal to the total internal energy directly plotted in *MATSUM.<\/p>\n<p>&nbsp;<\/p>\n<p>Best,<\/p>\n<p>Rong<\/p>\n","protected":false},"template":"","class_list":["post-356897","topic","type-topic","status-closed","hentry","topic-tag-energy-1","topic-tag-ls-dyna","topic-tag-sph","topic-tag-strain-energy-2"],"aioseo_notices":[],"acf":[],"custom_fields":[{"0":{"_bbp_author_ip":["23.52.43.16"],"_bbp_subscription":["22672","88"],"_btv_view_count":["991"],"_bbp_topic_status":["unanswered"],"_bbp_topic_id":["356897"],"_bbp_forum_id":["27814"],"_bbp_engagement":["88","22672"],"_bbp_voice_count":["2"],"_bbp_reply_count":["1"],"_bbp_last_reply_id":["357843"],"_bbp_last_active_id":["357843"],"_bbp_last_active_time":["2024-03-19 10:07:59"]},"test":"rji229uky-edu"}],"_links":{"self":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/356897","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":2,"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/356897\/revisions"}],"predecessor-version":[{"id":356937,"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/356897\/revisions\/356937"}],"wp:attachment":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/media?parent=356897"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}