Linear hay nonlinear trong SAP2000: chọn đúng mô hình cho giằng rod và cable

Một hệ giằng chữ X có thể cho chuyển vị rất nhỏ khi chạy linear, nhưng chuyển vị tăng rõ rệt khi chạy nonlinear. Nguyên nhân thường nằm ở cơ chế chịu lực: mô hình linear cho cả hai đường chéo tham gia, trong khi giằng rod thực tế được thiết kế chỉ làm việc khi chịu kéo.

Chọn đúng loại phân tích cần đi cùng chọn đúng phần tử, điều kiện liên kết và cách chất tải. Bài viết tập trung vào giằng rod trong nhà xưởng, đồng thời làm rõ khi nào cần P-Delta, large displacement hoặc mô hình cable để mô tả độ võng của chính dây.

1. Linear và nonlinear khác nhau ở cơ chế làm việc

Trong phân tích tuyến tính với độ cứng cố định:

$$\mathbf{K}\mathbf{u}=\mathbf{F}$$

• $\mathbf{K}$ — ma trận độ cứng của hệ.

• $\mathbf{u}$ — vectơ chuyển vị nút.

• $\mathbf{F}$ — vectơ tải trọng nút tương đương.

Khi giả thiết tuyến tính còn phù hợp, tải tăng hai lần thì đáp ứng tăng hai lần; có thể cộng kết quả của các trường hợp tải trên cùng nền độ cứng. Điều này giúp việc lập nhiều tổ hợp trở nên nhanh và dễ kiểm tra.

Trong phân tích phi tuyến, trạng thái làm việc có thể thay đổi theo tải. Phương trình cân bằng tổng quát được viết:

$$\mathbf{f}_{int}(\mathbf{u},\mathbf{s})-\mathbf{F}=\mathbf{0}$$

• $\mathbf{f}_{int}$ — vectơ nội lực quy về các bậc tự do của hệ.

• $\mathbf{s}$ — trạng thái phần tử, vật liệu và lịch sử biến dạng.

• $\mathbf{0}$ — vectơ phần dư bằng không tại trạng thái cân bằng lý tưởng.

Ví dụ, một rod ban đầu đang kéo có thể chùng khi gió đổi chiều; gối có thể nhấc lên; khớp dẻo có thể hình thành. Phần mềm phải tìm cân bằng đồng thời với việc cập nhật các trạng thái đó.

Ba nguồn phi tuyến cần phân biệt:

2. Vì sao giằng rod cần phân tích tension-only?

2.1. Một thanh chịu nén trong mô hình chưa chắc chịu nén được ngoài thực tế

Giằng thép góc, thép ống hoặc thanh tổ hợp có thể được thiết kế chịu cả kéo và nén. Ngược lại, rod mảnh với liên kết và cấu tạo chỉ phục vụ truyền lực kéo thường được lý tưởng hóa là tension-only. Cần dựa vào cấu tạo và mục tiêu thiết kế, không áp dụng giả thiết này cho mọi loại giằng.

Hình 1 cho thấy khác biệt ở một ô giằng: khi đỉnh ô dịch sang phải, đường chéo từ chân trái lên đỉnh phải dài ra; đường chéo còn lại ngắn lại. Mô hình thanh đàn hồi hai chiều cho đường chéo thứ hai chịu nén. Mô hình chỉ kéo cho đường chéo đó mất tác dụng dọc.

2.2. Gán Compression Limit bằng 0 chưa đủ

Theo CSI Help — Frame Tension/Compression Limits, giới hạn kéo/nén được kích hoạt trong nonlinear static hoặc nonlinear direct-integration time history. Linear dùng độ cứng tại trạng thái chưa ứng suất vẫn sử dụng độ cứng dọc đầy đủ như không có giới hạn. Linear dùng độ cứng cuối một case nonlinear sẽ giữ độ cứng dọc theo trạng thái cuối case đó.

Vì vậy, đối với rod chỉ kéo, cần đồng thời có khai báo giới hạn dọc phù hợp và Load Case nonlinear phù hợp. Chỉ sửa thuộc tính thanh rồi chạy lại case linear không giải quyết được vấn đề.

Hình 2 biểu diễn quy luật lý tưởng của thanh chỉ kéo không căng trước:

$$N=\max\left(0,\frac{EA}{L}\Delta L\right)$$

• $N$ — lực dọc, quy ước kéo dương (N).

• $E$ — mô đun đàn hồi (N/mm²).

• $A$ — diện tích tiết diện (mm²).

• $L$ — chiều dài ban đầu (mm).

• $\Delta L$ — độ thay đổi chiều dài (mm), dương khi dài ra.

