力学数学研究生-CSDN博客

原创非线性有限元

非线性有限元#include <deal.II/base/quadrature_lib.h>#include <deal.II/base/function.h>#include <deal.II/base/tensor.h>#include <deal.II/lac/vector.h>#include <deal.II/lac/full_matrix.h>#include <deal.II/lac/sparse_matrix.h&

2022-05-26 20:35:15 382

原创晶界研究程序

matlab获得蜂窝六边形程序clearclc%%%下面得到六边形的蜂窝形状。%这里需要输出一个文件,data=zeros(n,8),其数据格式如下；%n代表晶界的数量。%第一到二列：表示晶界第一个端点的横纵坐标，即初始点。%第三到四列：表示晶界第二个端点的横纵坐标，即被扩展点。%第五到六列：表示晶界的平移向量。%第七列：表示晶界的方向。%第八到九列：表示晶界的旋转正余弦值。%以三角形的方式向外扩展。%输入六边形边长；num=7;a=10^(-8);center=[0,0];

2022-02-17 23:40:08 326 1

原创安装pyPENELOPE

pip3 install Numpypip3 install Matplotlibpip3 install pyparsingpip install -U -f https://extras.wxpython.org/wxPython4/extras/linux/gtk3/centos-8 wxPython

2021-09-13 10:51:46 359 3

原创 2021-07-23

https://blog.csdn.net/qyb19970829/article/details/105451195?utm_source=app&app_version=4.11.0&code=app_1562916241&uLinkId=usr1mkqgl919blenlammps安装链接

2021-07-23 09:23:49 148

原创 dealii软件的链接库安装

