Effective Assessment of Inelastic Torsional Deformation of Plan-Asymmetric Shear Wall Systems

Torsional deformation may occur in plan-asymmetric wall structures during seismic events. In most current seismic design codes, the torsional deformation is handled using the concept of design eccentricity, in which design loads may be excessively amplified. This approach has limitations in accurately estimating the torsional deformation of the plan-asymmetric structures, mainly because it is based on linear elastic material behavior. In this paper, we propose a simple method that can accurately evaluate the inelastic lateral displacement and rotation of plan-asymmetric wall structures and is suitable for use in the displacement-based design method. The effectiveness of the proposed method is verified by comparing its predictions with the results of rigorous time history analyses for a model problem. The comparison shows that the proposed method is able to provide accurate estimates of the inelastic torsional deformation for the plan-asymmetric wall system, while requiring less computational cost than the time history analysis.


Introduction
Plan-asymmetric structures may undergo torsional vibration in addition to lateral oscillation if they are subjected to seismic ground motions. This induces an additional internal force, which may lead to excessive deformation in relatively flexible structural components. This excessive deformation often results in the unexpected brittle failure of structures, which has been reported during several earthquake events in the past such as the Mexico earthquake in 1985 [1], the Chile earthquake in 1985 [2], the Loma Prieta earthquake in 1989 [3], and the Kobe earthquake in 1995 [4]. The unsymmetrical distribution of mass and lateral load resisting stiffness of the plan-asymmetric structure is the primary source of the lateral torsional structural response. However, torsional deformation may occasionally be induced even in symmetric structures due to non-uniform ground motions along the foundation of the structure or the inherent torsional components of ground excitation. Details of the structural as well as non-structural members also may be sources of the torsional response.
For this reason, most current seismic design codes include criteria to minimize the torsional response of building structures as well as damage caused by it. The current seismic design codes such as the International Building Code (IBC) [5], New Zealand Standards (NZS) (Standards New Zealand, 2004) [6], the National Building Code of Canada (NBCC) [7] and Eurocode [8] suggest design eccentricity based on the linear elastic state of structural components. However, the inelastic torsional response of a structure is quite different from that in the elastic range. If the main lateral force resisting components begins to yield, the eccentricity, lateral force distribution, and amount of rotation may

Current Code Provisions on Torsional Design of Plan-Asymmetric Structures
In general, most of the current seismic design code provisions are based on the force-based design approach. In this approach, the static eccentricity (e) is defined as the distance between the centers of mass and rigidity and is used to compute the torsional moment acting on a plan-asymmetric structure. Based on this definition, the torsional moment (T) can be computed by where V is the lateral force applied to the structure. If more than one lateral force components are applied, the resulting torsional moment can be calculated by where e is the eccentricity of each lateral force components. From this relation, the effective eccentricity can be derived as follows: The current seismic design codes such as IBC, Eurocode, NZS and NBCC introduce the concept of the design eccentricity (e design ), which can generally be written in the form of where e is the eccentricity defined as the distance between the center of mass (CM) and the center of rigidity (CR), b is the floor plan length perpendicular to the excitation direction, and γ and η are the coefficients specified in the seismic design codes listed in Table 1. In Equation (4), the first term on the right-hand side accounts for the natural torsional effect, and the second term incorporates the accidental torsional effect. The design eccentricity proposed by IBC can be stated as e design = 1.0e ± 0.05b.
In this case, no dynamic amplification is considered, and the eccentricity by accidental torsion is defined as ±5% of the plan length in the direction perpendicular to the applied lateral force. If the diaphragm is not rigid, and the maximum deformation at the edge of the plan considered exceeds 1.2 times the average deformation, the entire design eccentricity expressed by Equation (5) is multiplied by a torsional amplification factor (A x ) given as Appl. Sci. 2019, 9, 2814 4 of 24 In Equation (6), δ max is the maximum deformation at the plan edge, and δ avg is the average value of the deformations measured at both of the plan edges.
NZS suggests a design eccentricity similar to that of IBC, which is expressed by e design = 1.0e ± 0.1b.
It can be noted from Equation (7) that no dynamic effect is considered in NZS and ±10% of the plan length in the direction perpendicular to the applied lateral force is used as the eccentricity by accidental torsion. The design eccentricity suggested in NBCC is similar to that of NZS, but the eccentricity amplification due to dynamic effects is also considered. This suggests that the design force should be computed with the maximum eccentricity value among the four possible combinations calculated by e design = (1.0 ± 0.5)e ± 0.1b. (8) In Eurocode 8, if the lateral stiffness and mass are symmetrically distributed on the plan and the accidental torsion is not considered by a more accurate method, the accidental torsional effect can be accounted for by multiplying the load applied to an individual lateral load resisting element with δ, which is given by where x is the distance between the individual element considered and the center of mass in the perpendicular direction to seismic excitation, and L e is the distance between the two-outermost lateral load resisting elements measured in the same direction as x.
An alternative approximate method to the one above is presented in the annex of Eurocode 8. To consider the dynamic effect of simultaneous lateral and torsional vibrations, two design eccentricity equations, i.e., (10) and (11), are suggested. e max = e 0 + e 1 + e 2 , (10) e min = e 0 − e 1 .
Here, e 0 is the static eccentricity defined by Equation (3), e 1 the accidental eccentricity corresponding to 5% of the plan length in the direction perpendicular to seismic excitation, and e 2 is the additional eccentricity reflecting the dynamic effect induced by simultaneous lateral and torsional vibrations. The dynamic eccentricity e 2 is determined as the smaller of those computed in Equations (12) and (13).
where L and B represent the plan lengths in the directions perpendicular and parallel to seismic excitation, respectively, r is the ratio of the story torsional and lateral stiffness, and l s indicates the radius of gyration approximated by The inelastic deformation and actual failure mechanism of plan-asymmetric structures are not considered in any of the seismic design codes discussed up to this point. As a consequence, they may require excessive strength and stiffness for lateral force resisting components. A three-dimensional modal or dynamic time-history analysis can be performed for a more accurate evaluation. However, this is highly complicated and too time-consuming for practical design purposes. Therefore, it is necessary to develop a simple and accurate method for assessing inelastic torsional deformation of plan-asymmetric wall structures.