Quy luật này chưa bao gồm độ chùng ban đầu, lực căng trước hoặc chảy dẻo. Giới hạn kéo/nén của Frame là mô hình đàn hồi phi tuyến; không nên coi nó là một mô hình phá hoại hoặc khớp dẻo của rod.

3. Ví dụ tính tay: chuyển vị của ô giằng chữ X

3.1. Số liệu và giả thiết

Xét một ô giằng lý tưởng để kiểm chứng mô hình. Hai nút đỉnh cùng dịch ngang một lượng $\Delta$, không dịch đứng; hai nút chân cố định vị trí. Bỏ qua độ cứng uốn ngang của khung, tự trọng rod và lực căng trước. Hai đường chéo không liên kết tại giao điểm. Phân tích biến dạng nhỏ, vật liệu đàn hồi.

Đây là ví dụ tính tay với số liệu giả định, không phải kết quả đã chạy SAP2000. Các giá trị dưới đây là mốc đối chiếu cho một mô hình có đúng các giả thiết trên.

Chiều dài rod và hệ số chiếu theo phương ngang:

$$L=\sqrt{4^2+3^2}=5\ \mathrm{m}=5.000\ \mathrm{mm},\qquad c=\frac{b}{L}=0{,}8$$

• $b,h$ — kích thước ô giằng.

• $c$ — cosin góc nghiêng của rod so với phương ngang.

3.2. Trường hợp chỉ một rod chịu kéo

Cân bằng ngang cho rod đang làm việc:

$$N=\frac{H}{c}=\frac{20}{0{,}8}=25\ \mathrm{kN}$$

Độ giãn dọc của rod:

$$\Delta L=\frac{NL}{EA}=\frac{25.000\times5.000}{200.000\times314{,}16}=1{,}989\ \mathrm{mm}$$

Do tương thích biến dạng nhỏ $\Delta L\approx c\Delta$, chuyển vị ngang là:

$$\Delta_{TO}=\frac{HL}{EA c^2}=\frac{20.000\times5.000}{200.000\times314{,}16\times0{,}8^2}=2{,}487\ \mathrm{mm}$$

• $\Delta_{TO}$ — chuyển vị ngang với giằng chỉ kéo (mm).

Rod còn lại có lực dọc bằng 0 theo lý tưởng hóa. Khi đảo chiều tải, hai rod đổi vai trò kéo và chùng.

3.3. Trường hợp linear huy động cả hai rod

Hai thanh giống nhau cùng đóng góp độ cứng ngang. Lực dọc có trị số bằng nhau, một kéo và một nén:

$$N_1=+\frac{H}{2c}=+12{,}5\ \mathrm{kN},\qquad N_2=-12{,}5\ \mathrm{kN}$$

Chuyển vị ngang:

$$\Delta_{LIN}=\frac{HL}{2EA c^2}=1{,}243\ \mathrm{mm}$$

• $N_1,N_2$ — lực dọc của hai đường chéo (kN).

• $\Delta_{LIN}$ — chuyển vị ngang của mô hình đàn hồi hai chiều (mm).

Bảng sau tổng hợp kết quả của cùng một ô giằng:

Trong ví dụ này, linear cho chuyển vị bằng một nửa và lực kéo trong rod bằng một nửa so với mô hình chỉ kéo. Kết quả có vẻ thuận lợi hơn vì mô hình đã huy động thêm một đường truyền lực nén không có trong giả thiết thực tế.

4. Chọn None, P-Delta hay Large Displacement?

Đặt Analysis Type là Nonlinear chưa xác định hết bài toán. Cần chọn riêng mức phi tuyến hình học.

CSI Help — Geometric Nonlinearity Parameters khuyến nghị kiểm tra mô hình với mức đơn giản trước, sau đó tăng mức phi tuyến nếu cần. Với ví dụ mục 3, case nonlinear có None đã đủ để khảo sát sự bật/tắt giằng; bật large displacement ngay sẽ làm phép đối chiếu tính tay kém trực tiếp hơn.

Ở một khung có tải đứng $P$ và chuyển vị ngang $\Delta$, mô men bổ sung có dạng:

$$M_{add}=P\Delta$$

• $P$ — tải đứng gây hiệu ứng bậc hai (kN).

• $\Delta$ — chuyển vị tương ứng (m).

• $M_{add}$ — mô men bổ sung (kN·m).

P-Delta và tension-only là hai cơ chế khác nhau. Bật P-Delta không tự biến một giằng thông thường thành thanh chỉ kéo. Tương tự, bài toán trị riêng buckling có thể dùng để sàng lọc ổn định; muốn khảo sát sau mất ổn định cần mô hình và phương pháp phi tuyến thích hợp, thường có hình học không hoàn hảo.

