您现在的位置:首页 >> 技术文章 >> MATLAB技术 >> 内容

MATLAB代做|FPGA代做|python代做-无网格mlpg方法解决集中荷载悬臂梁

时间:2020-10-13 20:36:53 点击:

  核心提示:MATLAB代做|FPGA代做|python代做-无网格mlpg方法解决集中荷载悬臂梁...
% MLPG-main.m 
clear 
tic  
Lx=48;                                                               %������x����ı߳� 
Ly=12;                                                               %������y����ı߳� 
E=3e4;                                                               %����ģ�� 
nu=0.3;                                                              %���ɱ� 
q=-1;                                                                %����ֲ��� 
nx=21;                                                               %x����Ľڵ���� 
ny=11;                                                               %y����Ľڵ���� 
as=2.8;                                                              %֧����뾶  2.8 
dmax=1.5;                                                            %�ֲ�������뾶 1.8 
nc=[2 2];                                                            %���ֲ������򻮷ֵ��ӻ�����ĸ��� 
quado=4;                                                             %һά�����˹��ĸ��� 
Dmat=(E/(1-nu^2))*[1 nu 0;nu 1 0;0 0 (1-nu)/2];                      %ƽ��Ӧ�����ԸնȾ��� 
[xc,numnod,ds]=coordinate(Lx,Ly,nx,ny);                              %���ڵ�����꣬�ڵ������ͽڵ��� 
[k,f]=liu(xc,numnod,nc,quado,ds,dmax,Dmat,q,E,nu,as);                %����k�����f���� 
u=k\f;                                                               %�õ����нڵ��λ�� 
u2(2,numnod)=0; 
for i=1:numnod 
    u2(1,i)=u(2*i-1); 
    u2(2,i)=u(2*i); 
end 
figure                % ������y����Ľڵ�λ�ƺ;�ȷλ�ƵıȽ�ͼ 
plot(xc(1,(ny+1)/2:ny:numnod),u2(2,(ny+1)/2:ny:numnod),'r.'); 
%plot(xc(1,(ny+1)/2:ny:numnod),u2(2,1:ny:numnod),'r.');    % �ϱ߽�y����ڵ�λ�� 
hold on 
x=0:2:48; 
y=-(198*x+144*x.^2-x.^3)/(18*144e4); 
%y=-(0.9*36*(48-x)+198*x+144*x.^2-x.^3)/(18*144e4);        % �ϱ߽�y����ڵ㾫ȷλ��  
plot(x,y,'b'); 
xlabel('xc/m','fontweight','bold'); 
ylabel('u2/m','fontweight','bold'); 
legend('Uynode','Exact Solution'); 
hold off 
 
%________����x=L/2�������Ӧ���ֲ�ͼ,���γɽ��� 
nx1=2; ny1=10; 
I = Ly^3/12; 
xjm=zeros(2,ny1+1); 
stressex=zeros(1,ny+1); 
stress12ex=zeros(1,ny+1); 
for i=1:(ny1+1) 
    xjm(1,i) =Lx/nx1;                 % ������x���� 
    xjm(2,i) =-(Ly/ny1)*(i-1)+Ly;     % ������y���� 
    t = xjm(2,i)-Ly/2; 
    stressex(i)=-q*(Lx-xjm(1,i))*t/I; 
    stress12ex(i)=q*(Ly^2/4-t^2)/(2*I); 
end 
stress=zeros(3,ny1+1); 
for s=1:ny1+1  
    gpos=xjm(1:2,s); 
    gg = gpos(1:2); 
    v=supportdomain(numnod,gg,xc,ds,as); 
    L = length(v); 
    en = zeros(1,2*L); 
    [phi,dphix,dphiy]=shape(gg,as,xc,v,ds); 
    Bmat=zeros(3,2*L); 
    for j=1:L 
        Bmat(1:3,(2*j-1):2*j) = [dphix(j) 0;0 dphiy(j);dphiy(j) dphix(j)]; 
    end  
    for i=1:L 
        en(2*i-1) = 2*v(i)-1; 
        en(2*i)   = 2*v(i); 
    end 
    stress(1:3,s) = Dmat*Bmat*u2(en)'; 
end 
figure       % x=Lx/2���棬�ڵ�Ӧ���;�ȷӦ���ıȽ�ͼ 
hold on 
plot( xjm(2,1:(ny1+1)),stress(1,1:(ny1+1)),'rs');    % x�������Ӧ������ֵ 
plot( xjm(2,1:(ny1+1)),stressex(1,1:(ny1+1)),'.-');  % x�������Ӧ����ȷֵ 
legend('MLPG Solution','Exact Solution'); 
xlabel('y/m','fontweight','bold'); 
ylabel('Stress ','fontweight','bold'); 
hold off   
 
figure       % x=Lx/2���棬�ڵ�Ӧ���;�ȷӦ���ıȽ�ͼ 
hold on 
plot( xjm(2,1:(ny1+1)),stress(3,1:(ny1+1)),'rs');    % ��Ӧ������ֵ 
plot( xjm(2,1:(ny1+1)),stress12ex(1,1:(ny1+1)),'.-');  % ��Ӧ����ȷֵ 
legend('MLPG Solution','Exact Solution'); 
xlabel('y/m','fontweight','bold'); 
ylabel('Stress ','fontweight','bold'); 
hold off    
toc
 

联系:highspeedlogic

QQ :1224848052

微信:HuangL1121

邮箱:1224848052@qq.com

网站:http://www.mat7lab.com/

网站:http://www.hslogic.com/

微信扫一扫:

作者:无网格mlpg方法 来源:无网格mlpg方法
  • 您是如何找到本站的?
  • 百度搜索
  • Google搜索
  • 查阅资料过程中
  • 论坛发现
  • 百度贴吧发现
  • 朋友介绍
本站最新成功开发工程项目案例
相关文章
  • 没有相关文章
相关评论
发表我的评论
  • 大名:
  • 内容:
  • matlab代做|matlab专业代做|matlab淘宝代做|python人工智能代做|FPGA项目合作(www.hslogic.com) © 2020 版权所有 All Rights Reserved.
  • Email:highspeed_logic@163.com 站长QQ: 1224848052

    专业代做/代写/承接、MATLAB、SIMULINK、FPGA项目、博士/硕士/本科毕业设计、课题设计、论文,毕业论文,Coursework、Eassy、Assignment