Inelastic Torsional Response of Plan-Asymmetric Structures
The torsional response of plan-asymmetric structures in the inelastic range becomes quite different from that in the linear elastic range. In the linear elastic range, the rotation of the asymmetric structure increases in proportion to the level of the lateral load applied, but if any of the lateral force resisting walls start to yield, there is change in the overall torsional stiffness of the structure, which leads to a change in the eccentricity value as well. This again affects the torsional moment induced in the asymmetric structure, and thus the rotation of the structure increases in a nonlinear fashion in the inelastic range.
To accurately assess the inelastic torsional response of the structure, we derived the generalized form of the inelastic lateral displacement-rotation relation of the structure by following a step by step procedure, and we represent this relation in the lateral displacement-rotation coordinate system in this section. According to the authors' knowledge, this coordinate system has never been introduced before and is very convenient to describe the inelastic torsional behavior of the asymmetric structure. In addition, we discuss the effects of transverse wall stiffness on the yielding procedure of asymmetric structures.

Torsional Response of Asymmetric Structures in the D-R Coordinate System
In this section, the generalized form of the displacement-rotation relation of plan-asymmetric structures is derived by assuming that the structure is in the linear elastic range. From this relation, we can derive the inelastic displacement-rotation relation of plan-asymmetric structures. Figure 1 illustrates a representative asymmetric wall structure subjected to a lateral force in the y-direction (V) as well as its deformed shape. In the figure, Wall X i and Wall Y i are the i-th lateral load resisting walls in the xand y-directions, respectively. The x-coordinate of Wall Y i and y-coordinate of Wall X i with respect to the center of rigidity (CR) in the undeformed configuration are denoted as r yi and r xi , respectively. Similarly, the same quantities with respect to the center of mass (CM) are expressed as d yi and d xi , respectively.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 5 of 25 require excessive strength and stiffness for lateral force resisting components. A three-dimensional modal or dynamic time-history analysis can be performed for a more accurate evaluation. However, this is highly complicated and too time-consuming for practical design purposes. Therefore, it is necessary to develop a simple and accurate method for assessing inelastic torsional deformation of plan-asymmetric wall structures.

Inelastic Torsional Response of Plan-Asymmetric Structures
The torsional response of plan-asymmetric structures in the inelastic range becomes quite different from that in the linear elastic range. In the linear elastic range, the rotation of the asymmetric structure increases in proportion to the level of the lateral load applied, but if any of the lateral force resisting walls start to yield, there is change in the overall torsional stiffness of the structure, which leads to a change in the eccentricity value as well. This again affects the torsional moment induced in the asymmetric structure, and thus the rotation of the structure increases in a nonlinear fashion in the inelastic range.
To accurately assess the inelastic torsional response of the structure, we derived the generalized form of the inelastic lateral displacement-rotation relation of the structure by following a step by step procedure, and we represent this relation in the lateral displacement-rotation coordinate system in this section. According to the authors' knowledge, this coordinate system has never been introduced before and is very convenient to describe the inelastic torsional behavior of the asymmetric structure. In addition, we discuss the effects of transverse wall stiffness on the yielding procedure of asymmetric structures.

Torsional Response of Asymmetric Structures in the D-R Coordinate System
In this section, the generalized form of the displacement-rotation relation of plan-asymmetric structures is derived by assuming that the structure is in the linear elastic range. From this relation, we can derive the inelastic displacement-rotation relation of plan-asymmetric structures. Figure 1 illustrates a representative asymmetric wall structure subjected to a lateral force in the y-direction (V) as well as its deformed shape. In the figure, Wall Xi and Wall Yi are the i-th lateral load resisting walls in the x-and y-directions, respectively. The x-coordinate of Wall Yi and y-coordinate of Wall Xi with respect to the center of rigidity (CR) in the undeformed configuration are denoted as ryi and rxi, respectively. Similarly, the same quantities with respect to the center of mass (CM) are expressed as dyi and dxi, respectively. The lateral force (V) and torsional moment (T) of the system can be calculated by Equations (15) and (16), respectively: The lateral force (V) and torsional moment (T) of the system can be calculated by Equations (15) and (16), respectively: Appl. Sci. 2019, 9, 2814 6 of 24 Here, K y and K T are the lateral stiffness in the y-direction and torsional stiffness of the asymmetric system, respectively. In addition, ∆ and ϕ are the lateral displacement and rotation, respectively, and they are estimated at the center of mass. The counterclockwise direction of ϕ is assumed to be positive. In the derivation of Equation (16), we assumed that the diaphragm has an infinite rigidity, and the out-of-plane deformation of the shear wall is ignored. If no external torsional moment is applied, the resulting torsional moment of the asymmetric system can be determined by multiplying the lateral force and the eccentricity in the x-direction between the centers of mass and rigidity (e x ), as given below: The eccentricity in the x-direction can be computed by where k yi is the lateral stiffness of the i-th wall in the y-direction (Wall Y i ). By substituting Equations (15) and (16) into Equation (17), the rotation can be expressed in terms of the lateral displacement as The lateral and torsional stiffness of the asymmetric structure can be estimated by where k xi represents the lateral stiffness of the i-th wall in the x-direction. By inserting Equations (20) and (21) into Equation (19), the relationship between the lateral displacement and rotation can be obtained as follows: As shown in Figure 2, this relationship can be represented by a solid line with an arrow on the lateral displacement (∆) and rotation (ϕ) coordinate system, which is denoted as the D-R coordinate system in this paper. This indicates that a linear relation holds between the lateral displacement and the rotation as long as no wall component yields. The yield condition for Wall Y 1 can be stated as where u y is the yield displacement of Wall Y 1 . The figure shows that the yielding of Wall Y 1 occurs at the intersection of Equations (22) and (23). As will be discussed later, this coordinate system is very convenient to describe the inelastic torsional behavior of plan-asymmetric structures. Appl. Sci. 2019, 9, x FOR PEER REVIEW 7 of 25

