Questions tagged [acegen]

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

Filter by
Sorted by
Tagged with
3 votes
1 answer
81 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 ...
user avatar
1 vote
1 answer
26 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. ...
user avatar
0 votes
0 answers
61 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 ...
user avatar
2 votes
1 answer
85 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 ...
user avatar
0 votes
1 answer
227 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 ...
user avatar
  • 612
0 votes
0 answers
50 views

Deformation in Cylindrical coordinates in AceGen

How is it possible to obtain the deformation gradient in cylindrical coordinates in AceGen? The usual procedure, which I understand, is given by the following: ...
user avatar
0 votes
0 answers
66 views

How to vectors of fibre direction tangent to element surface in AceGen

hope everybody is well and safe. So, I'm creating an AceGen element for the inflation of an anisotropic hyperelastic tube (2-fibre families). I have created the element in a strip formulation and it ...
user avatar
1 vote
1 answer
142 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 ...
user avatar
2 votes
1 answer
82 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 "...
user avatar
3 votes
1 answer
141 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 ...
user avatar
4 votes
1 answer
70 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 ...
user avatar
3 votes
2 answers
96 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 ...
user avatar
5 votes
1 answer
183 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 ...
user avatar
2 votes
2 answers
144 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 ...
user avatar
  • 612
0 votes
1 answer
80 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?
user avatar
5 votes
1 answer
133 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. ...
user avatar
1 vote
2 answers
188 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 ...
user avatar
0 votes
1 answer
95 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 ...
user avatar
10 votes
1 answer
170 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. ...
user avatar
  • 612
2 votes
0 answers
93 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 ...
user avatar
1 vote
0 answers
31 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 ,...
user avatar
3 votes
0 answers
87 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 ...
user avatar
16 votes
2 answers
528 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 ...
user avatar
  • 612
1 vote
0 answers
785 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 ...
user avatar
3 votes
1 answer
72 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 ...
user avatar
  • 51
2 votes
0 answers
44 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 ...
user avatar
  • 21
2 votes
2 answers
191 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,...
user avatar
  • 23
7 votes
1 answer
258 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. ...
user avatar
10 votes
1 answer
739 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 ...
user avatar
  • 1,159
2 votes
1 answer
102 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 ...
user avatar
1 vote
0 answers
60 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 ...
user avatar
  • 1,159
2 votes
1 answer
64 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 ...
user avatar
  • 444
1 vote
1 answer
81 views

Quick access to the numerical integration coordinates and weights?

How can one simply access the numerical integration coordinates and weights for a given integration code (e.g. 23) without generating the code itself? Is there some command for it? Something similar ...
user avatar
3 votes
1 answer
56 views

AceFEM Solver Pauses/Freezes

I have been running many simulations using AceGen/AceFEM and I usually leave them to run while I work on something else. However, every so often I notice that the solver seems to pause/freeze on a ...
user avatar
2 votes
1 answer
184 views

How to change boundary conditions during AceFEM analysis?

I'm trying to prestress a sample by a compressive strain of 10 percent and then stretch it by 30 percent. The way I intend to do this is apply the prestress in one step and then have the code run and ...
user avatar
1 vote
0 answers
38 views

Error in SMTRestart (Mathematica 12.0 + Ace 6.804)

I've recently upgraded Mathematica to version 12.0 and I'm using Ace 6.804 Windows (12 Jul 16). Now I have problems with command SMTRestart[]. I tested this using build-in example: AceGen -> Help -> ...
user avatar
  • 43
0 votes
1 answer
52 views

Can I convert my SMTShowMesh Object into a mathematica graphics object in AceFem?

My question is related to a string of questions Ive asked in the last few days. I'm trying to superimpose a vectorplot on a3D plot that I've produced on AceFem using the SMTShowMesh Command. I'm ...
user avatar
2 votes
1 answer
102 views

Plotting fibre direction at each node

I was wondering if there was a way to plot the direction vector of a fibre family at each node in AceGen/AceFem? Below is an example of exactly what I want to do but I do not know how to recreate it. ...
user avatar
4 votes
2 answers
148 views

Derivatives of function with 20 variables with respect to each variable

I want to find the derivative of a function with 20 variables and arrange it in a 20 by 1 vector. Si[a,b,c,d,.....20 variables] is a function. Is there a simple command to find the first derivative ...
user avatar
0 votes
0 answers
42 views

AceFem keeps crashing immediately after plotting

I hope everybody is having a great day/night. I have generated an element on AceGen in which now I am running an analysis on through AceFem. Upon the evaluation of the AceFem notebook, the plot runs ...
user avatar
1 vote
1 answer
71 views

Creating single AceGen element using the same constitutive model defined for 2 fibre families

I'm currently working with a single constitutive model that describes 2 collagen fibre families that have different element properties. The strain energy density function is W, where there would be 2 ...
user avatar
1 vote
1 answer
63 views

How to use SMSHookeToLame(E, v) when E is a tensor? (AceFem/Acegen Material properties)

I'm using a constitive model for an FEA analysis that is outlined as follows: In which Young's Modulus is an structure strain invariant defined by the following tensor: And I need to use the ...
user avatar
1 vote
0 answers
105 views

Newton iteration procedure without multiplier in AceFem [closed]

Is there any way in AceFem to use newton iterative scheme without using multiplier λ or time step t? For example, if one is trying to calculate solution to some nonlinear equation no multiplier is ...
user avatar
  • 11
4 votes
1 answer
98 views

AceGen: Equation solving fails (because of missing pivoting?)

I am looking to solve some nonlinear equations in a code generated by AceGen. To this end, i implemented a Newton-Raphson scheme making use of the SMSLUFactor[] and SMSLUSolve[] methods. This ...
user avatar
  • 43
3 votes
2 answers
231 views

Double shell wall contact with AceGen

I have a shell with double walls and some space between (Picture 1). The 2 walls are connected with three edges. If I apply force to outer wall and if the force is big enough, outer shell wall will ...
user avatar
  • 93
1 vote
1 answer
134 views

How to print the residual and tangent matrix? [closed]

Is there any command in AceGen to print the residual and tangent matrix in a simple and robust way? EDIT 2: It would be best to have the whole stiffness matrix and residual vector of the problem in ...
user avatar
0 votes
0 answers
74 views

SMSCorrect internal: Incorrect structure-2

Does anyone know what does the error: SMSCorrect internal: Incorrect structure-2. refer to? I mean, in general. Specifically, I get it from the following code. ...
user avatar
2 votes
1 answer
93 views

AceGen incremental values of degrees of freedom

I need to operate with incremental values of my rotational degrees of freedom in my FE code. But I wish this incremental rotations to be dependent on the actual rotations, that I use as degrees of ...
user avatar
  • 444
2 votes
0 answers
59 views

AceGen map result for new variable

I am using AceGen default material SED3H1DFHYH1DOgdenU1 and for this material I want to create resultant map which will show me map for the energy but only if it is above some treshold. I am using ...
user avatar
  • 33
1 vote
1 answer
108 views

SMTNodeData in BatchModeSteering

Does anyone know wether it is at all possible to read/modify node data or element data within the batch mode steering procedure of AceFEM? I'm looking for something like: ...
user avatar