24小时热门版块排行榜    

查看: 1498  |  回复: 7

JIM080

新虫 (初入文坛)

[求助] C++参数化建模程序出错,谁能帮忙看下错在哪里? 已有4人参与

如题。应该是程序格式错误,算法没有错误,可以忽略。错误是
fatal error C1083: Cannot open include file: 'iostream': No such file or directory
执行 cl.exe 时出错.程序如下:
CODE:
//程序头,不同版本的 c++软件,程序头有细微差别
#include<iostream>
using namespace std;
#include<fstream>
using namespace std;
#include<stdlib.h>
#include<math.h>
long double JDX[1000][1000],JDY[1000][1000],JDZ[1000][1000],JDDX[1000][1000],JDDY[1000][1000],JDDZ[1000][1000];             //定义数据类型
double PI=4*atan(1.0);
void main()
{
        ofstream fileout;                                           //定义输入输出文本
fileout.open("resultfile.dat");
        if(!fileout)
{
cerr<<"产生错误!无法打开testfile文件!"<<endl;
exit(1);
}
ofstream filein;
filein.open("testfile.dat");
if(!filein)
{
cerr<<"产生错误!无法打开testfile文件!"<<endl;
exit(1);
}
int m,n;
int i,j,p=1,k=1,e=1,g=1;
long M;
long double a,b,c,d,T,s,V;
double L,B,R,L1,L2,L12,q,H,Hc,Hba,Hcbh;
    filein>>L>>B>>R>>L1>>L2>>L12>>q>>Hba>>Hcbh;                 //输入蜂窝夹层壳参数
filein>>T;
    filein>>s;
    filein.close();
V=2*L*Hba*(1-Hcbh);
    H=2*L*Hba;
Hc=2*L*Hba*Hcbh;
L12=(L12/180)*PI;
c=L-L1/2-L2*sin(L12);
d=2*L1+2*L2*sin(L12);
m=c/d;
n=(B-L2*cos(L12))/(2*L2*cos(L12));
a=L-(L1/2+L2*sin(L12))-m*(2*L1+2*L2*sin(L12));
b=B-L2*cos(L12)-n*2*L2*cos(L12);
fileout<<"/PREP7"<<endl;                                    //定义材料属性
fileout<<"ET"<<","<<1<<","<<"shell63"<<endl;
fileout<<"T"<<","<<1<<","<<T<<", , , , , ,"<<endl;
    fileout<<"T"<<","<<2<<","<<2*T<<", , , , , ,"<<endl;                           //******注意生成命令流后,在命令流文本中把 T 改成 R ******
fileout<<"T"<<","<<3<<","<<V<<", , , , , ,"<<endl;
fileout<<"MP"<<","<<"EX"<<","<<1<<","<<68.98<<"e"<<9<<endl;
fileout<<"MP"<<","<<"PRXY"<<","<<1<<","<<0.33<<endl;
fileout<<"MP"<<","<<"DENS"<<","<<1<<","<<2768<<endl;
fileout<<"MP"<<","<<"EX"<<","<<2<<","<<68.98<<"e"<<9<<endl;
fileout<<"MP"<<","<<"PRXY"<<","<<2<<","<<0.33<<endl;
fileout<<"MP"<<","<<"DENS"<<","<<2<<","<<2768<<endl;
    if((b>0)&&(b<=L2*cos(L12))&&(a>0)&&(a<=L1))                  //在此种边界情况下的上蒙皮节点坐标
{
        for(i=0;i<2*m+3;i++)
{
if(i==0)
{
for(j=0;j<2*n+3;j++)
{
JDX[i][j]=0;
JDY[i][j]=R*sin(j*L2*cos(L12)/R);
JDZ[i][j]=sqrt(R*R-JDX[i][j]*JDX[i][j]-JDY[i][j]*JDY[i][j]);
if(j==2*n+2)
{
                        JDX[i][j]=0;
   JDY[i][j]=R*sin(((j-1)*L2*cos(L12)+b)/R);
JDZ[i][j]=sqrt(R*R-JDX[i][j]*JDX[i][j]-JDY[i][j]*JDY[i][j]);
}
}
}
else if(i==2*m+2)
{
for(j=0;j<2*n+3;j++)
{
JDX[i][j]=R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R);
JDY[i][j]=R*sin(j*L2*cos(L12)/R);
JDZ[i][j]=sqrt(R*R-JDX[i][j]*JDX[i][j]-JDY[i][j]*JDY[i][j]);
if(j==2*n+2)
{
                         JDX[i][j]=R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R);
    JDY[i][j]=R*sin(((2*n+1)*L2*cos(L12)+b)/R);
JDZ[i][j]=sqrt(R*R-JDX[i][j]*JDX[i][j]-JDY[i][j]*JDY[i][j]);
}
}
}
            else if(i%2==1)
{
for(j=0;j<2*n+3;j++)
{
                    if(j==2*n+2)
{
JDX[i][j]=R*sin(((i-1)/2*(2*L1+2*L2*sin(L12))+L1/2.0+b*tan(L12))/R);
                        JDY[i][j]=R*sin(((2*n+1)*L2*cos(L12)+b)/R);
JDZ[i][j]=sqrt(R*R-JDX[i][j]*JDX[i][j]-JDY[i][j]*JDY[i][j]);

}
   else if(j%2==0)
{
JDX[i][j]=R*sin(((i-1)/2*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12))/R);
JDY[i][j]=R*sin(j*L2*cos(L12)/R);
JDZ[i][j]=sqrt(R*R-JDX[i][j]*JDX[i][j]-JDY[i][j]*JDY[i][j]);
}
else if(j%2==1)
{
JDX[i][j]=R*sin(((i-1)/2*(2*L1+2*L2*sin(L12))+L1/2.0)/R);
JDY[i][j]=R*sin(j*L2*cos(L12)/R);
JDZ[i][j]=sqrt(R*R-JDX[i][j]*JDX[i][j]-JDY[i][j]*JDY[i][j]);
}
}
}
else if(i%2==0)
{
for(j=0;j<2*n+3;j++)
{
if(j==2*n+2)
{
JDX[i][j]=R*sin(((i-2)/2*(2*L1+2*L2*sin(L12))+3.0/2.0*L1+2*L2*sin(L12)-b*tan(L12))/R);
JDY[i][j]=R*sin(((2*n+1)*L2*cos(L12)+b)/R);
JDZ[i][j]=sqrt(R*R-JDX[i][j]*JDX[i][j]-JDY[i][j]*JDY[i][j]);
}
else if(j%2==0)
{
JDX[i][j]=R*sin(((i-2)/2*(2*L1+2*L2*sin(L12))+3.0/2.0*L1+L2*sin(L12))/R);
JDY[i][j]=R*sin(j*L2*cos(L12)/R);
JDZ[i][j]=sqrt(R*R-JDX[i][j]*JDX[i][j]-JDY[i][j]*JDY[i][j]);
}
else if(j%2==1)
{
JDX[i][j]=R*sin(((i-2)/2*(2*L1+2*L2*sin(L12))+3.0/2.0*L1+2*L2*sin(L12))/R);
JDY[i][j]=R*sin(j*L2*cos(L12)/R);
JDZ[i][j]=sqrt(R*R-JDX[i][j]*JDX[i][j]-JDY[i][j]*JDY[i][j]);
}
}
}
}
for(i=0;i<2*m+3;i++)                                         //输出节点号,节点坐标
{
if(i==0)
{
for(j=0;j<2*n+3;j++)
{
                    fileout<<"N"<<","<<"  "<<k<<","<<JDX[i][j]<<","<<JDY[i][j]<<","<<JDZ[i][j]<<endl;
k++;
}
}
else if(i==2*m+2)
{
for(j=0;j<2*n+3;j++)
{
                    fileout<<"N"<<","<<"  "<<k<<","<<JDX[i][j]<<","<<JDY[i][j]<<","<<JDZ[i][j]<<endl;
k++;
}
}
else
{
for(j=0;j<2*n+3;j++)
{
        fileout<<"N"<<","<<"  "<<k<<","<<JDX[i][j]<<","<<JDY[i][j]<<","<<JDZ[i][j]<<endl;
k++;
}
}
}
for(p=1;p<q+1;p++)                                            //输出下几层节点号以及节点坐标
  {
for(i=0;i<2*m+3;i++)
{
if(i==0)
{
for(j=0;j<2*n+3;j++)
{
                    JDDX[i][j]=((R-(p*Hc/q))/R)*JDX[i][j];
JDDY[i][j]=((R-(p*Hc/q))/R)*JDY[i][j];
                    JDDZ[i][j]=((R-(p*Hc/q))/R)*JDZ[i][j];
                    fileout<<"N"<<","<<"  "<<k<<","<<JDDX[i][j]<<","<<JDDY[i][j]<<","<<JDDZ[i][j]<<endl;
k++;
}
}
else if(i==2*m+2)
{
for(j=0;j<2*n+3;j++)
{
                    JDDX[i][j]=((R-(p*Hc/q))/R)*JDX[i][j];
JDDY[i][j]=((R-(p*Hc/q))/R)*JDY[i][j];
                    JDDZ[i][j]=((R-(p*Hc/q))/R)*JDZ[i][j];
                    fileout<<"N"<<","<<"  "<<k<<","<<JDDX[i][j]<<","<<JDDY[i][j]<<","<<JDDZ[i][j]<<endl;
k++;
}
}
else
{
for(j=0;j<2*n+3;j++)
{
   JDDX[i][j]=((R-(p*Hc/q))/R)*JDX[i][j];
JDDY[i][j]=((R-(p*Hc/q))/R)*JDY[i][j];
                    JDDZ[i][j]=((R-(p*Hc/q))/R)*JDZ[i][j];
                    fileout<<"N"<<","<<"  "<<k<<","<<JDDX[i][j]<<","<<JDDY[i][j]<<","<<JDDZ[i][j]<<endl;
k++;
}
}
}
}
M=(2*m+3)*(2*n+3);                                                           
fileout<<"TYPE"<<","<<1<<endl;                                          //定义蒙皮材料属性,生成蒙皮单元
fileout<<"MAT"<<","<<1<<endl;
fileout<<"REAL"<<","<<3<<endl;
for(g=1;g<=q+1;g++)
{
if(g==1)
   {
        for(i=1;i<=2*m+2;i++)
{
            for(j=(i-1)*(2*n+3)+1;j<i*(2*n+3);j++)
{
                fileout<<"E"<<","<<"  "<<(g-1)*M+j<<","<<(g-1)*M+j+2*n+3<<","<<(g-1)*M+j+2*n+4<<","<<(g-1)*M+j+1<<endl;
e++;
}
}
}
if(g==q+1)
   {
        for(i=1;i<=2*m+2;i++)
{
            for(j=(i-1)*(2*n+3)+1;j<i*(2*n+3);j++)
{
                fileout<<"E"<<","<<"  "<<(g-1)*M+j<<","<<(g-1)*M+j+2*n+3<<","<<(g-1)*M+j+2*n+4<<","<<(g-1)*M+j+1<<endl;
e++;
}
}
}
  }
        fileout<<"TYPE"<<","<<1<<endl;                                           //定义芯层材料属性,生成芯层单元(单壁厚单元)
