Heip Matlap kỹ thuật cơ khí

[jndvina]

Hệ thống quảng cáo SangNhuong.com

:adore:Tôi là dân kỹ thuật, nhưng nay mới đụng vào matlap nên mò hoài không ra. Nhờ Anh em help để được học hỏi thêm;
Tôi muốn vẽ đồ thị nội lực (lực cắt, độ võng, góc xoay, mô men) cho phần tử dầm
Ví dụ:
ví dụ Có dầm gồm 2 phần tử
- Tọa độ nút dầm, x=[0 100 300]
- Chiều dài các phần tử
L =[100 200]
- Mô men quán tính các phần tử
I =[5000000 75000000]
- Lực phân bố trên từng phần tử
q = 1.0e+006 *[0.7693 2.3079]
Tôi tính được
- Véc tơ chuyển vị và góc xoay:
d=[0 0.0007 0.0446 -0.0001 0 -0.0003]
- Véc tơ lực cắt
Qn =1.0e+008 *[-1.7284 0.8642 1.7284 -0.8642]
- Véc tơ mô men
Mn = 1.0e+010 *[ 0.0641 2.4977 -1.7925 -0.7693]

Tôi muốn vẽ đồ thị độ võng, lực cắt, góc xoay của dầm trên.

Tôi có tìm hiểu trên sách phần tử hữu hạn:
Hàm độ võng cho từng phần tử :
V(x)= [N]*{de}+q/24EIe(x^4-2Lex^3+Le^2x^2)
[N] là véc tơ hàm dạng
{de} là véc tơ chuyển vị nút phần tử
với [N]=[Nvi Npi Nvj Npj]
Nvi =1-(3x^2/Le^2)+2x^3/Le^3
Npi=x-(2x^2/Le)+x^3/Le^2
Nvj=(3x^2/Le^2)+2x^3/Le^3
Npj=-x^2/Le+x^3/Le^2
Hàm mô men là đạo hàm bậc 2 theo x của hàm độ võng
Mn=EI*(d^2V(x)/d(x^2))
Hàm lực cắt là đạo hàm bậc 3 theo x của hàm độ võng
Qn=EI*(d^3V(x)/d(x^3))

Tôi không biết vẽ đồ thị theo các hàm này hay vẽ theo kiểu nội suy theo các giá trị đã biết.
Do em còn lạ với mấy cái này, mà tìm tài liệu cũng không thấy nên nhờ anh em hướng dẫn.
Hoặc làm mẫu để học tập thêm.
Cảm ơn nhiều!
Có thể liên lạc:Bahung210@gmail.com