Inelastic Torsional Response Assessment Procedure
In this section, we discuss a procedure to assess the inelastic torsional response of plan-asymmetric wall structures. If any wall components yield during the application of lateral force, the linear relationship between the lateral displacement and rotation is no longer valid, and thus it must be modified by considering the reduced stiffness of the wall component that has yielded. Figure 3 illustrates a simple model problem of plan-asymmetric wall structures to explain the assessment procedure of the inelastic torsional response. This indicates that the structure has two lateral force resisting walls in the x-and y-directions, respectively, and the lateral force is applied in the y-direction. The two walls in the y-direction (Walls Y1 and Y2) have elastic lateral stiffnesses of k and αk (α > 1), respectively, and thus Wall Y2 is stiffer than Wall Y1. The two walls in the x-direction have the same elastic lateral stiffness βk. Here, β can be of any value, but it may affect the inelastic behavior of the asymmetric structure after yielding of Walls Y1 and Y2, as will be discussed later. For simplicity, the out-of-plane deformation of the walls is ignored. The post yielding behavior of the two lateral force resisting walls Y1 and Y2 is defined as a bilinear function, as shown in Figure 4. They have the same yield displacement ( y u ), and their post-yield stiffness is γ times the elastic stiffness ( E y K ). In contrast, we assumed that the other two walls in the transverse direction (Walls X1 and X2) remain in the linear elastic range even after the two lateral force resisting walls yield. The inelastic lateral-torsional response of this structure can be obtained by following the step-by-step procedure described below.

Inelastic Torsional Response Assessment Procedure
In this section, we discuss a procedure to assess the inelastic torsional response of plan-asymmetric wall structures. If any wall components yield during the application of lateral force, the linear relationship between the lateral displacement and rotation is no longer valid, and thus it must be modified by considering the reduced stiffness of the wall component that has yielded. Figure 3 illustrates a simple model problem of plan-asymmetric wall structures to explain the assessment procedure of the inelastic torsional response. This indicates that the structure has two lateral force resisting walls in the xand y-directions, respectively, and the lateral force is applied in the y-direction. The two walls in the y-direction (Walls Y 1 and Y 2 ) have elastic lateral stiffnesses of k and αk (α > 1), respectively, and thus Wall Y 2 is stiffer than Wall Y 1 . The two walls in the x-direction have the same elastic lateral stiffness βk. Here, β can be of any value, but it may affect the inelastic behavior of the asymmetric structure after yielding of Walls Y 1 and Y 2 , as will be discussed later. For simplicity, the out-of-plane deformation of the walls is ignored.

Inelastic Torsional Response Assessment Procedure
In this section, we discuss a procedure to assess the inelastic torsional response of plan-asymmetric wall structures. If any wall components yield during the application of lateral force, the linear relationship between the lateral displacement and rotation is no longer valid, and thus it must be modified by considering the reduced stiffness of the wall component that has yielded. Figure 3 illustrates a simple model problem of plan-asymmetric wall structures to explain the assessment procedure of the inelastic torsional response. This indicates that the structure has two lateral force resisting walls in the x-and y-directions, respectively, and the lateral force is applied in the y-direction. The two walls in the y-direction (Walls Y1 and Y2) have elastic lateral stiffnesses of k and αk (α > 1), respectively, and thus Wall Y2 is stiffer than Wall Y1. The two walls in the x-direction have the same elastic lateral stiffness βk. Here, β can be of any value, but it may affect the inelastic behavior of the asymmetric structure after yielding of Walls Y1 and Y2, as will be discussed later. For simplicity, the out-of-plane deformation of the walls is ignored. The post yielding behavior of the two lateral force resisting walls Y1 and Y2 is defined as a bilinear function, as shown in Figure 4. They have the same yield displacement ( y u ), and their post-yield stiffness is γ times the elastic stiffness ( E y K ). In contrast, we assumed that the other two walls in the transverse direction (Walls X1 and X2) remain in the linear elastic range even after the two lateral force resisting walls yield. The inelastic lateral-torsional response of this structure can be obtained by following the step-by-step procedure described below.
Yielding point Figure 3. Model problem of a plan-asymmetric structure.
The post yielding behavior of the two lateral force resisting walls Y 1 and Y 2 is defined as a bilinear function, as shown in Figure 4. They have the same yield displacement (u y ), and their post-yield stiffness is γ times the elastic stiffness (K E y ). In contrast, we assumed that the other two walls in the transverse direction (Walls X 1 and X 2 ) remain in the linear elastic range even after the two lateral force resisting walls yield. The inelastic lateral-torsional response of this structure can be obtained by following the step-by-step procedure described below. Appl. Sci. 2019, 9, x FOR PEER REVIEW 8 of 25 (1) Linear elastic stage In the linear elastic stage, the lateral displacement and rotation of the asymmetric structure can be computed by Equation (19), where the lateral and torsional stiffnesses and the static eccentricity are estimated by Equations (24)-(26): In the above equations, the superscript 0 indicates the linear elastic stage. Finally, Equation (19) can be rewritten as where 0 Δ and 0 ϕ denote the lateral displacement and rotation of the asymmetric structure in the linear elastic stage, respectively. The deformation history in the displacement-rotation coordinate system and location of the center of rigidity in the linear elastic range are provided in Figure 5. As shown in the figure, the displacement-rotation relationship in the linear elastic stage defined by Equation (27) is represented by a line with an arrow.  (1) Linear elastic stage In the linear elastic stage, the lateral displacement and rotation of the asymmetric structure can be computed by Equation (19), where the lateral and torsional stiffnesses and the static eccentricity are estimated by Equations (24)-(26): In the above equations, the superscript 0 indicates the linear elastic stage. Finally, Equation (19) can be rewritten as where ∆ 0 and ϕ 0 denote the lateral displacement and rotation of the asymmetric structure in the linear elastic stage, respectively. The deformation history in the displacement-rotation coordinate system and location of the center of rigidity in the linear elastic range are provided in Figure 5. As shown in the figure, the displacement-rotation relationship in the linear elastic stage defined by Equation (27)

