Mô hình hoá ứng xử của vật liệu bột kim loại biến dạng ở trạng thái nguội

bộ giáo dục đào tạo trường đại học bách khoa hà nội phạm văn cường mô hình hoá ứng xư cđa vËt liƯu bét kim lo¹i biÕn d¹ng ë trạng thái nguội luận án tiến sĩ kỹ thuật Hà nội - 2005 giáo dục đào tạo trường đại học bách khoa hà nội phạm văn cường mô hình hoá ứng xử vật liệu bột kim loại biến dạng trạng thái nguội Chuyên ngành : công nghệ máy gia công áp lực Mà số: 2.01.10 ln ¸n tiÕn sÜ kü tht Ng­êi h­íng dÉn khoa học: PGS, TS Nguyễn Trọng Giảng TS Đào Minh Ngừng Hà Nội - 2005 danh mục ký hiệu, chữ viết tắt A( i ), B( i ), K( i ), Các hệ số phụ thuộc vật liệu, hàm biến nội hệ số phụ thuộc vật liệu, hàm biến nội i biến trạng thái tỉ lệ hàm lượng cấu tử f hàm chất tải hàm nhiệt động hàm tiêu hao ten xơ ứng suất ' ten xơ ứng suất lệch ij thành phần ten xơ ứng suất 'ij thành phần ten xơ ứng suất lệch ten xơ biến dạng e ten xơ biến dạng đàn hồi ε VP biÕn d¹ng thĨ tÝch ε eq biÕn d¹ng tương đương eq ứng suất tương đương mật độ tương đối độ xốp S ứng suất chảy cđa vËt liƯu bét σ0 øng st ch¶y cđa vËt liệu đặc m ứng suất trung bình J1 bất biến bËc mét cđa ten x¬ øng st J2 bÊt biÕn bËc hai cđa ten x¬ øng st J3 bÊt biÕn bËc ba cđa ten x¬ øng st J '2 BÊt biÕn bËc hai cđa ten x¬ øng st lƯch λ hệ số nhân dẻo a i , b, c, C, m i , M, n i , p i , q i hệ số danh mục bảng TT Tên bảng Trang Bảng 1.1 Các hệ số hàm chất tải đề xuất tác giả 11 Bảng 2.1 Kích thước hợp lý phân tố thể tích đặc trưng 44 Bảng 4.1 Các tính chất đặc trưng bột Fe bột Cu 75 Bảng 4.2 Thành phần hỗn hợp bột 76 Bảng 4.3 Tỉ trọng lắc độ xốp ban đầu hỗn hợp bột 77 Bảng 4.4 Kết thí nghiệm ép thuỷ tĩnh hỗn hợp bột Fe- Cu ứng 81 với tỉ lệ hàm lượng cấu tử khác Bảng 4.5 Kết thí nghiệm trường hợp ép chiều 86 Bảng 4.6 Các giá trị ứng suất-biến dạng thiết lập từ đồ thị thực 88 nghiệm ép thuỷ tĩnh ép chiều danh mục hình vẽ, đồ thị TT Tên hình vẽ, đồ thị Trang Hình 1.1 Sơ đồ mô tả vai trò mô hình hoá ứng xử học vật liệu Hình 1.2 Mô tả hình học mặt dẻo Von- Mises không gian ứng suất Hình 1.3 Mô tả hình học mặt dẻo Tresca không gian ứng suất Hình 1.4 Mô tả hình học mặt dẻo Tresca mặt phẳng ten xơ ứng suất lệch 10 Hình 1.5 Minh hoạ hàm chất tải không gian ứng suất lệch ứng suất thuỷ tĩnh 15 Hình 1.6 Điều kiện dẻo Kuhn Downey không gian ứng suất lệch ứng suất thuỷ tĩnh 25 Hình 1.7 Quỹ chất tải đối xứng mặt dẻo tương ứng với mật độ khác không gian ứng suất lệch ứng suất thuỷ tĩnh 32 Hình 1.8 So sánh đường cong áp lực- mật độ nhận tính toán theo mô hình với kết thực nghiệm bột Cu có kích thước hạt 2,5à 100nm 48 Hình 2.1 Mặt đồng tiêu hao 50 10 Hình 2.2 Minh hoạ mặt dẻo không liên kết với hướng chảy dẻo 55 11 Hình 3.1 Sơ đồ mô tả phương pháp liên tục hoá môi trường xốp rời rạc 69 12 Hình 4.1 Quy trình thí nghiệm nhận dạng mô hình ứng xử hỗn hợp bột kim loại biến dạng trạng thái nguội 71 13 Hình 4.2 Minh hoạ mặt cắt dọc mẫu ép chiều 72 14 Hình 4.3 Minh hoạ ¶nh h­ëng cđa ¸p st khÝ c¸c khoang xèp liên thông đến biến dạng môi trường xốp 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 H×nh 4.4 Bng Ðp thủ tÜnh H×nh 4.5 Quan hệ ứng suất thuỷ tĩnh biến dạng thể tích hỗn hợp bột Fe- Cu Hình 4.6 ảnh HVĐT quét mẫu bột Fe(a), Cu(b) Hình 4.7 Máy trộn lệch tâm Hình 4.8 Thiết bị xác định tỉ trọng lắc hỗn hợp bột Hình 4.9 Thiết bị ép thuỷ tĩnh Hình 4.10 Buồng ép thuỷ tĩnh cïng víi mÉu H×nh 4.11 MÉu sau Ðp thủ tĩnh Hình 4.12 Quan hệ ứng suất - biến dạng trường hợp ép thuỷ tĩnh Hình 4.13 Sơ đồ nguyên lý thiết bị ép chiều Hình 4.14 Máy ép mẫu Hình 4.15 Mẫu ép chiều Hình 4.16 Các đường cong ứng suất- biến dạng thực nghiệm nhận trường hợp ép chiều Hình 4.17 Quan hệ hệ số B độ xốp hỗn hợp bột ứng với tỉ lệ hàm lượng cấu tử khác Hình 4.18 Quan hệ hệ số A độ xốp ứng với tỉ lệ hàm lượng cấu tử khác Hình 4.19 Quan hệ hệ số K' độ xốp ứng với tỉ lệ hàm lượng cấu tử khác Hình 4.20 Quan hệ biến thiên độ xốp với mức độ gia tăng biến dạng Hình 4.21 Quan hệ hệ số A( i ) biến nội Hình 4.22 Quan hệ hệ số B( i ) biến nội Hình 4.23 Quan hệ hệ số K'( i ) với độ xốp tỉ lệ hàm lượng cấu tử Hình 4.24 Đồ thị kiểm tra kết nhận dạng hệ số A( i ) Hình 4.25 Đồ thị kiểm tra kết nhận dạng hàm B( i ) Hình 4.26 Đồ thị kiểm tra kết nhận dạng hàm K( i ) 74 74 75 76 77 78 79 80 82 83 84 85 87 89 90 91 92 93 95 96 97 98 99 38 39 40 41 42 43 44 45 46 47 48 49 Hình 4.27 Quan hệ trở kháng biến dạng hỗn hợp bột Fe Cu (tỉ lệ hàm lượng cấu tử = 1,2) với độ xốp ứng suất thuỷ tĩnh để gây biến dạng Hình 4.28 a, b,c,d Đồ thị so sánh kết mô thực nghiệm trường hợp ép chiều Hình 4.29 a, b,c,d Đồ thị so sánh kết mô thực nghiệm trường hợp ép thuỷ tĩnh Hình 5.1 Sơ đồ biến dạng lớp tiếp xúc với thành khuôn Hình 5.2 Khả cạnh tranh ma sát Coulomb Ma sát Tresca trình ép tạo hình vật liệu bột (mũi tên chiều gia tăng) Hình 5.3 Sơ đồ thí nghiệm xác định hệ số ma sát thành khuôn Hình 5.4 Thiết bị đo độ nhám Hình 5.5 Quan hệ lực đẩy độ xốp mẫu ứng với cấp độ nhám thành khuôn khác Hình 5.6 Quan hệ ứng suất ma sát độ xốp ứng với cấp độ nhám thành khuôn khác Hình 5.7 Quan hệ hệ số ma sát m độ xốp ứng với cấp độ nhám thành khuôn khác Hình 5.8 Quan hệ ứng suất ma sát độ nhám thành khuôn Hình 5.9 Mặt biểu diễn biến đổi hệ số ma sát phụ thuộc vào độ xốp độ nhám thành khu«n 100 101 105 113 114 115 116 118 119 119 121 122 mục lục Trang Lời nói đầu Chương Tổng quan mô hình ứng xử học vật liệu bột kim loại trình biến dạng tạo hình 1.1 Vai trò mô hình hoá ứng xử học vật liệu 1.2 Các mô hình ứng xử học vật liệu bột kim loại trình biến dạng tạo hình 1.3 Kết luận Chương Cơ sở lý thuyết mô hình hoá thuộc tính học vật liệu 2.1 Phương pháp trạng thái cục 2.1.1 Biến trạng thái 2.1.2 Thế nhiệt động học- định luật trạng thái 2.1.3 Năng lượng tiêu hao - Các định luật bổ sung 2.2 Phương pháp thiết lập mô hình thuộc tính 2.2.1 Phương pháp tượng tổng thể xây dựng mô hình thuộc tính 2.2.2 Thiết lập mô hình dạng phương trình giải tích 2.2.3 Mô hình hoá đặc trưng phi tuyến 2.2.4 Sự phân tán hệ số ngẫu nhiên 2.2.5 Các phương pháp nhận dạng mô hình 2.3 Mặt chất tải tiêu hao 2.3.1.Mặt chất tải điều kiện dỡ tải-chất tải 2.3.2 Dẻo không liên kết 2.3.3 Biểu thức hệ số nhân dẻo 2.3.4 Dẻo liên kết 2.4 KÕt luËn Ch­¬ng 33 36 36 37 39 42 43 44 46 46 46 47 47 49 50 52 52 Xây dựng mô hình ứng xử hỗn hợp bột kim loại biến dạng trạng thái nguội 3.1 Xây dựng mô hình 3.1.1 Phương pháp tiếp cận 3.1.2 Mô hình ứng xử 3.2 Phương pháp nhận dạng mô hình 3.2.1 Quá trình ép chiều khuôn trụ 3.2.2 Quá trình ép thuỷ tĩnh 3.2.3 Xác định hệ số mô hình 54 54 56 61 62 65 68 3.3 Kết luận Chương áp dụng mô hình cho hỗn hợp bột Fe Cu 4.1 Phương pháp thực nghiệm nhận dạng mô hình 4.1.1 Đánh giá hiệu bôi trơn thành khuôn 4.1.2 Nghiên cứu ảnh hưởng áp suất khoang xốp liên thông đến ứng xử hỗn hợp bột thí nghiệm ép thuỷ tĩnh 4.2 áp dụng mô hình cho hỗn hợp bột Fe - Cu 4.2.1.Chuẩn bị hỗn hỵp bét 4.2.2 ThÝ nghiƯm Ðp thủ tÜnh 4.2.3 ThÝ nghiệm ép chiều 4.2.4 Nhận dạng hệ số 4.3 Kiểm tra kết nhận dạng mô hình 4.4 Kiểm tra mức độ tương hợp mô hình 4.5 Kết luận Chương Ma sát hỗn hợp bột kim loại thành Khuôn trình ép tạo hình 5.1.Vai trò ma sát tiếp xúc mô số trình tạo hình 5.2 Cơ sở lý thuyết 5.2.1 Các mô hình ma sát 5.2.2 Ma sát bột dụng cụ tạo hình 5.3 Phương pháp thực nghiệm 5.3.1.Mô tả phương pháp 5.3.2 Chuẩn bị mẫu 5.4 Kết thảo luận 5.5 Kết luận Kết luận chung Các công trình khoa học đà công bố Tài liệu tham khảo Phụ lục: Phần mềm nhận dạng hệ số mô hình 68 69 70 71 75 75 78 82 87 96 101 108 110 111 111 113 115 115 116 118 122 124 127 128 lời nói đầu Công nghệ chế tạo sản phẩm từ bột kim loại ngày phát triển mạnh mẽ Bằng công nghệ luyện kim bột (LKB) đà tạo nhiều sản phẩm ứng dụng ngành công nghiệp vũ trụ, hàng không, công nghiệp chế tạo ô tô vv Chế tạo chi tiết phương pháp LKB gồm bước sau: trộn bột, ép tạo hình , thiêu kết gia công sau thiêu kết Vật liệu bột có cấu trúc môi trường xốp rời rạc Do chất rời rạc bột tương tác hạt bột với chúng với dụng cụ ép mà nhiều tượng học phức tạp xảy trình tạo hình Các tượng tổ hợp biến dạng đàn dẻo, ma sát môi trường bột ma sát tiếp xúc bột với dụng cụ vv Các tượng tạo phân bố lỗ xốp không đồng gây vết nứt phôi sau ép Vấn đề đặt nhà công nghệ không chế tạo phôi không nứt, hình dạng kích thước xác, có độ bền học định để vận hành công đoạn mà phải khống chế phân bố độ xốp phù hợp với yêu cầu sử dụng Thiết kế tối ưu trình biến dạng tạo hình nói chung tạo hình vật liệu bột nói riêng thường khó khăn Thông thường người ta dựa vào kinh nghiệm sản phẩm đơn giản, sản phẩm có hình dạng phức tạp áp dụng phương pháp "thử sai" tiêu tốn nhiều thời gian, nhân lực kinh phí Mô số công cụ hữu hiệu cho phép dự báo vấn đề xảy trình tạo hình sản phẩm gợi ý phương án giải chúng Phân bố lỗ xốp yếu tố quan trọng định chất lượng sản phẩm sau tạo hình Vấn đề dự báo phân bố lỗ xốp thể tích vật liệu ép tạo hình, đặc biệt chi tiết có hình dạng phức tạp đòi hỏi phải giải toán tối ưu ảnh hưởng thông số công nghệ đến phân bố ứng suất 143 44 CONTINUE RN=REAL(N) S(M,M)=RN C C Vec-to phai C DO 46 I=1,M-1 DO 45 J=1,N SS(I)=SS(I)+XX(I,J)*YI(J) 45 46 CONTINUE CONTINUE DO 47 J=1,N SS(M)=SS(M)+YI(J) 47 CONTINUE C DO 49 I=1,M DO 48 J=1,M AS(I,J)=S(I,J) 48 49 CONTINUE CONTINUE DO 50 I=1,M BS(I)=SS(I) 50 CONTINUE C C Gan cac phan cuar ma tran A C -DO 52 I=1,M DO 51 J=3,4 144 AS(I,J)=-AS(I,J) 51 52 CONTINUE CONTINUE C WRITE (2,*)' MATRICE A !' WRITE(6,'(5(/,5F16.6))')AS WRITE(2,'(5(/,5F16.6))')AS WRITE (6,*)' MATRICE B !' WRITE (2,*)' MATRICE B !' WRITE(6,'(5F16.6)')BS WRITE(2,'(5F16.6)')BS C C Giai he phuong trinh dai so tuyen tinh bang phuong phap Gauss C C CALL GAUSS(AS,BS,XS,M,DET,RES) IF (.NOT.RES) STOP 77 WRITE(2,*)' matrices A,B,X :' WRITE(6,*)' matrices A,B,X :' WRITE(2,85)AS,BS,XS WRITE(6,85)AS,BS,XS 85 FORMAT (5(/,3X,5F14.6)) PUISQ=XS(1) PUISM=XS(2) PUISN=XS(3) PUISP=xs(4) CONSTA=10.**XS(5) C 145 C Tinh he so tuong hop C LOG10(Z(I)-3.0)=LOG10(a)+LOG10(x**q*(x+y)**m)/ C (y**p*(1-x)**n) C - C DO 90 I=1,N LOZ(I)=LOG10(X(I)**PUISQ*(x(i)+y(i))**PUISM/ & ((1.0-X(I))**PUISN*y(i)**PUISP)) LOS(I)=LOG10(Z(I)-3.0) WRITE (6,105)LOZ(I),LOS(I) WRITE (2,105)LOZ(I),LOS(I) 90 CONTINUE C SS1=0.0 E+00 SS2=0.0 E+00 DO 95 I=1,N SS1=SS1+LOZ(I) SS2=SS2+LOS(I) 95 CONTINUE XB=SS1/RN YB=SS2/RN SS11=0.0 E+00 SS22=0.0 E+00 SS33=0.0 E+00 DO 100 I=1,N XII=LOZ(I)-XB YII=LOS(I)-YB SS11=SS11+XII*XII 146 SS22=SS22+YII*YII SS33=SS33+XII*YII 100 CONTINUE KCOR=SS33/SQRT(SS11*SS22) WRITE(2,110)PUISQ,PUISN,PUISM,PUISP,CONSTA,KCOR 105 FORMAT (3X,2F12.6) 110 FORMAT (//,'Q=',F12.10,' ','N =',F12.10,3X,1X,' ','M=',F12.10,' ', & 'P=',F12.10,1X,' ','A=',F12.10,' ','COEF CORREL.=',F12.10,//) CLOSE(1) CLOSE(2) STOP END C C - C*************************************************** SUBROUTINE GAUSS (AC,BC,XC,M,D,RES) C*************************************************** REAL*8 AC(M,M),BC(M),XC(M) INTEGER H LOGICAL RES RES=.FALSE M1=M-1 DO 1000 I=1,M1 H=I+1 IF (AC(I,I).NE.0.) GO TO 150 C Tim mot Pivot khac C DO 120 K=H,M 147 IF (AC(K,I).NE.0.) GO TO 130 120 CONTINUE C PAS DE PIVOT NON NUL GO TO 1040 C Dao cho hai hang C 130 DO 140 J=I,M V=AC(I,J) AC(I,J)=AC(H,J) 140 AC(H,J)=V V=BC(I) BC(I)=BC(H) BC(H)=V C C Triet tieu mot phan AC(H,I) C 150 IF (AC(H,I).EQ.0.) GO TO 180 R=-AC(H,I)/AC(I,I) DO 160 J=I,M 160 AC(H,J)=AC(H,J)+R*AC(I,J) BC(H)=BC(H)+R*BC(I) 180 H=H+1 IF (H.LE.M) GO TO 150 1000 CONTINUE IF (AC(M,M).EQ.0.)GO TO 1040 C C Giai he phuong trinh voi ma tran he so dang tam giac C 148 XC(M)=BC(M)/AC(M,M) I=M1 1010 V=0 I1=I+1 DO 1020 H=I1,M 1020 V=V+AC(I,H)*XC(H) XC(I)=(BC(I)-V)/AC(I,I) I=I-1 IF (I.GE.1) GO TO 1010 D=1 DO 1030 I=1,M 1030 D=D*AC(I,I) RES=.TRUE 1040 RETURN END 149 C****************************************************** C CHUONG TRINH HOI QUI TUYEN TINH XAC DINH CAC HAM B, K’ C (dang: z=a*x**q*(x+y)**m/((1-x)**n*y**p) ) C****************************************************** C PARAMETER (N=20,M=5) REAL*8 X(N),Y(N),Z(N),XX(M,N),YI(N),XO(N),YO(N),ZO(N) 150 REAL*8 AS(M,M),XS(M),BS(M),LOZ(N),LOS(N) REAL*8 S(M,M),SS(M) REAL*8 KCOR LOGICAL RES CHARACTER*18 NOM1,NOM2 NOM1='resultb12' NOM2='ketquab' C OPEN (1,FILE=NOM1,STATUS='OLD',FORM='FORMATTED') OPEN (2,FILE=NOM2,STATUS='unknown',FORM='FORMATTED') C DO 20 I=1,N READ (1,*)XO(I),YO(I),ZO(I) C dao cot (ty le ham luong) cho cot (do xop) X(I)=YO(I) Y(I)=XO(I) Z(I)=ZO(I) C XX(1,I)=LOG10(X(I)) XX(2,I)=LOG10(X(i)+Y(I)) XX(3,I)=LOG10(1.0-X(I)) XX(4,i)=log10(Y(i)) YI(I)=LOG10(Z(I)) WRITE (6,30)XX(1,I),XX(2,I),XX(3,I),xx(4,i),YI(I) 20 CONTINUE 30 FORMAT (3X,5F14.6) C C Khoi thao cac ma tran 151 C DO 32 I=1,M DO 31 J=1,M S(I,J)=0.0 31 32 CONTINUE CONTINUE DO 33 I=1,M SS(I)=0.0 33 CONTINUE C C Matran trai C K=1 34 DO 40 I=1,M-1 DO 38 J=1,N S(K,I)=S(K,I)+XX(K,J)*XX(I,J) 38 40 CONTINUE CONTINUE K=K+1 IF(K.LT.M) THEN GOTO 34 ENDIF DO 44 I=1,M-1 DO 43 J=1,N S(M,I)=S(M,I)+XX(I,J) 43 CONTINUE S(I,M)=S(M,I) 44 CONTINUE 152 RN=REAL(N) S(M,M)=RN C C Vec-to phai C DO 46 I=1,M-1 DO 45 J=1,N SS(I)=SS(I)+XX(I,J)*YI(J) 45 46 CONTINUE CONTINUE DO 47 J=1,N SS(M)=SS(M)+YI(J) 47 CONTINUE C DO 49 I=1,M DO 48 J=1,M AS(I,J)=S(I,J) 48 49 CONTINUE CONTINUE DO 50 I=1,M BS(I)=SS(I) 50 CONTINUE C C Gan cac phan cua ma tran A C -DO 52 I=1,M DO 51 J=3,4 AS(I,J)=-AS(I,J) 153 51 52 CONTINUE CONTINUE C WRITE (6,*)' MATRICE A !' WRITE (2,*)' MATRICE A !' WRITE(6,'(5(/,5F16.6))')AS WRITE(2,'(5(/,5F16.6))')AS WRITE (6,*)' MATRICE B !' WRITE (2,*)' MATRICE B !' WRITE(6,'(5F16.6)')BS WRITE(2,'(5F16.6)')BS C C Giai he phuong trinh dai so tuyen tinh C bang phuong phap Gauss C C CALL GAUSS(AS,BS,XS,M,DET,RES) IF (.NOT.RES) STOP 90 WRITE(2,*)' matrices A,B,X :' WRITE(6,*)' matrices A,B,X :' WRITE(2,85)AS,BS,XS WRITE(6,85)AS,BS,XS 85 FORMAT (5(/,3X,5F14.6)) PUISQ=XS(1) PUISM=XS(2) PUISN=XS(3) PUISP=xs(4) CONSTA=10.**XS(5) 154 C C Tinh he so tuong hop C LOG10(Z(I))=LOG10(a)+LOG10(x**q*(x+y)**m)/ c (y**p*(1-x)**n) C - C DO 90 I=1,N LOZ(I)=LOG10(X(I)**PUISQ*(x(i)+y(i))**PUISM/ & ((1.0-X(I))**PUISN*y(i)**PUISP)) LOS(I)=LOG10(Z(I)) WRITE (6,105)LOZ(I),LOS(I) WRITE (2,105)LOZ(I),LOS(I) 90 CONTINUE C SS1=0.0 E+00 SS2=0.0 E+00 DO 95 I=1,N SS1=SS1+LOZ(I) SS2=SS2+LOS(I) 95 CONTINUE XB=SS1/RN YB=SS2/RN SS11=0.0 E+00 SS22=0.0 E+00 SS33=0.0 E+00 DO 100 I=1,N XII=LOZ(I)-XB YII=LOS(I)-YB 155 SS11=SS11+XII*XII SS22=SS22+YII*YII SS33=SS33+XII*YII 100 CONTINUE KCOR=SS33/SQRT(SS11*SS22) WRITE(2,110)PUISQ,PUISN,PUISM,PUISP,CONSTA,KCOR 105 FORMAT (3X,2F12.6) 110 FORMAT (//,'Q=',F12.10,' ','N =',F12.6,3X,1X,' ','M=',F12.10,' ', & 'P=',F12.10,1X,' ','A=',F12.6,' ','COEF CORREL.=',F12.10,//) CLOSE(1) CLOSE(2) STOP END C C============================================ C ************************************************ SUBROUTINE GAUSS (AC,BC,XC,M,D,RES) C ************************************************ REAL*8 AC(M,M),BC(M),XC(M) INTEGER H LOGICAL RES RES=.FALSE M1=M-1 DO 1000 I=1,M1 H=I+1 IF (AC(I,I).NE.0.) GO TO 150 C Tim mot pivot khac C - 156 DO 120 K=H,M IF (AC(K,I).NE.0.) GO TO 130 120 C CONTINUE Khong co pivot khac GO TO 1040 C Doi cho hai hang C 130 DO 140 J=I,M V=AC(I,J) AC(I,J)=AC(H,J) 140 AC(H,J)=V V=BC(I) BC(I)=BC(H) BC(H)=V C Triet tieu mot phan AC(H,I) C 150 IF (AC(H,I).EQ.0.) GO TO 180 R=-AC(H,I)/AC(I,I) DO 160 J=I,M 160 AC(H,J)=AC(H,J)+R*AC(I,J) BC(H)=BC(H)+R*BC(I) 180 H=H+1 IF (H.LE.M) GO TO 150 1000 CONTINUE IF (AC(M,M).EQ.0.)GO TO 1040 C C Giai he phuong trinh voi ma tran he so C dang tam giac 157 C -XC(M)=BC(M)/AC(M,M) I=M1 1010 V=0 I1=I+1 DO 1020 H=I1,M 1020 V=V+AC(I,H)*XC(H) XC(I)=(BC(I)-V)/AC(I,I) I=I-1 IF (I.GE.1) GO TO 1010 D=1 DO 1030 I=1,M 1030 D=D*AC(I,I) RES=.TRUE 1040 RETURN END ... có mô hình sử dụng thuận lợi cho mô số trình tạo hình vật liệu bột kim loại trạng thái nguội Như vậy, mô hình tổng quát có khả mô tả xác ứng xử vật liệu bột kim loại biến dạng trạng thái nguội. .. trình biến dạng tạo hình 1.1 Vai trò mô hình hoá ứng xử học vật liệu 1.2 Các mô hình ứng xử học vật liệu bột kim loại trình biến dạng tạo hình 1.3 Kết luận Chương Cơ sở lý thuyết mô hình hoá thuộc... phần lớn trình tạo hình vật liệu biến dạng dẻo 4 Chương tổng quan mô hình ứng xử học vật liệu bột kim loại trình biến dạng tạo hình 1.1.Vai trò mô hình hoá ứng xử học vật liệu Mô số xem công cụ