偏微分方程的有限元法求解.pdf

上传人:l**** 文档编号:71293358 上传时间:2023-02-01 格式:PDF 页数:2 大小:27.92KB
返回 下载 相关 举报
偏微分方程的有限元法求解.pdf_第1页
第1页 / 共2页
偏微分方程的有限元法求解.pdf_第2页
第2页 / 共2页
亲,该文档总共2页,全部预览完了,如果喜欢就下载吧!
资源描述

《偏微分方程的有限元法求解.pdf》由会员分享,可在线阅读,更多相关《偏微分方程的有限元法求解.pdf(2页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、-%对于 d u/d*=f 的 FEM 解算器,其中 f=*(1-*)%边界条件 u(0)=0,u(1)=0.%精确解用以比对*=linspace(0,1,101);%产生 0-1 之间的均分指令,101 为元素个数ue*=(1/6).*.3-(1/12).*.4-(1/12).*;%对力项设置高斯点的数目NGf=2;if(NGf=2)*iGf=-1/sqrt(3);1/sqrt(3);%1、2的值aGf=1 1;else,NGf=1;*iGf=0.0;aGf=2.0;end%单元数目Ne=5;%建立网格节点*=linspace(0,1,Ne+1);%零刚性矩阵K=zeros(Ne+1,Ne+

2、1);b=zeros(Ne+1,1);%对所有单元循环计算刚性和残差for ii=1:Ne,kn1=ii;kn2=ii+1;*1=*(kn1);*2=*(kn2);d*=*2-*1;%每一个单元的长度d*id*=2/d*;%d/d*d*d*i=1/d*id*;%d*/ddN1d*i=-1/2;%d1/ddN2d*i=1/2;%d2/ddN1d*=dN1d*i*d*id*;%-1/(*j-*j-1)dN2d*=dN2d*i*d*id*;%1/(*j-*j-1)K(kn1,kn1)=K(kn1,kn1)-2*dN1d*dN1d*d*d*i;%Rj的第二项K(kn1,kn2)=K(kn1,kn2)-

3、2*dN1d*dN2d*d*d*i;K(kn2,kn1)=K(kn2,kn1)-2*dN2d*dN1d*d*d*i;K(kn2,kn2)=K(kn2,kn2)-2*dN2d*dN2d*d*d*i;%用高斯积分估计力项的积分for nn=1:NGf%NGf=2*iG=*iGf(nn);%得到高斯点的N1=0.5*(1-*iG);%求 N1 和 N2(即在*iG 的权重/插值)形状函数在的值.z.22-N2=0.5*(1+*iG);%值fG=*iG*(1-*iG);%对点求 fgG1=N1*fG*d*d*i;%在节点处估计权函数在高斯点的被积函数gG2=N2*fG*d*d*i;%估计是个积分值b(

4、kn1)=b(kn1)+aGf(nn)*gG1;%aGf为 1b(kn2)=b(kn1)+aGf(nn)*gG2;endend%在*=0 处设置 Dirichlet 条件kn1=1;K(kn1,:)=zeros(size(1,Ne+1);K(kn1,kn1)=1;b(kn1)=0;%在*=1 处设置 Dirichlet 条件kn1=1;K(kn1,:)=zeros(size(1,Ne+1);K(kn1,kn1)=1;b(kn1)=0;%求解方程v=Kb;%v 为 K*=b 的解plot(*,v,*-);%画图并比较hold on;plot(*,ue*);hold off;*label(*);ylabel(u);.z.

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 应用文书 > 工作报告

本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知淘文阁网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

工信部备案号:黑ICP备15003705号© 2020-2023 www.taowenge.com 淘文阁