Topology Optimization - How to apply load at different nodes

20 次查看（过去 30 天）
Sagar Gupta Satrasala
评论： Sagar Gupta Satrasala 2021-3-14
I am working on a standard code for the topology optimization of a MBB beam as in the figure.The force has to be defined at the upper left corner. The design domain is descretized into finite elements, say 100 elements in the x-direction and 10 elements in the y-direction. Both node numbers and element numbers are numbered column wise from left to right.Now, I want to change the position of code the midpoint of right side end of the beam. How do I do that?

请先登录，再回答此问题。

采纳的回答

Anderson Pereira
Anderson Pereira 2021-1-27
Hi Sagar,
If you are refering to the 99-line code, available at https://www.topopt.mek.dtu.dk/Apps-and-software/A-99-line-topology-optimization-code-written-in-MATLAB, the corner nodes are labeled as
upper left corner = 1
bottom left corner = nely+1
upper right corner = (nelx)*(nely+1)+1
bottom left corner = (nelx+1)*(nely+1)
Finally, considering that nely is even, the middle node of the right edge is
midpoint = (nelx)*(nely+1)+1+nely/2
Best,
Anderson
  4 个评论
显示 2更早的评论
Anderson Pereira
Anderson Pereira 2021-2-3
Hi Sagar,
Lines 37–48 are related to the Optimality criteria based optimizer. Please take a look at the code article:
Sigmund, O. A 99 line topology optimization code written in Matlab. Struct Multidisc Optim 21, 120–127 (2001). https://doi.org/10.1007/s001580050176
Best,
Anderson
Sagar Gupta Satrasala
Sagar Gupta Satrasala 2021-3-14
Hi,
I am trying to calculate stress at each element of the beam in the same code and included the code to calculate that inside the while loop(lines 46-48). But I am getting very high stress values. Can anyone please help me to find the reason for that,

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Get Started with Optimization Toolbox 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by