(1) Linear elastic stage
In the linear elastic stage, the lateral displacement and rotation of the asymmetric structure can be computed by Equation (19), where the lateral and torsional stiffnesses and the static eccentricity are estimated by Equations (24)-(26): In the above equations, the superscript 0 indicates the linear elastic stage. Finally, Equation (19) can be rewritten as where 0 Δ and 0 ϕ denote the lateral displacement and rotation of the asymmetric structure in the linear elastic stage, respectively. The deformation history in the displacement-rotation coordinate system and location of the center of rigidity in the linear elastic range are provided in Figure 5. As shown in the figure, the displacement-rotation relationship in the linear elastic stage defined by Equation (27) is represented by a line with an arrow.  (2) First yielding point As the lateral force increases, the wall on the flexible side yields first since (i) it has greater displacement than the one on the stiff side and (ii) the two walls have the same yield displacement. The wall on the flexible side starts to yield if its y-displacement reaches the yield displacement, and this yield limit can be stated by Consequently, the first yielding point can be defined as the intersection of the lines represented by Equations (27) and (28), as shown in Figure 5. This first yielding point (YP 1 ) can be obtained by inserting Equation (28) into Equation (27), and its coordinates are given by (

3) First inelastic stage
After the yielding of the wall component on the flexible side, its lateral stiffness is reduced to γ. From this, the lateral and torsional stiffnesses and the eccentricity of the asymmetric structure can be re-estimated as follows.
Here, superscript 1 indicates the first inelastic stage. Similarly, the relationship between the increments of the lateral displacement (δ∆ 1 ) and rotation (δϕ 1 ) after the first yielding can be expressed by

(4) Second yielding point
A second yielding occurs if the lateral displacement of the wall on the stiff side also reaches the yield displacement with increasing lateral load. This limit state can be expressed by Like the first yielding, the second yielding point can be obtained by finding the intersection of Equations (33) and (34). This procedure is also illustrated in Figure 6. Thus, δ∆ 1 and δϕ 1 can be expressed by Appl. Sci. 2019, 9, x FOR PEER REVIEW 10 of 25 Finally, the coordinates of the second yielding point (YP2) can be written as (37) (5) Second inelastic stage (final stage) As in the case of the wall component on the flexible side, the lateral stiffness of the wall on the stiff side is also reduced to γk. After yielding of the two wall components in the y-direction, the modified lateral and torsional stiffnesses and the eccentricity of the asymmetric system can be respectively expressed by.
where the superscript 2 indicates the second inelastic stage. Note that the eccentricity in these equations is the same as the one in the linear elastic stage. This is mainly because the lateral stiffnesses of the two walls in the y-direction are assumed to be reduced by the same rate after yielding and the transverse wall components in the x-direction do not contribute to the x-directional stiffness eccentricity. The transverse walls in the x-direction still remain in the elastic stage and now offer relatively large torsional stiffness compared to the wall components in the y-direction. Finally, the relationship between the increase in lateral displacement ( 2 δΔ ) and rotation ( 2 δϕ ) at this stage can be given by and this is plotted in Figure 7. Figure 8 schematically shows the typical lateral-torsional responses of asymmetric structures obtained using three different approaches on the D-R coordinate system based on the discussion up to this point. In the figure, the dashed line represents the linear elastic torsional response given by Equation (22), and the dotted line is the amplified linear elastic torsional response, which is adopted in most of the current design codes in terms of the amplified design eccentricity, as already discussed in Section 2. The solid line represents the inelastic torsional response based on the procedure presented in this section. A comparison of these three curves shows that the amplified Yielding procedure on D-R coordinate system (b) Locations of centers of mass and rigidity Finally, the coordinates of the second yielding point (YP 2 ) can be written as (5) Second inelastic stage (final stage) As in the case of the wall component on the flexible side, the lateral stiffness of the wall on the stiff side is also reduced to γk. After yielding of the two wall components in the y-direction, the modified lateral and torsional stiffnesses and the eccentricity of the asymmetric system can be respectively expressed by.
where the superscript 2 indicates the second inelastic stage. Note that the eccentricity in these equations is the same as the one in the linear elastic stage. This is mainly because the lateral stiffnesses of the two walls in the y-direction are assumed to be reduced by the same rate after yielding and the transverse wall components in the x-direction do not contribute to the x-directional stiffness eccentricity. The transverse walls in the x-direction still remain in the elastic stage and now offer relatively large torsional stiffness compared to the wall components in the y-direction. Finally, the relationship between the increase in lateral displacement (δ∆ 2 ) and rotation (δϕ 2 ) at this stage can be given by and this is plotted in Figure 7.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 11 of 25 linear-elastic response, which is adopted in most of the current design codes, may both underestimate the torsional rotation at the initial yielding stage and overestimate it at the post yielding stage, resulting in an inaccurate estimation of the torsional response over the entire inelastic range.  Figure 8 schematically shows the typical lateral-torsional responses of asymmetric structures obtained using three different approaches on the D-R coordinate system based on the discussion up to this point. In the figure, the dashed line represents the linear elastic torsional response given by Equation (22), and the dotted line is the amplified linear elastic torsional response, which is adopted in most of the current design codes in terms of the amplified design eccentricity, as already discussed in Section 2. The solid line represents the inelastic torsional response based on the procedure presented in this section. A comparison of these three curves shows that the amplified linear-elastic response, which is adopted in most of the current design codes, may both underestimate the torsional rotation at the initial yielding stage and overestimate it at the post yielding stage, resulting in an inaccurate estimation of the torsional response over the entire inelastic range.
linear-elastic response, which is adopted in most of the current design codes, may both underestimate the torsional rotation at the initial yielding stage and overestimate it at the post yielding stage, resulting in an inaccurate estimation of the torsional response over the entire inelastic range.

Effect of Transverse Wall Stiffness on Yielding Procedure
In the derivation introduced in Section 3.2, we assumed for simplicity that the two lateral force resisting walls have the same yield displacement. As a result, the wall on the flexible side always yields first, and then the wall on the stiff side. However, wall stiffness and yield displacement are determined mainly by the length of the wall. According to [24,25], the stiffness (k) and yield displacement (uy) of the lateral force resisting wall can be determined by the following two relationships: where lw denotes the wall length. In this case, the stiffness of the transverse wall may affect the overall yielding procedure of asymmetric structures, and this issue is further discussed in this section.