5. Quy trình khai báo và chất tải trong SAP2000

5.1. Từ mô hình hình học đến case kiểm chứng

1. Kiểm tra nút, kết nối và điều kiện biên. Hai rod cắt nhau trên hình không đồng nghĩa có liên kết thực tại giao điểm. Khai báo theo đúng cấu tạo.

2. Khai báo tiết diện, vật liệu và liên kết đầu thanh. Với mô hình rod truyền lực dọc, kiểm tra Assign → Frame → Releases/Partial Fixity theo trục local; tránh vô tình để rod đóng góp độ cứng uốn không có thực. Đồng thời tránh release gây bậc tự do không được khống chế.

3. Chọn rod, vào Assign → Frame → Tension/Compression Limits; bật Compression Limit, nhập 0. Không đặt Tension Limit bằng 0 khi muốn rod vẫn chịu kéo; giữ không giới hạn kéo nếu chưa chủ định mô hình giới hạn đó.

4. Vào Define → Load Cases, tạo case Static → Nonlinear. Chạy case tải ngang nhỏ với Geometric Nonlinearity = None để kiểm tra đường truyền lực trước.

5. Tạo tải hai chiều độc lập. Sau khi cơ chế đã đúng, thêm tải đứng, P-Delta hoặc lực căng trước theo mục tiêu công trình.

Tên mục và bố cục hộp thoại có thể khác giữa các phiên bản; cần đối chiếu Help của phiên bản đang sử dụng.

5.2. Không cộng tuyến tính các đáp ứng phi tuyến độc lập

Với hệ đổi trạng thái, cộng kết quả nonlinear của $D$ với kết quả nonlinear của $W$ thường không bằng kết quả khi cả hai tải cùng tác dụng. Để xét một tổ hợp, đưa các Load Pattern và hệ số của tổ hợp đó vào một case nonlinear, hoặc xây dựng chuỗi chất tải có chủ đích.

Ví dụ minh họa cách đặt case, không phải tổ hợp thiết kế theo tiêu chuẩn:

$D,W$ là các Load Pattern; $a,b$ là hệ số do người thiết kế xác định theo tiêu chuẩn và mục đích kiểm tra. Mỗi tổ hợp yêu cầu hệ số tải đứng khác nhau cần case nền tương ứng hoặc một case trực tiếp phù hợp.

Hình 3 mô tả cách chia nhánh từ cùng trạng thái tải đứng. CSI Help — Initial Conditions quy định tải khai báo trong case tiếp nối là tải tăng thêm vào trạng thái trước đó. Vì vậy, nếu case nền đã có $D$, khai báo thêm $D+W$ ở case tiếp nối sẽ tạo tổng tải $2D+W$.

Có thể lập Envelope của các case nonlinear đã giải đầy đủ để tìm trị số bao. Các cực trị của envelope có thể đến từ các case khác nhau; không ghép chúng thành một trạng thái cân bằng đồng thời.

5.3. Linear dùng độ cứng cuối nonlinear có thay thế được không?

CSI Help — Stiffness to Use cho phép case linear dùng độ cứng cuối case nonlinear. Cách này hữu ích khi đáp ứng sau đó chỉ là nhiễu nhỏ quanh một trạng thái ổn định.

Tuy nhiên, độ cứng được giữ cố định trong case linear tiếp theo. Nếu đổi chiều gió làm rod đang kéo chuyển sang chùng, cần giải nonlinear để cập nhật trạng thái. Khi tổng hợp kết quả, phải hiểu rõ đâu là đáp ứng tăng thêm của case linear và đâu là đáp ứng của case nền, tránh bỏ sót hoặc cộng lặp.

6. Chuyển vị khung khác độ võng của rod/cable

Một đường chéo trên hình biến dạng vẫn nối hai nút đã di chuyển, kể cả khi lực dọc bằng 0. Nó không trở thành “không có chuyển vị” vì ngừng chịu nén. Hình vẽ cũng có thể dùng hệ số phóng đại lớn.

Với Frame, giới hạn dọc không tự loại bỏ khả năng chịu uốn, cắt và xoắn; xem CSI Analysis Reference Manual, Chapter VII — Nonlinear Properties, trang 133. Nếu rod vẫn xuất hiện mô men, kiểm tra release, tải ngang thân thanh và loại phần tử trước khi kết luận.

Hình 4 tách ba câu hỏi mô hình thường bị gộp làm một:

Theo CSI Knowledge Base — Cable tips, lực trong cable phụ thuộc hình dạng nên cần phân tích nonlinear và xác định đúng hình học ban đầu. Muốn xem độ võng thực, cần thêm tự trọng/tải phân bố, chiều dài không ứng suất hoặc thông tin căng trước, gối và điều kiện liên kết.

