只有注册用户登录后才能阅读该文。 阅读全文
摘要：
以下面的问题为例，讨论FEM进行梁的模态分析：单跨梁受约束如图所示，梁的密度ρ=3e3kg/m³，试求梁的自振频率及振型。设梁的弹性模量为E=200GPa,横截面积A=1.0e-2㎡，长度L=2，界面惯性矩为I=1E-4m4。从题意，不考虑轴向拉伸，因此，刚度矩阵和质量矩阵为4*4的矩阵（对于6*6矩阵，见文章有限元法中，求解2D梁问题（http://www.cnblogs.com/heaventian/articles/2766227.html）下面给出其计算的程序：首先，给定单元划分，单元编号，节点的单元及整体编号的子程序coord2DBeam.mView Code % THE 阅读全文
摘要：
以下面的问题为例：对于如图所示的平面传热问题，若上端有给定的热流-2W/m2，即从下往上传输热量，结构下端有确定的温度100°，周围介质温度为20°，在两侧有换热，换热系数为α=100W/㎡/K，热导率200W/m/K,试分析其稳定温度场。使用下面的程序包进行求解：首先：将区域划分为4个单元，各单元包含的节点在element3.dat中显示：1 2 42 5 42 6 52 3 6每个节点的横纵坐标在coordinates.dat文件中显示：-10 00 010 0-5 100 105 10边界包含三个：恒温边界的节点编号，在dirichelet.dat文件中显示：1 22 阅读全文
摘要：
如果假定地表的温度由下述公式决定：Ts=T0+beta*h;其中，T0是0海拔处地表的平均温度，h为海拔，beta为空气的温度梯度，一般为-6.5K/km。则在地表各处温度时不相同的。这也导致了地表热流的不同。下面的程序用来计算地表地形有起伏，上表面温度由上述公式决定，下表面温度为1365°的等温面时，地表热流值随着地形的变化。使用http://www.cnblogs.com/heaventian/archive/2012/11/23/2784488.html文中的程序，来进行稳态温度场的计算。首先，代码coord3.m进行网格划分，单元，节点的编号与连接。View Code % t 阅读全文
摘要：
本文程序，使用有限元计算2D瞬态热传导方程。本程序只使用2D三角形单元求解，没有四边形单元求解。不过，也非常容易就可做到。主程序为:Heat_Conduction_transient.m大部分子程序和这篇文章相同，但是多了下面的初始温度场的子程序。u_init.m现将各个程序列于下：%function Heat_Conduction_transient( )%*****************************************************************************80%%% Applies the finite element method to 阅读全文
摘要：
加法中的误差传递：X=u±v则X的均方差为：σX=sqrt(σu^2+σv^2);乘法中的误差传递：除法中的误差传递：有限次幂的误差的传播：可以使用蒙特卡罗法来验证其误差：如下面的程序用来验证出发的误差：N=1e6;x=10+randn(N,1);y=5+randn(N,1)*2;std(x./y)mean(x./y) 阅读全文
摘要：
下面的有限差分程序计算2D稳态热传导方程。其模型和文章有限元计算2D稳态热传导方程的模型一样。从后面也可以看到，其结果也和该文章的结果一样。即有限差分和有限元法，解决同一问题，得到了相同的解。程序见下面：-----------------------------------------------------------------------------------% Solution of 2D steady heat conduction equation(Poisson equation): d2FI/dx2+d2FI/dy2=f % with finite differences o 阅读全文
摘要：
本文程序，使用有限元计算2D稳态热传导方程。本程序只使用2D三角形单元求解，没有四边形单元求解。不过，也非常容易就可做到。主程序为:Heat_Conduction_steady.m需要用到的子程序有：coord3.m 用来生成节点坐标，节点坐标与单元的链接。边界条件等等。生成文件coordinates.dat 包含两列，第一列节点x坐标，第二列，节点y坐标elements3.dat 包含三列：分别是各个单元的i,j,k节点的总体编号。dirichlet.dat 包含两列：分别是第一类边界条件的两端点节点的总体编号。另外，coord4.m与coord3相同，但是，生成的是四边形单元。即：... 阅读全文
摘要：
cited from wiki:http://en.wikipedia.org/wiki/Picrite_basaltPicrite basalt, picrobasalt or oceanite[citation needed] is a variety of high-magnesium olivine basalt that is very rich in the mineral olivine. It is dark with yellow-green olivine phenocrysts (20 to 50%) and black to dark brown pyroxene, m 阅读全文
摘要：
cited from:http://blog.sciencenet.cn/blog-4600-13085.html大多数期刊 均向作者免费发出校样(PROOFS)，供作者检查。经作者检查过的校样返回出版社后，经过排版公司对校样的更正和必要的技术处理，文章才可以在网上发 表。通过电子信箱传递的电子校样(PDF)不受地理和天气状况的 影响，只需几秒钟便可传递到联系作者(CORRESPONDINGAUTHOR)的手中。首先需要注意的是联系作者电子信箱的变动。一旦有变，要及时通知出版社。这种变动既包括联系作者因个人单位调动而导致的地址改变，也包括单位内计算机网络系统的更新换代和单位更名而造成的电邮地址 阅读全文