应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » C/C++ » C++参数化建模程序出错，谁能帮忙看下错在哪里？
北京石油化工学院2026年研究生招生接收调剂公告
81/1返回列表
查看: 1625  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 ]
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-02-21 12:58:52
已阅   回复此楼   关注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 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +5 一样YWY 2026-04-06 5/250 2026-04-06 13:08 by 无际的草原
[考研] 308求调剂 +12 倘若起风了呢 2026-04-05 12/600 2026-04-06 10:42 by lqwchd
[考研] 求调剂 +19 111623 2026-04-04 21/1050 2026-04-06 09:35 by guanxin1001
[考研] 一志愿南昌大学，085600，344分求调剂 +9 调剂上岸玘 2026-04-05 10/500 2026-04-06 09:30 by dongzh2009
[考研] 考研调剂生寻找导师 +3 顾瞻考研啊 2026-04-05 3/150 2026-04-05 18:18 by 啵啵啵0119
[考研] 270求调剂 +9 小杰pp 2026-03-31 11/550 2026-04-05 11:02 by 风雨无晴
[考研] 283分求调剂 +9 试试看呗 2026-04-04 9/450 2026-04-05 10:27 by 果冻大王
[考研] 考研求调剂 +3 木心想继续深造 2026-04-03 3/150 2026-04-03 21:56 by 啵啵啵0119
[考研] 320调剂 +4 农业工程与信息� 2026-04-03 4/200 2026-04-03 21:40 by lbsjt
[考研] 材料专硕322分 +13 哈哈哈吼吼吼哈 2026-04-01 13/650 2026-04-03 16:08 by 哦哦123
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +9 哇呼哼呼哼 2026-04-02 9/450 2026-04-03 12:05 by 1753564080
[考研] 重庆大学材料与化工085600，初试370+，求求调剂建议 +8 shzhou_ 2026-04-01 9/450 2026-04-03 09:31 by 蓝云思雨
[考研] 295求调剂 +7 愿旅途永远坦然 2026-04-02 7/350 2026-04-03 08:22 by fangshan711
[考研] 26考研调剂 +4 Wnz.20030617 2026-04-01 5/250 2026-04-02 16:11 by 1939136013狗壮
[考研] 318求调剂 +3 笃行致远. 2026-03-31 4/200 2026-04-02 15:56 by Jaylen.
[考研] 0856初试324分求调剂 +6 想上学求调 2026-04-01 6/300 2026-04-02 11:42 by 星空星月
[考研] 一志愿厦门大学化学工程（专硕）-数二英二406分-求调剂 +5 厦大化工 2026-04-01 5/250 2026-04-02 10:03 by jp9609
[考研] 284求调剂 +12 小熊～～ 2026-03-31 12/600 2026-04-01 20:23 by 花??
[考研] 材料调剂 +10 Eujd1 2026-03-31 11/550 2026-04-01 11:23 by ivanqyq
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
信息提示
关闭
请填处理意见
关闭