Questions tagged [acegen]

AceGen is a Mathematica package for automatic code generation, focused on finite element method subroutines.

Filter by
Sorted by
Tagged with
2 votes
1 answer
65 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 ...
Frida's user avatar
  • 45
0 votes
0 answers
59 views

Acegen and wsl missing window decoration

I'm working with Linux on Windows with wslg (Windows Subsystem Linux). Thus, the linux version of Acegen and Mathematica installed in Ubuntu running in a wsl under Windows 11. The two problems I'm ...
rcontinuum's user avatar
0 votes
0 answers
60 views

How to install AceFEM in MacOS?

I downloaded the Evaluation/student version and ran Install.nb according to the official installation instructions^1, it returned some errors as follows. Cannot load installation package. Problematic ...
lumw's user avatar
  • 583
1 vote
0 answers
69 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: ...
rcontinuum's user avatar
0 votes
2 answers
125 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 ...
VarunM's user avatar
  • 101
0 votes
1 answer
79 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 ...
user1234321's user avatar
0 votes
1 answer
77 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 ...
Frida's user avatar
  • 45
0 votes
1 answer
75 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? ...
user1234321's user avatar
0 votes
1 answer
57 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, ...
Julia Sunten's user avatar
1 vote
1 answer
88 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 ...
AceRox's user avatar
  • 13
2 votes
1 answer
80 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 ...
user1234321's user avatar
3 votes
1 answer
73 views

AceGen/ AceFEM: Wrong number of nodes?

I am trying to create an element T1 (linear triangle) to solve the steady, incompressible Navier-Stokes equations. Therefore, my SMSTemplate command reads as follows ...
user1234321's user avatar
1 vote
1 answer
79 views

Numbering of auxiliary symbols in AceGen generated code

