


{"id":283669,"date":"2023-05-12T13:35:51","date_gmt":"2023-05-12T13:35:51","guid":{"rendered":"\/forum\/forums\/topic\/calculate-by-hand-frequency-responses-including-residual-vectors\/"},"modified":"2023-05-12T13:35:51","modified_gmt":"2023-05-12T13:35:51","slug":"calculate-by-hand-frequency-responses-including-residual-vectors","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/calculate-by-hand-frequency-responses-including-residual-vectors\/","title":{"rendered":"Calculate by hand frequency responses including residual vectors"},"content":{"rendered":"<p>Hello,&nbsp;<\/p>\n<p>I performed a modal analysis and exported in a text file X,Y,Z,RX,RY,RZ displacement (in that case, mode shape) of some interesting nodes: here an example<\/p>\n<p><img decoding=\"async\" src=\"\/forum\/wp-content\/uploads\/sites\/2\/2023\/05\/12-05-2023-1683897845-mceclip3.png\"><\/p>\n<p>Line : node (6 nodes here)<\/p>\n<p>column : mode (10 modes)<\/p>\n<p>Then, I performed a harmonic analysis with the MSUP method (mode-superposition) linked with the previous modal analysis, and write the MCF file related to it (HROPT,MSUP, , ,YES).<\/p>\n<p>I know that the idea of the MSUP method is to calculate the displacement like this :<\/p>\n<p><img decoding=\"async\" style=\"display: block;margin-left: auto;margin-right: auto\" src=\"\/forum\/wp-content\/uploads\/sites\/2\/2023\/05\/12-05-2023-1683897433-mceclip0.png\"><\/p>\n<p>That is to say, express the displacement vector in the modal basis formed with the mode shapes calculated during the eigenproblem in the previous modal analysis.<\/p>\n<p>Then, i want to take into account the residual vector, which is the highest modes contribution. When i enable including the residual mode (RESVEC, ON), or the harmonic analysis settings <img decoding=\"async\" src=\"\/forum\/wp-content\/uploads\/sites\/2\/2023\/05\/12-05-2023-1683897672-mceclip1.png\">, the MCF file contains one more mode (last column) that is the residual mode :<\/p>\n<p><img loading=\"lazy\" decoding=\"async\" src=\"\/forum\/wp-content\/uploads\/sites\/2\/2023\/05\/12-05-2023-1683897724-mceclip2.png\" width=\"736\" height=\"408\"><\/p>\n<p>In my modal analysis, I only chose the 10 first modes, and the 11th one refers to the residual mode.<\/p>\n<p>Now, I have all i need to calculate the frequency response.<\/p>\n<p>I guess I just have to make &#8220;sum(phi_i * y_i) + residual vector&#8221; at each frequency excitation to calculate that, but it doesn&#8217;t give me the good result. In fact, i&#8217;m comparing result to ANSYS and it doesnt correspond. Here is the comparison response at one particular node :<\/p>\n<p><img decoding=\"async\" src=\"\/forum\/wp-content\/uploads\/sites\/2\/2023\/05\/12-05-2023-1683898048-mceclip4.png\"><\/p>\n<p>Blue : &#8220;by hand&#8221; calculation in MATLAB with MSUP method, <strong>with <\/strong>residual mode (RM)<\/p>\n<p>Red : &#8220;by hand&#8221; calculation in MATLAB with MSUP method, <strong>without <\/strong>residual mode<\/p>\n<p>Yellow : Ansys frequency response with MSUP method, <strong>without <\/strong>residual mode<\/p>\n<p>Purple : Ansys frequency response with MSUP method, <strong>with <\/strong>residual mode<\/p>\n<p>The yellow and red curves are corresponding, which is normal because i just made by hand the same calculation than ANSYS.<\/p>\n<p>Yet, the red and purple curves does not correspond, and it should. Is the operation &#8220;sum(phi_i * y_i) + residual vector&#8221; wrong ?<\/p>\n<p>To resume : I obtain the good results when i&#8217;m performing the MSUP method whithout RM, but not the good results in MSUP method with RM.<\/p>\n<p>Thanks,<\/p>\n<p>Adrien<\/p>\n","protected":false},"template":"","class_list":["post-283669","topic","type-topic","status-closed","hentry","topic-tag-Modal_Analysis-1","topic-tag-eigenmode","topic-tag-eigensolver-analysis-2","topic-tag-frequency-response","topic-tag-harmonic","topic-tag-harmonicmodalanalysis-1","topic-tag-modal","topic-tag-modal-analysis","topic-tag-mode-superposition","topic-tag-residual"],"aioseo_notices":[],"acf":[],"custom_fields":[{"0":{"_bbp_subscription":["284164","88","20905"],"_bbp_author_ip":["23.39.172.6"]," _bbp_last_reply_id":["0"]," _bbp_likes_count":["0"],"_bbp_likes_count":["1"],"_btv_view_count":["1721"],"_bbp_topic_status":["unanswered"],"_bbp_status":["publish"],"_bbp_topic_id":["283669"],"_bbp_forum_id":["27791"],"_bbp_engagement":["88","20905","284164"],"_bbp_voice_count":["3"],"_bbp_reply_count":["7"],"_bbp_last_reply_id":["284525"],"_bbp_last_active_id":["284525"],"_bbp_last_active_time":["2023-05-17 17:23:13"]},"test":"adrien-vetsafrangroup-com"}],"_links":{"self":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/283669","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\/283669\/revisions"}],"wp:attachment":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/media?parent=283669"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}