Không nên hiểu “chia nhỏ rod” là thao tác đủ để biến Frame thành cable. Với mô hình nhiều đoạn Frame để khảo sát uốn/mất ổn định, phải lựa chọn độ cứng uốn, hình học ban đầu và large displacement phù hợp. Với Cable object, số đoạn cần theo vị trí gối, điểm đặt tải và yêu cầu mô hình; kiểm tra hội tụ theo cách rời rạc hóa đang dùng.

Rod có lực căng trước có thể vẫn còn lực kéo khi tải ngang làm nó ngắn lại. Ví dụ “một thanh kéo, một thanh bằng 0” ở mục 3 không áp dụng tự động cho trạng thái đã căng trước.

7. Kiểm tra kết quả và hội tụ

Kết quả nonlinear chỉ có ý nghĩa khi mô hình đạt trạng thái cân bằng cần khảo sát. Kiểm tra theo thứ tự:

1. Mức tải đạt được: case đã hoàn thành toàn bộ tải mục tiêu hay dừng giữa chừng? Đọc log và cảnh báo, không chỉ mở hình bước cuối.

2. Cân bằng: tổng phản lực có phù hợp tổng tải tại trạng thái đang xem không? Với case tiếp nối, xét tổng trạng thái đã tích lũy.

3. Lực và trạng thái rod: đường chéo nào kéo, đường chéo nào chùng? Đảo tải có đổi đúng vai trò không? Lực nén rất nhỏ có thể do dung sai; lực nén đáng kể cần kiểm tra case và khai báo.

4. Mô men và chuyển vị: rod có vô tình tham gia chịu uốn? Hình biến dạng có dùng tỷ lệ phóng đại gây hiểu nhầm?

5. Độ nhạy: giảm bước tải hoặc siết dung sai có làm kết quả thay đổi đáng kể không? Với mô hình cable chi tiết, kiểm tra thêm độ nhạy lưới và hình học ban đầu.

CSI Help — Nonlinear Iteration giải thích phần mềm dùng lặp để kiểm tra cân bằng, có thể giảm bước tải khi không hội tụ; Newton–Raphson thường hiệu quả hơn với cable và phi tuyến hình học. Thay đổi tham số giải phải đi cùng kiểm chứng cơ học.

8. Chọn cách phân tích theo câu hỏi cần trả lời

• Dùng linear khi hệ làm việc đàn hồi trên nền độ cứng cố định phù hợp; đây vẫn là công cụ hiệu quả cho các bài toán thiết kế thông thường.

• Với rod chỉ kéo, bắt đầu bằng nonlinear static + giới hạn nén bằng 0, kiểm tra liên kết và giải cả hai chiều tải.

• Thêm P-Delta khi cần xét ảnh hưởng tải đứng lên đáp ứng ngang; dùng large displacement khi hình học biến dạng là cơ chế quan trọng.

• Muốn biết độ võng/chùng thực của cable, bổ sung mô hình hình học, tự trọng, lực căng ban đầu và kiểm tra hội tụ; không kết luận từ một đường Frame trên hình biến dạng.

• Giải từng trạng thái tải phi tuyến cần thiết rồi lấy bao; kiểm tra cân bằng, mức tải đạt được và trạng thái rod trước khi dùng nội lực để thiết kế.

Tài liệu tham khảo

• Tài liệu gốc do người dùng cung cấp: Phân tích kỹ thuật — Khi nào dùng linear và khi nào dùng nonlinear trong SAP, ngày 03/04/2026, tệp phan_tich_linear_vs_nonlinear_SAP_giang_rod.docx. Bài đã biên tập lại, bổ sung ví dụ tính tay và sơ đồ.

• CSI — Frame Tension/Compression Limits: điều kiện kích hoạt giới hạn kéo/nén.

• CSI — Load Case Data, Nonlinear Static Form: các ứng dụng và khai báo tải nonlinear static.

• CSI — Geometric Nonlinearity Parameters: None, P-Delta và large displacement.

• CSI — Initial Conditions: case tiếp nối, chia nhánh và tải tăng thêm.

• CSI — Stiffness to Use: phân tích linear trên độ cứng cuối case nonlinear.

• CSI Analysis Reference Manual, Rev. 15, July 2016, Chapter VII — The Frame Element, Nonlinear Properties, trang 133.

• CSI Knowledge Base — Cable tips: hình học và phân tích cable.

• CSI SAP2000 v26 Help — Nonlinear Iteration: kiểm soát lặp và hội tụ.