Kiểm tra ứng suất tại chân răng khi ăn khớp bằng Catia các Pro check giúp !!!

Mình đang làm bài có liên quan tới tính ứng suất tại chân răng khi ăn khớp.Mình muốn dùng các phần mềm cơ khí hiện để kiểm tra giá trị này.Bước đâu mình tính giá trị ứng suất theo lí thuyết,vào sau đó dùng fần CAE của các fần mềm,catia mình dùng modul Analysis để thức hiện.Nhưng kết quả tính từ fần mềm mình thấy khác nhiều so với tính lý thuyết.Các Pro kiểm tra giúp mình có sai sót ở điểm nào.Thanks!!!
Đây là quá trình mình làm:
Thông số của các bánh răng :
m = 2,5 mm, anpha= 20 độ ,
Z1 = 19 , Z2 = 103 , P = 5.7 kW , n = 960 v/ph ,
X1 = 0, X2 = 0, b=49 mm,
vật liệu thép 45 HB<350.



T=9.55*10^6 *P/n =9.55*10^6 *5.7/960 =56703.125
Ứng suất = 2*56703.125*4*1.1*1/(49*2.5*49*2.5) = 85.8 MPa

Tính các lực thành fần để đưa vào Catia tính :
Ft =2T/d =2*56703/(2.5*19) =2387.5 N
Fr = Ft*tg(anpha) =2387.5* tg(20 độ) =868.99 N
Mình đã vẽ bánh răng Z1 theo cách vẽ dùng các tham số của nó có trong diễn đàn và quá trình tính toán như các hình dưới:



Nhìn vào thấy ứng suất thu đ từ fần mềm lên tới 105 MPa.
Từ tr tới j mình chỉ tự tìm hiểu về Catia nên chắc có lỗi j đó,mong mọi người chỉ giúp !!!
(Nãy type bài chỉ link ảnh theo photobuket mào post rồi k nhìn thấy,mình up lại ảnh ở dưới)


Tham khảo Topic này để biết cách up hình lên Forum .
http://www.meslab.org/mes/threads/22475-tong-hop-huong-dan-ho-tro-cho-cac-thanh-vien-moi#3
Nguyenthanh.

Ðề: Kiểm tra ứng suất tại chân răng khi ăn khớp bằng Catia các Pro check giúp !!!

Trước hết cám ơn bạn đã quan tâm.Còn fần tính mình giải thích thêm thế này.Fần tính theo lí thuyết,thì hầu như tất cả các công thức tính đều là lấy gần đúng,và nó là ứng suất tổng của cả ứng suất uốn và nén rồi.Công thức mình sd trên là lấy theo công thức này.

Còn công thức tính Ft,Fr ở trên là mình tính lực từ điều kiện đề bài để cho vào Catia.Khi ăn khớp thì chỉ có lực từ bánh răng chủ động tác dụng lên bánh bị động tại đường ăn khớp trên 1 răng ăn khớp đó,mình chỉ đưa lực vào 1 răng đó và Catia cũng chỉ tính ứng suất trên chiếc răng đó(Ngàm giữ thì ở trục bánh răng).Catia tính thì nó là ứng suất tổng cộng rồi.Mình nghĩ là các bước của mình đã hợp lí.
Còn 1 vấn đề nữa là,kết quả tính ứng suất trên catia có khác nhau với cùng thông số lực tác dụng đưa vào khi ta thay đổi kích thước của lưới chia như hình dưới :
Cái này mình đọc thì hướng dẫn nói kích thước càng nhỏ thì kết quả càng chính xác,nhưng mình thử thì tấy kết quả cũng chưa gần tới kết quả tính theo lí thuyết.
Mong bạn đóng góp! Thân!