Effect of Transverse Wall Stiffness on Yielding Procedure
In the derivation introduced in Section 3.2, we assumed for simplicity that the two lateral force resisting walls have the same yield displacement. As a result, the wall on the flexible side always yields first, and then the wall on the stiff side. However, wall stiffness and yield displacement are determined mainly by the length of the wall. According to [24,25], the stiffness (k) and yield displacement (u y ) of the lateral force resisting wall can be determined by the following two relationships: where l w denotes the wall length. In this case, the stiffness of the transverse wall may affect the overall yielding procedure of asymmetric structures, and this issue is further discussed in this section. An example problem illustrated in Figure 9 is considered to explore this issue. This is similar to the model problem shown in Figure 3 except that its floor plan has a square shape, not a rectangular shape. The lateral stiffness (k w1 ), wall length (l w1 ) and yield displacement (u y,w1 ) of the flexible side wall Y 1 are denoted by k, l w and u y , respectively. The same quantities for the stiff side wall Y 2 (k w2 , l w2 and u y,w2 , respectively) are given by where α is greater than 1 and is the stiffness ratio of the stiff and flexible side walls. Note that the parameters given by Equations (44)-(46) satisfy the conditions stated by (42) and (43). We assumed that the two transverse walls in the x-direction (Walls X 1 and X 2 ) have the same elastic lateral stiffness as βk and never yield.
where α is greater than 1 and is the stiffness ratio of the stiff and flexible side walls. Note that the parameters given by Equations (44)-(46) satisfy the conditions stated by (42) and (43). We assumed that the two transverse walls in the x-direction (Walls X1 and X2) have the same elastic lateral stiffness as βk and never yield. By using Equations (18)-(21), the relationship between lateral displacement and rotation angle of the example problem can be expressed by From this equation, the displacement of the two lateral force resisting walls can be derived as The flexible side wall yields first if the following condition is met: By substituting Equations (48)  By using Equations (18)-(21), the relationship between lateral displacement and rotation angle of the example problem can be expressed by From this equation, the displacement of the two lateral force resisting walls can be derived as The flexible side wall yields first if the following condition is met: By substituting Equations (48) and (49) into Equation (50), we obtain The condition stated by Equation (51) is plotted in Figure 10. The horizontal and vertical axes of the figure represent the total wall stiffness in the excitation and transverse directions, respectively. The results in the figure indicate that first yielding may occur even in the stiff side wall if the total transverse wall stiffness is located in the region above the curve plotted in the figure. For this reason, some researchers have classified plan-asymmetric structures into torsionally unrestrained (TU) and torsionally restrained (TR) structures to consider this effect [11][12][13]. The inelastic torsional response assessment procedure proposed in this study can properly handle this effect without introducing the classification as in the previous research, but none of the current code provisions like IBC and Eurocode consider this effect. some researchers have classified plan-asymmetric structures into torsionally unrestrained (TU) and torsionally restrained (TR) structures to consider this effect [11][12][13]. The inelastic torsional response assessment procedure proposed in this study can properly handle this effect without introducing the classification as in the previous research, but none of the current code provisions like IBC and Eurocode consider this effect.

Comparison between Predictions by the Proposed Method and Time History Analysis Results
The effectiveness of the inelastic torsional response assessment procedure proposed in Section 3 was verified by comparing its results with those of time history analyses. The proposed method was implemented using the well-known script language MATLAB [26]. The time history analyses were performed by using well-known seismic analysis software, OPENSEES [27], and their solutions are taken as the reference solution.

Model Problem and Analysis Method
The model problem for verification of the proposed method is the same as the one illustrated in Figure 3. In the model problem, the size of the plan is given as a = 5 m, and the lateral stiffness of Wall Y1 on the flexible side is taken as k = 1000 kN/m. We assumed that a mass of 500 kg is lumped at the center of mass and the yield displacement of the two lateral force resisting walls in the y-direction is defined as 0.05 m y u = . In contrast, we also assumed that the two x-directional walls never yield. A Rayleigh damping ratio of 5% was adopted.
For the time history analysis, each structural wall in the asymmetric structure was modeled using the equivalent truss members, as illustrated in Figure 11. In the modeling, we assumed that the horizontal and vertical trusses are rigid, and the diagonal truss members have zero compressive

Comparison between Predictions by the Proposed Method and Time History Analysis Results
The effectiveness of the inelastic torsional response assessment procedure proposed in Section 3 was verified by comparing its results with those of time history analyses. The proposed method was implemented using the well-known script language MATLAB [26]. The time history analyses were performed by using well-known seismic analysis software, OPENSEES [27], and their solutions are taken as the reference solution.

