The governing Eqs. (11), (14), (16)

The governing Eqs. (11), (14), (16) and (18) are discretized using
the finite difference method and solved simultaneously to obtain
the pressure and the film profiles. The finite difference formulation
for the Reynolds equation is provided in Appendix A. The dimensionless
input parameters are the load W, speed U, material G,
surface roughness σ̅, surface hardness V, and ellipticity parameter
κ. It should be noted that for the sake of minimizing the number of
input parameters, the asperity radius β and the asperity density n
are omitted from the input by the assumption of nβσ¼0.05
[31,33,44] and σ=β ¼ 0:01 [33]. Also, when both surfaces are
rough, the combined roughness of σ ¼
ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
σ2
1 þσ2
2
q
should be used.
For the line-contact EHL problem, the widely-used solution
method is the Newton–Raphson algorithm. The convergence rate
of this method is fast (quadratic), but the procedure requires
decomposition of the Jacobian matrix to solve the equations For
1-D problem of line-contact, applying Newton–Raphson algorithm
was found to be the best choice, while the use of this approach for
2-D problem of point contact is not efficient. In other words, for N
nodes along the line contact domain, the number of elements in
the Jacobian matrix is equal to N2 while for point contact problem
with N N number of nodes, the Jacobian matrix has N4 elements,
which is difficult to handle.
The widely-used method for solving the point-contact EHL
problem is relaxation method which is also utilized in this study.
In fact, Successive Over-relaxation Method (S.O.R) is used to solve
the discretized Reynolds equation in an inner loop, while the
hydrodynamic and asperity pressures are updated by an underrelaxation
factor (starts from 0.01 and decreases gradually
throughout the process) in an outer loop. The other unknown
parameter H00 is also under-relaxed by the load balance in the
outer loop as well. Appendix B shows a flowchart of the numerical
procedure.
To ensure that the fully flooded condition prevails, the minimum
amount of X is set to 4 in all cases simulated. Appendix C
shows the effect of inlet variation on the results for both smooth
aand rough surfaces. Also, while the hydrodynamic pressure
becomes nil at a point close to X¼1, the asperity pressure may
have a value far beyond that point and thus the maximum amount
of X is set to 2. Choosing the limits of Y as 72, the computational
domain is a rectangle with 4rXr2 and 2rYr2. The
hydrodynamic pressure is zero along all boundaries. Also, the
hydrodynamic pressure is set to zero whenever it becomes
negative (commonly known as Reynolds condition [1]). It should
be mentioned that since the problem is symmetric with respect to
the X axis (Y¼0 plane), only half of the problem is modeled and
solved. Therefore, the actual computational domain is 4rXr2
and 0rYr2 with symmetric boundary conditions on Y¼0 plane.
The convergence of the results is ensured by refining the mesh
until the change in the results becomes negligible. In general,
within the range of the operating conditions studied, 30 nodes in
one unit along the X direction and 20 nodes in one unit along the Y
direction (which results in a 180 40 mesh for 4rXr2 and
0rYr2) is proven to be sufficient. Appendix D presents the
results of a mesh refinement study for both the smooth and rough
surfaces. As shown, the error uniformly decreases as the computational
mesh is refined. However, further refinement beyond
180 40 nodes yields less than 1% error at the experience of
much greater computational demand.
The run time of the simulations is dependent on the input
parameters. In general, solution for the higher loads takes longer
because more iteration is needed to converge. Also, the solution
for rough surfaces is generally slower than the smooth solution,
because the asperity pressure must be calculated in each loop. The
typical run time on an Intel i7 3.47 GHz machine for the chosen
domain and grid setup is generally between 20 and 100 minutes
depending on the input.
The final results consist of the hydrodynamic, asperity, and
total pressure distributions as well as the film profile. Based on
these results, the central and the minimum film thickness are
easily determined. Also predicted is the asperity load ratio which
is the ratio of the load carried by the asperities to the total
normal load.
4.

ผลลัพธ์ (ไทย) 1:

Eqs ควบคุม (11), (14), (16) และ (18) discretized ใช้การจำกัดความแตกต่างวิธีการ และแก้ไขพร้อมรับความดันและค่าฟิล์ม กำหนดความแตกต่างแน่นอนสำหรับเรย์โนลด์สที่ ให้สมการในภาคผนวกเอ ที่ dimensionlessอินพุตพารามิเตอร์ไม่โหลด W, G วัสดุ ความเร็ว Uσ̅ความหยาบผิว ความแข็งของผิว V และพารามิเตอร์ ellipticityΚมันควรจะสังเกตที่เพื่อลดจำนวนพารามิเตอร์ βรัศมี asperity และ n ความหนาแน่น asperityละเว้นจากการป้อนข้อมูล โดยสมมติฐานของ nβσ¼0.05[31,33,44] และσ =β¼ 0:01 [33] นอกจากนี้ เมื่อทั้งพื้นผิวจะหยาบ ความหยาบรวมของσ¼ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiΣ21 þσ22qควรใช้ติดต่อบรรทัด EHL ปัญหา การแก้ปัญหาที่ใช้กันอย่างแพร่หลายวิธีคือ อัลกอริทึมนิวตัน – Raphson อัตราการลู่เข้าวิธีนี้จะรวดเร็ว (กำลังสอง), แต่กระบวนการต้องการแยกส่วนประกอบของ Jacobian เมทริกซ์การแก้สมการปัญหา 1 D บรรทัดติดต่อ ใช้อัลกอริทึมนิวตัน – Raphsonพบจะ ดีที่สุด ในขณะที่การใช้วิธีการนี้สำหรับปัญหา 2 D จุดไม่มีประสิทธิภาพ ในคำอื่น ๆ สำหรับ Nโหนตามแนวติดต่อโดเมน หมายเลขขององค์ประกอบในเมตริกซ์ Jacobian มีค่าเท่ากับ N2 สำหรับปัญหาจุดติดต่อมี N N โหน เมตริกซ์ Jacobian มีองค์ประกอบ N4ซึ่งเป็นเรื่องยากในการจัดการวิธีการใช้กันอย่างแพร่หลายสำหรับแก้ EHL จุดติดต่อปัญหาคือ เป็นวิธีที่ยังใช้ประโยชน์ในการศึกษานี้ในความเป็นจริง ใช้วิธีการผ่อนคลายมากกว่าต่อเนื่อง (S.O.R) แก้สมการของเรย์โนลด์ส discretized ในวงการภายใน ในขณะhydrodynamic และ asperity ความดันจะถูกปรับปรุง โดยการ underrelaxationปัจจัย (เริ่มต้นจาก$ 0.01 และลดลงทีละน้อยตลอดการในวงการภายนอก ไม่รู้จักอื่น ๆพารามิเตอร์ H00 ยังอยู่ภายใต้ผ่อนคลาย โดยโหลดยอดดุลในการวนภายนอกเช่นการ แผนผังลำดับงานของตัวเลขแสดงภาคผนวก Bขั้นตอนการเพื่อให้แน่ใจว่าสภาพน้ำท่วมเต็มชัย ต่ำได้ 4 ในทุกกรณีจำลองจำนวน X ภาคผนวก Cแสดงผลของการเปลี่ยนแปลงทางเข้าของผลสำหรับทั้งเรียบพื้นผิวหยาบ aand ในขณะความดัน hydrodynamic ยังกลายเป็น nil ที่จุดใกล้กับ X¼1, asperity ความดันอาจมีค่ามากกว่าจุดนั้น และทำเงินสูงสุดของ X ถูกกำหนดเป็น 2 เลือกขีดจำกัดของ Y เป็น 72 การคำนวณโดเมนคือ สี่เหลี่ยมกับ 4rXr2 และ 2rYr2 ที่ความดัน hydrodynamic เป็นศูนย์ตามขอบเขตทั้งหมด ยังความดัน hydrodynamic ถูกตั้งค่าเป็นศูนย์เมื่อใดก็ ตามจะกลายเป็นลบ (โดยทั่วไปเรียกว่าเรย์โนลด์สเงื่อนไข [1]) ควรจะกล่าวถึงตั้งแต่ปัญหาคือสมมาตรกับ respect เพื่อแกน X (เครื่องบิน Y¼0), จำลองครึ่งเดียวของปัญหา และหาทางแก้ไข ดังนั้น โดเมนคำนวณจริงอยู่ 4rXr2และ 0rYr2 มีเงื่อนไขขอบเขตสมมาตรบนเครื่องบิน Y¼0การบรรจบกันของผลลัพธ์จะมั่นใจ โดยปรับตาข่ายจนกระทั่งการเปลี่ยนแปลงในผลลัพธ์กลายเป็นระยะ ทั่วไปภายในช่วงของเงื่อนไขปฏิบัติศึกษา โหน 30 ในหนึ่งหน่วยในทิศทาง X และโหน 20 ในหนึ่งหน่วยตาม Yทิศทาง (ซึ่งผลลัพธ์ใน 180 40 ตาข่ายสำหรับ 4rXr2 และ0rYr2) คือพิสูจน์ให้เพียงพอ ภาคผนวก D แสดงการผลของการศึกษารีไฟน์เมนท์ตาข่ายสำหรับทั้งเรียบและหยาบพื้นผิว แสดง ข้อผิดพลาดสม่ำเสมอเมื่อเทียบเคียงลดลงเป็นการคำนวณตาข่ายได้รับการปรับปรุง อย่างไรก็ตาม เพิ่มเติมรีไฟน์เมนท์เลย180 โหน 40 อัตราผลตอบแทนน้อยกว่า 1% ข้อผิดพลาดที่ประสบการณ์ของคำนวณอุปสงค์มากมากเวลารันของแบบจำลองจะขึ้นอยู่กับการป้อนข้อมูลพารามิเตอร์ ทั่วไป โซลูชั่นสำหรับโหลดสูงใช้เวลาเนื่องจากเกิดซ้ำเพิ่มเติมจะต้องมาบรรจบกัน ยัง โซลูชั่นสำหรับพื้นผิวหยาบจะโดยทั่วไปจะช้ากว่าโซลูชันเรียบเนื่องจากต้องคำนวณแรงดัน asperity ในแต่ละวงการ ที่โดยทั่วไปเวลาเครื่องมี Intel i7 3.47 GHz สำหรับเลือกการตั้งค่าโดเมนและตารางการคิดระหว่าง 20 ถึง 100 นาทีขึ้นอยู่กับการป้อนข้อมูลผลสุดท้ายประกอบด้วยการ hydrodynamic, asperity และการกระจายความดันรวมตลอดจนประวัติฟิล์ม ขึ้นอยู่กับผลลัพธ์เหล่านี้ ใจกลาง และความหนาของฟิล์มต่ำสุดได้ง่ายขึ้น คาดว่า ยัง มีอัตราส่วน asperity โหลดซึ่งเป็นอัตราส่วนของงานที่ดำเนินการ โดย asperities รวมโหลดปกติ4

ผลลัพธ์ (ไทย) 2:

ปกครอง EQS (11) (14) (16) และ (18) จะ discretized โดยใช้
วิธีการที่แตกต่างกันแน่นอนและแก้ไขไปพร้อม ๆ กันที่จะได้รับ
แรงกดดันและรูปแบบภาพยนตร์ สูตรที่แตกต่างกันแน่นอน
สำหรับสมนาดส์จะอยู่ในภาคผนวกมิติ
ป้อนพารามิเตอร์เป็น W โหลดความเร็ว U, วัสดุ G,
พื้นผิวที่ขรุขระσ̅พื้นผิวแข็ง V และ ellipticity พารามิเตอร์
κ มันควรจะตั้งข้อสังเกตว่าเพื่อประโยชน์ในการลดจำนวนของการ
ป้อนพารามิเตอร์, βรัศมีความรุนแรงและความหนาแน่นของความรุนแรง n
จะถูกตัดออกจากการป้อนข้อมูลโดยข้อสันนิษฐานของnβσ¼0.05
[31,33,44] และσ = β¼ 0 01 [33] นอกจากนี้เมื่อทั้งพื้นผิวที่มีความ
หยาบขรุขระรวมของσ¼
ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
σ2
1 þσ2
2
คิว
ควรใช้.
สำหรับสายติดต่อ EHL ปัญหาการแก้ปัญหากันอย่างแพร่หลาย
วิธีการขั้นตอนวิธีนิวตันราฟสัน อัตราการบรรจบกัน
ของวิธีการนี้เป็นไปอย่างรวดเร็ว (สม) แต่ขั้นตอนต้องมี
การสลายตัวของเมทริกซ์จาโคเบียนการแก้สมการสำหรับ
ปัญหา 1-D ของสายการติดต่อใช้นิวตันราฟสันอัลกอริทึม
ถูกพบว่าเป็นทางเลือกที่ดีที่สุดในขณะที่ การใช้วิธีการนี้สำหรับ
ปัญหา 2 มิติของการติดต่อจุดไม่ได้มีประสิทธิภาพ ในคำอื่น ๆ หา N
โหนดพร้อมโดเมนติดต่อสายจำนวนขององค์ประกอบใน
เมทริกซ์จาโคเบียนเท่ากับ N2 ในขณะที่สำหรับปัญหาจุดติดต่อ
ที่มีจำนวนของโหนด NN, เมทริกซ์จาโคเบียนมีองค์ประกอบ N4,
ซึ่งเป็นเรื่องยากที่จะจัดการ.
วิธีการใช้กันอย่างแพร่หลายสำหรับการแก้จุดติดต่อ EHL
ปัญหาคือวิธีการผ่อนคลายซึ่งถูกนำมาใช้ในการศึกษาครั้งนี้.
ในความเป็นจริงต่อเนื่องกว่าวิธีการผ่อนคลาย (SOR) ถูกนำมาใช้ในการแก้
สมการนาดส์ discretized ในวงด้านในขณะที่
อุทกพลศาสตร์ ความรุนแรงและความกดดันที่มีการปรับปรุงโดย underrelaxation
ปัจจัย (เริ่มต้นที่ 0.01 และลดลงทีละน้อย
ตลอดกระบวนการ) ในวงด้านนอก อื่น ๆ ที่ไม่รู้จัก
พารามิเตอร์ H00 นี้ยังอยู่ภายใต้การผ่อนคลายจากความสมดุลภาระใน
วงด้านนอกเช่นกัน ภาคผนวก B แสดงผังของตัวเลข
ขั้นตอน.
เพื่อให้แน่ใจว่าสภาพน้ำท่วมอย่างเต็มที่ชัยขั้นต่ำ
จำนวนเงินที่มีการตั้งค่า X 4 ในทุกกรณีจำลอง ภาคผนวก C
แสดงให้เห็นถึงผลกระทบของการเปลี่ยนแปลงทางเข้าผลทั้งเรียบ
พื้นผิวหยาบ aand นอกจากนี้ในขณะที่ความดันอุทกพลศาสตร์
จะกลายเป็นศูนย์ที่จุดใกล้กับX¼1ดันความรุนแรงอาจ
มีค่าไกลเกินกว่าจุดนั้นและทำให้จำนวนเงินสูงสุด
ของ X ถูกตั้งค่าให้เลือก 2 ข้อ จำกัด ของ Y 72, การคำนวณ
โดเมน สี่เหลี่ยมผืนผ้าที่มี 4rXr2 และ 2rYr2
ความดันอุทกพลศาสตร์เป็นศูนย์ตามแนวรอยต่อทั้งหมด นอกจากนี้
ความดันอุทกพลศาสตร์ถูกกำหนดเป็นศูนย์เมื่อใดก็ตามที่มันจะกลายเป็น
ลบ (ที่รู้จักกันทั่วไปว่าเป็นสภาพที่นาดส์ [1]) มันควร
จะกล่าวว่าเนื่องจากปัญหาคือส่วนที่เกี่ยวกับ
แกน X (เครื่องบินY¼0) เพียงครึ่งหนึ่งของปัญหาที่เกิดขึ้นเป็นแบบจำลองและ
การแก้ไข ดังนั้นการคำนวณโดเมนจริง 4rXr2
และ 0rYr2 กับเงื่อนไขขอบเขตสมมาตรบนเครื่องบินY¼0.
บรรจบกันของผลคือมั่นใจโดยการปรับแต่งตาข่าย
จนการเปลี่ยนแปลงในผลจะกลายเป็นเล็กน้อย โดยทั่วไป
อยู่ในช่วงของสภาพการดำเนินงานศึกษา 30 โหนดใน
หนึ่งหน่วยตามทิศทาง X และ 20 โหนดในหนึ่งหน่วยพร้อม Y
ทิศทาง (ซึ่งส่งผลให้ 180 40 ตาข่ายสำหรับ 4rXr2 และ
0rYr2) จะพิสูจน์ให้เพียงพอ ภาคผนวก D นำเสนอ
ผลการศึกษาการปรับแต่งตาข่ายทั้งเรียบและขรุขระ
พื้นผิว ดังแสดงข้อผิดพลาดเหมือนกันลดลงเมื่อคำนวณ
ตาข่ายคือการกลั่น อย่างไรก็ตามการปรับแต่งเพิ่มเติมเกิน
180 40 โหนดอัตราผลตอบแทนน้อยกว่า 1% ข้อผิดพลาดที่ประสบการณ์ในการ
คำนวณความต้องการที่มากขึ้น.
เวลาทำงานของแบบจำลองจะขึ้นอยู่กับการป้อนข้อมูล
พารามิเตอร์ โดยทั่วไปวิธีการโหลดสูงใช้เวลานาน
เพราะซ้ำขึ้นเป็นสิ่งจำเป็นที่จะมาบรรจบกัน นอกจากนี้ยังมีวิธีการแก้ปัญหา
สำหรับพื้นผิวที่หยาบกร้านโดยทั่วไปจะช้ากว่าการแก้ปัญหาได้อย่างราบรื่น
เพราะดันความรุนแรงต้องได้รับการคำนวณในแต่ละวง
เวลาทำงานปกติบน Intel i7 เครื่อง 3.47 GHz สำหรับเลือก
โดเมนและการตั้งค่าตารางโดยทั่วไประหว่าง 20 และ 100 นาที
ขึ้นอยู่กับการป้อนข้อมูล.
ผลสุดท้ายประกอบด้วยอุทกพลศาสตร์ความรุนแรงและ
การกระจายความดันรวมเช่นเดียวกับรายละเอียดภาพยนตร์ . ขึ้นอยู่กับ
ผลลัพธ์เหล่านี้ภาคกลางและความหนาของฟิล์มขั้นต่ำ
ที่กำหนดได้อย่างง่ายดาย ที่คาดการณ์ไว้นอกจากนี้ยังมีอัตราการโหลดความรุนแรงซึ่ง
เป็นอัตราส่วนของการโหลดที่ดำเนินการโดย asperities เพื่อรวม
โหลดปกติ.
4

ผลลัพธ์ (ไทย) 3:

การปกครอง EQS . ( 11 ) ( 14 ) ( 16 ) และ ( 18 ) โดยใช้วิธีผลต่างสืบเนื่องแบบจุด

และแก้ไขพร้อมกันเพื่อให้ได้แรงดันและ ฟิล์ม โพรไฟล์ มีการกำหนดให้สมการเรย์โนลด์
ความแตกต่างที่ระบุไว้ในภาคผนวก A
ป้อนพารามิเตอร์ไร้มิติคือความเร็วโหลด w , U , วัสดุพื้นผิวขรุขระ g ,
σ̅ความแข็งผิว ellipticity
V และพารามิเตอร์κ . มันควรจะสังเกตว่าเพื่อลดจำนวนของ
พารามิเตอร์การป้อนข้อมูล , asperity รัศมีบีตาและ asperity ความหนาแน่น n
จะละเว้นจากการป้อนข้อมูลโดยสมมติฐานของβσ¼ 0.05
[ 31,33,44 ] และσ = บีตา¼ 0:01 [ 33 ] นอกจากนี้ เมื่อทั้งสองมีพื้นผิวขรุขระผิว
, รวมσ¼

σ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2
1 þσ 2
2
q

ควรใช้สำหรับสายติดต่อ ehl ปัญหา ใช้กันอย่างแพร่หลายโซลูชั่น
วิธีนิวตัน - ราฟสันขั้นตอนวิธี อัตราการลู่เข้า
วิธีนี้รวดเร็ว ( กำลังสอง ) แต่ขั้นตอนต้อง
การเน่าเปื่อยของจาโคเบียนเมทริกซ์เพื่อแก้สมการ
ปัญหาภายในของเส้นการติดต่อของนิวตันซึ่งใช้ขั้นตอนวิธี
ถูกพบว่าเป็นทางเลือกที่ดีที่สุด ในขณะที่การใช้วิธีการนี้สำหรับ
2 ปัญหาของจุดติดต่อไม่ได้มีประสิทธิภาพ ในคำอื่น ๆสำหรับ n
จุดตามเส้นติดต่อโดเมน จำนวนขององค์ประกอบใน
จาโคเบียนเมทริกซ์เท่ากับ n2 ในขณะที่จุดปัญหาติดต่อ
n n จำนวนของโหนด , จาโคเบียนเมทริกซ์มีองค์ประกอบ N4
, ซึ่งยากที่จะจัดการ
ใช้กันอย่างแพร่หลาย วิธีการแก้จุดติดต่อ ehl
ปัญหาคือการผ่อนคลายวิธีการที่นำมาใช้ในการศึกษานี้
ในความเป็นจริง ต่อเนื่องกว่าวิธีคลายเครียด ( s.o.r ) จะใช้ในการแก้สมการเรย์โนลด์
แบบจุดในลูปด้านใน ในขณะที่แรงกดดันดัชนี asperity
และมีการปรับปรุงโดย underrelaxation
ปัจจัย ( เริ่มต้นที่ 0.01 และค่อยๆลดลง
ตลอดกระบวนการ ) ลูปนอก
ไม่รู้จักอื่น ๆh00 พารามิเตอร์ยังอยู่ภายใต้ผ่อนคลายโดยโหลดสมดุลใน
ห่วงด้านนอกเช่นกัน ภาคผนวก B แสดงผังของกระบวนการเชิงตัวเลข
.
เพื่อให้แน่ใจว่าน้ำท่วมเต็มสภาพ prevails , จำนวนเงินขั้นต่ำ
x คือ ชุดที่ 4 ในทุกกรณีจำลอง ภาคผนวก C
แสดงให้เห็นถึงผลกระทบของการเปลี่ยนแปลงในผลทั้งเรียบและพื้นผิวขรุขระ
. นอกจากนี้ในขณะที่
ความดันอุทกพลศาสตร์กลายเป็นศูนย์ที่จุดใกล้ x ¼ 1 , asperity ความดันอาจ
มีมูลค่าเกินกว่าจุดที่และดังนั้นเงินสูงสุด
x เป็นชุด 2 เลือกขอบเขตของ Y เป็น 72 , โดเมนคอมพิวเตอร์
เป็นรูปสี่เหลี่ยมผืนผ้าและมี 4rxr2 2ryr2 .
ความดันอุทกพลศาสตร์เป็นศูนย์ ด้วยขอบเขต นอกจากนี้ ดัชนี
ความดันตั้งศูนย์เมื่อมันกลายเป็น
ลบ ( หรือที่เรียกกันทั่วไปว่า เรย์โนลด์ส เงื่อนไข [ 1 ] ) มันควรที่จะกล่าวว่าเนื่องจากปัญหา

จะสมมาตรตามแกน X ( Y ¼ 0 เครื่องบิน ) , เพียงครึ่งหนึ่งของปัญหาและแก้ไขแบบ
. ดังนั้นโดเมนคอมพิวเตอร์ที่แท้จริงคือ 4rxr2
0ryr2 กับสมมาตรและเงื่อนไขขอบบนเครื่องบิน¼ 0 Y .
บรรจบผลคือมั่นใจโดยการกลั่นตาข่าย
จนกว่าจะมีการเปลี่ยนแปลงในผลกลายเป็นกระจอก โดยทั่วไป
ภายในช่วงของเงื่อนไขที่ศึกษา 30 จุด
หนึ่งหน่วยตามทิศทาง x และ 20 จุดในหนึ่งหน่วยตาม Y
ทิศทาง ( ซึ่งผลลัพธ์ใน 180 40 ตาข่ายสำหรับ 4rxr2 และ
0ryr2 ) พิสูจน์ให้เพียงพอ ภาคผนวก D เสนอผลศึกษาการปรับแต่งตาข่าย
ทั้งเรียบและพื้นผิวขรุขระ

เป็นแสดงข้อผิดพลาดโดยการลดลงเป็นคอมพิวเตอร์
ตาข่ายคือการกลั่น อย่างไรก็ตาม การปรับแต่งเพิ่มเติมนอกเหนือจาก
180 40 โหนดผลผลิตน้อยกว่า 1% ของข้อผิดพลาดที่พบมากในการคำนวณความต้องการ
.
วิ่งเวลาจำลองขึ้นอยู่กับการป้อนข้อมูล
พารามิเตอร์ ในทั่วไป , โซลูชั่นสำหรับโหลดสูงกว่า ใช้เวลานานกว่า
เพราะซ้ำมากขึ้นเป็นสิ่งจำเป็นที่จะมาบรรจบกัน นอกจากนี้ โซลูชั่น
สำหรับพื้นผิวที่ขรุขระโดยทั่วไปจะช้ากว่าวิธีเรียบ
เพราะ asperity ความดันจะต้องคำนวณในแต่ละวง
ปกติใช้เวลาบน Intel i7 3.47 GHz เครื่องเพื่อเลือกโดเมนและการตั้งค่าตาราง
โดยทั่วไประหว่าง 20 และ 100 นาที

ขึ้นอยู่กับการป้อนข้อมูล ผลลัพธ์สุดท้าย ประกอบด้วย asperity
ดัชนีและการกระจายความดันทั้งหมดรวมทั้งภาพยนตร์ประวัติ โดย
ผลลัพธ์เหล่านี้กลางและความหนาขั้นต่ำ
กําหนดได้อย่างง่ายดาย . ยังคาดการณ์ว่าเป็น asperity โหลดอัตราส่วนซึ่ง
คือสัดส่วนของภาระการดําเนินการโดย asperities ที่จะโหลดปกติรวม
.
4