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
14 juin 2011, 16:33 UTC−4
Hi
First of all Point loads, as well as Edge loads in 3D are singularities, hence are NOT good practice, but sometimes you can use the results if you do NOT analyse in detail the local stress concentration.
If you do not have a point I'm not sure how to make one, what I normally do is to use a surface load (edge load in 2D) and distribute the force over the area/edge length i.e. with a Gaussian. This gives you a concentrated load, but no pure singularities, if the Gaussian is larger than a few mesh elements. THe CoG of the Gaussian can then be moved over the surface/edge with a parametric sweep
Hi
First of all Point loads, as well as Edge loads in 3D are singularities, hence are NOT good practice, but sometimes you can use the results if you do NOT analyse in detail the local stress concentration.
If you do not have a point I'm not sure how to make one, what I normally do is to use a surface load (edge load in 2D) and distribute the force over the area/edge length i.e. with a Gaussian. This gives you a concentrated load, but no pure singularities, if the Gaussian is larger than a few mesh elements. THe CoG of the Gaussian can then be moved over the surface/edge with a parametric sweep
Please login with a confirmed email address before reporting spam
Posted:
1 decade ago
27 juin 2011, 10:17 UTC−4
Thank you for your answer.
I tried to do what you said but I don´t know how to use the function (interpolation) in the edge load label. I tried with: 1000[N/m]*int1(x[1/m]) but it don´t works. How should I do?
Other question: Is it posible to get a point graph without defing a node in that point?
Thank you!
Bruno.
Thank you for your answer.
I tried to do what you said but I don´t know how to use the function (interpolation) in the edge load label. I tried with: 1000[N/m]*int1(x[1/m]) but it don´t works. How should I do?
Other question: Is it posible to get a point graph without defing a node in that point?
Thank you!
Bruno.
Please login with a confirmed email address before reporting spam
Posted:
1 decade ago
5 juil. 2011, 01:38 UTC−4
Hi Ivar,
This Gaussian distribution trick to simulate point load is great!
I have however a question:
How can I really control the total load acting on my beam?
A simple example is a beam of meter fixed at both ends with a point load actin in the middle (let’s say 1000N).
I define a normalized gaussien by
FCT= (1/sqrt(2*pi*c^2))*exp(-(x-b)^2/(2*c^2))
X = curvilinear coordinate
b = point where the peak of distribution is located
c = std deviation of the distribution
The integral of FCT = 1. However if the std deviation is small enough, the peak is larger than one. Meaning that the “point load” acting on the beam center have force larger than 1000N (if F =1000*FCT). So the sum of all the “point loads” created by the Gaussian distribution is then substantially higher than the wished value.
I have plotted the reaction force vs c: the curve increase up to maximum (1000N per fixed point!) and then decrease….
I’m totally lost with this outcome. I’m sure the element size has smth to do with this. But my principal problem is still how to define properly my Gaussian in function of the element size to obtain a displacement field similar to the one obtained by COMSOL with a point load.
How to you proceed to circumvent this issue?
I have attached my model so that you can have a closer look.
Thank you for your feedback,
Yannick
Hi Ivar,
This Gaussian distribution trick to simulate point load is great!
I have however a question:
How can I really control the total load acting on my beam?
A simple example is a beam of meter fixed at both ends with a point load actin in the middle (let’s say 1000N).
I define a normalized gaussien by
FCT= (1/sqrt(2*pi*c^2))*exp(-(x-b)^2/(2*c^2))
X = curvilinear coordinate
b = point where the peak of distribution is located
c = std deviation of the distribution
The integral of FCT = 1. However if the std deviation is small enough, the peak is larger than one. Meaning that the “point load” acting on the beam center have force larger than 1000N (if F =1000*FCT). So the sum of all the “point loads” created by the Gaussian distribution is then substantially higher than the wished value.
I have plotted the reaction force vs c: the curve increase up to maximum (1000N per fixed point!) and then decrease….
I’m totally lost with this outcome. I’m sure the element size has smth to do with this. But my principal problem is still how to define properly my Gaussian in function of the element size to obtain a displacement field similar to the one obtained by COMSOL with a point load.
How to you proceed to circumvent this issue?
I have attached my model so that you can have a closer look.
Thank you for your feedback,
Yannick