iMechanica - Comments for "One time calculation in UEL"
https://imechanica.org/node/16831
Comments for "One time calculation in UEL"enHi Brunda,As you know I'm
https://imechanica.org/comment/26135#comment-26135
<a id="comment-26135"></a>
<p><em>In reply to <a href="https://imechanica.org/node/16831">One time calculation in UEL</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>Hi Brunda,<br />As you know I'm using UEL. I have two beam elements in a row (elements is defined in UEL) and every element has two nodes. So, this structure will have three nodes. The node #1 is fixed completely (ENCASTRE) and node #3 is the tip of this structure. There is a load on node#3 (let's say a static one) with a user defined amplitude (defined in UAMP subroutine). In HISTORY OUTPUT part of the input file, I'm asking ABAQUS to give me the displacement for the node #2 (the node in the middle of the structure). I will refer to it as DIS2 after this (n*1 dimension). During UEL calculations I'm getting a big matrix (let's call it AA), n*n dimension. I must multiply AA by DIS2 (AA*DIS2). The point is this, I don't know how I can have access to history output results for the node #2.<br />One of the methods which came to my mind is defining another step. But I'm not sure is that going to work? How can I write results for a step and call them during another step? How can I have access to AA matrix in the later step?<br />I would appreciate if you can help through this out as usual.<br />Best Regards,<br />P.S.: I've tried to make the simplest example by using two elements in a row. What I have is a lot more than two.</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/16831%23comment-form">Log in</a> or <a href="/user/register?destination=node/16831%23comment-form">register</a> to post comments</span></li>
</ul>Thu, 31 Jul 2014 20:06:30 +0000ashkan khalilicomment 26135 at https://imechanica.orgHi Brunda,
https://imechanica.org/comment/26134#comment-26134
<a id="comment-26134"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/26124#comment-26124">Hi ashkan,</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>Hi Brunda,</p>
<p>Thanks for your great comments. They are really helpful.</p>
<p>Best Regards,</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/16831%23comment-form">Log in</a> or <a href="/user/register?destination=node/16831%23comment-form">register</a> to post comments</span></li>
</ul>Thu, 31 Jul 2014 18:51:50 +0000ashkan khalilicomment 26134 at https://imechanica.orgHi ashkan,
https://imechanica.org/comment/26124#comment-26124
<a id="comment-26124"></a>
<p><em>In reply to <a href="https://imechanica.org/node/16831">One time calculation in UEL</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>Hi <a class="username" title="View user profile." href="user/15872">ashkan</a>,</p>
<p>a) you have both matrices, multiply them .. write the code for in your uel file.</p>
<p>b) you can store the values on SVARS. assuming you have three displacements after step 1, then you will need three SVARS.</p>
<p> </p>
<p>as far as I remember you can use SVARS to store the solution dependent variables.</p>
<p> </p>
<p>hope this helps</p>
<p>brunda</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/16831%23comment-form">Log in</a> or <a href="/user/register?destination=node/16831%23comment-form">register</a> to post comments</span></li>
</ul>Mon, 28 Jul 2014 16:56:43 +0000brundacomment 26124 at https://imechanica.orgHi Brunda,
https://imechanica.org/comment/26123#comment-26123
<a id="comment-26123"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/26091#comment-26091">Hi Ashkan</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>Hi Brunda,</p>
<p>Thanks for your response.</p>
<p>I'm using UEL. Displacements that I'm calculating by UEL are in frequency domain, U(f). I want to change them to time domain, U(t). For doing that I must multiply a trnsfer matrix (n*n) by U (n*1). I'm calculating transfer matrix inside of the UEL and I have access to displacement (U). In history output part of ABAQUS, what I have as displacement is U(f), displacement in frequency domain. Is there anyway that I can add a user-defined history output, let's name it UT, and ask ABAQUS to calculate it for me?</p>
<p>Let's say there is no way to do that. Now I have to write the U(f) and transfer matrix to two files and then in another fortran code call them and multiply them and plot them. Writing transfer matrix to a file is not a big deal but, writing displacement, U, is difficult because it's been calculating in every time incerement. When I'm trying to write it to a file, because for ever time incerement it has a new value, it is been over writing on the value of last incerement. But I need the value of every time incerement to be written in the file. How can I do that?</p>
<p>I really appreciate your consideration. Thanks in advance.</p>
<p>Best Regards,</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/16831%23comment-form">Log in</a> or <a href="/user/register?destination=node/16831%23comment-form">register</a> to post comments</span></li>
</ul>Mon, 28 Jul 2014 16:12:28 +0000ashkan khalilicomment 26123 at https://imechanica.orgHi Ashkan
https://imechanica.org/comment/26091#comment-26091
<a id="comment-26091"></a>
<p><em>In reply to <a href="https://imechanica.org/node/16831">One time calculation in UEL</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>Hello Ashkan,</p>
<p>I do not follow your question completely.</p>
<p>Will you be able to access the displacements and transfer matrix in the subroutine? If yes, then you can create some loops and multiply them. </p>
<p>Also, what do you mean by " It should have the same ability, at the results part, as other outputs have.". </p>
<p> </p>
<p>B</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/16831%23comment-form">Log in</a> or <a href="/user/register?destination=node/16831%23comment-form">register</a> to post comments</span></li>
</ul>Sat, 26 Jul 2014 03:54:12 +0000brundacomment 26091 at https://imechanica.orgHi Brunda,
https://imechanica.org/comment/26090#comment-26090
<a id="comment-26090"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/25992#comment-25992">hi akshan</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>Hi Brunda,</p>
<p> </p>
<p>I really appreciate your helpful and valuable comments.</p>
<p> </p>
<p>When I'm running the analysis in ABAQUS I'm getting displacements as output (history output). Let's say at the end of analysis, for displacement in x direction I have an n*1 array (displacements in every n time increment). At the other hand I have an n*n transfer matrix which must be multiplied by n*1 displacement array. Because it's a matrix multiplication, it can't be done in every time increment and it must be done at the end of analysis. I need to define a new output (I have no idea how I can do that), lets name it Transformed-Displacement, and have it on history output. It should have the same ability, at the results part, as other outputs have.</p>
<p> </p>
<p>It would be very kind of you to help me with this.</p>
<p> </p>
<p>Best Regards,</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/16831%23comment-form">Log in</a> or <a href="/user/register?destination=node/16831%23comment-form">register</a> to post comments</span></li>
</ul>Fri, 25 Jul 2014 21:08:08 +0000ashkan khalilicomment 26090 at https://imechanica.org
https://imechanica.org/comment/26086#comment-26086
<a id="comment-26086"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/25992#comment-25992">hi akshan</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p> </p>
<p>Hi Brunda,</p>
<p> </p>
<p>I really appreciate your helpful and valuable comments.</p>
<p> </p>
<p>When I'm running the analysis in ABAQUS I'm getting displacements as output (history output). Let's say at the end of analysis, for displacement in x direction I have an n*1 array (displacements in every n time increment). At the other hand I have an n*n transfer matrix which must be multiplied by n*1 displacement array. Because it's a matrix multiplication, it can't be done in every time increment and it must be done at the end of analysis. I need to define a new output (I have no idea how I can do that), lets name it Transformed-Displacement, and have it on history output. It should have the same ability, at the results part, as other outputs have.</p>
<p> </p>
<p>It would be very kind of you to help me with this.</p>
<p> </p>
<p>Best Regards,</p>
<p> </p>
<p> Normal<br />
0</p>
<p> false<br />
false<br />
false</p>
<p> EN-US<br />
X-NONE<br />
AR-SA</p>
<p> /* Style Definitions */<br />
table.MsoNormalTable<br />
{mso-style-name:"Table Normal";<br />
mso-tstyle-rowband-size:0;<br />
mso-tstyle-colband-size:0;<br />
mso-style-noshow:yes;<br />
mso-style-priority:99;<br />
mso-style-parent:"";<br />
mso-padding-alt:0in 5.4pt 0in 5.4pt;<br />
mso-para-margin-top:0in;<br />
mso-para-margin-right:0in;<br />
mso-para-margin-bottom:10.0pt;<br />
mso-para-margin-left:0in;<br />
line-height:115%;<br />
mso-pagination:widow-orphan;<br />
font-size:11.0pt;<br />
font-family:"Calibri","sans-serif";<br />
mso-ascii-font-family:Calibri;<br />
mso-ascii-theme-font:minor-latin;<br />
mso-hansi-font-family:Calibri;<br />
mso-hansi-theme-font:minor-latin;}</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/16831%23comment-form">Log in</a> or <a href="/user/register?destination=node/16831%23comment-form">register</a> to post comments</span></li>
</ul>Thu, 24 Jul 2014 16:27:32 +0000ashkan khalilicomment 26086 at https://imechanica.orgHi Brunda,
https://imechanica.org/comment/26005#comment-26005
<a id="comment-26005"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/25992#comment-25992">hi akshan</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>Hi Brunda,</p>
<p>Thanks again for your valuable response.</p>
<p>It helped me so much.</p>
<p>If you don't mind I would be in touch with you in the future in case of having any problem.</p>
<p>Best Regards,</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/16831%23comment-form">Log in</a> or <a href="/user/register?destination=node/16831%23comment-form">register</a> to post comments</span></li>
</ul>Thu, 10 Jul 2014 21:47:18 +0000ashkan khalilicomment 26005 at https://imechanica.orghi akshan
https://imechanica.org/comment/25992#comment-25992
<a id="comment-25992"></a>
<p><em>In reply to <a href="https://imechanica.org/node/16831">One time calculation in UEL</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>akshan,</p>
<p>i would do the following.</p>
<p>1) at the begining of increment 1 or at time t=0, you compute the eigen values in your uel. use an if loop to achieve this.</p>
<p>2) at the end of increment 1, use uexternaldb to compute the eigen values and store the values in a variable.</p>
<p>3) from increment 2, use the stored eigen values from step 2.</p>
<p>4) the datatype of the variable would be common. this kind of variable can be accessed across multiple subroutines and further the stored data is not deleted during increments. </p>
<p>5) further in uexternaldb, you can use an if loop to control when to calculate the eigen values.</p>
<p>ex: if (kinc .eq. 1) then</p>
<p>calculate eigen values</p>
<p>end if</p>
<p>using the above loop, the eigen values are calculated only once at the end of increment 1</p>
<p>i hope this helps</p>
<p> </p>
<p>brunda</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/16831%23comment-form">Log in</a> or <a href="/user/register?destination=node/16831%23comment-form">register</a> to post comments</span></li>
</ul>Wed, 02 Jul 2014 16:46:21 +0000brundacomment 25992 at https://imechanica.orgHi Brunda,
https://imechanica.org/comment/25987#comment-25987
<a id="comment-25987"></a>
<p><em>In reply to <a href="https://imechanica.org/comment/25979#comment-25979">One time calculation in UEL</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>Hi Brunda,</p>
<p>Thanks for your response and I'm sorry for my late reply on your comment.</p>
<p>I checked uexternaldb. I couldn't realize that if there is any way to make this subroutine to do some calculations just in the beginning of analysis or not. If there is how can I do that?</p>
<p>You mentioned COMMON TYPE variables. I didn't get what do you mean. would you please let me know more about it?</p>
<p>Let's say I have UEL and UEXTERNALDB subroutines. UEXTENALDB is performing some computations at the beginning of analysis and the results (which is a matrix) is needed in all of the time increments (UEL need them). How can I call the UEXTERNALDB's output inside of the UEL in every time increment?</p>
<p>It would be very kind of you to help with this issue.</p>
<p>Best Regards,</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/16831%23comment-form">Log in</a> or <a href="/user/register?destination=node/16831%23comment-form">register</a> to post comments</span></li>
</ul>Tue, 01 Jul 2014 15:12:27 +0000ashkan khalilicomment 25987 at https://imechanica.orgOne time calculation in UEL
https://imechanica.org/comment/25979#comment-25979
<a id="comment-25979"></a>
<p><em>In reply to <a href="https://imechanica.org/node/16831">One time calculation in UEL</a></em></p>
<div class="field field-name-comment-body field-type-text-long field-label-hidden"><div class="field-items"><div class="field-item even"><p>hi ashkan,</p>
<p> </p>
<p>check out the user subroutine uexternaldb. this subroutine can be used to perfom the computations (sprinc?) at the end of the increment. further you can use common type variable to store the values and recall them during the next increment.</p>
<p> </p>
<p>brunda</p>
</div></div></div><ul class="links inline"><li class="comment_forbidden first last"><span><a href="/user/login?destination=node/16831%23comment-form">Log in</a> or <a href="/user/register?destination=node/16831%23comment-form">register</a> to post comments</span></li>
</ul>Sat, 28 Jun 2014 03:18:26 +0000brundacomment 25979 at https://imechanica.orgError | iMechanica