How to Apply Boundary Conditions in Finite Element Analysis Using C++?

AI Thread Summary
The discussion focuses on applying boundary conditions in a finite element analysis program using C++. The user successfully assembled the global stiffness matrix and load vector and is now seeking guidance on implementing boundary conditions through the elimination approach. They have found a book that outlines an algorithm for transient heat conduction but are unclear about the step involving the correction of the nodal heat vector for boundary conditions after each time iteration. The user emphasizes their need for hints or guidance on this specific aspect, indicating a desire to learn more about finite element methods. Clarification on boundary condition corrections during time-stepping is the primary request for assistance.
dinaharchery
Messages
22
Reaction score
0
Hello all,

I am trying to create a simple finite element program in C++ and have been able to build and assemble the GLOBAL stiffness matrix as well as the load vector. I want to apply the boundary conditions using the elimination approach, can anyone tell me how to do this?

I am sorry for the simplicity of the problem but I have been unable to find any information on this particular issue. Thanks in advance.
 
Engineering news on Phys.org
Well, I figured out the boundary "elimination" approach to the GLOBAL stiffness matrix as well as the vector Load/Heat-Rate:smile: (I finally got a book on FEM)

The book has a section on transient heat conduction using a LUMPED capacity matrix - complete with algorithm. However, there is one step in the algorithm that I am not sure about.

The last step before iterating the time-loop again is to "correct" the just determined nodal heat vector (i.e. solution vector at current time) for boundary conditions. Can anyone please explain what is meant by correcting the solution vector for boundary conditions, I thought the bounds were set-up before entering the Time-Loop?

BTW, the algorithm is explicit forward difference for LUMPED capacity matrix. I really need some help if anyone has a hint?

Thank you
 
Anyone?

I just need some kind of hint in the direction to go for the application of boundary corrections at each time step. I am not an engineer but am trying to learn finite element for a program.

Once again, thank you.
 

Thread 'Race car suspension Class'

Posted June 2024 - 15 years after starting this class. I have learned a whole lot. To get to the short course on making your stock car, late model, hobby stock E-mod handle, look at the index below. Read all posts on Roll Center, Jacking effect and Why does car drive straight to the wall when I gas it? Also read You really have two race cars. This will cover 90% of problems you have. Simply put, the car pushes going in and is loose coming out. You do not have enuff downforce on the right...
View full post »
Thread 'Physics of Stretch: What pressure does a band apply on a cylinder?'

Thread 'Physics of Stretch: What pressure does a band apply on a cylinder?'

Scenario 1 (figure 1) A continuous loop of elastic material is stretched around two metal bars. The top bar is attached to a load cell that reads force. The lower bar can be moved downwards to stretch the elastic material. The lower bar is moved downwards until the two bars are 1190mm apart, stretching the elastic material. The bars are 5mm thick, so the total internal loop length is 1200mm (1190mm + 5mm + 5mm). At this level of stretch, the load cell reads 45N tensile force. Key numbers...
View full post »

Thread 'Cantilever help'

I'm trying to decide what size and type of galvanized steel I need for 2 cantilever extensions. The cantilever is 5 ft. The space between the two cantilever arms is a 17 ft Gap the center 7 ft of the 17 ft Gap we'll need to Bear approximately 17,000 lb spread evenly from the front of the cantilever to the back of the cantilever over 5 ft. I will put support beams across these cantilever arms to support the load evenly
View full post »

Similar threads

Abuse of the Finite Element Method
Replies
9
Views
3K
How to apply divergence free (∇.v=0) in nodal finite element method?
Replies
3
Views
1K
B Incorporating boundary conditions in the Finite Element Method (FEM)
Replies
4
Views
3K
Corrugated tube (bellow) under axial load
Replies
2
Views
2K
Writing Finite Element Code for Structural Analysis
Replies
2
Views
4K
Modal mass and kinetic energy in FEM modal analysis
Replies
3
Views
3K
Stiffness matrix in finite element method
Replies
3
Views
10K
I Finding the weak form of an eq. with DG elements (finite elements)
Replies
6
Views
3K
Dynamic Structural Analysis Using Finite Element Method
Replies
7
Views
4K
Boundary conditions in Ansys
Replies
3
Views
3K

Hot Threads

Recent Insights

Back
Top