Search found 778 matches

Go to advanced search

by HarryvL
Sat Feb 16, 2019 8:06 am
Forum: FEM
Topic: fcFEM - FEA from start to finish
Replies: 57
Views: 1431

Re: fcFEM - FEA from start to finish

I plan to provide basic finite element functionality and bolt on features if and when needed for a specific design challenge. As this will be FOSS at some point, others can do the same.
by HarryvL
Fri Feb 15, 2019 11:47 pm
Forum: FEM
Topic: fcFEM - FEA from start to finish
Replies: 57
Views: 1431

Re: fcFEM - FEA from start to finish

Hi Thomas, General/motivation ; there are quite a few practical and interesting engineering problems that are hard to do or even impossible with the current FreeCAD/Calculix combination and I don’t think adding more external solvers or waiting for their developer(s) is the answer. For example: mixed...
by HarryvL
Fri Feb 15, 2019 2:30 pm
Forum: FEM
Topic: fcFEM - FEA from start to finish
Replies: 57
Views: 1431

Re: fcFEM - FEA from start to finish

And it exports well to VTK / ParaView:


Sliced_Cube_1.png
Sliced_Cube_1.png (110.47 KiB) Viewed 118 times
by HarryvL
Fri Feb 15, 2019 12:56 am
Forum: FEM
Topic: fcFEM - FEA from start to finish
Replies: 57
Views: 1431

Re: fcFEM - FEA from start to finish

Got the interface elements working. Identifying the interface elements from FC mesh data and extracting the relevant information was the greatest challenge. The actual finite element implementation is a matter of a few hours. Example of a sliced block with interface normal stiffness high and shear s...
by HarryvL
Wed Feb 13, 2019 9:24 pm
Forum: FEM
Topic: fcFEM - FEA from start to finish
Replies: 57
Views: 1431

Re: fcFEM - FEA from start to finish

I am planning to profile the code and report out on where time is spent. Clearly, solving the matrix equations is the most computationally intensive, but here is where NumPy/SciPy comes the rescue. Including smart use of preconditioning with classical techniques, like Cuthill McKee band optimizatio...
by HarryvL
Tue Feb 12, 2019 9:28 pm
Forum: FEM
Topic: fcFEM - FEA from start to finish
Replies: 57
Views: 1431

Re: fcFEM - FEA from start to finish

CORRECTED TO RUN STAND ALONE: elfaces=[] numelfaces=0 mesh = App.ActiveDocument.getObject("FEMMeshGmsh").FemMesh for obj in App.ActiveDocument.Compound.Links: print("obj: {}\n".format(obj.Name)) bfcount=0 for bfface in obj.Shape.Faces: bfcount+=1 print(">> object_face {}: {}\n".format(bfcount,bfface...
by HarryvL
Tue Feb 12, 2019 9:25 pm
Forum: FEM
Topic: fcFEM - FEA from start to finish
Replies: 57
Views: 1431

Re: fcFEM - FEA from start to finish

Yes, that is actually what the attached example is. However, in general I would need Boolean Fragments for more complex bodies than the solid primitives in part WB?
by HarryvL
Tue Feb 12, 2019 6:13 pm
Forum: FEM
Topic: fcFEM - FEA from start to finish
Replies: 57
Views: 1431

Re: fcFEM - FEA from start to finish

Interface.FCStd elfaces=[] numelfaces=0 for obj in App.ActiveDocument.Compound.Links: print("obj: {}\n".format(obj.Name)) bfcount=0 for bfface in obj.Shape.Faces: bfcount+=1 print(">> object_face {}: {}\n".format(bfcount,bfface)) for elface in mesh.getFacesByFace(bfface): numelfaces+=1 elfaces.appe...
by HarryvL
Tue Feb 12, 2019 6:09 pm
Forum: FEM
Topic: fcFEM - FEA from start to finish
Replies: 57
Views: 1431

Re: fcFEM - FEA from start to finish

Have you tried? There is no method in part wb to check for equal geometry. I did some basic one for FEM in mesh tools. See https://github.com/FreeCAD/FreeCAD/blob/51b32d75b411eaf5e75761be99e41a1b16f85943/src/Mod/Fem/femmesh/meshtools.py#L1422-L1454 It might not work on a compound because of CenterO...
by HarryvL
Tue Feb 12, 2019 5:13 am
Forum: FEM
Topic: fcFEM - FEA from start to finish
Replies: 57
Views: 1431

Re: fcFEM - FEA from start to finish

FreeCAD has some methods isEqual, isPartner, isSame but they all only work on exact the same geomety, but not on equal geometry. Means two boxes could have equal geometry but they are not the same. It is like two women could wear an equal dress (normally they never would, but they could :mrgreen: )...

Go to advanced search