嵇康的专栏

1.complex.h#ifndef __COMPLEX_H #define __COMPLEX_Hclass CComplex{public: CComplex(); CComplex(double x,double y);  CComplex(const CComplex& other);  bool operator==(const CComplex

*Heading** Job name: doubleCrack Model name: Model-1** Generated by: Abaqus/CAE 6.12-1*Preprint, echo=NO, model=NO, history=NO, contact=NO**** PARTS***Part, name=crack1*End Part**

reference:http://wenku.baidu.com/link?url=KE3pz7qs2J5gxE2h-CR7rl0uIg3NyozFZVdOMTnCcOkamQxfIoHhATt6yEbWvR1qDbUm_FonuQh5CBSmtS9jhX85-WmbRuvknz3t2LHatI3网上很多关于一个节点的一维动态有限差分matlab程序，但是如果要做多个节点呢（也即是一维），

heat.mfunction heat_neil% Numerical Modeling in the Geosciences,  Humphrey% homework problem in heat flow% Mainly to get you working with Matrices% **************************************

function heat_neil_varK% program to solve temperature in the earth% note this does not use realistic values, it is mainly to % demonstrate correct usage of variable conductivity (K) layers% su

function heat_transit_1d_varKq=0.065;Z=20;K=20;               % conductivityTT=4;times=100;dt=TT/times;n=200;delz=Z/(n-1);z=0:delz:Z;A=zeros(n,n);    K2=0.3*K;for ii=2:n-1   

function heat_2d% Geol 5200, Humphrey 2013% 2D, steady state temperature.  Simplest approach% full (not sparse) matrices, with BCnodes included in matrixclear% solve temperature field in cru

function heat_transit_1dq=0.065;Z=20;K=20;               % conductivityTT=4;times=100;dt=TT/times;n=200;delz=Z/(n-1);z=0:delz:Z;A=zeros(n,n);    K2=0.3*K;    alpha=K*dt/(delz

% Geol 5200, Humphrey 2013% 2D, steady state temperature.  Simplest approach% full (not sparse) matrices, with BCnodes included in matrixclear% solve temperature field in crust with a constant

边界条件： |--------------------------------------------------------------|x=0m，0 摄氏度                                                                   x=100000m，q=K*dT/dx=0.014（表示地底下有稳定热源供给）参数：% h

mash显示用：https://www.pygimli.org/_examples_auto/1_meshing/plot_gmsh-example.html网格划分：https://pypi.org/project/pygmsh/fem后台计算：solidspy整合项目见：https://github.com/jimingkang/SolidPyGmsh...

1.生成网格：C:\Users\Yuan\Downloads\SolidsPy\tests\testpygmsh.py注意pygmsh.generate_mesh的方法默认是mesh4的，我们要生成mesh2的给Solidspy计算用，修改C:\Users\Yuan\Anaconda2\envs\env-freecad\Lib\site-packages\pygmsh\helpers....

staticFEM.javapackage scu.edu.fem;import java.math.BigInteger;class ETNode{                            //单元结点结构体    double x=0,y=0;                            //单元结点坐标    int number=0;

一受均布荷载的简支梁，跨度为3000mm,截面为100x200mm，配有两根钢筋，承受的均布载荷为0.04N/mm*mm。命令中采用了1/4模型，材料参数详见命令中。由于选择时采用的是实体号而不是坐标，可能在有些系统上会受到影响，各位注意。 !钢筋混凝土简支梁分析 /COM, Structural ---------定义单元及材料等--------------------- 

for i=1:10    a=2;    b=1;    c=2+0.1*(i-1);%left crack line%sqr_alpha_1=-(c-b)/(b+c);sqr_alpha=(c-b)/(b+c);sqr_k=(a-b)*(c-b)/((a+b)*(b+c));g=2/sqrt((a+b)*(b+c));v0=mfun('EllipticK',

int  rguass(int n,double a[],double b[]){ int *js,l,k,i,j,is,p,q; double d,t; js=(int *)malloc(n*sizeof(int)); l=1; for(k=0;k  d=0.0;  for(i=k;i   for(j=k;j    t=fabs(a[i*n+j]);

char tmp[100];memset(tmp,0,100);CString ss; for(int m=0;m{  sprintf(&tmp[10*m],"%10.5f,",pb[m]);//sprinf将double转成char* //ss+=}::MessageBox(NULL,tmp,"result is x=",2);}

1. 原书写作：int Matrix::GetColVector(int Col,double* pvector)const{ if(pvector==NULL)  delete pvector; pvector=new double[m_NumRows];//在里面分配 for (int j=0;j {  pvector[j]=GetElement(j,Col);

pv=[0,0;3,0;3,6;0,6;0,3.1;1.5,3;0,2.9;0,0][p,t]=distmesh2d(@dpoly,@huniform,0.02,[0,0; 3,6],pv,pv);

文件1：fem.h// Download by http://www.codefans.net#ifndef FEM__H#define FEM__H#include #include #include int Gauss(double a[],double b[],int n);    //全选主元高斯消去法/*double GaussIntegral(i

>>fh=@(p) 0.05+0.3*dellipse(p,[0.5,0.2]);>>fd=@(p) ddiff(drectangle(p,-10,10,-10,10),dellipse(p,[0.5,0.2]));>>[p,t]=distmesh2d(fd,fh,0.05,[-10,-10;10,10],[-10,-10;-10,10;10,-10;10,10]);

1.节点类若用模板则实现只能写在h里头    template     class MyNode{public:    T* data;    MyNode *pre;    MyNode *next;public :    MyNode(){            this->data=NULL;    }        MyNode(T *&

1.distmesh是一个较好的网格划分程序，具体可以参考：http://persson.berkeley.edu/distmesh/ 2.matlab有限元可以参考徐荣桥的书3.这里本人打算画一个园内包含一个椭圆的模型： 具体程序如下：a.网格划分： >>fh=@(p) 0.05+0.3*dellipse(p,[0.5,0.2]);>> fd=@(p) d

DynamicFDM.javapackage scu.edu.fem; class FDMETNode_2D{                            //单元结点结构体    double x=0,y=0;                            //单元结点坐标    int number=0;

Crack_DynamicFEM.javapackage scu.edu.fem; class Crack_DynamicETNode{                            //单元结点结构体    double x=0,y=0;                            //单元结点坐标    int number=0;

DynamicFEM.javapackage scu.edu.fem; class DynamicETNode{                            //单元结点结构体    double x=0,y=0;                            //单元结点坐标    int number=0;

mpi安装dowload./configure --prefix=/opt/mpichmakemake installsudo ln -s /opt/mpich/include /usr/local/include/mpipetsc 安装：./configure --with-mpi-dir=/opt/mpichmakemake install DES...