復(fù)合資料力學(xué)上機(jī)編程作業(yè)(計(jì)算層合板剛度)重點(diǎn)復(fù)合資料力學(xué)上機(jī)編程作業(yè)(計(jì)算層合板剛度)重點(diǎn)復(fù)合資料力學(xué)上機(jī)編程作業(yè)(計(jì)算層合板剛度)重點(diǎn)復(fù)合資料力學(xué)上機(jī)編程作業(yè)學(xué)院：SchoolofCivilEngineering專業(yè)：EngineeringMechanics小構(gòu)成員信息：JamesWilson（）、TauYoung（）復(fù)合資料力學(xué)學(xué)了五個(gè)禮拜，這是這門課的第一次編程作業(yè)。我和楊濤結(jié)成一個(gè)小組，我用的是Fortran編制的程序，TauYoung用的是matlab編制。此中的算例以我的Fortran計(jì)算結(jié)果為準(zhǔn)。Matlab作為可視化界面有其獨(dú)到之處，在附錄2中將會(huì)有所展現(xiàn)。作業(yè)的內(nèi)容是層合板的剛度的計(jì)算和驗(yàn)算，包含拉伸剛度A、曲折剛度D以及耦合剛度B。第一要給定層合板的各個(gè)參數(shù)，詳盡有：層合板的層數(shù)N；各單層的彈性常數(shù)E1、E2、υ21、G12；各單層對(duì)應(yīng)的厚度；各單層對(duì)應(yīng)的主方向夾角θ。而后就要計(jì)算每個(gè)單層板的二維剛度矩陣Q，詳盡公式以下：12=υ21E2E1；Q11=E11-υ12υ；21Q22=E21-υ12υ；21Q12=υ12E1；1-υ12υ21Q66=G12獲取Q矩陣后，依據(jù)課本上講到的Q=(T-1)TQ(T-1)獲取Q。而后依據(jù)z坐標(biāo)的定義求出z0到zn，接下來(lái)，最重要的一步，依據(jù)下式計(jì)算、B、D。n??Aij=∑(Qij)k(zk-zk-1)k=1??1n22?Bij=∑(Qij)k(zk-zk-1)2k=1??1n33?Dij=∑(Qij)k(zk-zk-1)3k=1?一、書上P110的幾個(gè)問(wèn)題可以歸納為以下幾個(gè)種類。第1頁(yè)共1頁(yè)（4）6層反對(duì)稱角鋪設(shè)層合板（T5-10）第2頁(yè)共2頁(yè)（5）我還想考據(jù)一個(gè)書上的例題，在課本P114。三層層合板，外層厚度t1，內(nèi)層10t1，正交鋪設(shè)比m=0.2,。二、考據(jù)Verchery的論文里給出的數(shù)值算例。這里向來(lái)到Table5的數(shù)據(jù)都是從Verchery的論文中截獲。Verchery論文中的18層序列，第（21）式【laminateswithoutbending-extensioncoupling】的擺列有兩種資料，一種是Boron-Epoxy，一種是Glass-Epoxy。并且都給出了最后的計(jì)算結(jié)果Q，A*，D*。下邊是我的Fortran計(jì)算數(shù)據(jù)文檔和結(jié)果文檔。（1）Boron-Epoxy資料。第3頁(yè)共3頁(yè)（2）Glass-Epoxy資料。第4頁(yè)共4頁(yè)（3）自然我也考據(jù)了第（22）【laminateswithequalelasticpropertiesinbendingandextension】、（23）【quasi-homogeneouslaminates】的排序，資料是Boron-Epoxy，下邊給出計(jì)算的結(jié)果。從下邊的兩個(gè)結(jié)果表中可以知道，（22）排序的確是C=0，（23）的排序的確是B=0且C=0?？紦?jù)成功。第5頁(yè)共5頁(yè)附件1：計(jì)算所用的程序代碼。PROGRAMCOMPOSITE!CodedbyJamesWilsonIMPLICITNONEREAL(8)::A(3,3),B(3,3),D(3,3),MC(5),TEMP,ROT(3,3)!A拉伸剛度；B耦合剛度；D曲折剛度；!MC讀入資料常數(shù)；ROT旋轉(zhuǎn)矩陣REAL(8)::TOTAL_TH,HALF_TH!總厚度；半厚度REAL(8),ALLOCATABLE::Q(:,:,:),AL(:),T(:),Z(:),Z1(:),Z2(:),Z3(:)!Q每層板相應(yīng)剛度；AL轉(zhuǎn)角；T每層板的厚度；Z坐標(biāo)量INTEGER(4)::N,I,J,K,SEQ,L!____IJK循環(huán)變量；N板的層數(shù)；SEQ序數(shù)CHARACTER(50)::CHR(8),TEMPC,filename1,filename2!CHR、TEMPC:charactervariablesWRITE(*,*)"PleaseinserttheINPfilename(a.txtforexample):"READ(*,*)filename1OPEN(8,file=filename1)!Opendatafile!ReadindataREAD(8,*)TEMPC,NALLOCATE(Q(3,3,N),AL(N),T(N),Z(N+1),Z1(N),Z2(N),Z3(N))READ(8,*)CHR(1:8)DOI=1,NREAD(8,*)SEQ,T(I),MC(1:5),AL(I)ROT(1,1)=(cos(AL(I)))2!WorkoutRotMatrixROT(2,2)=ROT(1,1)ROT(3,3)=cos(2*AL(I))ROT(2,1)=1-ROT(1,1)ROT(1,2)=ROT(2,1)ROT(3,1)=0.5*sin(2*AL(I))ROT(3,2)=-ROT(3,1)ROT(1,3)=-2*ROT(3,1)ROT(2,3)=-2*ROT(3,2)Q(:,:,I)=MATMUL(MATMUL(ROT,Q(:,:,I)),TRANSPOSE(ROT))ENDDOTOTAL_TH=sum(T)HALF_TH=TOTAL_TH/2Z(1)=-HALF_TH!WorkoutZDOI=1,NZ(I+1)=Z(I)+T(I)ENDDO!calculatetensorA、BandDDOK=1,NZ1(K)=(Z(K+1)-Z(K))Z2(K)=(Z(K+1)-Z(K))*(Z(K+1)+Z(K))/2Z3(K)=(Z(K+1)3-Z(K)3)/3ENDDOA=0;B=0;D=0WRITE(*,*)"PleaseinserttheOUPfilename(b.txtforexample):"READ(*,*)filename2OPEN(9,file=filename2)!WriteinstiffnesstensorforeachsingleplyDOK=1,NWRITE(9,100)K100FORMAT("Thestiffnessofnumber",1X,I2,2X,"plyis:")DOI=1,3WRITE(9,200)Q(I,:,K)200FORMAT(ES12.4,6X,ES12.4,6X,ES12.4)ENDDOWRITE(9,"(/)")A=A+Q(:,:,K)*Z1(K)B=B+Q(:,:,K)*Z2(K)D=D+Q(:,:,K)*Z3(K)ENDDO!OutputtheactualstiffnessofthelaminateWRITE(9,"(/)");WRITE(9,"(/)")WRITE(9,*)"TheACTUALstiffnesstensorofthelaminate:"WRITE(9,"(/)")WRITE(9,*)"TheextensionstiffnessAequals:"DOI=1,3WRITE(9,200)A(I,1:3)ENDDOWRITE(9,"(/)")WRITE(9,*)"ThecouplingstiffnessBequals:"DOI=1,3WRITE(9,200)B(I,1:3)ENDDOWRITE(9,"(/)")WRITE(9,*)"ThebendingstiffnessDequals:"DOI=1,3WRITE(9,200)D(I,1:3)ENDDO!NormalisedtensoroutputWRITE(9,"(/)");WRITE(9,"(/)")WRITE(9,*)"TheNORMALISEDstiffnesstensorofthelaminate:"WRITE(9,"(/)")WRITE(9,*)"TheNORMALISEDextensionstiffnessA*equals:"DOI=1,3WRITE(9,200)A(I,1:3)/TOTAL_THENDDOWRITE(9,"(/)")WRITE(9,*)"TheNORMALISEDcouplingtensorCequals:"DOI=1,3WRITE(9,200)A(I,1:3)/TOTAL_TH-12*D(I,1:3)/TOTAL_TH3ENDDOWRITE(9,"(/)")WRITE(9,*)"TheNORMALISEDbendingstiffnessD*equals:"DOI=1,3WRITE(9,200)12*D(I,1:3)/TOTAL_TH3ENDDOWRITE(*,*)"OUTPUTsuccessfully,pleasepressanykeytoendprogram!"READ(*,*)ENDPROGRAMCOMPOSITE附2楊濤同學(xué)的MATLAB（GUI）計(jì)算程序。主要程序：（編了個(gè)小界面，程序略長(zhǎng)，刪掉一些程序自帶解說(shuō)語(yǔ)句，增加了一些對(duì)重點(diǎn)語(yǔ)句的解說(shuō)。）界面是：第6頁(yè)共6頁(yè)作的一個(gè)算比方下：該算例結(jié)果與組內(nèi)伙伴JamesWilson同學(xué)基本一致，其他算例結(jié)果也基本一致，不過(guò)在趨近于零時(shí)有稍微差異，在此不贅于。后邊附上源代碼：functionvarargout=composit_plate(varargin)gui_Singleton=1;gui_State=struct('gui_Name',mfilename,...'gui_Singleton',gui_Singleton,...'gui_OpeningFcn',@composit_plate_OpeningFcn,...'gui_OutputFcn',@composit_plate_OutputFcn,...'gui_LayoutFcn',[],...'gui_Callback',[]);ifnargin&&ischar(varargin{1})gui_State.gui_Callback=str2func(varargin{1});endifnargout[varargout{1:nargout}]=gui_mainfcn(gui_State,varargin{:});elsegui_mainfcn(gui_State,varargin{:});endfunctioncomposit_plate_OpeningFcn(hObject,eventdata,handles,varargin)handles.output=hObject;guidata(hObject,handles);ha=axes('units','normalized','position',[0011]);%嵌入座標(biāo)，為嵌入背景圖片準(zhǔn)備uistack(ha,'down')%作為背景II=imread('武漢大學(xué).jpg');%讀入圖片信息image(II)第7頁(yè)共7頁(yè)colormaphsvset(ha,'handlevisibility','off','visible','off')functionvarargout=composit_plate_OutputFcn(hObject,eventdata,handles)varargout{1}=handles.output;functionedit1_Callback(hObject,eventdata,handles)functionedit1_CreateFcn(hObject,eventdata,handles)ifispc&&isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunctionedit2_Callback(hObject,eventdata,handles)functionedit2_CreateFcn(hObject,eventdata,handles)ifispc&&isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunctionedit3_Callback(hObject,eventdata,handles)functionedit3_CreateFcn(hObject,eventdata,handles)ifispc&&isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunctionedit4_Callback(hObject,eventdata,~)functionedit4_CreateFcn(hObject,eventdata,handles)ifispc&&isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunctionpushbutton1_Callback(hObject,eventdata,handles)symse1e2v21g12ang%本程序采納符號(hào)運(yùn)算v12=v21*e2/e1;q=[e1/(1-v12*v21),v21*e2/(1-v12*v21),0v21*e2/(1-v12*v21),e2/(1-v12*v21),00,0,g12];tran=[cos(ang)^2,sin(ang)^2,-sin(2*ang)sin(ang)^2,cos(ang)^2,sin(2*ang)sin(2*ang)/2,-sin(2*ang)/2,cos(2*ang)];q1=tran*q*tran';%獲取Qn=str2num(get(handles.edit2,'string'));%讀入層數(shù)n=floor(n);nn=0;A=0;B=0;D=0;t=str2num(get(handles.edit3,'string'));%讀入每層厚度t1=zeros(1,n+1);t1(1)=0;fornn=1:nt1(nn+1)=t1(nn)+t(nn);endang1=str2num(get(handles.edit4,'string'));%讀入每層角度const=str2num(get(handles.edit1,'string'));%讀入資料系數(shù)t0=sum(t)/2;t1=t1-t0;e1=const(1);e2=const(2);v21=const(3);g12=const(4);q11=subs(q1);fornn=1:nang=ang1(nn);nn=nn+1;A=A+subs(q11*(t1(nn)-t1(nn-1)));B=B+subs(0.5*q11*(t1(nn)^2-t1(nn-1)^2));D=D+subs(1/3*q11*(t1(nn)^3-t1(nn-1)^3));end%累加計(jì)算set(handles.edit5,'string',num2str(A(1,:)));%以下為輸出結(jié)果set(handles.edit6,'string',num2str(A(2,:)));set(handles.edit7,'string',num2str(A(3,:)));set(handles.edit8,'string',num2str(B(1,:)));set(handles.edit9,'string',num2str(B(2,:)));set(handles.edit10,'string',num2str(B(3,:)));set(handles.edit11,'string',num2str(D(1,:)));set(handles.edit12,'string',num2str(D(2,:)));第8頁(yè)共8頁(yè)set(handles.edit13,'string',num2str(D(3,:)));functionedit5_Callback(hObject,eventdata,handles)functionedit5_CreateFcn(hObject,eventdata,handles)ifispc&&isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunctionedit6_Callback(hObject,eventdata,handles)functionedit6_CreateFcn(hObject,eventdata,handles)ifispc&&isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunctionedit7_Callback(hObject,eventdata,handles)functionedit7_CreateFcn(hObject,eventdata,handles)ifispc&&isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunctionedit8_Callback(hObject,eventdata,handles)functionedit8_CreateFcn(hObject,eventdata,handles)ifispc&&isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunctionedit9_Callback(hObject,eventdata,handles)functionedit9_CreateFcn(hObject,eventdata,handles)ifispc&&isequal(get(hObject,'BackgroundColor'),get(0,'defaultUicontrolBackgroundColor'))set(hObject,'BackgroundColor','white');endfunctionedit10_Callback(hObject,eventdata,handles)f