搜索
返回列表 发新帖
查看: 853|回复: 8
打印 上一主题 下一主题

C++参数化建模程序出错，谁能帮忙看下错在哪里？

[复制链接]
acer080
跳转到指定楼层
楼主
发表于 2014-2-21 12:53:25 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
5啊哈币
如题。应该是程序格式错误，算法没有错误，可以忽略。错误是
fatal error C1083: Cannot open include file: 'iostream': No such file or directory
执行 cl.exe 时出错.程序如下：
//程序头，不同版本的 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();
}
}

回复

举报

rosynirvana
沙发
发表于 2014-2-21 15:43:53 | 只看该作者
fatal error C1083: Cannot open include file: 'iostream': No such file or directory

这明显就是编译环境的问题吧
回复

举报

acer080
板凳
 楼主| 发表于 2014-2-22 17:33:20 | 只看该作者
该怎么修改呢?
回复

举报

acer080
地板
 楼主| 发表于 2014-2-22 17:33:56 | 只看该作者
rosynirvana 发表于 2014-2-21 15:43
fatal error C1083: Cannot open include file: 'iostream': No such file or directory

这明显就是编译 ...

怎么修改呢？
回复

举报

rosynirvana
5#
发表于 2014-2-22 22:54:40 | 只看该作者
acer080 发表于 2014-2-22 17:33
怎么修改呢？

你在用vc6？
删掉然后去下个visual c++ 2010 express或者更新的版本
回复

举报

cad20020601
6#
发表于 2014-2-23 08:49:34 | 只看该作者
rosynirvana 发表于 2014-2-22 22:54
你在用vc6？
删掉然后去下个visual c++ 2010 express或者更新的版本

大神？？           
回复

举报

张晟
7#
发表于 2014-2-23 15:13:35 | 只看该作者
差不多吧，我没看出来
回复

举报

1531590867
8#
发表于 2014-2-25 21:11:34 | 只看该作者
呵呵 小白兼即将中考初三党路过~
回复

举报

linwl1234
9#
发表于 2014-2-25 22:27:30 | 只看该作者
唉    天书呀    看不懂！！！！！
回复

举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

广播台
特别关注
阿哈王国（武汉）教育科技有限公司©2018      鄂ICP备17029633号-1
AhaKingdom (Wuhan) Education & Technology Co., Ltd.
快速回复 返回顶部 返回列表