fileout<<"MAT"<<","<<2<<endl;
fileout<<"REAL"<<","<<1<<endl;
for(g=1;g<=q;g++)
{
for(i=1;i<=2*m+2;i++)
{
for(j=(i-1)*(2*n+3)+1;j<i*(2*n+3);j++)
{
fileout<<"E"<<","<<"  "<<(g-1)*M+j<<","<<(g-1)*M+j+1<<","<<j+g*M+1<<","<<j+g*M<<endl;
e++;
}
}
}
        fileout<<"TYPE"<<","<<1<<endl;                                               //定义芯层材料属性,生成芯层单元(双壁厚单元)
fileout<<"MAT"<<","<<2<<endl;
fileout<<"REAL"<<","<<2<<endl;
for(g=1;g<=q;g++)
{
for(i=1;i<=2*m+2;i++)
{
if(i==1)
{
for(j=2;j<=2*n+2;j++)
{
if(j%2==0)
{
fileout<<"E"<<","<<"  "<<(g-1)*M+j<<","<<(g-1)*M+j+2*n+3<<","<<j+g*M+2*n+3<<","<<j+g*M<<endl;
j++;
e++;
}
}
}
else if(i==2*m+2)
{   
                for(j=(i-1)*(2*n+3)+1;j<i*(2*n+3);j++)
{
if((j-(i-1)*(2*n+3))%2==1)
{
   
fileout<<"E"<<","<<"  "<<(g-1)*M+j<<","<<(g-1)*M+j+2*n+3<<","<<j+g*M+2*n+3<<","<<j+g*M<<endl;
e++;
}
}
}
else if(i%2==0)
{
for(j=(i-1)*(2*n+3)+1;j<i*(2*n+3);)
{
                    
fileout<<"E"<<","<<"  "<<(g-1)*M+j<<","<<(g-1)*M+j+2*n+3<<","<<j+g*M+2*n+3<<","<<j+g*M<<endl;
j=j+2;
e++;
}
}
else if(i%2==1)
{
for(j=(i-1)*(2*n+3)+1;j<i*(2*n+3)-1;)
{
fileout<<"E"<<","<<"  "<<(g-1)*M+j+1<<","<<(g-1)*M+j+2*n+4<<","<<j+g*M+2*n+4<<","<<j+g*M+1<<endl;
j=j+2;
e++;
}
}
}
}

   fileout<<"NSYM"<<","<<"Y"<<","<<k-1<<","<<"ALL"<<endl;
