Questions tagged [acegen]
AceGen is a Mathematica package for automatic code generation, focused on finite element method subroutines.
122 questions
1 vote
0 answers
44 views
How to set up line search in AceFEM for plasticity with non-associative asymmetric yield?
I would need some help when introducing line-search inside Newton-Raphson method inside AceGEN. I'm working with a custom plasticity model implemented in AceGEN/AceFEM using standard Tangent and ...
1 vote
0 answers
158 views
Pure bending boundary conditions in AceGEN/AceFEM
I have a column beam (x = 10, y = 1, z = 1) and would like to prescribe pure bending conditions around Y axis. How can I achieve that in AceGEN/AceFEM? Edit: by pure bending I mean applying equal and ...
2 votes
1 answer
171 views
Negative Jacobian for small finite element in AceFEM
I want to make simulation for 1 hexahedral finite element of dimensions 1mm x 1mm x 1mm. I need to use SI units (metres) as below in the code ...
1 vote
1 answer
111 views
SMTShowMesh with additional graphics in SMTRestart session
Some time ago I posted a question concerning additional graphics for SMTShowMesh[] SMTShowMesh with additional graphics Now I want to include additional graphics in SMTRestart[] session for ...
5 votes
1 answer
362 views
AceFEM-generated mesh versus Mathematica mesh: How to extrude 2D mesh to create a 3D mesh in AceFEM
The packages MeshTools (@Pinti) and FEMAddons seem to feature similar functions. I sometimes had strange (i.e. inconsistent results when meshing would work/not work without changing anything). Which ...
2 votes
1 answer
161 views
How to enforce DOFs to have the same values in AceFEM?
I want to make simulation of shearing of one cubic finite element as in the code below. How can I constrain all horizontal displacements (in X direction) of nodes "X"==L to have the same ...
2 votes
1 answer
141 views
AceGen: Assemble a user specific global array (similar to "Residual")
I have an issue which I know is definitely solvable but I don't know how to implement this in AceGen. Basically, I want to assemble a matrix $\mathbf{L} \in \mathbb{R}^{n_{dof} \times 6}$ for a global ...
1 vote
0 answers
112 views
SMSDo Error when try to use counter variable i as vector index in loop
I try to evaluate the following code, but I get "Error during numerical evaluation of expression with random numbers." Expression contains parts that canot be evaluated. Messages: ...
0 votes
2 answers
201 views
Why are my AceGen results not matching with the commercial FE software?
I conducted a linear FE analysis of a 3D cantilever beam using AceGen/Acefem and Ansys. I find that results are a little bit off from the ansys results. For e.g., the deflection reported at the ...
0 votes
1 answer
131 views
Question regarding initial guesses for Newton's method in AceGen/ AceFEM
I am trying to solve a system of equations with the command SMTNewtonIteration[]. This method needs initial guesses for the degrees of freedom I am looking for. Is ...
0 votes
1 answer
135 views
Split module "Tangent and Residual" to two modules "Tangent" and "Residual"
I am struggling with an issue and need advice. In general: My objective is to implement a Quasi-Newton (QN) method using a small strain elasto-plastic element from the main library and a user defined ...
0 votes
1 answer
117 views
Meaning of divergence message?
I am using the command SMTNewtonIteration[] to solve a boundary value problem in AceFEM and get the following error message: What does this message mean exactly? Am ...
0 votes
1 answer
107 views
AceGen/AceFEM: Divergence in a term, that should be equal to zero
I have a problem concerning the order of the evaluation of the AceFEM parameters. My AceGen Element contains, among many others, an expression of the form B (aa - bb/(cc Exp[dd w^2] + ff)) with B, aa, ...
1 vote
1 answer
137 views
Integrate doesn't give a result for non-linear functions
I use Acegen for Finitie Element formulation that gives out residuals and tangents by taking in the field values as inputs. For this in my residual formulation I need to integrate a nonlinear term ...
2 votes
1 answer
169 views
How to extract gauss coordinates in AceGen/ AceFEM?
I am trying to implement a 2D turbulence problem with the model of Spalart-Allmaras. Therefore, I need to extract the global gauss coordinates to define parameter 'd' which stands for the distance to ...