FreeCad FEM exploitation des résultats
Forum rules
Be nice to others! Read the FreeCAD code of conduct!
Be nice to others! Read the FreeCAD code of conduct!
Re: FreeCad FEM exploitation des résultats
Merci Alex pour tous ces points complémentaires, mais je ne suis pas doué en informatique ...
Je n'arrive même pas à inclure une image ou un fichier dans ce débat comme tu le fais si bien
Pour info j'ai 30 plots mais passer par la modif du fichier.inp ne me cause pas de problème
Ce qui me pose problème sont la fiabilité des résultats des forces équivalentes aux déplacements nuls.
Je pense devoir créer un nouveau sujet sur ce point.
Mais déjà si je pouvais inclure une image ou fichier cela serait plus lisible ...
Cordialement
46Cpi
Je n'arrive même pas à inclure une image ou un fichier dans ce débat comme tu le fais si bien
Pour info j'ai 30 plots mais passer par la modif du fichier.inp ne me cause pas de problème
Ce qui me pose problème sont la fiabilité des résultats des forces équivalentes aux déplacements nuls.
Je pense devoir créer un nouveau sujet sur ce point.
Mais déjà si je pouvais inclure une image ou fichier cela serait plus lisible ...
Cordialement
46Cpi
Re: FreeCad FEM exploitation des résultats
je ne suis pas spécialement doué non plus, j'en sais juste assez pour comprendre que je ne sais pas grand chose
regardes l'onglet "pièces jointes" en bas à gauche lorsque tu édites ton messageJe n'arrive même pas à inclure une image ou un fichier dans ce débat
ci-dessous un fichier de miroir ep.10mm D400 sur 16 appuis, pour exemple (pas tout à fait réaliste concernant les appuis, notamment le central fixe)Pour info j'ai 30 plots mais passer par la modif du fichier.inp ne me cause pas de problème
Le poids du miroir: P=mG=2²xPixep=(2x2x3.14x0.1x2.52)x9.81= 31N à comparer à la somme des réactions en Z du fichier .dat. J'ai comparé vite fait et le total a l'air proche.
De mon point de vu chaque appui reprend le poids de la zone de matière qu'il supporte, au prorata des déformations environnantes (merci la FEM).Ce qui me pose problème sont la fiabilité des résultats des forces équivalentes aux déplacements nuls.
Code: Select all
** outputs --> dat file
** reaction forces for Constraint fixed
*NODE PRINT, NSET=ConstraintFixed_centrale, TOTALS=ONLY
RF
*NODE PRINT, NSET=ConstraintDisplacement_midd1, TOTALS=ONLY
RF
*NODE PRINT, NSET=ConstraintDisplacement_midd2, TOTALS=ONLY
RF
*NODE PRINT, NSET=ConstraintDisplacement_midd3, TOTALS=ONLY
RF
*NODE PRINT, NSET=ConstraintDisplacement_midd4, TOTALS=ONLY
RF
*NODE PRINT, NSET=ConstraintDisplacement_midd5, TOTALS=ONLY
RF
*NODE PRINT, NSET=ConstraintDisplacement_ext1, TOTALS=ONLY
RF
*NODE PRINT, NSET=ConstraintDisplacement_ext2, TOTALS=ONLY
RF
*NODE PRINT, NSET=ConstraintDisplacement_ext3, TOTALS=ONLY
RF
*NODE PRINT, NSET=ConstraintDisplacement_ext4, TOTALS=ONLY
RF
*NODE PRINT, NSET=ConstraintDisplacement_ext5, TOTALS=ONLY
RF
*NODE PRINT, NSET=ConstraintDisplacement_ext6, TOTALS=ONLY
RF
*NODE PRINT, NSET=ConstraintDisplacement_ext7, TOTALS=ONLY
RF
RF
*NODE PRINT, NSET=ConstraintDisplacement_ext8, TOTALS=ONLY
RF
*NODE PRINT, NSET=ConstraintDisplacement_ext9, TOTALS=ONLY
RF
*NODE PRINT, NSET=ConstraintDisplacement_ext10, TOTALS=ONLY
RF
OS: Ubuntu 20.04 LTS (XFCE/xubuntu)
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.19.21125 (Git)
Build type: Unknown
Branch: master
Hash: 40600a55c2fe71ff589be677f6e427ccc937d003
Python version: 3.8.2
Qt version: 5.12.8
Coin version: 4.0.0
OCC version: 7.3.0
Locale: English/United States (en_US)
Re: FreeCad FEM exploitation des résultats
Alex, ton exemple ressemble beaucoup à ce que j'essaye de faire avec en plus une vérification c'est à dire remplacer les forces otenues par des contraintes de force et voir si la pièce retrouve la déformée initiale, baser sur le principe action/réaction de la mécanique.
Voici un test en 3 étapes pour exemple, une simple plaque de 100x100mmx1mm en acier.
La plaque de 1mm est chanfreinée avec des chanfreins de 0,45mm ce qui donne des arêtes vericales de 0,1mm, donc des zones d’appuis très petites qui ne devrait pas influencer la déformation de la plaque : 1- Test Z0 initial
La plaque est portée par contraintes de déplacement nul en ses 4 coins
- coin1 : X= 0 Y = 0 Z = 0
- coin2 : Y = 0 Z = 0
- coin3 : Z = 0
- coin4 : Z = 0
Tout semble normal la plaque se déforme en son centre avec une flèche d’environ 10,16 microns
Et cette simulation donne les forces correspondantes FZ (les FX et FY étant logiquement très faibles)
- coin1 : FZ = 0,1922729
- coin2 : FZ = 0,1895192
- coin3 : FZ = 0,18849
- coin4 : FZ = 0,1887194
On observe de petites variations, la moyenne est de 0,18975 N, et la somme de ces forces correspond assez bien au poids de la pièce. 2- Test Z1
Les contraintes de déplacement nul, seulement en Z, sont simplement remplacées par une contrainte de force fournis par les valeurs de l’essai Z0 aux 3 coins (de coin 2 à 4).
Il ne reste donc les contraintes de déplacement nul pour le coin1 en X, Y et Z (comme en Z0) et pour le coin2 en Y
Là, la pièce veut s’envoler et le flèche avoisine les 51 microns (voir fichiers joints).
Voici un test en 3 étapes pour exemple, une simple plaque de 100x100mmx1mm en acier.
La plaque de 1mm est chanfreinée avec des chanfreins de 0,45mm ce qui donne des arêtes vericales de 0,1mm, donc des zones d’appuis très petites qui ne devrait pas influencer la déformation de la plaque : 1- Test Z0 initial
La plaque est portée par contraintes de déplacement nul en ses 4 coins
- coin1 : X= 0 Y = 0 Z = 0
- coin2 : Y = 0 Z = 0
- coin3 : Z = 0
- coin4 : Z = 0
Tout semble normal la plaque se déforme en son centre avec une flèche d’environ 10,16 microns
Et cette simulation donne les forces correspondantes FZ (les FX et FY étant logiquement très faibles)
- coin1 : FZ = 0,1922729
- coin2 : FZ = 0,1895192
- coin3 : FZ = 0,18849
- coin4 : FZ = 0,1887194
On observe de petites variations, la moyenne est de 0,18975 N, et la somme de ces forces correspond assez bien au poids de la pièce. 2- Test Z1
Les contraintes de déplacement nul, seulement en Z, sont simplement remplacées par une contrainte de force fournis par les valeurs de l’essai Z0 aux 3 coins (de coin 2 à 4).
Il ne reste donc les contraintes de déplacement nul pour le coin1 en X, Y et Z (comme en Z0) et pour le coin2 en Y
Là, la pièce veut s’envoler et le flèche avoisine les 51 microns (voir fichiers joints).
Re: FreeCad FEM exploitation des résultats
3- Test Z2
Réduction des 3 forces identiques à 0,189 N. Cela donne encore une flèche de 19,44 micron (le double de Z0) et montre surtout une énorme dissymétrie de la déformation que je ne peux m’expliquer. De plus la force équivalente au coin1 est passé à 0,192 N.
En résumer les forces calculés correspondant aux déplacements ne semblent pas cohérentes, de plus une force identique dans les coins donne une déformée dissymétrique.
Tout cela me semble illogique, Calculix calculerait faussement ? Ou, où fais-je des erreurs ?
J'ai remarqué dans ton dernier post une petite différence de script dans la modif du fichier.inp, au niveau de :
*NODE PRINT , NSET=Nall
U Cela pourrait venir aussi de la version 0.18 par rapport à la version 0.19 ?
Réduction des 3 forces identiques à 0,189 N. Cela donne encore une flèche de 19,44 micron (le double de Z0) et montre surtout une énorme dissymétrie de la déformation que je ne peux m’expliquer. De plus la force équivalente au coin1 est passé à 0,192 N.
En résumer les forces calculés correspondant aux déplacements ne semblent pas cohérentes, de plus une force identique dans les coins donne une déformée dissymétrique.
Tout cela me semble illogique, Calculix calculerait faussement ? Ou, où fais-je des erreurs ?
J'ai remarqué dans ton dernier post une petite différence de script dans la modif du fichier.inp, au niveau de :
*NODE PRINT , NSET=Nall
U Cela pourrait venir aussi de la version 0.18 par rapport à la version 0.19 ?
Re: FreeCad FEM exploitation des résultats
Test avec des pressions équivalentes aux contraintes de déplacement nul :
sur une simple plaque carrée de 100mm x 100mm x 1mm avec 4 appuis de rayon 2mm placés à 0,5mm des coins
les forces calculées sont respectivement en X Y Z
3.126444E-12 -1.695736E-11 1.925470E-01
2.333391E-15 1.570553E-15 1.904788E-01
-3.107160E-15 1.258867E-15 1.911145E-01
-9.094631E-16 2.409700E-15 1.904979E-01
soit 0,19115955 N en Z avec une flèche de 3,027 microns Pour retrouver une déformée équivalente et symétrique tout en gardant le coin1 fixe doit être de
14,122 kPa pour chaque des 3 autres coins.
14,122 kPa appliquer sur les 3 disques de diamètre 4mm représente une force de 0,17746 N ce qui constitue une belle différence avec la force moyenne calculée.
De plus avec cette pression la flèche passe à 8 microns soit plus du double qu’avec les contraintes de déplacement nul en Z.
Comment expliquer ces différences de résultats ?
sur une simple plaque carrée de 100mm x 100mm x 1mm avec 4 appuis de rayon 2mm placés à 0,5mm des coins
les forces calculées sont respectivement en X Y Z
3.126444E-12 -1.695736E-11 1.925470E-01
2.333391E-15 1.570553E-15 1.904788E-01
-3.107160E-15 1.258867E-15 1.911145E-01
-9.094631E-16 2.409700E-15 1.904979E-01
soit 0,19115955 N en Z avec une flèche de 3,027 microns Pour retrouver une déformée équivalente et symétrique tout en gardant le coin1 fixe doit être de
14,122 kPa pour chaque des 3 autres coins.
14,122 kPa appliquer sur les 3 disques de diamètre 4mm représente une force de 0,17746 N ce qui constitue une belle différence avec la force moyenne calculée.
De plus avec cette pression la flèche passe à 8 microns soit plus du double qu’avec les contraintes de déplacement nul en Z.
Comment expliquer ces différences de résultats ?
Re: FreeCad FEM exploitation des résultats
Bonjour,
J'appliquerai les quatre réactions aux quatre coins et je bloquerai un sommet de la plaque pour qu'elle ne "s'envole" pas (ou mathématiquement de permettre la diagonalisation de la matrice de rigidité).
Malheureusement chez moi, je n'arrive pas à faire refonctionner la FEM.
En appliquant les trois forces, je pense -mais je n'en suis pas sûr- que tu appliques un moment à la plaque qu'il faudrait contrecarrer par un moment de même valeur en sens inverse au point fixe. Je ne pense pas que ce soit possible de le faire.
J'appliquerai les quatre réactions aux quatre coins et je bloquerai un sommet de la plaque pour qu'elle ne "s'envole" pas (ou mathématiquement de permettre la diagonalisation de la matrice de rigidité).
Malheureusement chez moi, je n'arrive pas à faire refonctionner la FEM.
Fada de 3D.
Linux Mint
Linux Mint
Re: FreeCad FEM exploitation des résultats
Sorry for asking in english. I came accross and stumbled on this file. Just a question. Is there a special reason you use multiple displacement objects? They all have the same attributes. You could have used one displacement object and put all faces in there as references.
Re: FreeCad FEM exploitation des résultats
Bernd, The aim is to find the equivalent force in each of the supports and to check if the values given by CalcuiX are correct.
Re: FreeCad FEM exploitation des résultats
Hi Bernd, 46cpi told you his need. In this regard another aim is to ease the match of each reaction force with the relevant support, so I applied one (and only one) face per constraint in this regards.bernd wrote: ↑Fri May 22, 2020 1:31 pmSorry for asking in english. I came accross and stumbled on this file. Just a question. Is there a special reason you use multiple displacement objects? They all have the same attributes. You could have used one displacement object and put all faces in there as references.
Plus the constraints naming (not default naming) in order to locate them easily in .dat file. This last point due to the fact the NODEPRINT card does not use the constraint label but the initial constraint name in .dat file. Something could be improved maybe about that, I'll open a topic in FEM forum maybe.
Anyway, for now FEM WB is great to investigate and make some progress about the 46cpi's project.