Say I'd like to compute values within my element, related to the centre of mass(COM) of my domain. Therefore I need to know the coordinates of the COM in advance. This can obviously done in Mathematica as well and then given as an input into the element. However Id like to compute this directly within the element. For this purpose the Area of the domain must be known in advance. Access to the AceFEM data fields is needed, as this computations involve for example the nodal coordinates at element level (numerical integration) and because I need to summarise my results of all elements over the whole domain, as indicated by the code:
SMSExport[w detJ , ed$$["Data", 2], "AddIn" -> True];
Thus the procedure is as follows:
- Call a Task to compute the Area of the domain
- Call a Task to compute the COM of the domain
- Call SKR for the actually computation (part of SMTNewtonIteration)
As yous see I've already developed a procedure which uses the SMSStandardSubroutine["Tasks"], but they have to be called manually after SMTAnalysis. Now I wonder wether it is possible to automatize this procedure, such that it is more prone to errors, as if I forget to call one of these tasks or choose a different order, the results are wrong.
Thanks for your replies. in advance :)