fileout<<"ESYM"<<","<<" "<<","<<k-1<<","<<"ALL"<<endl;
         fileout<<"NSYM"<<","<<"X"<<","<<2*(k-1)<<","<<"ALL"<<endl;
fileout<<"ESYM"<<","<<" "<<","<<2*(k-1)<<","<<"ALL"<<endl;
        fileout<<"NUMMRG,NODE, , , ,LOW"<<endl;                                     //压缩节点
fileout<<"NUMCMP,NODE"<<endl;
fileout<<"FINISH"<<endl;
if(s==1)                                                                    //对边间支,对边固支边界
{
fileout<<"/SOL"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<((R-(Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<((R-(2*Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<((R-(3*Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<-R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<-(((R-(Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<-(((R-(2*Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<-(((R-(3*Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
        fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<R*sin(((2*n+1)*L2*cos(L12)+b)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<((R-(Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<((R-(2*Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<((R-(3*Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<-R*sin(((2*n+1)*L2*cos(L12)+b)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<-(((R-(Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<-(((R-(2*Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<-(((R-(3*Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"ALLSEL,ALL"<<endl;
}
if(s==2)                                                                   //对边间支边界
{
fileout<<"/SOL"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<((R-(Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<((R-(2*Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<((R-(3*Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<-R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<-(((R-(Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<-(((R-(2*Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<-(((R-(3*Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ, , , ,"<<endl;
        fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<R*sin(((2*n+1)*L2*cos(L12)+b)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<((R-(Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<((R-(2*Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<((R-(3*Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<-R*sin(((2*n+1)*L2*cos(L12)+b)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<-(((R-(Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<-(((R-(2*Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<-(((R-(3*Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ, , , ,"<<endl;
fileout<<"ALLSEL,ALL"<<endl;
}
if(s==3)                                                                          //四边固支边界
{
fileout<<"/SOL"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<((R-(Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<((R-(2*Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<((R-(3*Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<-R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<-(((R-(Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<-(((R-(2*Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"x"<<","<<-(((R-(3*Hc/q))/R)*R*sin((m*(2*L1+2*L2*sin(L12))+L1/2.0+L2*sin(L12)+a)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UX,UZ,ROTY, , ,"<<endl;
        fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<R*sin(((2*n+1)*L2*cos(L12)+b)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ,ROTX, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<((R-(Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ,ROTX, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<((R-(2*Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ,ROTX, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<((R-(3*Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ,ROTX, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<-R*sin(((2*n+1)*L2*cos(L12)+b)/R)<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ,ROTX, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<-(((R-(Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ,ROTX, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<-(((R-(2*Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ,ROTX, , ,"<<endl;
fileout<<"nsel"<<","<<"s"<<","<<"loc"<<","<<"y"<<","<<-(((R-(3*Hc/q))/R)*R*sin(((2*n+1)*L2*cos(L12)+b)/R))<<endl;
fileout<<"D"<<","<<"ALL"<<", , , , , ,UY,UZ,ROTX, , ,"<<endl;
fileout<<"ALLSEL,ALL"<<endl;
}
fileout<<"EPLOT"<<endl;
fileout<<"ANTYPE,2"<<endl;
fileout<<"MSAVE,0"<<endl;
fileout<<"MODOPT,LANB,30"<<endl<<"EQSLV,SPAR"<<endl<<"MXPAND,30, , ,0"<<endl;
fileout<<"LUMPM,0"<<endl<<"PSTRES,0"<<endl<<"MODOPT,LANB,30,0,10000000, ,OFF"<<endl;
fileout<<"/STATUS,SOLU"<<endl<<"SOLVE"<<endl;                                       //ansys开始分析计算
fileout.close();
}
}



[ Last edited by jjdg on 2014-2-23 at 00:46 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色


jjdg: 金币+1, 感谢参与 2014-02-23 00:46:48
   把第一句换成 iostream.h 试试
我同意叔本华的观点,人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活,逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
2楼2014-02-21 15:48:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖


感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与 2014-02-23 00:46:57
项目配置问题,连iostream都找不到,你是用的visual studio吗?还是其他编辑器
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
3楼2014-02-21 18:32:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

JIM080

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by baobiao007 at 2014-02-21 15:48:51
   把第一句换成 iostream.h 试试

换了 好像不行啊
4楼2014-02-22 17:18:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

JIM080

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by libralibra at 2014-02-21 18:32:39
项目配置问题,连iostream都找不到,你是用的visual studio吗?还是其他编辑器

就是Visual C++6.0我不太懂 这是我毕设里的一个部分 我不是学编程的 只是大一学过C语言 老师给我让我把错误改了 难为我了。
5楼2014-02-22 17:23:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yaozhq

金虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与 2014-02-23 00:47:06
提示你找不到头文件 如果你调整一下顺序可能会发现后面的也找不到 那应该是你的编译环境有问题
6楼2014-02-23 00:25:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rbs

木虫 (小有名气)

【答案】应助回帖

这是编译器开发环境的问题,把VC6重装一遍。不过推荐你使用VS2008以后版本,或者Dev c++或者code blocks。总之不要用VC6。
7楼2014-03-25 17:58:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jmppok

新虫 (初入文坛)

【答案】应助回帖

新建一个VC6的 控制台程序, 将该代码粘贴进去,编译试试看。

找不到iostream,要么是vc没装好,要么是工程设置不对。
8楼2014-03-26 09:56:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 JIM080 的主题更新
信息提示
请填处理意见