Multi - control optimization and volume constraint implementation
4 months ago by
I wish to perform topology optimization so that I have more than one control parameter (P1 to P7). I have a function 'volfrac' that defines the volume fraction given these control parameters and so I intend integrate this to prescribe my volume constraint. The problem is the examples for structural optimization on the repository all have one parameter as the control (usually density). How do I define the volume constraint jacobian in a way that will accomodate a vector of control parameters rather than just one control parameter?
def jacobian(self, m):
Viii = volfrac(m,m,m,m,m,m,m)
self.tmpvec.vector()[:] = Viii
dintegral = assemble(derivative(self.tmpvec*dx, ?))
What should be at '?' since my m is a vector ie. m = (Control(P1), Control(P2)........Control(P7))?
Community: FEniCS Project
Please login to add an answer/comment or follow this question.