安装MPICH./configure --prefix=/softer/mpich --with-device=ch4:ofimake make install 添加环境变量export PATH=/softer/mpi/bin:$PATHexport INCLUDE=/softer/mpi/include:$INCLUDEexport LD_LIBRARY_PATH=/softer/mpi/lib:$LD_LIBRARY_PATH配置pestcexport PETSC_DIR=`pw

2021-05-29 13:18:12 1279 4

原创 deal-II计算刃位错周围弹性场程序与数据处理程序

#include <deal.II/base/quadrature_lib.h>#include <deal.II/base/function.h>#include <deal.II/base/tensor.h>#include <deal.II/lac/vector.h>#include <deal.II/lac/full_matrix.h>#include <deal.II/lac/sparse_matrix.h>#in

2021-03-18 18:40:51 308

原创本征应变问题有限元弱形式推导（Maple）

二维情况下> restart:> with(linalg):> C:=array(1..3,1..3,[[2*mu+lambda,lambda,0],[lambda,2*mu+lambda,0],[0,0,mu]]);> B:=array(1..3,1..8,sparse):> NN1:=N1(x1,x2):> NN2:=N2(x1,x2):> NN3:=N3(x1,x2):> NN4:=N4(x1,x2):> for i from 1

2021-03-18 18:39:11 821

原创推导Eshelby夹杂解析解的maple程序

计算内点的弹性场#计算内点:restart:> with(linalg):> assume(a>0):> a1:=a:> a2:=a:> a3:=a:> engstrain:=array(1..3,1..3,sparse):> for i from 1 to 3 do:> engstrain[i,i]:=10^(-6);> od: > deta:=array(1..3,1..3,sparse):> for i f

2021-03-18 18:36:49 521

原创 deal-II软件计算Eshelby夹杂程序，以及matlab处理数据程序

deal-ii程序#include <deal.II/base/quadrature_lib.h>#include <deal.II/base/function.h>#include <deal.II/base/tensor.h>#include <deal.II/lac/vector.h>#include <deal.II/lac/full_matrix.h>#include <deal.II/lac/sparse_matrix

2021-03-18 18:34:39 484

原创有限差分法解第二类边界条件的椭圆形偏微分方程

有限差分法解第二类边界条件的椭圆形偏微分方程clearclcax=[0,1];ay=[0,1];x1=0.1;xh=ax(1):x1:ax(2);%保持步长相同yh=ay(1):x1:ay(2);n1=length(xh);%下面输入第二类边界条件，即neumann边界条件,为了保证唯一解，下面赋予函数在原点的初值为0；syms x yun=x+y;A1=zeros(n1,n1);A2=zeros(n1,n1);I=zeros(n1,n1);B=zeros(n1,n1);fo

2020-11-18 15:30:48 3747 3

原创用三角函数做最佳平方逼近（C++实现）

用三角函数做最佳平方逼近（C++实现）#include<iostream>#include<cstdio> #include<cmath>#define pi 3.1415926535using namespace std;double lim(int num,double x[],double y[],double t,int aa){ //这里判断num和aa为奇数还是偶数； int a[2]; int a1=num/2; int a2=aa/2

2020-10-28 17:18:30 1224 1

原创 LU矩阵分解法解线性方程组（C++实现）

LU矩阵分解法解线性方程组（C++实现）#include<cstdio>#include<iostream>using namespace std;void fenjie(double a[],double b[],int n){ double aa[n][n]; for(int i=0;i<n;i++){ for(int j=0;j<n;j++){ aa[i][j]=a[n*i+j]; } } double l[n][n]; do

2020-10-28 17:14:55 1627

原创列主元消去法解线性方程组（C++实现）

列主元消去法解线性方程组（C++实现）#include<cstdio>#include<iostream>using namespace std;int quen(double a[],int n){ double x; int k; for(int i=0;i<n;i++){ if(a[i]>=0){ a[i]=a[i]; }else{ a[i]=-a[i]; } } x=a[0]; int p=0;

2020-10-28 17:13:48 1668 2

原创高斯消元法解线性方程组（C++实现）

高斯消元法解线性方程组#include<cstdio>#include<iostream>using namespace std;void gaosi(double a[],double b[],int n){ double aa[n][n]; for(int i=0;i<n;i++){ for(int j=0;j<n;j++){ aa[i][j]=a[n*i+j]; } } double A[n][n],B[n]; for (int i=

2020-10-28 17:12:14 2531 4

原创加权六点格式求解抛物型偏微分方程（Matlab实现）

加权六点格式求解抛物型偏微分方程clearclcat=[0,1];ax=[0,1];t1=0.01;x1=0.01;theta=0.5;th=at(1):t1:at(2);xh=ax(1):x1:ax(2);n1=length(xh);n2=length(th);u=zeros(n2,n1);syms x t;bot=20;%时间边界boxr=0;%空间左边界boxl=100;%空间右边界%下面求边界值bt=zeros(1,n1);for i=1:n1 bt(i)

2020-10-28 16:37:49 1725 1

原创古典显示格式求解抛物型偏微分方程（Matlab实现）

古典显示格式求解抛物型偏微分方程clearclcat=[0,1];ax=[0,1];t1=0.01/6;x1=0.1;th=at(1):t1:at(2);xh=ax(1):x1:ax(2);n1=length(xh);n2=length(th);u=zeros(n2,n1);syms x t;bot=sin(pi*x);boxr=0;boxl=0;%下面求边界值bt=zeros(1,n1);for i=1:n1 bt(i)=subs(bot,x,xh(i));

2020-10-28 16:35:39 4117 1

原创解常微分方程的多步法（有admas外插法，admas内插法，一般多步法（Milne法））

解常微分方程的多步法（有admas外插法，admas内插法，一般多步法（Milne法））admas外插法%admas外插法，具有k-1阶的精度,下面编写k等于3的外插公式clearclcsyms x yf=3*x^2;a=0;b=10;y0=0;y1=0.001;y2=0.008;y3=0.027;h=0.1;n=(b-a)/h+1;xx=zeros(n,1);yy=zeros(n,1);yy=zeros(n,1);yy(1)=y0;yy(2)=y1;yy(3)=y2

2020-09-27 20:32:15 4716 1

原创微分方程数值解法的matlab程序

微分方程数值解法的matlab程序欧拉法clearclcsyms x yf=2*x;a=0;b=10;y0=0;h=0.1;n=(b-a)/h+1;xx=zeros(n,1);yy=zeros(n,1);for i=1:n xx(i)=a+h*(i-1); endyy=zeros(n,1);yy(1)=y0;for i=2:n yy(i)=yy(i-1)+h*subs(f,{x,y},{xx(i-1),yy(i-1)}); endff1=[xx,yy];

2020-09-24 16:11:45 1493

原创三次样条插值法的C++程序实现（追赶法求解线性方程组）

三次样条插值法的C++程序实现（追赶法求解线性方程组）求解线性方程组追赶法的程序为：double capway(int n,double nu[],double d[],double lan[],double h[],double x[],double y[], double t){ double beita[n-1]; beita[0]=lan[0]/2; for (int i=1;i<n-1;i++){ beita[i]=lan[i]/(2-nu[i-1]*b

2020-09-20 09:49:03 3614 4

原创秦九邵算法计算多项式（C语言实现）

秦九邵算法计算多项式程序（C语言实现）#include<iostream>#include<cstdio>using namespace std;double algr(int n,double a[],double x){ double b[n]; b[0]=a[0]; for (int i=1;i<=n-1;i++){ b[i]=b[i-1]*x+a[i]; } printf("%f",b[n-1]);}int main(){

2020-09-19 22:48:08 1727

原创牛顿插值法的C++程序实现

标题使用C++编写牛顿插值法程序（可以输出差商表）#include<iostream>#include<cstdio> using namespace std;double aver(int num,double x[],double y[],double aa){ double a[num]; a[0]=y[0]; cout<<"下面输出差商表"<<"\n"; for (int i=1;i<=num-1;i++){

2020-09-19 22:31:10 5521 1

原创拉格朗日插值法的C++程序实现

使用C++编写拉格朗日插值法程序程序如下：在这里#include<cstdio> #include<iostream> using namespace std; double function(int num,double x[],double y[],double a){ double fund[num]; for (int i=0;i<num;i++){ double sum1=1,sum2=1; for (int j=0;j<num;j++){ if(

2020-09-19 22:25:26 5986 3

py753159123456的博客