Questions tagged [acegen]
AceGen is a Mathematica package for automatic code generation, focused on finite element method subroutines.
76
questions
2
votes
0answers
26 views
Deformation following loads in 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. ...
10
votes
1answer
292 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 ...
2
votes
0answers
35 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 ...
1
vote
0answers
30 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 ...
2
votes
1answer
39 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 ...
1
vote
1answer
53 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 ...
3
votes
1answer
41 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 ...
2
votes
1answer
97 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 ...
0
votes
0answers
60 views
slow calculations for Mathematica 12.0 + Ace 7.006
I've upgraded Matehmatica to version 12.0 and Ace to version 7.006. I'm surprised that now calculations last much longer than using older version of Ace 6.804. For example using build-in example:
...
0
votes
0answers
48 views
How to perform adaptive mesh refinement on a 3D element consisting of 4 superimposed meshes on AceFem?
I'm trying to perform adaptive mesh refinement on a body that is composed of 4 separate meshes joined together. The body is composed of the following individual meshes:
...
1
vote
0answers
34 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 -> ...
0
votes
1answer
39 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 ...
2
votes
1answer
77 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.
...
4
votes
2answers
125 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 ...
0
votes
0answers
39 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 ...
1
vote
1answer
41 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 ...
1
vote
1answer
48 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 ...
1
vote
0answers
79 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 ...
4
votes
1answer
75 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 ...
3
votes
2answers
193 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 ...
1
vote
1answer
73 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 ...
0
votes
0answers
65 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.
...
2
votes
1answer
78 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 ...
2
votes
0answers
51 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 ...
1
vote
1answer
67 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:
...
3
votes
2answers
151 views
AceGen/AceFEM Explicit Meshless
Hello fellow Ace Users.
Currently I'm working on a project to implement Peridynamics.
This is a discretization technique in the fashion of a meshless particle method.
AceGen/AceFEM provides the ...
5
votes
2answers
188 views
Contact with rigid obstacle in AceFEM
Does anyone of you know how to specify the simulation when a deformable body is affected by contact with some rigid body? I think it should be possible, but I am unable to find it in the manual.
3
votes
1answer
84 views
PostProcessing the SMTNodeData values in the deformed configuration
I have been struggling with this problem for a long time. I know that this question is not specifically related to AceGen or AceFEM, but I would like to know if anyone had a similar experience in ...
0
votes
0answers
61 views
Derive a tangent using the AceGen package
I want to derive the algorithmic tangent of a material model, but all I have given is the potential and the flow function. Are there any exmaple codes or tutorials for the implementation with AceGen? ...
2
votes
1answer
72 views
Extract the integer part from a real valued variable in AceGen
I would like to extract the integer part of a real variable. To be more specific, given a vector alpha I would like to extract the integer part of ...
1
vote
0answers
56 views
Implementing different integration rules for different fields
I am back with another issue that I couldn't solve, although trying different ways. The problem is in two-dimensional plane strain setup with hyperelasticity. The only degrees of freedom that I have ...
4
votes
1answer
209 views
Running Abaqus with AceGen
I am an AceGen/AceFEM user, but I think codes generated with AceGen could be more useful if people accustomed to Abaqus could also use them. The AceGen help claims that generation of Abaqus codes is ...
0
votes
0answers
60 views
AceGen list parameter
I am working with AceGen and I try to define my own materiaÅ and in this material I have a parameter K which is calculated at every incement. My problem is that ...
2
votes
1answer
76 views
Using My own colorfunction in SMTShowMesh
When plotting the data using SMTShowMesh in AceFEM package I would like to use my own color function not the one that SMTShowMesh uses by default. For this to be achieved, I created my own color ...
0
votes
0answers
65 views
Time dependent parameter in AceGen
The AceGen is a third party Mathematica package (website) used for the automatic derivation of formulae needed in numerical procedures.
Is there possible to define in AceGen code, parameter which ...
2
votes
0answers
104 views
Using ArcLength method for my simulation
Me again with another problem :-)
Since I couldn't get the convergence in one of my simulations due to some problems which I couldn't identify, I decided to try the Arc length procedure in order to ...
1
vote
0answers
54 views
Is there any limitations for SMTPut?
Recently, I figured out that there is a problem in separating the simulation and post-processing phases in my analysis. Since my simulation is quite big so I decided to separate the simulation phase ...
4
votes
1answer
164 views
What could be the reason for for the residual error in post-iteration to high…?
Recently I started having this error in my simulations in AceFEM. For me it is quite strange because it happens in a particular step when the convergence has ...
4
votes
4answers
218 views
Prescribe the rotation as the boundary condition and the plane section constraint in AceFEM
I am going to analyze a rectangular beam under bending moment in AceFEM. To prescribe the bending I would like to prescribe the rotation angle of the section at the point $z=L_z$ (see the AceFEM file ...
3
votes
1answer
144 views
Using a fast solver to speed up the simulation in AceFEM
My question this time is more general and I would like to know if you have some suggestions.
I am going to simulate a very big 3D problem in AceFEM (109368 elements, 1419432 nodes, and 3594150 ...
6
votes
1answer
396 views
AceGen/AceFEM Material Subroutine
I'm back with an unusual issue.
I`d like to generate a finite element for AceFEM, that calls an external (AceGen generated) subroutine, for the material.
Upfront, as the code is a little broad, I ...
1
vote
1answer
137 views
Define two elements, e.g. H2 and S1 elements, in SMSTemplate
Suppose that I have developed a code for a 3D heat problem (without considering heat convection of the surface) with H2 elements. Now I am going to enhance the code and add the heat convection of the ...
3
votes
1answer
131 views
How to generate such a complicated 3D object in AceFEM?
I would like to simulate a specimen with the following shape in AceFEM. I know that I have to create each part separately and call SMTMesh multiple times but I really don't know how to create the ...
2
votes
2answers
47 views
User Task evaluated after SMTAnalysis[]
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 ...
0
votes
0answers
73 views
AceGen: Transfer the code for H1 element to a code for H2S
I have recently generated a code based on 8-node Hexahedral (tri-linear) element that gives the response of a hyperelastic material. However, now I have found out that with 20-node serendipity (bi-...
3
votes
2answers
105 views
AceFEM: Continue despite error
I am testing different FE formulations on a series of simple benchmark tests. Not all elements converge to the final solution for all the examples, so it occurs that a following (or a similar) message ...
2
votes
0answers
69 views
AceFEM: Matrix condensation and elimination of local unknowns
Lately I have been working with non-linear mixed hybrid elements. The basic article that I took for my research is one by T.H.H. Pian dn K. Sumihara: Rational approach for assumed stress finite ...
2
votes
1answer
50 views
AceFEM: Session Identifier
I wonder wether there is a way to identify an AceFEM session. Some variable with a Hash or the startup time of the current session, which is reset after an SMTInputData[] command.
Thanks in advance !
1
vote
0answers
52 views
AceFEM: Manipulation of SMTData[“MatrixGlobal”] after building but before solving
I want to manipulate the global stiffness matrix before I solve the system of equations. Unfortunately, the command SMTNewtonIteration[] builds up the global matrix and the global vector and solve the ...
2
votes
0answers
85 views
AceFEM: Tie - Option Issue
I try to use the "Tie" Option of the SMTAnalysis[] function.
In general, I want to have some additional Nodes (two kinds of them) at certain (not all) coordinates tied and at some not.
What I have is ...
