Dear All,<div><br></div><div>I am PhD student of Northeastern University, specializing in computational material science. My current work is to determine the equilibrium position of dislocation in the metal bar under tensile load. This is achieved by minimizing interaction energy of dislocation. The optimization model, for a single degree of freedom, is </div>
<div><br></div><div>find: x</div><div>min: W = - ln ( x^2 + a^2) - b*x</div><div>st: lb &lt; x &lt; ub</div><div><br></div><div>where optimization variable, x, is the position of dislocation; objective function, W, is the interaction energy; a and b are constant, and lb and ub are the upper and lower bound. The attached figure shows the x versus W. The boundary is also considered. </div>
<div><br></div><div>The initial value x0 is the equilibrium position in last loading step, which is pretty close to the optimal x*</div><div><br></div><div>The ipopt works perfectly well when the number of degree of freedom is low, like 50. However, when the no. of optimization variable, x, is greater than 100, ipopt takes many search around x*, up to 500 times. The gradient is analytically available for this problem and the minimum is always there. But I cannot figure out how this happens. </div>
<div><br></div><div>If you have any idea to this situation. Please feel free to contact me so that I can provide more details.</div><div><br></div><div>Thanks</div><div><br></div><div>Cheers</div><div><br></div><div>Run Zhu</div>
<div><br></div><div><br></div><div>-- <br><div>Run Zhu, PhD student</div>
<div>Department of Mechanical and Industrial Engineering</div>
<div>Northeastern University</div>
<div>Boston, MA 02115</div>
</div>