Search found 238 matches
- 30 Nov 2023, 12:59
- Forum: General
- Topic: 1D beam modelling
- Replies: 38
- Views: 3600
Re: 1D beam modelling
I made a small change to the devel branch so that mass-proportional damping (Rayleigh Damping Alpha) for translational motions can now be given for the beam solver, whereas giving Rayleigh Damping Beta = ... causes a warning message.
- 29 Nov 2023, 18:18
- Forum: General
- Topic: 1D beam modelling
- Replies: 38
- Views: 3600
Re: 1D beam modelling
The beam solver doesn't yet understand the Rayleigh damping parameters, but adding the Rayleigh damping would be an easy modification
- 29 Nov 2023, 15:27
- Forum: General
- Topic: 1D beam modelling
- Replies: 38
- Views: 3600
Re: 1D beam modelling
Both these cases are free vibrations and were supposed to vibrate forever with a magnitude of the initial displacement as there was no damping First, it might be good to mention that while this is expected for an analytic solution, computational time integration methods often have numerical damping...
- 22 Nov 2023, 13:16
- Forum: General
- Topic: 1D beam modelling
- Replies: 38
- Views: 3600
Re: 1D beam modelling
As another test, I now checked that if the initial deflection is such that the displacement component U 2 depends linearly on the coordinate X in the absense of external loads, U 2 at the free end starts to oscillate like this end-u2.png Probably a refined mesh would be needed, so that all eigenmode...
- 22 Nov 2023, 11:20
- Forum: General
- Topic: 1D beam modelling
- Replies: 38
- Views: 3600
Re: 1D beam modelling
There seems to be a problem that the beam solver has not automatically set the order of time derivatives. Therefore one has needed an explicit command Time Derivative Order = 2 in the sif file. I however updated the code as https://github.com/ElmerCSC/elmerfem/commit/7fe6d9b916d272db36f02019b202aa2c...
- 10 Nov 2023, 16:01
- Forum: ElmerPost
- Topic: Compute or store the CRS matrices
- Replies: 11
- Views: 2679
Re: Compute or store the CRS matrices
I would like to learn how to extract the CRS matrices You might try to use the ready subroutine SaveLinearSystem contained in the file SolverUtils.F90. For example, you could try to add the following lines to the solver section linear system save = true linear system save slot = assembly so that th...
- 11 Sep 2023, 11:05
- Forum: Software development
- Topic: Edge Element Matrices / Microwave Examples
- Replies: 2
- Views: 988
Re: Edge Element Matrices / Microwave Examples
A rather simple example of coding a solver which uses H(curl)-conforming finite elements could be https://github.com/ElmerCSC/elmerfem/blob/devel/fem/tests/SecondOrderEdgeElement2D_BCs/EdgeFEInterpolationTest.F90 In a special case (Material Param = Real 0.0) it assembles just the mass matrix. Howeve...
- 02 Sep 2023, 11:21
- Forum: General
- Topic: Boundary conditions [Solved]
- Replies: 6
- Views: 827
Re: Boundary conditions
in general? This may depend on the type of a boundary condition. Prescribing the value of the primary variable of a solver (a Dirichlet-type constraint) should generally work for 1-dimensional boundaries (provided the Elmer mesh file mesh.boundary contains information about the 1-dimensional bounda...
- 16 Aug 2023, 11:10
- Forum: ElmerSolver
- Topic: Eddy Current Power Unit
- Replies: 4
- Views: 479
Re: Eddy Current Power Unit
It is the result of integrating the energy density over volume and thus expressed in watts (W).
-- Mika
-- Mika
- 05 Jul 2023, 12:07
- Forum: General
- Topic: Simulation of Metallic Sphere in Uniform Magnetic Field
- Replies: 27
- Views: 1998
Re: Simulation of Metallic Sphere in Uniform Magnetic Field
I'm wondering if we can create a desired magnetic field, given the formula for B (Bx, By, Bz) for all x, y, z, using only the boundary conditions A mathematical feature of the vector potential formulation is that a vector potential producing the desired B as B = curl A is not unique (if A is a suit...