ใช้อัลกอริทึม Newton-Raphson (NR) แบบเต็มรูปแบบมาตรฐานเพื่อหาคำตอบของปัญหา FEM แบบไม่เชิงเส้น
โดยทั่วไป อัลกอริทึม NR มักไม่ลู่เข้าเมื่อให้แรงกระทำเต็มค่าในขั้นตอนเดียว วิธีที่นิยมใช้ ซึ่งนำมาใช้ที่นี่เช่นกัน คือให้แรงกระทำตามลำดับเป็นหลายส่วนเพิ่ม และใช้ผลลัพธ์จากส่วนเพิ่มของแรงกระทำก่อนหน้าเป็นจุดเริ่มต้นของการแก้ด้วยวิธี Newton ในขั้นถัดไป เพื่อวัตถุประสงค์นี้ จึงมีการพัฒนาอัลกอริทึมการควบคุมแรงกระทำเพิ่มเติมบน Newton-Raphson ในกรณีที่การวนซ้ำของ NR ไม่ลู่เข้า ส่วนเพิ่มของแรงกระทำปัจจุบันจะถูกลดลงเหลือครึ่งหนึ่งของค่าเดิม แล้วทำการวนซ้ำ NR อีกครั้ง
วัตถุประสงค์ประการที่สองของอัลกอริทึมการควบคุมแรงกระทำคือการหาแรงกระทำวิกฤต ซึ่งสอดคล้องกับ “เกณฑ์การหยุด” บางประการ ได้แก่ ความเครียดสูงสุดใน Concrete การเลื่อนไถลสูงสุดในองค์ประกอบแรงยึดเหนี่ยว การเคลื่อนที่สูงสุดในองค์ประกอบยึดเหนี่ยว และความเครียดสูงสุดในเหล็กเสริม แรงกระทำวิกฤตหาได้โดยใช้วิธีแบ่งครึ่งช่วง (bisection method) ในกรณีที่เกินเกณฑ์การหยุดที่จุดใดก็ตามในแบบจำลอง ผลลัพธ์ของส่วนเพิ่มของแรงกระทำล่าสุดจะถูกละทิ้ง และคำนวณส่วนเพิ่มใหม่ที่มีขนาดครึ่งหนึ่งของส่วนเพิ่มก่อนหน้า กระบวนการนี้ทำซ้ำจนกว่าจะพบแรงกระทำวิกฤตภายในค่าความคลาดเคลื่อนที่กำหนด
สำหรับ Concrete เกณฑ์การหยุดถูกกำหนดไว้ที่ความเครียด 5% ในแรงอัด (กล่าวคือ ใหญ่กว่าความเครียดที่เกิดการวิบัติจริงของ Concrete ประมาณหนึ่งอันดับ) และ 7% ในแรงดึง ที่จุดอินทิเกรตขององค์ประกอบเชลล์ ในแรงดึง ค่านี้ถูกกำหนดเพื่อให้ความเครียดจำกัดในเหล็กเสริม ซึ่งโดยปกติอยู่ที่ประมาณ 5% โดยไม่คำนึงถึงการเสริมความแข็งจากแรงดึง ถึงค่าจำกัดก่อน ในแรงอัด ค่านี้ถูกเลือกจากหลายทางเลือกให้เป็นค่าที่ใหญ่พอที่จะมองเห็นผลของการบดอัดเสียหายในผลลัพธ์ได้ แต่เล็กพอที่จะไม่ก่อให้เกิดปัญหาด้านเสถียรภาพเชิงตัวเลขมากเกินไป
สำหรับเหล็กเสริม เกณฑ์การหยุดถูกกำหนดในรูปของความเค้น เนื่องจากแบบจำลองความเค้นที่รอยแตกร้าว เกณฑ์ในแรงดึงจึงสอดคล้องกับกำลังรับแรงดึงของเหล็กเสริมโดยคำนึงถึงสัมประสิทธิ์ความปลอดภัย ใช้ค่าเดียวกันสำหรับเกณฑ์ในแรงอัด
เกณฑ์การหยุดในองค์ประกอบแรงยึดเหนี่ยวและ Spring ยึดเหนี่ยวคือ α·δumax โดยที่ δumax คือการเลื่อนไถลสูงสุดที่ใช้ในการตรวจสอบตามมาตรฐาน และ α = 10