Model Problem and Analysis Method
The model problem for verification of the proposed method is the same as the one illustrated in Figure 3. In the model problem, the size of the plan is given as a = 5 m, and the lateral stiffness of Wall Y 1 on the flexible side is taken as k = 1000 kN/m. We assumed that a mass of 500 kg is lumped at the center of mass and the yield displacement of the two lateral force resisting walls in the y-direction is defined as u y = 0.05 m. In contrast, we also assumed that the two x-directional walls never yield. A Rayleigh damping ratio of 5% was adopted.
For the time history analysis, each structural wall in the asymmetric structure was modeled using the equivalent truss members, as illustrated in Figure 11. In the modeling, we assumed that the horizontal and vertical trusses are rigid, and the diagonal truss members have zero compressive stiffness (tension only element) to effectively define the hysteretic behavior of the wall component. The elastic stiffness (k diag ) and yield force (P y,diag ) of the diagonal truss elements are defined as follows: where k wall is the elastic lateral stiffness of the wall considered, P y,wall is its yield force. In addition, P sys and ∆ sys are the applied lateral force and the lateral displacement at the top of the wall, respectively. Each diagonal truss member consists of two parallel elements, as shown in Figure 12. The primary element among the two elements has an elastic-fully plastic behavior, and the values corresponding to 94% of the total wall yield force and elastic stiffness were set to its yield force and elastic lateral stiffness, respectively. The material behavior of the secondary element is modeled by a fully elastic relationship, of which the elastic stiffness is the 6% of the total wall stiffness. As a result, the structural wall with this diagonal truss member can exhibit the combined bilinear behavior shown in the figure, where the yield force and tangent stiffness after yielding are set to P y,diag and 0.06k diag , respectively.
and Δsys are the applied lateral force and the lateral displacement at the top of the wall, respectively. Each diagonal truss member consists of two parallel elements, as shown in Figure 12. The primary element among the two elements has an elastic-fully plastic behavior, and the values corresponding to 94% of the total wall yield force and elastic stiffness were set to its yield force and elastic lateral stiffness, respectively. The material behavior of the secondary element is modeled by a fully elastic relationship, of which the elastic stiffness is the 6% of the total wall stiffness. As a result, the structural wall with this diagonal truss member can exhibit the combined bilinear behavior shown in the figure, where the yield force and tangent stiffness after yielding are set to Py,diag and 0.06kdiag, respectively.  The modified Takeda model was used for the nonlinear hysteretic relationship between the time history analysis because its effectiveness has been verified by numerous experimental results. [28,29] Figure 13 shows the hysteretic response of the model problem when subjected to static cyclic  and Δsys are the applied lateral force and the lateral displacement at the top of the wall, respectively. Each diagonal truss member consists of two parallel elements, as shown in Figure 12. The primary element among the two elements has an elastic-fully plastic behavior, and the values corresponding to 94% of the total wall yield force and elastic stiffness were set to its yield force and elastic lateral stiffness, respectively. The material behavior of the secondary element is modeled by a fully elastic relationship, of which the elastic stiffness is the 6% of the total wall stiffness. As a result, the structural wall with this diagonal truss member can exhibit the combined bilinear behavior shown in the figure, where the yield force and tangent stiffness after yielding are set to Py,diag and 0.06kdiag, respectively.  The modified Takeda model was used for the nonlinear hysteretic relationship between the time history analysis because its effectiveness has been verified by numerous experimental results. [28,29] Figure 13 shows the hysteretic response of the model problem when subjected to static cyclic  The modified Takeda model was used for the nonlinear hysteretic relationship between the time history analysis because its effectiveness has been verified by numerous experimental results. [28,29] Figure 13 shows the hysteretic response of the model problem when subjected to static cyclic loading. Figure 13a plots the lateral force-displacement relationships for the flexible side wall Y 1 and stiff side wall Y 2 , which are distinguished by dotted and solid lines, respectively. Similarly, the relationship between the lateral displacement and rotation angle of the model problem estimated at its center of mass is plotted in Figure 13b. loading. Figure 13a plots the lateral force-displacement relationships for the flexible side wall Y1 and stiff side wall Y2, which are distinguished by dotted and solid lines, respectively. Similarly, the relationship between the lateral displacement and rotation angle of the model problem estimated at its center of mass is plotted in Figure 13b.

Analysis Parameters
This section discusses various parameters that are considered in the time history analysis. As given in Table 2, the time history analyses were performed by changing three different types of parameters, which are the relative stiffness of the stiff side wall (α), the relative stiffness of the x-directional (or transverse directional) walls (β), and seismic excitation data applied to the model problem. Three different values (1.3, 1.6 and 2.0) were considered for α, and result in three elastic eccentricity values of the plan-asymmetric model problem, which are 0.65 m, 1.15 m and 1.67 m, respectively. Three values (i.e., 0.5, 1.0 and 2.0) were selected for β and can be categorized as torsionally unrestrained, normal, and torsionally restrained cases, respectively. Seven historical seismic records were selected as input ground excitations for the model problem. Time history analyses were performed for 63 total combinations of the test parameters.

Analysis Parameters
This section discusses various parameters that are considered in the time history analysis. As given in Table 2, the time history analyses were performed by changing three different types of parameters, which are the relative stiffness of the stiff side wall (α), the relative stiffness of the x-directional (or transverse directional) walls (β), and seismic excitation data applied to the model problem. Three different values (1.3, 1.6 and 2.0) were considered for α, and result in three elastic eccentricity values of the plan-asymmetric model problem, which are 0.65 m, 1.15 m and 1.67 m, respectively. Three values (i.e., 0.5, 1.0 and 2.0) were selected for β and can be categorized as torsionally unrestrained, normal, and torsionally restrained cases, respectively. Seven historical seismic records were selected as input ground excitations for the model problem. Time history analyses were performed for 63 total combinations of the test parameters. The detailed characteristics of the seven excitation data including the peak ground acceleration (PGA), the peak ground velocity (PGV), the peak ground displacement (PGD), differential time (DT) and total seismic excitation duration are summarized in Table 3. As can be seen from the table, they have different peak ground acceleration and displacement values. For example, the Parkfield excitation has the highest peak ground acceleration value, 0.476 g, and the Turkey excitation showed the largest peak ground displacement of 0.705 m. Figure 14 shows the acceleration and displacement spectra plotted using the seven excitation data.  The detailed characteristics of the seven excitation data including the peak ground acceleration (PGA), the peak ground velocity (PGV), the peak ground displacement (PGD), differential time (DT) and total seismic excitation duration are summarized in Table 3. As can be seen from the table, they have different peak ground acceleration and displacement values. For example, the Parkfield excitation has the highest peak ground acceleration value, 0.476 g, and the Turkey excitation showed the largest peak ground displacement of 0.705 m. Figure 14 shows the acceleration and displacement spectra plotted using the seven excitation data.

