TAILIEUCHUNG - The boundary element method with programming for engineers and scientists - phần 5

Để thực hiện các chương trình dễ đọc hơn và dễ dàng hơn để chỉnh sửa, đọc của đầu vào đã được giao cho chương trình con. Điều này cũng cung cấp cho người đọc một số tự do để xác định định dạng đầu vào và thực hiện các phương tiện thế hệ đơn giản lưới. | 192 The Boundary Element Method with Programming IF NCode Ldest nel nd 0 THEN Elres_u nel nd Elres_u nel nd u1 Ldest nel nd ELSE Elres_t nel nd Elres_t nel nd u1 Ldest nel nd END IF END IF END DO D_o_F1 Elres_u nel Elres_u nel Scad Elres_t nel Elres_t nel Scat WRITE 2 Elres_u nel m m 1 Ndofe WRITE 2 Elres_t nel m m 1 Ndofe END DO Elements_2 END PROGRAM To make the program more readable and easier to modify the reading of the input has been delegated to subroutines. This also gives the reader some freedom to determine the input FORMAT and implement simple mesh-generation facilities. SUBROUTINE Jobin Title Cdim Ndof Toa Nreg Ltyp Con E ny Isym nodel nodes maxe ------------------------------------------------- Subroutine to read in basic job information ------------------------------------------------- CHARACTER LEN 80 INTENT OUT Title INTEGER INTENT OUT Cdim Ndof Toa Nreg Ltyp Isym nodel INTEGER INTENT OUT Nodes Maxe REAL INTENT OUT Con E ny READ 1 A80 Title WRITE 2 Project Title READ 1 Cdim WRITE 2 Cartesian dimension Cdim READ 1 Ndof Degrees of freedom per node IF NDof 1 THEN WRITE 2 Potential Problem ELSE WRITE 2 Elasticity Problem END IF IF Ndof 2 THEN READ 1 Toa Analysis type plane strain 1 plane stress 2 IF Toa 1 THEN WRITE 2 Type of Analysis Solid Plane Strain ELSE WRITE 2 Type of Analysis Solid Plane Stress END IF END IF READ 1 Nreg Type of region IF NReg 1 THEN WRITE 2 Finite Region ASSEMBLY AND SOLUTION 193 ELSE WRITE 2 Infinite Region END IF READ 1 Isym Symmetry code SELECT CASE isym CASE 0 WRITE 2 No symmetry CASE 1 WRITE 2 CASE 2 Symmetry about y-z plane WRITE 2 CASE 3 Symmetry about y-z and x-z planes WRITE 2 Symmetry about all planes END SELECT READ 1 Ltyp Element type IF Ltyp 1 THEN WRITE 2 Linear Elements ELSE WRITE 2 Quadratic Elements END IF Determine number of nodes per element IF Cdim 2 THEN Line elements IF Ltyp 1 THEN Nodel 2 ELSE Nodel 3 END IF ELSE S urface elements IF Ltyp 1 THEN Nodel 4 ELSE Nodel 8 END IF END IF Read .

TỪ KHÓA LIÊN QUAN
Đã phát hiện trình chặn quảng cáo AdBlock
Trang web này phụ thuộc vào doanh thu từ số lần hiển thị quảng cáo để tồn tại. Vui lòng tắt trình chặn quảng cáo của bạn hoặc tạm dừng tính năng chặn quảng cáo cho trang web này.