I have a general question regarding the numbering of auxiliary symbols in AceGen. Consider this example (taken from the AceGen manual, only the language is switched to ...
isteinbrecher's user avatar
3 votes
1 answer
130 views

AceFEM higher order elements fail basic test

I was trying to perform a simple test, as some results that we obtained from AceFEM were not matching with the results we obtained from other methods. Here, I have a simple biaxial test on a cube of ...
Mike's user avatar
  • 33
0 votes
0 answers
124 views

How to modify global stiffness matrix before solving it by SMTNewtonIteration[]?

I'm trying to use the Model Order Reduction technique to reduce the stiffness matrix in AceFEM/AceGen. But there is a problem with the manipulation of the global stiffness before solving it by ...
Tsinghua_Zhang's user avatar
3 votes
1 answer
101 views

Dynamic allocation of memory within the given array

Let say a have an array of known dimension array={-1, 3, 6, -5, -6, 12} and indexes where the negative numbers appear is not known in advance. I want to create an array for which elements satisfy ...
Maciej's user avatar
  • 33
4 votes
1 answer
141 views

Is it possible to remove an element during computation in AceFEM?

The title says it all. I know there have been a lot of enhancements in the recent versions of AceGen/AceFEM. I would like to know if there is a possibility of removing an element during computation (...
KratosMath's user avatar
  • 1,243
1 vote
1 answer
94 views

Different boundary conditions

I have a problem regarding the definition of multiple boundary conditions in AceFEM. I am dealing with coupled problem simulations (displacement and magnetic potential) and therefore I have to apply ...
Max's user avatar
  • 634
1 vote
2 answers
125 views

Strange visualization for H2S element

I am experiencing a strange behavior in the postprocessing of H2S element. I have provided a MWE below. Reference: H2 element ...
KratosMath's user avatar
  • 1,243
4 votes
1 answer
95 views

Package AceFEM error: cannot be automatically installed

I have a problem to perform the package AceFEM in Mathematica 12.0.0.0 for Mac (software: macOS Monterey Version 12.4 | processor: 1,8 GHz Dual-Core Intel Core i5). I downloaded the newest AceGen and ...
user1234321's user avatar
3 votes
1 answer
114 views

How to stop simulation after a given model variable achieves certain value?

I would like to stop my simulation after a given model variable (say, plastic strain averaged over entire domain) achieves a given threshold value. Do you know how to do it? I'm supplying the MWE. How ...
Karol Frydrych's user avatar
1 vote
1 answer
33 views

Number of input parameters (nargin equivalent) for AceGen

Is there an equivalent to the MATLAB parameter nargin that can be used in AceGen? For example, consider a function that takes in two vectors and an optional scalar. ...
DvanHuyssteen's user avatar
0 votes
0 answers
78 views

Clear AceFEM memory

I am doing multiple AceFEM simulations by invoking the mathematica notebooks using NotebookEvaluate. After a number of simulations everything breaks down, I suppose the problem is too much memory ...
Karol Frydrych's user avatar
2 votes
1 answer
135 views

Calculate ArcTan in AceGen

I am working on a material model SMSStandardModule["Tangent and residual"] in AceGen and want to use command ToSphericalCoordinates for but AceGen is ...
Luka Roblek's user avatar
0 votes
1 answer
313 views

AceGen/ AceFem on new M1 Mac

Apple has updated its processors some time ago and I am curious if AceGen/ AceFem are compatible with this latest generation of the Apple M1 processors? And maybe someone can already report on the ...
Max's user avatar
  • 634
1 vote
1 answer
202 views

AceFEM GlobalTask

I wonder whether there is any way to add data directly to the linear equation system in AceFEM. I´m having external data that I want to superimpose to the linear system before solving it. The targeted ...
Sascha Maassen's user avatar
2 votes
1 answer
92 views

AceFEM: MIEL example aborts for many load steps

I am currently using AceFEM version 7.303 and Mathematica 12.2 but had the same problem with other versions before. My laptop as well as the server I use work with Linux. When I run the example "...
Simon M.'s user avatar
4 votes
1 answer
208 views

3D meshing of cyllindrical domain into hexahedral elements

I need a solution for generating a hexahedral mesh for a cyllindrical domain in Mathematica/AceFEM. Ideally, it should be close to cyllindrical symmetry not very far from the center. The best option ...
Karol Frydrych's user avatar
4 votes
1 answer
116 views

AceGen SMSFreeze[] function

I'm quite new to AceGen and I'm not getting the function of SMSFreeze[]. The manual writes the following: SMSFreeze[symbol,{exp1 ,exp2 ,{exp3…}..},allOptions] create data objects that represent ...
Duc Hoang Nguyen's user avatar
3 votes
2 answers
123 views

AceGen: Keep full symmetric matrix

I am using AceGen to compile a subroutine for MATLAB. The subroutine involves several steps of linear algebra that yields a number of symmetric matrices. However, when using AceGen there seems to be ...
DvanHuyssteen's user avatar
5 votes
1 answer
219 views

AceFEM: Divergence in iterative procedure (Newton-Raphson) for fine meshes

I am trying to model a problem of a nearly incompressible $10~\rm{m} \times 2~\rm{m}$ beam with a uniformly distributed end load. The beam has a Young's modulus of $200~\rm{Pa}$ and a Poisson's ratio ...
DvanHuyssteen's user avatar
2 votes
2 answers
154 views

AceFem Error in SMTShowMesh

I consider for a calculation a sphere embedded within a matrix. To plot the results and especially the behavior around the spherical inclusion I want to "cut" through the matrix. That seems ...
Max's user avatar
  • 634
0 votes
1 answer
150 views

Running AceGen / AceFEM on a cluster

Did anyone here ever sucessfully installed AceGen/AceFEM on a cluster? Does anyone has experience with performing AceFEM simulations on a cluster system? Is it possible at all?
Karol Frydrych's user avatar
5 votes
1 answer
170 views

Parallel/Multi-thread element code generation using AceGen

I am using AceGen to generate source code for finite element simulations in AceFEM. I have elements with a very large number of vertices/degrees of freedom which can take a very long time to compile. ...
DvanHuyssteen's user avatar
1 vote
2 answers
269 views

Print global stiffness matrix in AceFEM

I have a 2D linear elastic simulation in AceFEM and for post-processing purposes I would like to print the global stiffness matrix and global displacement vector. I am able to get the global x and y ...
DvanHuyssteen's user avatar
0 votes
1 answer
124 views

initial condition in finite element method

I want to apply the initial condition below to the linear element c(x; 0) = sin(10πx); if x in [0; 0:1] and c(x; 0) = 0; if x in [0:1; 1]. I tried like ...
user75507's user avatar
10 votes
1 answer
191 views

Meshing an inclusion in an inclusion

I do have FEM mesh related question. I am trying to mesh some inclusion from a png file using ToElementMesh. Thats working quite nicely even for more complex structures. ...
Max's user avatar
  • 634
2 votes
0 answers
131 views

Help with script for running AceGen generated routine in ABAQUS

I am using the below configuration to run a AceGen-generated viscoelastic element in Abaqus. Abaqus version: ABAQUS 6.14-6, OS: Windows 10, Visual studio: 2013, Fortan compiler: Intel Fortran compiler ...
Manu Hegde's user avatar
1 vote
0 answers
36 views

AceGen set up for Linux [closed]

I would like to run Abaqus User Element (UEL) subroutine generated by AceGen. I have Abaqus 6.14 coupled with gfortran compiler in Linux Ubuntu operating system. How do I configure ABAQUSCompile.bat ,...
Manu Hegde's user avatar
3 votes
0 answers
124 views

Coupling AceGen with Abaqus

I am using Abaqus for nonlinear FEA and, I have generated a UEL (User defined Element Subroutine)using AceGen. But the .log file generated by Abaqus shows a syntax error when I try to run a job using ...
Manu Hegde's user avatar
16 votes
2 answers
600 views

3D Inclusion with structured mesh and coarse and arbitrary matrix

I am wondering if there is a simple way to define a structured 3D mesh (inclusion) like e.g. that and surround it with a corse and unstructured Matrix. It should be possible to refine both more or ...
Max's user avatar
  • 634
2 votes
0 answers
1k views

General:: Exp[-717.401] is too small to represent as a normalized machine number; precision may be lost

I am solving pde. For the post processing bvp, there is solution given in output which is the solution from mathematica version 10.0.1. Output from above in mathematica 10.0.1 is, But when the similar ...
user75507's user avatar
3 votes
1 answer
102 views

How to obtain eigenvectors in AceGen

I'm trying to obtain eigenvectors of a 3D plastic stretch tensor $C_p$. There is no SMSEigenvectors function, so I tried with ...
LukaP's user avatar
  • 51
2 votes
0 answers
48 views

AceFEM MacOS not responding

I just upgraded to v7.114 of AceFEM running on MacOS 10.15.5 and Mathematica 12.0.0.0. I was running v7.101 without a problem. AceGEN seems to work but all AceFEM code starts but does nothing. The ...
Andrew's user avatar
  • 21
2 votes
2 answers
235 views

LinkOpen::linke: Specified file is not a MathLink executable. in AceGen

I tried to run AceGen/AceFEM on my Mac(10.15.5). I installed Mathematica (12.1.0.0) successfully and also the latest version of the AcePackage (7.113). The installation of the AcePackage is successful,...
Carina's user avatar
  • 23
7 votes
1 answer
302 views

Deformation following loads in AceFEM/AceGen

I have a simple 2D finite element problem comprising a unit domain that is fully constrained on the left, vertically constrained on the bottom and subject to a uniformly distributed load at the top. ...
DvanHuyssteen's user avatar
11 votes
1 answer
999 views

The performance of AceFEM in Intel Processor VS. AMD Processor

I am going to buy a new powerful desktop that I will use to run finite element simulations using AceFEM. I found out that the new AMD Ryzen 9 3950X obtained more points in the CPU Passmark benchmark ...
KratosMath's user avatar
  • 1,243
2 votes
1 answer
113 views

Application of spatially varying boundary conditions and source term in AceFEM

I am trying to model a simple 2D steady state heat conduction problem using automatically generated Q2 elements. The problem consists of a square domain with prescribed temperatures along the bottom ...
DvanHuyssteen's user avatar
1 vote
0 answers
76 views

Generating the tetrahedral element with nodal quadrature rule

By default, we do not have the possibility of using linear Tetrahedral elements with nodal quadrature in AceGen/AceFEM. But I suppose that there should be a way to construct it through the ...
KratosMath's user avatar
  • 1,243
2 votes
1 answer
78 views

SMSPrint causes error in FE code

I am using Mathematica 12.0, AceGen version 7.101. I have not detected this problem when I was using older versions of both, but I've had no need to use the ...
Marko's user avatar
  • 464