Analysis Results and Interpretation
In this section, the effectiveness of the inelastic torsional deformation assessment procedure proposed in Section 3 is verified by comparing its results with those of the time history analyses of the model problem. Figure 15 shows the lateral displacement history for α = 1.3, β = 1.0 under El Centro excitation. The lateral displacements of the two y-directional walls are plotted in Figure 15a. In the figure, the solid line represents the displacement of the flexible side wall Y 1 , and that of the stiff side wall Y 2 (dotted line). The difference between the lateral displacements of the two walls can be interpreted as the rotation of the asymmetric system. Figure 15b shows the lateral displacement history of the two x-directional walls, in which the lateral displacements of the two walls are of almost the same magnitude, but in the opposite direction. This indicates that the displacement history of the x-directional walls is mainly determined by the rotational behavior of the asymmetric system. In Figure 15c, the displacement of the flexible side wall Y 1 is decomposed into the component induced by lateral shear force and that by torsion, which are differentiated by a dark line and a light one, respectively. The sum of the two displacement components is equivalent to the solid line in Figure 15a.

Analysis Results and Interpretation
In this section, the effectiveness of the inelastic torsional deformation assessment procedure proposed in Section 3 is verified by comparing its results with those of the time history analyses of the model problem. Figure 15 shows the lateral displacement history for α = 1.3, β = 1.0 under El Centro excitation. The lateral displacements of the two y-directional walls are plotted in Figure 15a. In the figure, the solid line represents the displacement of the flexible side wall Y1, and that of the stiff side wall Y2 (dotted line). The difference between the lateral displacements of the two walls can be interpreted as the rotation of the asymmetric system. Figure 15b shows the lateral displacement history of the two x-directional walls, in which the lateral displacements of the two walls are of almost the same magnitude, but in the opposite direction. This indicates that the displacement history of the x-directional walls is mainly determined by the rotational behavior of the asymmetric system. In Figure 15c, the displacement of the flexible side wall Y1 is decomposed into the component induced by lateral shear force and that by torsion, which are differentiated by a dark line and a light one, respectively. The sum of the two displacement components is equivalent to the solid line in Figure 15a.  The results of Figure 15 are replotted in Figure 16, which presents the lateral displacement and rotation history estimated at the mass center of the asymmetric model problem. The maximum lateral displacement of 0.279 m occurs at 4.61 s, while the maximum rotation of −0.0119 rad occurs at 9.30 s. This clearly indicates that the two peak values do not occur simultaneously due to the dynamic effect on the structure by the excitation loading. The results of Figure 15 are replotted in Figure 16, which presents the lateral displacement and rotation history estimated at the mass center of the asymmetric model problem. The maximum lateral displacement of 0.279 m occurs at 4.61 s, while the maximum rotation of −0.0119 rad occurs at 9.30 s. This clearly indicates that the two peak values do not occur simultaneously due to the dynamic effect on the structure by the excitation loading.
(c) Displacement components of Wall Y1 induced by the lateral force and torsional moment. The results of Figure 15 are replotted in Figure 16, which presents the lateral displacement and rotation history estimated at the mass center of the asymmetric model problem. The maximum lateral displacement of 0.279 m occurs at 4.61 s, while the maximum rotation of −0.0119 rad occurs at 9.30 s. This clearly indicates that the two peak values do not occur simultaneously due to the dynamic effect on the structure by the excitation loading.  The maximum rotations of the 63 time history analysis simulations are plotted with respect to the relative eccentricity (e x /L) and relative stiffness in the transverse direction (β) in Figure 17a,b, respectively. It can be noted from Figure 17a that the maximum rotations of all 63 simulations show an increasing trend with increasing relative eccentricity. As already discussed in Section 4.2, β can be regarded as the torsional constraint ratio, and thus Figure 17b clearly indicates that the maximum rotations of the time history analyses decrease as β is increased. Table 4 provides the results of the time history analysis, the proposed method and the linear elastic torsional response computed for the 9 different combinations of α and β discussed in Section 4.2. The ground motion applied in the time history analysis is the El Centro excitation. In the time history analysis results of the table, the maximum displacement in the y-direction at the center of mass (∆ cen_TH ) was given, and then the rotation at the center of mass (ϕ cen_TH ) and the displacement on the flexible side wall Y 1 (∆ flex_TH ) corresponding to this displacement were listed. The maximum displacement ∆ cen_TH was taken as the target displacement in the results of the proposed method and the linear elastic response, and it was used to compute the rotation at the center of mass (ϕ cen_pro and ϕ cen_elas ) and the displacement on the flexible side wall Y 1 (∆ flex_pro and ∆ flex_elas ) by each of the two methods, respectively. In addition, the relative error of the displacement on the wall Y 1 was calculated by taking the corresponding time history analysis value ∆ flex_TH as the reference value. The average of the absolute relative errors for the other two methods are also provided at the end of the table.   Table 4 provides the results of the time history analysis, the proposed method and the linear elastic torsional response computed for the 9 different combinations of α and β discussed in Section 4.2. The ground motion applied in the time history analysis is the El Centro excitation. In the time history analysis results of the table, the maximum displacement in the y-direction at the center of mass (Δcen_TH) was given, and then the rotation at the center of mass (ϕcen_TH) and the displacement on the flexible side wall Y1 (Δflex_TH) corresponding to this displacement were listed. The maximum displacement Δcen_TH was taken as the target displacement in the results of the proposed method and the linear elastic response, and it was used to compute the rotation at the center of mass (ϕcen_pro and ϕcen_elas) and the displacement on the flexible side wall Y1 (Δflex_pro and Δflex_elas) by each of the two methods, respectively. In addition, the relative error of the displacement on the wall Y1 was calculated by taking the corresponding time history analysis value Δflex_TH as the reference value. The average of the absolute relative errors for the other two methods are also provided at the end of the table.   The results of Table 4 show that the proposed method slightly underestimated the displacement on the flexible side wall in general, and its average of the absolute relative errors (relative error average) is smaller than that of the linear elastic response. The reason why both the proposed and linear elastic methods underestimate the flexible side wall displacement seems to be that they are all based on simple static analysis procedure and have some limitations in accurately capturing the effect of complex dynamic behavior. Therefore, some safety margin needs to be introduced, if the proposed approach is applied to the design of an actual plan-asymmetric structure. The ratio between the two average values is 80.7% (=10.03/12.43 × 100%). A similar trend can be found from Table 5; Table 6, which provide the results obtained by applying Kobe and Northridge excitations, respectively. In this case, the relative error averages of the proposed method are even much smaller than those of the linear elastic response.
In the case of the Kobe excitation, the relative error average of the proposed method is only 55.4% of the linear elastic response value (=13.65/24.66 × 100%). Similarly, the former is only 71.4% of the latter (=10.28/14.39 × 100%) in the case of the Northridge excitation.    Tables 4-6 with respect to the target displacement (∆ cen_TH ) for the proposed method and the linear elastic response. A linear regression analysis is performed on each set of the results computed with the two methods by having the target displacement as an independent variable, and their regression lines are drawn in the figure. The coefficients of determination obtained by the regression results by the proposed and linear elastic methods are 0.0002 and 0.4198, respectively. It can be seen from the results of the figure that the linear regression line of the proposed method is located closer to the horizontal axis than that of the linear elastic response in all ranges of the target displacement, confirming that the proposed method can provide more accurate solutions than the linear elastic response. their regression lines are drawn in the figure. The coefficients of determination obtained by the regression results by the proposed and linear elastic methods are 0.0002 and 0.4198, respectively. It can be seen from the results of the figure that the linear regression line of the proposed method is located closer to the horizontal axis than that of the linear elastic response in all ranges of the target displacement, confirming that the proposed method can provide more accurate solutions than the linear elastic response. In Figure 19, the two sets of the results presented in Figure 18 are plotted with respect to the ratio of the stiff wall stiffness to that of the flexible wall (α), and the linear regression analysis is performed again. The coefficients of determination obtained by the regression results by the proposed and linear elastic methods are 0.0473 and 0.063, respectively. It can be seen from the figure that the linear regression line of the proposed method gets closer to the horizontal axis with increasing α. This means that the proposed method is able to provide more accurate solutions when the level of plan-asymmetry gets higher. Similarly, Figure 20 plots the same data as in Figure 18; Figure 19 with respect to β. The coefficients of determination obtained by the regression results by the proposed and linear elastic methods are 0.2191 and 0.0961, respectively. Since the regression line of the proposed method gets closer to the horizontal axis with decreasing β, the proposed method can provide more accurate solutions as the level of torsional constraint gets smaller, and thus the same conclusion can be made as in Figure 19. In Figure 19, the two sets of the results presented in Figure 18 are plotted with respect to the ratio of the stiff wall stiffness to that of the flexible wall (α), and the linear regression analysis is performed again. The coefficients of determination obtained by the regression results by the proposed and linear elastic methods are 0.0473 and 0.063, respectively. It can be seen from the figure that the linear regression line of the proposed method gets closer to the horizontal axis with increasing α. This means that the proposed method is able to provide more accurate solutions when the level of plan-asymmetry gets higher. Similarly, Figure 20 plots the same data as in Figure 18; Figure 19 with respect to β. The coefficients of determination obtained by the regression results by the proposed and linear elastic methods are 0.2191 and 0.0961, respectively. Since the regression line of the proposed method gets closer to the horizontal axis with decreasing β, the proposed method can provide more accurate solutions as the level of torsional constraint gets smaller, and thus the same conclusion can be made as in Figure 19.

