Questions tagged [acegen]
AceGen is a Mathematica package for automatic code generation, focused on finite element method subroutines.
1
vote
0answers
36 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 ...
3
votes
1answer
79 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
42 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
57 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
46 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
84 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
47 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
118 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 ...
3
votes
4answers
149 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 ...
2
votes
1answer
106 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
112 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
97 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 ...
2
votes
1answer
104 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
41 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
55 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-...
2
votes
2answers
75 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
47 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
44 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
40 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
62 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 ...
4
votes
1answer
188 views
How to draw a generic 3D closed shape to be analyzed in AceFEM?
I would like to draw a 3D shape to test my AceGen procedure. I will insert the AceFEM part of my code. In this example, I have performed the analysis on a hexahedron element.
...
2
votes
1answer
70 views
Index out of bound error in AceFEM
I have written an AceGEN code for a 3D element. The degrees of freedom at each point are 4, 3 spatial displacements and one additional degree of freedom that acts like temperature, i.e. for the ...
4
votes
0answers
197 views
AceGen: User element generation for FEAP or ABAQUS to solve coupled problems
Is possible to generate elements (SKR module) to use in FEAP or Abaqus for problems involving additional degrees of freedom?
I'm trying to generate an 3D element to establish the mechanical ...
3
votes
1answer
73 views
SMTShowMesh with additional graphics
I make some simulations in AceGen and AceFEM using H1 finite elements. After computations I want to present results, in particular, deformed (or undeformed) mesh with some map (e.g. stress component) ...
0
votes
0answers
63 views
Does anyone have experience in simulating reinforced concrete member in AceGen and AceFEM?
I am going to start generating an acegen code in which the flexural behaviour of a reinforced concrete beam can be simulated. Since there is no documentary about reinforced concrete members in AceGEN ...
4
votes
1answer
75 views
Assumed units of material parameters in AceFEM
I want to be sure of the values I am obtaining using AceFEM, so wanted to know that the units of elastic modulus that the AceShare assumes is it MPa or N/m^2 and for density of solid is it Kg/m^3 ...
2
votes
0answers
112 views
AceGen: Solving the dynamic equation in the middle of the time step
Note: This question is about AceFEM and AceGen package for automatic code generation for finite element analysis.
If the question is not specific enough, I will try to refrase it, but I am a fairly ...
4
votes
1answer
144 views
Example of dynamic analysis in AceFEM
I am interested to see how to create a simple example of dynamic FEM analysis implemented with AceFEM package. Dynamic analysis means that inertial effects are taken into account. The package ...
0
votes
0answers
62 views
How to change the element connectivity?
After each load step, I check if the stress has reached a critical value. If this is true, I would like to change the elements connectivity as shown in the figure. Is it possible to do it with AceFem ...
5
votes
1answer
183 views
Application of natural boundary conditions
I am trying to apply a time dependent load in AceFEM, for which i have used adaptive time method for analysis, also i need it to be applied at the nodal points from 1 to 15. So i have used natural ...
4
votes
1answer
77 views
Mesh import - Contra-intuitive results of simple AceFEM analysis
I am doing AceFEM analysis with mesh imported from external software. I use ELFEN to create tetrahedral mesh, write its properties to a text file and then read this file into Mathematica, extracting ...
5
votes
1answer
92 views
GUI button with AceFEM functions
I would like to have a GUI with a button that starts AceFEM analysis and then some results appear. I have taken one example from the documentation on boundary condition specification and I wrapped it ...
5
votes
2answers
107 views
How to collect results even if AceFEM analysis diverges?
Sometimes (implicit) finite element analysis with AceFEM diverges and solution procedure is aborted before reaching the final time/load step. I have followed the method shown in documentation for ...
7
votes
1answer
158 views
How to include AceFEM analysis in a package
I would like to wrap my AceFEM analysis in one function and include it in a Mathematica package (.wl). I have followed the instructions on how to create packages from documentation and saved the ...
7
votes
1answer
148 views
AceFEM: Can I change boundary conditions in the middle of analysis session?
I'd like to compute an example of a clamped column that is subjected to mixed types of load. Namely, the load consists of the horizontal force λF0 and the vertical prescribed displacement v0 that ...
3
votes
2answers
202 views
AceFEM: Show mesh with average element stress
How can I get value of stress (or some other field) in each integration point of finite element with AceFEM package? I need this information to show mesh with values of stress averaged over each ...
6
votes
1answer
202 views
AceFEM: Arc-Length Method and Load Element
I've got another issue regarding the arc length method. In the analysis below, I introduce a load using a finite element.
I then try to solve the problem using arc-length procedure.
First issue is: ...
9
votes
1answer
190 views
AceFEM: Efficient calculation of volume of axisymmetric elements
I am looking for the most efficient calculation of volume of axisymmetric elements in finite element analysis with AceFEM package. Ideally I would like to have this implemented in AceGen element code ...
6
votes
1answer
98 views
AceFEM: Manipulating Elements Domain
I'm trying to change the Domain of an Element.
For that purpose usage of
makes sense for me.
Anyway, it applies I've done a syntax error here.
Can anyone tell how to do it right ?
...
9
votes
1answer
295 views
AceFEM: Conversion of triangular to quadrilateral mesh
I am using AceFEM command SMTTriangularToQuad to convert 2D triangular mesh (created by Mathematica) to unstructured quadrilateral mesh. This mesh has internal "...
5
votes
1answer
149 views
AceGen: Deformation Dependent Load Element in Arc Length Procedure
for a structural analysis, I'd like to use an finite element that introduces a deformation dependent load within an arc-length procedure.
The Element gets the load value via an input parameter, which ...
4
votes
2answers
94 views
AceFEM: Why SMTShowMesh doesn't show changes in initial configuration?
I am changing initial positions of nodes after the analysis has been initialized (after the SMTAnalysis[...]) and this procedure works fine, but the problem is that ...
2
votes
1answer
114 views
AceGen: How to use MinGW in Compile[]? [duplicate]
Note: This question is about AceGen package for automatic code generation for finite element analysis.
MinGW compiler is recommended for compiling finite element subroutines generated by AceGen and ...
4
votes
2answers
162 views
AceFEM: How to improve legend in SMTShowMesh?
Note: This question is about AceFEM package for finite element analysis in Mathematica.
I would like to improve styling of color legend produced by default use of ...
6
votes
1answer
148 views
AceFEM: How to resolve a pre constrained DOF?
today I came up with the idea to build an element in AceGen which has a degree of freedom that can be "turned on and off". By this I want to control via an input parameter wether the element should ...
7
votes
1answer
159 views
AceFEM: How to mirror a domain in postprocessing?
Note: This question is about AceFEM package for finite element analysis in Mathematica.
In finite element analyses we usually take advantage of symmetry of geometry and loading to model only a ...
