Ivar KJELBERG
COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)
Please login with a confirmed email address before reporting spam
Posted:
1 decade ago
Jun 21, 2010, 1:35 p.m. EDT
Hi
in 2D or 3D there is no easy way to add a "point" rigid body mass, in beam or truss mode they are included as part of the BCs, apart that you can add a small surface/volume, that you give an artificially high density and stiffness. The only other way is to add it by equation, in the weak form and check that all your boundary conditions considered it, including in time and eigenfrequency mode, so that your global equations remain correct.
Point masses are anyhow slightly cumbersome as they make a easily singularities, you should better distribute them over an edge or an area.
Have fun Comsoling
Ivar
Hi
in 2D or 3D there is no easy way to add a "point" rigid body mass, in beam or truss mode they are included as part of the BCs, apart that you can add a small surface/volume, that you give an artificially high density and stiffness. The only other way is to add it by equation, in the weak form and check that all your boundary conditions considered it, including in time and eigenfrequency mode, so that your global equations remain correct.
Point masses are anyhow slightly cumbersome as they make a easily singularities, you should better distribute them over an edge or an area.
Have fun Comsoling
Ivar
Please login with a confirmed email address before reporting spam
Posted:
1 decade ago
Jun 22, 2010, 5:32 a.m. EDT
Hi,
Many thanks for such fast reply. I have tried to do many things :
1. instead of force excerted by the earth on the mass of 10mg, I have specified in terms of force/area on one of the faces.
2. I have tried force/length option on one of the edge faces.
3. I have tried specifying the initial normal stress component on along the length
4. I have tried specifying the initial shear stress component on along the cross section
5. I have tried to create a composite object - cantilever and a small cube of extremely high density value of such a volume so as to match 10mg load
All these work for static analysis with no problem whatsoever
I have used the following parameters for simulation
density : 2300kg/m3
L = 1000um, H = 50um, W = 10um rectangular cross section
But the problem with eigen frequency analysis is as follows. Irrespective of any force applied, since I am unable to attach the 10mg mass(let us forget about the mass for the time being.), the total mass of the cantilever beam is = density*volume = (2300)*(1000u*50u*10u) = 1.15e-9.
The value of k calculated theoretically is k = 3EI/(L^3) = 59.39 N/m . where I = ((10u)*(50u)^3)/12
Now if I use f = (1/2*pi)sqrt(k/m) i get around 36.16k Hz. But the simulation shows 14.713 kHz. Can you please explain this discrepancy?
Hi,
Many thanks for such fast reply. I have tried to do many things :
1. instead of force excerted by the earth on the mass of 10mg, I have specified in terms of force/area on one of the faces.
2. I have tried force/length option on one of the edge faces.
3. I have tried specifying the initial normal stress component on along the length
4. I have tried specifying the initial shear stress component on along the cross section
5. I have tried to create a composite object - cantilever and a small cube of extremely high density value of such a volume so as to match 10mg load
All these work for static analysis with no problem whatsoever
I have used the following parameters for simulation
density : 2300kg/m3
L = 1000um, H = 50um, W = 10um rectangular cross section
But the problem with eigen frequency analysis is as follows. Irrespective of any force applied, since I am unable to attach the 10mg mass(let us forget about the mass for the time being.), the total mass of the cantilever beam is = density*volume = (2300)*(1000u*50u*10u) = 1.15e-9.
The value of k calculated theoretically is k = 3EI/(L^3) = 59.39 N/m . where I = ((10u)*(50u)^3)/12
Now if I use f = (1/2*pi)sqrt(k/m) i get around 36.16k Hz. But the simulation shows 14.713 kHz. Can you please explain this discrepancy?
Ivar KJELBERG
COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)
Please login with a confirmed email address before reporting spam
Posted:
1 decade ago
Jun 22, 2010, 8:23 a.m. EDT
Hi
Yes as most programmes, eigen mode analysis is done without external forces, and if you define the "mass" through a force it will not be introduced into the equations for eigenfrequency or buckling analysis, you must add the mass as amultiplier to the j*omega term.
To learn how COMSOL does it, use the 2D beam physics and add a mass on a point and go the the equations (Physics - Equation Settings - Point Variables in 3.5, turn on the equations in the Preferences in V4 and click on the node equation icons).
Have fun Comsoling
Ivar
Hi
Yes as most programmes, eigen mode analysis is done without external forces, and if you define the "mass" through a force it will not be introduced into the equations for eigenfrequency or buckling analysis, you must add the mass as amultiplier to the j*omega term.
To learn how COMSOL does it, use the 2D beam physics and add a mass on a point and go the the equations (Physics - Equation Settings - Point Variables in 3.5, turn on the equations in the Preferences in V4 and click on the node equation icons).
Have fun Comsoling
Ivar