Conclusions
In this paper, we proposed a simple method that can accurately assess the inelastic lateral displacement and rotation of plan-asymmetric wall structures. The relationship between the lateral displacement and rotation was obtained by considering the inelastic behavior of asymmetric structures. This relationship was presented and interpreted in the displacement-rotation (D-R) coordinate system. This coordinate system was first introduced in this paper and very conveniently describes the inelastic torsional behavior of the asymmetric structure. The validity of the proposed method was verified by comparing its predictions and inelastic time-history analysis results for a model problem. The main conclusions of this paper are as follows:

Conclusions
In this paper, we proposed a simple method that can accurately assess the inelastic lateral displacement and rotation of plan-asymmetric wall structures. The relationship between the lateral displacement and rotation was obtained by considering the inelastic behavior of asymmetric structures. This relationship was presented and interpreted in the displacement-rotation (D-R) coordinate system. This coordinate system was first introduced in this paper and very conveniently describes the inelastic torsional behavior of the asymmetric structure. The validity of the proposed method was verified by comparing its predictions and inelastic time-history analysis results for a model problem. The main conclusions of this paper are as follows: 1.
Most of the current seismic design codes such as IBC, NZS, NBCC and Eurocode are based on the force-based design approach and introduce the concept of the design eccentricity consisting of static and accidental eccentricities. They do not properly consider the inelastic deformation and actual failure mechanism of plan-asymmetric structures. As a result, they may require excessive strength and stiffness for lateral force resisting components.

2.
The inelastic torsional response assessment procedure proposed in this paper shows that the linear elastic approach may both underestimate the torsional rotation at the initial yielding stage and overestimate it at the post yielding stage, resulting in inaccurate estimation of the torsional response over the entire inelastic range.

3.
The stiffness of the transverse wall may affect the overall yielding procedure of asymmetric structures. As a result, first yielding may occur even in the stiff side wall if the transverse wall stiffness and yield displacement are determined mainly by its length. This effect is not considered by any of the current code provisions such as IBC and Eurocode, while the procedure proposed in this study can handle it properly.

4.
A comparison between the results of the time history analysis and the predictions of the proposed method for the model problem under three different seismic excitations (El Centro, Kobe and Northridge) shows that the relative error averages of the proposed method range approximately from 10.0% to 13.7% and are considerably smaller than those of the linear elastic response. It also indicates that the proposed method is able to provide more accurate solutions if the level of plan-asymmetry is higher.

5.
In general, the proposed method underestimates the flexible side wall displacement as it is based on simple static analytical procedure and has some limitations in accurately capturing the effect of complex dynamic behavior. Therefore, some safety margin needs to be introduced, if it is applied to the design of an actual plan-asymmetric structure.
Currently, we are extending the inelastic torsional response assessment approach proposed in this paper and developing an iterative method utilizing the secant stiffness of each wall for the effective estimation of the inelastic torsional behavior of complex plan-asymmetric structures. This approach also can be incorporated into the adaptive push-over analysis discussed in [31][32][33][34] that explicitly considers the nonlinear torsional and higher-mode effects for the analysis of multi-story structures.